Moth species richness and diversity decline in a 30-year time series in Norway, irrespective of species’ latitudinal range extent and habitat

Insects are reported to be in decline around the globe, but long-term datasets are rare. The causes of these trends are elusive, with changes in land use and climate among the top candidates. Yet if species traits can predict rates of population change, this can help identify underlying mechanisms. If climate change is important, for example, high-latitude species may decline as temperate species expand. Land use changes, however, may impact species that rely on certain habitats. We present 30 years of moth captures (comprising 97,032 individuals of 808 species) from a site in southeast Norway to test for population trends that are correlated with species traits. We use time series analyses and joint species distribution models combined with local climate and habitat data. Species richness declined by 8.2% per decade and total abundance appeared to decline as well (−9.4%, p = 0.14) but inter-annual variability was high. One-fifth of species declined, although 6% increased. Winter and summer weather were correlated with annual rates of abundance change for many species. Opposite to general expectation, many species responded negatively to higher summer and winter temperatures. Surprisingly, species’ northern range limits and the habitat in which their food plants grew were not strong predictors of their time trends or their responses to climatic variation. Complex and indirect effects of both land use and climate change may play a role in these declines. Our results provide additional evidence for long-term declines in insect abundance. The multifaceted causes of population changes may limit the ability of species traits to reveal which species are most at risk.


