The relative importance of local and regional processes to metapopulation dynamics

Metapopulation dynamics - patch occupancy, colonization and extinction - are the result of complex processes at both local (e.g. environmental conditions) and regional (e.g. spatial arrangement of habitat patches) scales. A large body of work has focused on habitat patch area and connectivity (area-isolation paradigm). However, these approaches often do not incorporate local environmental conditions or fully address how the spatial arrangement of habitat patches (and resulting connectivity) can influence metapopulation dynamics. Here, we utilize long-term data on a classic metapopulation system - the Glanville fritillary butterfly occupying a set of dry meadows and pastures in the angstrom land islands - to investigate the relative roles of local environmental conditions, geographic space and connectivity in capturing patch occupancy, colonization and extinction. We defined connectivity using traditional measures as well as graph-theoretic measures of centrality. Using boosted regression tree models, we find roughly comparable model performance among models trained on environmental conditions, geographic space or patch centrality. In models containing all of the covariates, we find strong and consistent evidence for the roles of resource abundance, longitude and centrality (i.e. connectivity) in predicting habitat patch occupancy and colonization, while patch centrality (connectivity) was relatively unimportant for predicting extinction. Relative variable importance did not change when geographic coordinates were not considered and models underwent spatially stratified cross-validation. Together, this suggests that the combination of regional-scale connectivity measures and local-scale environmental conditions is important for predicting metapopulation dynamics and that a stronger integration of ideas from network theory may provide insight into metapopulation processes.Peer reviewe

Combining range and phenology shifts offers a winning strategy for boreal Lepidoptera

Species can adapt to climate change by adjusting in situ or by dispersing to new areas, and these strategies may complement or enhance each other. Here, we investigate temporal shifts in phenology and spatial shifts in northern range boundaries for 289 Lepidoptera species by using long-term data sampled over two decades. While 40% of the species neither advanced phenology nor moved northward, nearly half (45%) used one of the two strategies. The strongest positive population trends were observed for the minority of species (15%) that both advanced flight phenology and shifted their northern range boundaries northward. We show that, for boreal Lepidoptera, a combination of phenology and range shifts is the most viable strategy under a changing climate. Effectively, this may divide species into winners and losers based on their propensity to capitalize on this combination, with potentially large consequences on future community composition.Peer reviewe

A regionally informed abundance index for supporting integrative analyses across butterfly monitoring schemes

1. The rapid expansion of systematic monitoring schemes necessitates robust methods to reliably assess species' status and trends. Insect monitoring poses a challenge where there are strong seasonal patterns, requiring repeated counts to reliably assess abundance. Butterfly monitoring schemes (BMSs) operate in an increasing number of countries with broadly the same methodology, yet they differ in their observation frequency and in the methods used to compute annual abundance indices. 2. Using simulated and observed data, we performed an extensive comparison of two approaches used to derive abundance indices from count data collected via BMS, under a range of sampling frequencies. Linear interpolation is most commonly used to estimate abundance indices from seasonal count series. A second method, hereafter the regional generalized additive model (GAM), fits a GAM to repeated counts within sites across a climatic region. For the two methods, we estimated bias in abundance indices and the statistical power for detecting trends, given different proportions of missing counts. We also compared the accuracy of trend estimates using systematically degraded observed counts of the Gatekeeper Pyronia tithonus (Linnaeus 1767). 3. The regional GAM method generally outperforms the linear interpolation method. When the proportion of missing counts increased beyond 50%, indices derived via the linear interpolation method showed substantially higher estimation error as well as clear biases, in comparison to the regional GAM method. The regional GAM method also showed higher power to detect trends when the proportion of missing counts was substantial. 4. Synthesis and applications. Monitoring offers invaluable data to support conservation policy and management, but requires robust analysis approaches and guidance for new and expanding schemes. Based on our findings, we recommend the regional generalized additive model approach when conducting integrative analyses across schemes, or when analysing scheme data with reduced sampling efforts. This method enables existing schemes to be expanded or new schemes to be developed with reduced within-year sampling frequency, as well as affording options to adapt protocols to more efficiently assess species status and trends across large geographical scales

Widespread forest vertebrate extinctions induced by a mega hydroelectric dam in lowland Amazonia

Mega hydropower projects in tropical forests pose a major emergent threat to terrestrial and freshwater biodiversity worldwide. Despite the unprecedented number of existing, underconstruction and planned hydroelectric dams in lowland tropical forests, long-term effects on biodiversity have yet to be evaluated. We examine how medium and large-bodied assemblages of terrestrial and arboreal vertebrates (including 35 mammal, bird and tortoise species) responded to the drastic 26-year post-isolation history of archipelagic alteration in landscape structure and habitat quality in a major hydroelectric reservoir of Central Amazonia. The Balbina Hydroelectric Dam inundated 3,129 km2 of primary forests, simultaneously isolating 3,546 land-bridge islands. We conducted intensive biodiversity surveys at 37 of those islands and three adjacent continuous forests using a combination of four survey techniques, and detected strong forest habitat area effects in explaining patterns of vertebrate extinction. Beyond clear area effects, edge-mediated surface fire disturbance was the most important additional driver of species loss, particularly in islands smaller than 10 ha. Based on species-area models, we predict that only 0.7% of all islands now harbor a species-rich vertebrate assemblage consisting of ≥80% of all species. We highlight the colossal erosion in vertebrate diversity driven by a man-made dam and show that the biodiversity impacts of mega dams in lowland tropical forest regions have been severely overlooked. The geopolitical strategy to deploy many more large hydropower infrastructure projects in regions like lowland Amazonia should be urgently reassessed, and we strongly advise that long-term biodiversity impacts should be explicitly included in pre-approval environmental impact assessments

