Changes in nocturnal insect communities in forest-dominated landscape relevant to artificial light intensity

Artificial light at night has recently been identified as a major factor adversely affecting global insect diversity. Here, we compared the insect diversity in Gwangneung Forest Biosphere Reserve, specifically in the Korea National Arboretum (with no artificial light at night), with that of three nearby urban sites with a gradient of artificial light at night (five locations at each site). We analyzed the effects of the artificial night lighting index, mean annual temperature, and field light intensity (lux) at night on the insect community structure. The urban sites generally exhibited higher species richness and abundance as well as clear indicator species compared with the control site. The size distribution of the collected insects markedly differed between the control and the three urban sites. The abundance of herbivorous and omnivorous insects increased and decreased, respectively, with the increase in light intensity. Species richness of herbivorous and omnivorous insects was likely correlated with the field light intensity at night and artificial night lighting index, respectively. This study demonstrates the association between nighttime environment and marked changes in insect community structure and revealed consequent transition of ecosystem services by changes in trophic group composition.


Background
Insect species diversity and abundance are showing rapidly declining trends in all habitats (Hallmann et al. 2017). Although habitat loss, pollution, invasive species, and climate change have recently been identified as major drivers of the decline, the effect of artificial light at night (ALAN), which has never occurred before in evolutionary time scales, is relatively overlooked (Owens et al. 2020). ALAN, from streetlights, car headlights, and signs, is correlated with urban expansion (Elvidge et al. 2014) and is now widespread (Falchi et al. 2016;Davies and Smyth 2018). Approximately 90% of Europe and 50% of the USA are exposed to various light sources at night (Falchi et al. 2016). In Korea, the ratio of land area exposed to light pollution is the second largest in the world, and the geographical and social characteristics of small land area, urbanization, and high-density population are the main causes of light pollution (Government of Korea 2018). Furthermore, approximately 60% of invertebrate species worldwide are active at nighttime (Hölker et al. 2010). Thus, ALAN can affect not only the biological traits of these species, such as physiology (Eisenbeis and Hänel 2009), behavior (Yoon et al. 2010;Russo et al. 2019), and reproduction (Witherington 1992;Mbugua et al. 2020) but also their ecological communities (Davies et al. 2012). However, because of a bias toward studies on diurnal species, the effect of ALAN on whole ecological communities remains relatively less extensively studied Owens et al. 2020).
Despite the lack of strong evidence demonstrating the ecological effect of ALAN on invertebrates, some recent studies have reported changes in insect communities due to ALAN (Davies et al. 2012;Holzhauer et al. 2015;Spoelstra et al. 2015;Plummer et al. 2016;De Medeiros et al. 2016;Firebaugh and Haynes 2020). For example, ALAN has affected the composition of arthropod communities, leading to increases in the abundance of predatory and scavenging arthropods (Davies et al. 2012). Holzhauer et al. (2015) reported that specific groups of insects [e.g., beetles (Coleoptera), flies (Diptera), and bugs (Hemiptera)] were more abundant at sites with artificial light than in those without artificial light and that such sites were characterized by high insect abundance compared with those without artificial light. Plummer et al. (2016) andDe Medeiros et al. (2016) demonstrated that radiations from ultraviolet light at night attracted large numbers of moths (Lepidoptera) and beetles (Coleoptera), respectively, likely reducing insect diversity. Furthermore,  showed that ALAN increased the density of aphids, changing the food web structure.
In contrast to studies reporting changes in insect communities due to ALAN, there were also studies that reported no effective change due to ALAN. For example, moth communities sensitive to light were not influenced by ALAN (Spoelstra et al. 2015). Firebaugh and Haynes (2020) monitored changes in an arthropod community under ALAN for 2 years, and they did not observe any effect of ALAN on the abundance of herbivores or predators. Studies reporting the effect of ALAN on insect communities were generally conducted by controlled experiment by installation of artificial light sources; thus, field experiments in real-world situation were needed.
The present study was conducted in Gwangneung Forest Biosphere Reserve, one of the major areas for biodiversity protection in South Korea. The reserve consists of an area without artificial lighting (Korea National Arboretum) and several urbanized and light-polluted areas. The Korea National Arboretum has not contained ALAN since the 1990s, although the night light conditions in areas adjacent to the arboretum have changed considerably over the past few decades. This environmental gradient has been maintained across Gwangneung Forest Biosphere Reserve for a long time, and therefore, the insect communities present in these areas reflect the ecological effects of ALAN as well as urbanization on insect communities over the last few decades.
In the present study, we identified the ecological effect of ALAN on insect communities in Gwangneung Forest Biosphere Reserve including a control site (no ALAN) and three nearby sites along a gradient of development situations (DS). We aimed to determine the following: (1) how insect species richness and abundance and indicator species differ along the gradient of sites; (2) which of the potential environmental factors [e.g., temperature, artificial night lighting index (ANLI), and field light intensity] cause differences in the community structure of nocturnal insects in the study area; (3) how distribution of insect body size differs along the gradient of sites; (4) how the composition of trophic groups differs along the gradient of sites; and (5) which factors significantly affect changes in the composition of trophic groups.