Introduction
Insects are reported to be in decline around the globe, based on evidence spanning many different ecosystems and at both local and regional scales (Sánchez-Bayo and Wyckhuys 2019, Didham et al. 2020, van Klink et al. 2020, Wagner 2020). These declines may be due to a host of factors, including land use changes (especially agricultural intensification), pesticides and pollution (including light pollution), introduced species, and climate change (Habel et al. 2019a, Seibold et al. 2019, Raven and Wagner 2021.
Detecting and measuring declines, however, requires long-term datasets, and these are rare (Didham et al. 2020, Montgomery et al. 2020. Continuous long-term datasets, rather than comparisons between several time points (Seibold et al. 2019), are especially important for taxa (like many insects) that exhibit population cycles and high inter-annual variability (Welti et al. 2020). Unfortunately, there are few such datasets from Fennoscandia, even though high-latitude systems are experiencing some of the strongest effects of climate change (Hunter et al. 2014, Loboda et al. 2017. Insects are a hyper-diverse class, and no single sampling scheme can adequately monitor all groups (Montgomery et al. 2020). Lepidopterans, however, may be good early indicators of broader changes to insect communities (Sánchez-Bayo and Wyckhuys 2019, Montgomery et al. 2020), in part because they may be particularly susceptible to environmental change (Hunter et al. 2014). They are also important for pollination, a valuable ecosystem service (Macgregor et al. 2015). Moths, which are much more diverse than butterflies, serve as an important food source for birds, bats and other insectivorous taxa (Sánchez-Bayo and Wyckhuys 2019).
Lepidopterans may also be among those insect groups that are declining most rapidly (Sánchez-Bayo and Wyckhuys 2019). For example, moths and butterflies in the Netherlands (Van Dyck et al. 2009, Groenendijk and Ellis 2011, Langevelde et al. 2018) and the United Kingdom (Conrad et al. 2004, Roy et al. 2015, Dennis et al. 2019 appear to have experienced considerable declines in both abundance and species richness in the last 30 years. Historic baselines are inevitably relative, meaning that longer time series are more informative than shorter ones (Macgregor et al. 2019b(Macgregor et al. , 2021 and critical analyses and interpretation are required . Farther north, there is some evidence from Finland that total moth abundance has declined over the last 20 years, although species richness appears to have increased due to range expansions of southern species (Antão et al. 2020). Swedish butterflies inhabiting meadows and grazed open forests also appear to have declined in abundance during this same period Johannesson 2007, Nilsson et al. 2013).
Even as these trends in richness and abundance have been documented, their causes remain elusive; correlations with proposed causal factors provide only weak evidence. However, if trends of the individual species that make up these overall abundance patterns can be predicted using their traits (Mattila et al. 2006, Habel et al. 2019a, this can provide important clues to the underlying processes. Traits can help link species responses to land use changes , Coulthard et al. 2019, and can also help predict responses to climatic variation in the short term (Roy et al. 2015), reducing the noise in inter-annual population variability and providing clues to the role of climate in long-term trends in abundance.
Range limits are a trait that can be used to better understand a species' response to climate change (Kirkpatrick and Barton 1997, Sexton et al. 2009. In moths, for example, species with more northerly ranges (Itämies et al. 2011) are more likely to be declining in some regions. Many of the most common species in northern Finland, however, appear to have been stable or even to have increased over the last several decades (Hunter et al. 2014). If climate change is causing long-term poleward shifts in moth species' ranges, as has been demonstrated in some other insects (Parmesan et al. 1999, Pöyry et al. 2009, Boggs 2016, then many species in the Northern Hemisphere may be becoming more abundant near their northern range limits, or declining at their southern range edges (Warren et al. 2001). Range limits may also shift as a function of changing abundance of a species (Mair et al. 2014, Macgregor et al. 2019a.
If, on the other hand, land use change is an important contributor to population changes, these changes might be predicted by the habitat requirements of a species' preferred food plant (Habel et al. 2019b). Moth species with higher host plant specificity are more likely to be declining than other species in Finland (Mattila et al. 2006) and Sweden (Franzén and Johannesson 2007), and biotic homogenization appears to be occurring in Hungarian moths as species that require grasslands or specialized diets are disappearing (Valtonen et al. 2017). Land use and climate change are not independent, and they have been shown to interact in affecting species range limits and shifts (Mair et al. 2014, Platts et al. 2019. Here we examine a 30-year moth capture dataset, in which 808 species were collected using consistent methods in southeast Norway. We use time series analyses and Bayesian community modeling to ask: (a) How do moth community richness, diversity, and abundance vary through time? (b) Do rates of abundance change differ among species with different range extents or food plant habitats? (c) On an annual basis, are species' range limits good predictors of their responses to climatic variation?

Study area
Moths were trapped each summer from 1984 to 2013 at a site (59.7477° N, 10.5925° E, 70 m) in Nesodden municipality, southeast Norway (Kobro 1991) (Fig. 1). The trap site was located on a peninsula in the Oslo Fjord, 20 km south of Oslo (Fig. 1). It was at the edge of a garden, shaded from most moonlight, and surrounded by deciduous, mixed, and coniferous forest. The locations of trees and buildings in the immediate area did not change through the years. The landscape within five kilometers of the trap site is largely forested (70%), and much of the remaining land (22%) was farmland and open areas in 1984 based on our comparison of Landsat images using the rpart R-package (Therneau and Atkinson 2019) (see Supplementary Methods). Forest cover remained relatively constant throughout the trapping period, but by 2014 half of the fields and farmland had been developed, largely for housing. These changes are similar to those that have occurred within a larger 25 km radius (Supplementary Table S1; Supplementary Fig. S1). Within the forests of surrounding Akershus county, total wood volume and tree species composition were similar in 1983 and 2013, although there was a slight shift towards more mature forests (Statistics Norway; www. ssb. no).

Sampling methods
A funnel light trap (160 W/235 V mixed light bulb) was hung 1 m above the ground in the same location through the night for the first three nights of each week from early June to mid-October, 1984-2013 as described by Kobro (1991).On nights with rain, wind, or otherwise poor trapping conditions, that night was dropped and the trap was deployed an additional night. Moths were identified using genitalia and external characters by SK. Taxonomy follows the gbif taxonomic backbone (GBIF Secretariat 2021). To summarize the abundance of each species for a given year, we filtered the dataset to include only nights within the range of dates that were sampled every year (02 July-16 October). This left an average of 42.7 sampling nights per year (range 36-46). We then divided total counts of each species per year by the number of sampling nights in that year to get mean number of individuals per night, and multiplied this value by the mean number of sampling nights (42.7), rounded to the nearest integer, to get comparable estimates of abundance among years.
The filtered dataset included 808 species after excluding those species captured only outside of the core summer trapping season (n = 74).

Species traits
To determine if species' trends and responses to weather were correlated with their geographic ranges, we assigned 336 species (those captured in 10 or more years) to one of three groups based on their ranges in Norway, using information from Bengtsson et al. (2008), Johansson (2011), Svensson (2006), and Artsdatabanken (www. artsd ataba nken. no). These included: 'southern' species (n = 90), which occur only in the south of Norway (below 60.5° N), 'central' species (n = 113), which occur as far north as central Norway (64° N), and 'widespread' species (n = 133) which occur the entire length of the country (Fig. 3).
To test for relationships between species trends and land use changes, a list of preferred food plants was also obtained for each species from Goater et al. (1986), Palm (1989), Svensson (2006) (2011), Sterling and Parsons (2012), and the authors' observations. Using our knowledge of Norwegian plants, we then classified these groups of host plants as occurring primarily in: (a) forests (n = 162), (b) gardens (n = 24), (c) meadows (including agricultural field edges; n = 67), (d) meadows and forests (n = 51), or (e) unknown/other (n = 32).

Testing for trends
We summed the total numbers of species in the dataset by year (using the full dataset, with all 808 species from the core sampling period included). For total individuals, we used the number of individuals corrected for effort as described above. We also estimated true species richness in each year using the iChao1 metric in the SpadeR R-package (Chao et al. 2016) in R (R Core Team 2021). This metric uses the number of singletons and doubletons in a dataset to estimate the number of present but unobserved species. Finally, we calculated the Shannon diversity index for each year using the vegan R-package (Oksanen et al. 2017).
For each of these community metrics we tested for evidence of a trend through time using a modified Mann-Kendall test in the modifiedmk R package (Patakamuri and O'Brien 2020). This test is robust for highly variable nonnormal count data and estimates the slope of the trend line, if significant (p < 0.05). The test requires that data are not highly autocorrelated, and we found that serial autocorrelation in our richness and diversity data was significant only at a time interval of one year (p < 0.05), making this an appropriate test. For plotting purposes the intercept of the trend line (the estimate for year zero, 1984), which isn't provided by the Mann-Kendall test, was calculated such that the line would pass through the median value of the response variable in year 15 of 30. We expressed rates of increase or decline as the percent change across a 10-year period relative to the intercept (Sánchez-Bayo and Wyckhuys 2019). Similarly, we also tested for a trend in each species that was detected in ten or more years (n = 336), because detecting a significant trend is difficult in species that are present at fewer time steps. Species trends were expressed as change per decade as a % of the median annual count value for that species, because many species were rare or absent early or late in the study and so a rate based on intercept values could be misleading. In cases where the median count was less than three (due to the presence of many zeros), we used the median value when present to avoid inflating rates of change. Moving averages for plotting were calculated using the 'sma' function in the smooth R package (Svetunkov 2020), which uses BICc to determine the number of years to include in the moving average.

Climate covariates
We summarized climate covariates for the area surrounding the trap site for each summer (June-August) and winter (December-March) of the study period using estimates from the ERA5 climate reanalysis (Copernicus Climate Change Service [C3S] 2017). Values were taken from the grid cell (0.5° latitude x 0.5° longitude) surrounding the trap site. These included mean summer air temperature (at 2 m), cumulative summer precipitation, mean winter soil temperature (at 1 cm depth), and mean winter snow depth. We chose mean winter soil temperature, rather than air temperature, to better represent the temperature that overwintering eggs, larvae and pupae were subjected to.

Responses to covariates
To test the hypothesis that population growth rates respond to inter-annual variation in weather, with southern species benefiting from more benign weather, we fit Bayesian hierarchical joint species distribution models (JSDMs) using the Hmsc R-package (Ovaskainen et al. 2017a(Ovaskainen et al. , 2017b. We modeled population rate of change (capture rates relative to previous year) as a function of summer and winter temperature and precipitation, based on estimates from the ERA5 climate reanalysis (C3S 2017). Determining a population rate of change for a given year requires capture records from that year as well as from the previous year; we included all species for which rates of change could be calculated in ten or more years (n = 184). To calculate rate of population change r, we used log-transformed abundance data for each count C of species s in year t: For species that were not captured either in a given year or in the preceding year, the rate of population change response variable in that year was treated as missing (NA). These JSDMs jointly estimate the response of each species to the environmental covariates, incorporating trait (species range) and phylogenetic information as a hierarchical level to estimate whether functionally or phylogenetically related species respond to the environment in similar ways. We included a phylogeny based on the species-level insect tree from Chesters (2017) to test for phylogenetic signal in the responses. Missing species were added to the correct genus (when present) or family using the ape R package (Paradis and Schliep 2019). Year was included in models as a temporally structured random effect. We fit a series of six community models (Table 1) with Gaussian error distributions and an identity link function. Models differed in their combinations of environmental covariates (including covariates from the previous year, to test for lag effects). Models were fitted using five MCMC chains with 6800 iterations each. The first 5000 iterations of each chain were discarded as burn-in, and the remaining iterations were thinned by 6 to produce 300 samples per chain (1500 in total). Model fitting was conducted with high performance computational resources provided by Louisiana State University (http:// www. hpc. lsu. edu), and models were ranked using the widely-applicable information criterion (WAIC) from Watanabe (2010).

Results
We recorded 808 moth species and 97,032 individuals in 30 years of sampling during the core summer trapping period. These species came from 43 families in 14 superfamilies (Supplementary Table S2). and were all used to calculate the community metrics. An additional 74 species (11,308 individuals) were detected only outside the core summer trapping period and were excluded from all analyses. Annual observed species richness ranged from 217 to 429 (mean = 284.9, sd = 48.5). Of 808 total species, 336 were captured in ten or more years and were used to calculate species-specific metrics.

Trends in community metrics
Observed and estimated species richness, as well as diversity, decreased through the study period (Fig. 2), based on Mann-Kendall tests (alpha = 0.05). Observed species richness declined by 8.2% per decade (p < 0.01), and abundance of individuals appeared to decline as well (−9.4% per decade; p = 0.14); most seasons in the later years of the study included 50-100 fewer species than most seasons earlier in the study.
Moth species were assigned to three groups (southern, central, and widespread species) based on their ranges in Norway (Fig. 3), because of the expectation that more southerly species may respond more positively to climate change (Warren et al. 2001). When these groups with different ranges were considered separately, species richness declined for each group (range: −5.7% to −7.3% per decade; Supplementary Fig. S2). The number of individuals captured declined (−17.0% per decade) for widespread species (driving the decreasing trend in overall abundance), whereas we did not detect overall abundance changes in the other groups ( Supplementary Fig. S3).

Trends in individual species
We also tested for time trends in the 336 species that were detected in ten or more years (Supplementary  Table S3). Of these, 18% (n = 59) show a negative trend Fig. 2 Trends in total observed richness, estimated species richness (iChao1), total individuals, and Shannon diversity from a 30-year moth capture dataset from southeastern Norway (10year rates of decline: −8.2%, −8.1%, −9.4%, and −5.0%, respectively). Plots show raw values (blue), moving average (black), and trend (red) lines. Slopes were estimated using a modified Mann-Kendall test and plotted such that the line passes through the median value in year fifteen of the study. Counts of total individuals were standardized by number of trapping nights. Significance of slopes is indicated at the 0.05 (*), 0.01 (**), and 0.001 (***) levels. The number of time points included in each moving average (10, 5, 15, and 11, respectively) was selected using BIC. (Color figure online) (Mann-Kendall test; alpha = 0.05) and 6% a positive trend (n = 19; Table 2). The 10-year rate of change for the declining species varied from − 14% to − 111% (median = − 50%) relative to the median annual count, and for the increasing species the rate of change varied from 7% to 155% (median = 44%). There was some evidence that mean rates of population change differed among species with different ranges (df = 2, F = 2.46, p = 0.087; Fig.   S4). Southern species were somewhat less likely to be decreasing and widespread species somewhat more likely to be decreasing than expected by chance (Chi-squared test, df = 4, X 2 = 8.559, p value = 0.073; Table 2), based on standardized residuals.
We also classified species based on the habitat (forest, garden, meadow, or forest/meadow) of their preferred food plants, to test the hypothesis that land use change was contributing to population trends (Habel et al. 2019b). Trends in these species did not differ by habitat group (df = 4, F = 1.58, p = 0.178), nor did the proportion of species that were increasing or declining (Table 3; Chi-squared test, df = 8, X 2 = 6.415, p value = 0.601). Trend estimates for each species, along with range and habitat designations and total numbers of years and individuals detected, are shown in Supplementary Table S3.

Responses to environmental covariates
Model convergence and mixing were good, based on effective sample sizes and potential scale reduction factors. The two best-performing models for predicting species' rates of change included climatic covariates from winter (soil temperature and snow depth; Table 1) and summer (temperature and Fig. 3 Three moth species representing those with northern range limits in northern (top; 'widespread' species), central (middle), and southern (bottom) Norway. Plots (center column) show raw values (blue), moving average (black), and trend lines (red) based on captures from a site in southeastern Norway. Maps of detections (right) are from Artsdatabanken (artsd ataba nken. no). Slopes were estimated using a modified Mann-Kendall test and plotted such that the line passed through the median value in year fifteen of the study. Significance of the slope is indicated at the 0.05 (*) level. The number of time points included in each moving average (12, 2, and 12, respectively) was selected using BIC. (Color figure online) precipitation) of year t. Models including all covariates, or lag covariates from the previous year (t-1) were ranked lower. Many species responded negatively to increased winter soil temperature, summer temperature, and summer precipitation (Table 4). Responses to winter snow depth were mixed and less common. Winter soil temperature and summer precipitation explained the most variation in their respective models ( Supplementary Fig. S5). Responses to covariates did not differ among species with different ranges (95% CIs). Instead, responses were strongly correlated with phylogeny (95% CI: 0.93-0.99, where 1.0 indicates the strongest possible phylogenetic effect). Models explained a relatively small proportion of total differences in moth abundance between years, as indicated by the mean R 2 values (Table 1).

Discussion
We analyzed a 30-year moth capture dataset from southeastern Norway and found that the richness and diversity of the community have declined considerably. Abundance of species with the most northerly range limits declined as well. Nearly one-fifth of species declined, whereas only 6% increased.
Inter-annual variability in community metrics and individual species was high, even when trends were strong. Southern species appear less likely to be in decline, and widespread species more likely, although statistical support for this result was weak. This is compatible with the hypothesis that Scandinavia, which is becoming warmer and wetter (Ljungqvist et al. 2019), is becoming more hospitable to southerly species (and possibly less hospitable to more northerly species). However, opposite to a general expectation that warmer seasons could have an enriching effect on local richness and abundance of poikilothermal species like moths, rates of population change for many species in our study were negatively correlated with summer temperature and winter soil temperature. Many species also responded negatively to wetter summers. Rates of decline did not appear to differ among species whose preferred food plants grew in different types of habitat, despite the conversion of many meadows and farmlands in the surrounding area to residential or commercial property over the sampling period. Our declines in moth richness and abundance mirror many of those found around the globe, where insect biomass and species richness appear to be falling in many ecosystems (Sánchez-Bayo and Wyckhuys 2019, Didham et al. 2020, Wagner 2020. Some communities, however, don't follow this pattern, highlighting that biodiversity loss is a complex phenomenon (Antão et al. 2020). Moths in subarctic Finland, for example, are as a whole stable or increasing (Hunter et al. 2014, Antão et al. 2020, and biomass of British moths may be stable (Macgregor et al. 2019b(Macgregor et al. , 2021 although over one third of species are declining . Understanding the processes that cause declines in some groups but not in others remains a challenge. Traits can be helpful in some cases (Habel et al. 2019a), but (as in our study) they often do not explain much variation in species responses (Mattila et al. 2006). This could mean that stochastic processes play a role, but given our relative ignorance about the life histories, diets, and indirect effects of climate on many species it seems likely that additional study could reveal important traits. High phylogenetic signal in species' responses to climate covariates in our JSDMs implies that unmeasured but phylogenetically structured traits do play a role (Abrego et al. 2017). Two potential causes of the declines we report-climate and land use changes-have occurred simultaneously and are difficult to unravel (Mair et al. 2014, Platts et al. 2019. Models based on annual variation in weather had some predictive value in our system, but rates of decline varied only slightly among species with different range limits (and, presumably, different climatic optima), although southern species appear less likely to be declining.
Direct effects of climate change may not be driving these community changes in straightforward ways, but climate may be important for other reasons. Changes in phenology can cause mismatches between lepidoptera and their food sources, for example, and more frequent severe weather events (including droughts) could increase strain on populations (Van Dyck et al. 2015, Boggs 2016. Climate change may cause ranges to contract for species not successful at climate tracking, as has been demonstrated for several other insect groups, including bumblebees (Soroye et al. 2020). A long-term study of Hungarian macro-moths also found that southern species were no more likely to appear in new sites than were northern species (Valtonen et al. 2017).
Forest cover has remained constant in southeastern Norway while much farm and field habitat has been lost to housing and other development, but these land use changes were poor predictors of which species were in decline. This contrasts with studies of moths in Germany, where grassland-dependent species appear to have some of the highest rates of decline (Habel et al. 2019b, c). Land use change might be less influential for moths in Fennoscandia or southeast Norway specifically, or our limited knowledge of host plant specificity for many species may reduce our ability to detect such effects. For instance, moth host plants related to field edges might still be present along roadsides in developed areas. Use of pesticides and herbicides in Norway appears to have remained relatively stable since the mid 1990's (Stenrød 2015). That many sub-Arctic moth species in Finland appear to be stable, despite the speed at which high-latitude climates are warming, may imply that these far northern habitats that have been subject to less intensive land use can act as a buffer to stabilize populations (Hunter et al. 2014, Zellweger et al. 2020. Regardless of the causes, the negative trends in moth populations that we report can impact entire ecosystems; for example, moths are important for pollination (Macgregor et al. 2015), and as food for insectivorous birds, many of which are also in sharp decline across Europe (Bowler et al. 2019), and for bats. Population trends for moths could also serve as an early-warning sign for other insects (Hunter et al. 2014).
One limitation of our study is that our capture records are from a single site (Hallmann et al. 2017), which limits our ability to generalize our results; rather, they should be used to form hypotheses for larger-scale monitoring projects. However, our time series is unbroken, consistent methods were used throughout, and the immediate surroundings of the trap site were unchanged. Records from long temporal and large spatial scales are rare, yet they are vital as we try to characterize insect population trends, understand their causes, and prevent precipitous declines (Montgomery et al. 2020. Increased long-term insect monitoring is needed. Funders should make provisions for these long-term studies that have traditionally been difficult to sustain, and study designs should allow the roles of alternative causal mechanisms to be disentangled.