Scaling up the effects of inbreeding depression from individuals to metapopulations

Abstract Inbreeding is common in nature, and many laboratory studies have documented that inbreeding depression can reduce the fitness of individuals. Demonstrating the consequences of inbreeding depression on the growth and persistence of populations is more challenging because populations are often regulated by density- or frequency-dependent selection and influenced by demographic and environmental stochasticity. A few empirical studies have shown that inbreeding depression can increase extinction risk of local populations. The importance of inbreeding depression at the metapopulation level has been conjectured based on population-level studies but has not been evaluated. We quantified the impact of inbreeding depression affecting the fitness of individuals on metapopulation persistence in heterogeneous habitat networks of different sizes and habitat configuration in a context of natural butterfly metapopulations. We developed a spatial individual-based simulation model of metapopulations with explicit genetics. We used Approximate Bayesian Computation to fit the model to extensive demographic, genetic, and life-history data available for the well-studied Glanville fritillary butterfly (Melitaea cinxia) metapopulations in the Åland islands in SW Finland. We compared 18 semi-independent habitat networks differing in size and fragmentation. The results show that inbreeding is more frequent in small habitat networks, and consequently, inbreeding depression elevates extinction risks in small metapopulations. Metapopulation persistence and neutral genetic diversity maintained in the metapopulations increase with the total habitat amount in and mean patch size of habitat networks. Dispersal and mating behavior interact with landscape structure to determine how likely it is to encounter kin while looking for mates. Inbreeding depression can decrease the viability of small metapopulations even when they are strongly influenced by stochastic extinction-colonization dynamics and density-dependent selection. The findings from this study support that genetic factors, in addition to demographic factors, can contribute to extinctions of small local populations and also of metapopulations. This article is protected by copyright. All rights reserved.Peer reviewe

Oviposition Cues for a Specialist Butterfly–Plant Chemistry and Size

The oviposition choice of an insect herbivore is based on a complex set of stimuli and responses. In this study, we examined the effect of plant secondary chemistry (the iridoid glycosides aucubin and catalpol) and aspects of size of the plant Plantago lanceolata, on the oviposition behavior of the specialist butterfly Melitaea cinxia. Iridoid glycosides are known to deter feeding or decrease the growth rate of generalist insect herbivores, but can act as oviposition cues and feeding stimulants for specialized herbivores. In a previous observational study of M. cinxia in the field, oviposition was associated with high levels of aucubin. However, this association could have been the cause (butterfly choice) or consequence (plant induction) of oviposition. We conducted a set of dual- and multiple-choice experiments in cages and in the field. In the cages, we found a positive association between the pre-oviposition level of aucubin and the number of ovipositions. The association reflects the butterfly oviposition selection rather than plant induction that follows oviposition. Our results also suggest a threshold concentration below which females do not distinguish between levels of iridoid glycosides. In the field, the size of the plant appeared to be a more important stimulus than iridoid glycoside content, with bigger plants receiving more oviposition than smaller plants, regardless of their secondary chemistry. Our results illustrate that the rank of a cue used for oviposition may be dependent on environmental context

More than 75 percent decline over 27 years in total flying insect biomass in protected areas

Global declines in insects have sparked wide interest among scientists, politicians, and the general public. Loss of insect diversity and abundance is expected to provoke cascading effects on food webs and to jeopardize ecosystem services. Our understanding of the extent and underlying causes of this decline is based on the abundance of single species or taxonomic groups only, rather than changes in insect biomass which is more relevant for ecological functioning. Here, we used a standardized protocol to measure total insect biomass using Malaise traps, deployed over 27 years in 63 nature protection areas in Germany (96 unique location-year combinations) to infer on the status and trend of local entomofauna. Our analysis estimates a seasonal decline of 76%, and mid-summer decline of 82% in flying insect biomass over the 27 years of study. We show that this decline is apparent regardless of habitat type, while changes in weather, land use, and habitat characteristics cannot explain this overall decline. This yet unrecognized loss of insect biomass must be taken into account in evaluating declines in abundance of species depending on insects as a food source, and ecosystem functioning in the European landscape

Bioclimatic context of species' populations determines community stability

Aim: It is important to understand the factors affecting community stability because ecosystem function is increasingly at risk from biodiversity loss. Here, we evaluate how a key factor, the position of local environmental conditions within the thermal range of the species, influences the stability of butterfly communities at a continental scale. Location: Spain, UK and Finland. Time period: 1999–2017. Major taxa studied: Butterflies. Methods: We tested the following hypotheses about how species responses to temperature anomalies aggregate to influence stability: Hypothesis 1, species have contrasting responses to local temperature anomalies at opposing edges of their thermal range; hypothesis 2, communities with central thermal range positions have higher community stability; and the impacts of thermal range position on community stability are driven by hypothesis 3, population asynchrony, or hypothesis 4, additive population stability. Data were analysed at 876 sites for 157 species. Results: We found some support for hypothesis 1, because there were interactions between thermal range and response to temperature anomalies such that species at different range edges could provide weak compensatory dynamics. However, responses were nonlinear, suggesting strong declines with extreme anomalies, particularly at the hot range edge. Hypothesis 2 was supported in part, because community stability increased with central thermal range positions and declined at the edges, after accounting for species richness and community abundance. Thermal range position was weakly correlated with asynchrony (hypothesis 3) and population stability (hypothesis 4), although species richness and population abundance had larger impacts. Main conclusions: Future extreme heat events will be likely to impact species negatively across their thermal range, but might be particularly impactful on populations at the hottest end of the thermal range. Thermal range position influenced community stability because range edge communities were stable. However, the prediction of community stability from thermal range position is challenging because of nonlinear responses to temperature, with small temperature anomalies producing weak compensatory dynamics, but large extreme events synchronizing dynamics