Study area
We conducted field surveys in Gwangneung Forest Biosphere Reserve in June and August 2017. The mean annual temperature and precipitation for this area were 12.2°C and 1387.9 mm, respectively (Cho et al. 2020).

Site description
We selected the exhibition garden in the Korea National Arboretum as a control site, and three urban sites with different DS, namely, Neungnae Village (DS1), Jikdong Samgeori (DS2), and Gwangneungnae (DS3) (Fig. 1). The area of each site was 100 ha. The distance from the center of the control site to the nearest light source (convenience store) is 1 km, and the streetlight on the road where cars sometimes travel is 600 m. The other three sites received diverse levels of ALAN from various sources, such as streetlights, shops, residential areas, and vehicular traffic. The DS1 site was surrounded by a forest and included commercial and residential buildings with some two-way, two-lane roads and streetlights. There was no artificial lighting on the roads outside the DS1 site. The DS2 site was surrounded by residential buildings as well as many lodging establishments and shops that were well lit at night. The DS3 site was located in an urban area with a complex road network and densely packed apartment complexes and commercial facilities.
We measured the ANLI and mean annual air temperature (°C) at each site. The ANLI was calculated by dividing the area occupied by shops and residential buildings, including streetlights and roads, by the total area, based on high-resolution satellite images from Arc-GIS 10.3. ANLI ranged from 0 to 100, where 0 indicated no ALAN and 100 indicated the maximum levels of ALAN. For example, if there is a 100-ha area composed of 30% of vegetation with no light and 70% of commercial areas and roads, ANLI is calculated as 30. ANLI as simple index is equal to vegetation area in the specific area. We applied simple approach because it is very difficult to quantify ALAN (e.g., type and amount of artificial light at night) in complex urban situation. To incorporate temperature changes during industrialization of this area, we determined the mean annual air temperature in the southern Korean peninsula between 1970 and 2014, using the map of mean annual air temperatures (Arc/ Info binary grid), which is stored at the Korea National Arboretum. This was applied to each site. We also measured the field light intensity (lux) at night using a light meter (LX-1332, Custom Corp., Tokyo, Japan) at 50 locations at each site (100 ha) divided into 100 grid cells. We randomly selected 50 grid cells to measure lux using the random-number generator in Excel 2016.

Insect collection
At each study site, insects were collected along roads and at the edges of commercial areas where plants such as Rhododendron spp. and Buxus microphylla var. koreana were found. These sites have different levels of human activity and consumption of plants by insects. Two nighttime surveys of insects were conducted in 2017, i.e., before (June) and after (August) the rainy season. Each survey was conducted on the same day of the month (8 th ) in the lunar calendar between 21:00 and 23:00, to account for natural nighttime light conditions and insect activity. For each site, we conducted three round trips of five belt transects (5 × 20 m) and collected insects using a sweep net (diameter 45 cm) to assess species richness and abundance. Insect abundance data across the two seasons (before and after the rainy season) were pooled for analysis.

Insect species richness and abundance and indicator species along the gradient of sites
To identify differences in insect species richness and abundance along the gradient of sites, we identified the proportional richness and abundance of insect orders at each site. Then, we used Venn diagrams (R package "VennDiagram") (Chen and Boutros 2011) to illustrate the proportion of shared and unique insect species at the sites. Finally, we quantitatively compared insect species richness and abundance at the sites.

Differences in insect community composition
Based on the insect abundance data at each site, we determined the patterns of turnover in insect communities and determined the factors that contributed most to the changes. Meta-non-parametric multidimensional scaling was performed using the R package "vegan" (Oksanen et al. 2019) to explore the turnover patterns. Of the environmental factors (i.e., ANLI, mean annual air temperature, and field light intensity), those that best explained the changes with high predictability (p < 0.001) were exclusively displayed.

Size class distribution by site
We referred to information from online databases (e.g., Korea Biodiversity Information System) for data on the size of the collected insects. Except for earwigs (Dermaptera), we defined insect body size as the distance from the front of the head to the tip of the abdomen. The distance from the front of the head to the base of the forceps was measured for determining the body size of earwigs. We also analyzed patterns in the collected insects for abundance distribution by size class and checked the differences in the abundance distribution compared to the control site, using the two-sample Kolmogorov-Smirnov test.

Insect trophic structure by site
We classified the collected insects according to the trophic groups (i.e., herbivore, carnivore, omnivore, and detritivore). The classification was based on the investigators' knowledge and discussions, plus information from online databases (i.e., Korea Biodiversity Information System, daum.net, naver.com, and google.com).

Statistical analysis
The mean annual air temperature and field light intensity for each site and the species richness and abundance data for the collected insects were not normally distributed, so the differences in data between the sites were analyzed using the rank-based non-parametric Kruskal-Wallis test followed by a Bonferroni post hoc test. Bonferroni correction (α = 0.05/ number of tests) was considered to avoid type I errors. All statistical analyses were conducted in R software (R Core team 2019). Data are expressed as means ± standard deviations.

ANLI, temperature, and light intensity
Overall, the ANLI and field light intensity at the study sites increased with increasing development. The ANLI was the highest in DS3 (56.5), followed by DS2 (30.5), DS1 (17.1), and the control had an ANLI of zero (Table 1). The urban site DS2 (median: 11.80°C) had the highest mean annual air temperature, followed by DS3 (11.58°C), the control (11.36°C), and DS1 (11.28°C). However, these differences were not statistically significant. Field light intensity differed significantly among the sites (χ 2 = 125.48, df = 3, p = 0.001) and was the highest in DS2 (12.2 lux), followed by DS3 (8.7 lux), DS1 (5.2 lux), and the control (0 lux). Field light intensity showed considerable variance depending on the distance from the light sources, such as vehicular traffic and especially streetlights.

Proportional insect species richness and abundance of insect orders by site
The Gwangneung Forest Biosphere Reserve has long maintained the gradient of DS. The control site and three nearby sites (DS1, DS2, and DS3) showed notable differences in proportional insect species richness and abundance ( Fig. 2; Supplementary Table S1). This was especially evident with the proportion of hemipteran species increasing along the gradient of DS (Fig. 2). Neuropteran species were collected only from urban sites (DS1, DS2, and DS3), while dermapteran species were unique to the control site.
The number of insect species common to the control and urban sites was markedly lower than that common to the urban sites only (3 vs. 16, Fig. 3). The only three species common to the control and urban sites were Metcalfa pruinosa, Monolepta shirozui, and Pristomyrmex pungens, and the number of species unique to the control, DS1, DS2, and DS3 sites was 33, 22, 48, and 15, respectively (Fig. 3).

Insect community composition and size distribution
When the differences in species composition at each study site were analyzed by rank, the NMDS1 differentiated the control site and three urban sites (DS1-3), and their spatial distribution was significantly correlated with ANLI (r 2 = 0.582, p = 0.002) (Fig. 5).
By applying the two-sample Kolmogorov-Smirnov test between the control site and each of the three urban sites (DS1, DS2, and DS3), we confirmed significant differences in the size class distribution of insects (control-DS1: Kolmogorov-Smirnov (KS) = 0.12, p < 0.0001; control-DS2: KS = 0.17, p < 0.0001; control-DS3: KS = 0.12, p < 0.0001) (Fig. 6). Compared to the control site, Fig. 2 A Proportional richness and B abundance of insect orders at the control site and the three urban sites (DS1, DS2, and DS3). Raw data is presented in Supplementary Table S1. Refer to Fig. 1 for site descriptions and locations. He, Hemiptera; N, Neuroptera Fig. 3 Venn diagram showing groups of insect species shared between study sites (control, DS1, DS2, and DS3). Refer to Fig. 1 for site descriptions and locations the urban sites showed considerably increased numbers of small insects (≤ 15 mm). In particular, the control site showed a reverse-J-shaped distribution of insect sizes, but the three urban sites showed right-skewed unimodal distributions (Fig. 6).

Insect trophic structure along the sites
The collected insects were divided into four trophic groups (herbivores, omnivores, carnivores, detritivores), and the median and proportional species richness and abundance of these groups were compared (Fig. 7). Of the four trophic groups, only herbivores (richness: H = 7.944, df = 3, p = 0.047, abundance: H = 13.316, df = 3, p = 0.004) and omnivores (richness: H = 13.428, df = 3, p = 0.004, abundance: H = 10.895, df = 3, p = 0.012) showed significant differences in species richness and abundance between study sites (Fig. 7A). However, pairwise comparisons between study sites were not significant for the two trophic groups. The proportional species richness and abundance of herbivores were approximately 54% and higher at the urban sites with ALAN, compared with those at the control site. The proportion of omnivores at the urban sites was lower than that of the control site, and in DS3 specifically, which was experiencing the highest development pressure, no omnivores were observed (Fig. 7).

Discussion
The direct relationship between the effect of ALAN on insects and urban expansion (Elvidge et al. 2014) is implicated in a reduction in insect diversity (Owens et al., 2020); however, community-based studies in real-world ecosystems are limited. We demonstrated that ALAN affected the insect community structure (species diversity, Fig. 4 A Differences in insect species richness and B abundance among the four study sites (control, DS1, DS2, and DS3). Refer to Fig. 1 for site descriptions and locations. Bars and error bars indicate mean ± standard deviations, and pairwise comparisons of richness and abundance between the sites are not statistically significant by Bonferroni correction (p > 0.008) Fig. 5 Results of meta-NMDS ordination using abundance data for communities of insects collected from each study site (control, DS1, DS2, and DS3). Refer to Fig. 1 for site descriptions and locations. A significant correlation between insect communities and environmental factors (ANLI, temperature, lux) was only observed for ANLI (r 2 = 0.582, p = 0.002) Fig. 6 Size class distribution of insects collected from each site. The distribution in DS1, DS2, and DS3 all showed significant differences to that of the control site using the two-sample Kolmogorov-Smirnov test. Refer to Fig. 1 for site descriptions and locations size distribution, and trophic level) in Gwangneung Forest Biosphere Reserve, where a constant light regime has been maintained at night for the last few decades. We considered one dark site and the nearby three sites with a gradient of ALAN. Specifically, the insect communities showed an over-dominance of small insects (≤15 mm) and herbivore groups, related to an increasing gradient of DS and ALAN intensity.
Unlike previous studies that manipulated ALAN (Davies et al. 2012;Holzhauer et al. 2015;Spoelstra et al. 2015;Plummer et al. 2016;De Medeiros et al. 2016;Firebaugh and Haynes 2020), the present study provides evidence that the change of night light environment results in the compositional change of nighttime collected insect communities in real-world ecosystems. The three urban sites with ALAN (DS1-3) showed higher heterogeneity in nocturnal insect species composition than the control site, and this was correlated with ANLI, an index of artificial light.
Research on insect size class distribution is limited. In the present study, the distribution of the size classes under natural nighttime conditions showed a reverse-J shape. Unlike the control site, a marked increase in the number of small insects at nighttime in developed environments resulted in a unimodal distribution. Biometric variables, such as dry mass and length, are closely related to ecological functions (Hallmann et al. 2017). The change in size classes revealed in the present study may potentially be related to changes in species composition and function in groups. In studies on insect communities, it is important to determine the values of variables such as organism size (e.g., biomass) and evaluate their correlations with other structures and functional elements (Shortall et al. 2009). In addition, more attention should be paid to the body size distribution in order to understand the relationships between structural attributes and habitat types, such as urban, agricultural, and remote areas (Chown and Gaston 2010). Fig. 7 A Differences in insect species richness and B the abundance between the study sites for each trophic group and C proportional richness and D the abundance by trophic group between study sites. Asterisks indicate significant differences (p < 0.05). The results of multiple pairwise comparisons between the study sites were not significant (A and C). The error bars represent standard deviation. Site descriptions and locations are presented in Fig. 1 The large increase in the herbivore population at the three urban sites was the result of reduced overall predation pressure due to lower numbers of omnivores. This increase in herbivore numbers could be attributed to their attraction to nitrogen-and phosphorus-rich food (leaves), which is a commonly reported limiting factor for insects in urban environments (Moreira et al. 2019). The reduction in omnivores was the main factor responsible for the simplification of the food web. Omnivores, which comprise approximately 40 families in 12 orders, have generally adapted to supplemental nutrition that is lacking in a plant-based diet and are an important element in natural and agricultural ecosystems (Krimmel 2011). Because omnivores perform the roles of both herbivores and carnivores within the food web, they can consistently alter the intensity and direction (top-down or bottom-up) of trophic cascades. In the present study, it is likely that the three urban sites contained highquality edible leaves. In these areas, the increased abundance and richness of herbivores may lead to increased spatial competition between herbivores and omnivores and thus affect the extent to which omnivorous insects could enter the area. Moreover, with regard to omnivore-herbivore-plant interactions, the presence of high-quality plant-based food reduces predatory behavior, rather than increasing the attraction for omnivores. This reduction, in turn, is beneficial to herbivores (Eubanks and Denno 2000). Omnivores can also enhance top-down control by pressuring predator groups, but this pattern was not observed in the present study. Future studies that assess the effects of high-quality plant resources on competition between herbivores and omnivores may provide important insights into the food web changes associated with urbanization.
Although many studies have reported nighttime urban light pollution using satellite images Sutton et al. 2007;Gaston et al. 2015), the impact on biodiversity at the community level has not been fully clarified. However, it is an extremely challenging task to understand the effects of nighttime light on biological communities in urban environments because of the diversity and interconnectedness of factors within urban nighttime environments. In the present study, we used ANLI, mean annual air temperature, and field light intensity at the survey sites as the factors characterizing the urban nighttime environment. However, this approach and these factors are not sufficient to explain the causal relationships between ALAN and the observed changes in the insect communities clearly.