Urbanization extends flight phenology and leads to local adaptation of seasonal plasticity in Lepidoptera.

Urbanization is gaining force globally, which challenges biodiversity, and it has recently also emerged as an agent of evolutionary change. Seasonal phenology and life cycle regulation are essential processes that urbanization is likely to alter through both the urban heat island effect (UHI) and artificial light at night (ALAN). However, how UHI and ALAN affect the evolution of seasonal adaptations has received little attention. Here, we test for the urban evolution of seasonal life-history plasticity, specifically changes in the photoperiodic induction of diapause in two lepidopterans, Pieris napi (Pieridae) and Chiasmia clathrata (Geometridae). We used long-term data from standardized monitoring and citizen science observation schemes to compare yearly phenological flight curves in six cities in Finland and Sweden to those of adjacent rural populations. This analysis showed for both species that flight seasons are longer and end later in most cities, suggesting a difference in the timing of diapause induction. Then, we used common garden experiments to test whether the evolution of the photoperiodic reaction norm for diapause could explain these phenological changes for a subset of these cities. These experiments demonstrated a genetic shift for both species in urban areas toward a lower daylength threshold for direct development, consistent with predictions based on the UHI but not ALAN. The correspondence of this genetic change to the results of our larger-scale observational analysis of in situ flight phenology indicates that it may be widespread. These findings suggest that seasonal life cycle regulation evolves in urban ectotherms and may contribute to ecoevolutionary dynamics in cities.

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.

Spatial synchrony is related to environmental change in Finnish moth communities.

Spatially distinct pairs of sites may have similarly fluctuating population dynamics across large geographical distances, a phenomenon called spatial synchrony. However, species rarely exist in isolation, but rather as members of interactive communities, linked with other communities through dispersal (i.e. a metacommunity). Using data on Finnish moth communities sampled across 65 sites for 20 years, we examine the complex synchronous/anti-synchronous relationships among sites using the geography of synchrony framework. We relate site-level synchrony to mean and temporal variation in climatic data, finding that colder and drier sites-and those with the most drastic temperature increases-are important for spatial synchrony. This suggests that faster-warming sites contribute most strongly to site-level estimates of synchrony, highlighting the role of a changing climate to spatial synchrony. Considering the spatial variability in climate change rates is therefore important to understand metacommunity dynamics and identify habitats which contribute most strongly to spatial synchrony.

Landscape simplification weakens the association between terrestrial producer and consumer diversity in Europe.

Land-use change is one of the primary drivers of species loss, yet little is known about its effect on other components of biodiversity that may be at risk. Here, we ask whether, and to what extent, landscape simplification, measured as the percentage of arable land in the landscape, disrupts the functional and phylogenetic association between primary producers and consumers. Across seven European regions, we inferred the potential associations (functional and phylogenetic) between host plants and butterflies in 561 seminatural grasslands. Local plant diversity showed a strong bottom-up effect on butterfly diversity in the most complex landscapes, but this effect disappeared in simple landscapes. The functional associations between plant and butterflies are, therefore, the results of processes that act not only locally but are also dependent on the surrounding landscape context. Similarly, landscape simplification reduced the phylogenetic congruence among host plants and butterflies indicating that closely related butterflies become more generalist in the resources used. These processes occurred without any detectable change in species richness of plants or butterflies along the gradient of arable land. The structural properties of ecosystems are experiencing substantial erosion, with potentially pervasive effects on ecosystem functions and future evolutionary trajectories. Loss of interacting species might trigger cascading extinction events and reduce the stability of trophic interactions, as well as influence the longer term resilience of ecosystem functions. This underscores a growing realization that species richness is a crude and insensitive metric and that both functional and phylogenetic associations, measured across multiple trophic levels, are likely to provide additional and deeper insights into the resilience of ecosystems and the functions they provide.

Density of insect-pollinated grassland plants decreases with increasing surrounding land-use intensity.

Pollinator declines have raised concerns about the persistence of plant species that depend on insect pollination, in particular by bees, for their reproduction. The impact of pollinator declines remains unknown for species-rich plant communities found in temperate seminatural grasslands. We investigated effects of land-use intensity in the surrounding landscape on the distribution of plant traits related to insect pollination in 239 European seminatural grasslands. Increasing arable land use in the surrounding landscape consistently reduced the density of plants depending on bee and insect pollination. Similarly, the relative abundance of bee-pollination-dependent plants increased with higher proportions of non-arable agricultural land (e.g. permanent grassland). This was paralleled by an overall increase in bee abundance and diversity. By isolating the impact of the surrounding landscape from effects of local habitat quality, we show for the first time that grassland plants dependent on insect pollination are particularly susceptible to increasing land-use intensity in the landscape.

Is climate warming more consequential towards poles? The phenology of Lepidoptera in Finland.

The magnitude and direction of phenological shifts from climate warming could be predictably variable across the planet depending upon the nature of physiological controls on phenology, the thermal sensitivity of the developmental processes and global patterns in the climate warming. We tested this with respect to the flight phenology of adult nocturnal moths (3.33 million captures of 334 species) that were sampled at sites in southern and northern Finland during 1993-2012 (with years 2005-2012 treated as an independent model validation data set). We compared eight competing models of physiological controls on flight phenology to each species and found strong support for thermal controls of phenology in 66% of the species generations. Among species with strong thermal control of phenology in both the south and north, the average development rate was higher in northern vs. southern populations at 10 °C, but about the same at 15 and 20 °C. With a 3 °C increase in temperature (approximating A2 scenario of IPPC for 2090-2099 relative to 1980-1999) these species were predicted to advance their phenology on average by 17 (SE ± 0.3) days in the south vs. 13 (±0.4) days in the north. The higher development rates at low temperatures of poleward populations makes them less sensitive to climate warming, which opposes the tendency for stronger phenological advances in the north from greater increases in temperature.

Recent range shifts of moths, butterflies, and birds are driven by the breadth of their climatic niche.

Species are altering their ranges as a response to climate change, but the magnitude and direction of observed range shifts vary considerably among species. The ability to persist in current areas and colonize new areas plays a crucial role in determining which species will thrive and which decline as climate change progresses. Several studies have sought to identify characteristics, such as morphological and life-history traits, that could explain differences in the capability of species to shift their ranges together with a changing climate. These characteristics have explained variation in range shifts only sporadically, thus offering an uncertain tool for discerning responses among species. As long-term selection to past climates have shaped species' tolerances, metrics describing species' contemporary climatic niches may provide an alternative means for understanding responses to on-going climate change. Species that occur in a broader range of climatic conditions may hold greater tolerance to climatic variability and could therefore more readily maintain their historical ranges, while species with more narrow tolerances may only persist if they are able to shift in space to track their climatic niche. Here, we provide a first-filter test of the effect of climatic niche dimensions on shifts in the leading range edges in three relatively well-dispersing species groups. Based on the realized changes in the northern range edges of 383 moth, butterfly, and bird species across a boreal 1,100 km latitudinal gradient over c. 20 years, we show that while most morphological or life-history traits were not strongly connected with range shifts, moths and birds occupying a narrower thermal niche and butterflies occupying a broader moisture niche across their European distribution show stronger shifts towards the north. Our results indicate that the climatic niche may be important for predicting responses under climate change and as such warrants further investigation of potential mechanistic underpinnings.

Century-long butterfly range expansions in northern Europe depend on climate, land use and species traits.

Climate change is an important driver of range shifts and community composition changes. Still, little is known about how the responses are influenced by the combination of land use, species interactions and species traits. We integrate climate and distributional data for 131 butterfly species in Sweden and Finland and show that cumulative species richness has increased with increasing temperature over the past 120 years. Average provincial species richness increased by 64% (range 15-229%), from 46 to 70. The rate and direction of range expansions have not matched the temperature changes, in part because colonisations have been modified by other climatic variables, land use and vary according to species characteristics representing ecological generalisation and species interactions. Results emphasise the role of a broad ecological filtering, whereby a mismatch between environmental conditions and species preferences limit the ability to disperse and establish populations in emerging climates and novel areas, with potentially widespread implications for ecosystem functioning.

Environmental controls on the phenology of moths: predicting plasticity and constraint under climate change.

Ecological systems have naturally high interannual variance in phenology. Component species have presumably evolved to maintain appropriate phenologies under historical climates, but cases of inappropriate phenology can be expected with climate change. Understanding controls on phenology permits predictions of ecological responses to climate change. We studied phenological control systems in Lepidoptera by analyzing flight times recorded at a network of sites in Finland. We evaluated the strength and form of controls from temperature and photoperiod, and tested for geographic variation within species. Temperature controls on phenology were evident in 51% of 112 study species and for a third of those thermal controls appear to be modified by photoperiodic cues. For 24% of the total, photoperiod by itself emerged as the most likely control system. Species with thermal control alone should be most immediately responsive in phenology to climate warming, but variably so depending upon the minimum temperature at which appreciable development occurs and the thermal responsiveness of development rate. Photoperiodic modification of thermal controls constrains phenotypic responses in phenologies to climate change, but can evolve to permit local adaptation. Our results suggest that climate change will alter the phenological structure of the Finnish Lepidoptera community in ways that are predictable with knowledge of the proximate physiological controls. Understanding how phenological controls in Lepidoptera compare to that of their host plants and enemies could permit general inferences regarding climatic effects on mid- to high-latitude ecosystems.

Life-history traits predict species responses to habitat area and isolation: a cross-continental synthesis.

There is a lack of quantitative syntheses of fragmentation effects across species and biogeographic regions, especially with respect to species life-history traits. We used data from 24 independent studies of butterflies and moths from a wide range of habitats and landscapes in Europe and North America to test whether traits associated with dispersal capacity, niche breadth and reproductive rate modify the effect of habitat fragmentation on species richness. Overall, species richness increased with habitat patch area and connectivity. Life-history traits improved the explanatory power of the statistical models considerably and modified the butterfly species-area relationship. Species with low mobility, a narrow feeding niche and low reproduction were most strongly affected by habitat loss. This demonstrates the importance of considering life-history traits in fragmentation studies and implies that both species richness and composition change in a predictable manner with habitat loss and fragmentation.

Habitat fragmentation causes immediate and time-delayed biodiversity loss at different trophic levels.

Intensification or abandonment of agricultural land use has led to a severe decline of semi-natural habitats across Europe. This can cause immediate loss of species but also time-delayed extinctions, known as the extinction debt. In a pan-European study of 147 fragmented grassland remnants, we found differences in the extinction debt of species from different trophic levels. Present-day species richness of long-lived vascular plant specialists was better explained by past than current landscape patterns, indicating an extinction debt. In contrast, short-lived butterfly specialists showed no evidence for an extinction debt at a time scale of c. 40 years. Our results indicate that management strategies maintaining the status quo of fragmented habitats are insufficient, as time-delayed extinctions and associated co-extinctions will lead to further biodiversity loss in the future.

Dispersal capacity and diet breadth modify the response of wild bees to habitat loss.

Habitat loss poses a major threat to biodiversity, and species-specific extinction risks are inextricably linked to life-history characteristics. This relationship is still poorly documented for many functionally important taxa, and at larger continental scales. With data from five replicated field studies from three countries, we examined how species richness of wild bees varies with habitat patch size. We hypothesized that the form of this relationship is affected by body size, degree of host plant specialization and sociality. Across all species, we found a positive species-area slope (z = 0.19), and species traits modified this relationship. Large-bodied generalists had a lower z value than small generalists. Contrary to predictions, small specialists had similar or slightly lower z value compared with large specialists, and small generalists also tended to be more strongly affected by habitat loss as compared with small specialists. Social bees were negatively affected by habitat loss (z = 0.11) irrespective of body size. We conclude that habitat loss leads to clear shifts in the species composition of wild bee communities.

Meta-analysis of multidecadal biodiversity trends in Europe.

Local biodiversity trends over time are likely to be decoupled from global trends, as local processes may compensate or counteract global change. We analyze 161 long-term biological time series (15-91 years) collected across Europe, using a comprehensive dataset comprising ~6,200 marine, freshwater and terrestrial taxa. We test whether (i) local long-term biodiversity trends are consistent among biogeoregions, realms and taxonomic groups, and (ii) changes in biodiversity correlate with regional climate and local conditions. Our results reveal that local trends of abundance, richness and diversity differ among biogeoregions, realms and taxonomic groups, demonstrating that biodiversity changes at local scale are often complex and cannot be easily generalized. However, we find increases in richness and abundance with increasing temperature and naturalness as well as a clear spatial pattern in changes in community composition (i.e. temporal taxonomic turnover) in most biogeoregions of Northern and Eastern Europe.

Relative contributions of local and regional factors to species richness and total density of butterflies and moths in semi-natural grasslands.

Metapopulation theory predicts that species richness and total population density of habitat specialists increase with increasing area and regional connectivity of the habitat. To test these predictions, we examined the relative contributions of habitat patch area, connectivity of the regional habitat network and local habitat quality to species richness and total density of butterflies and day-active moths inhabiting semi-natural grasslands. We studied butterflies and moths in 48 replicate landscapes situated in southwest Finland, including a focal patch and the surrounding network of other semi-natural grasslands within a radius of 1.5 km from the focal patch. By applying the method of hierarchical partitioning, which can distinguish between independent and joint contributions of individual explanatory variables, we observed that variables of the local habitat quality (e.g. mean vegetation height and nectar plant abundance) generally showed the highest independent effect on species richness and total density of butterflies and moths. Habitat area did not show a significant independent contribution to species richness and total density of butterflies and moths. The effect of habitat connectivity was observed only for total density of the declining butterflies and moths. These observations indicate that the local habitat quality is of foremost importance in explaining variation in species richness and total density of butterflies and moths. In addition, declining butterflies and moths have larger populations in well-connected networks of semi-natural grasslands. Our results suggest that, while it is crucial to maintain high-quality habitats by management, with limited resources it would be appropriate to concentrate grassland management and restoration to areas with well-connected grassland networks in which the declining species currently have their strongest populations.

New insights into butterfly-environment relationships using partitioning methods.

Variation partitioning and hierarchical partitioning are novel statistical approaches that provide deeper understanding of the importance of different explanatory variables for biodiversity patterns than traditional regression methods. Using these methods, the variation in occupancy and abundance of the clouded apollo butterfly (Parnassius mnemosyne L.) was decomposed into independent and joint effects of larval and adult food resources, microclimate and habitat quantity. The independent effect of habitat quantity variables (habitat area and connectivity) captured the largest fraction of the variation in the clouded apollo patterns, but habitat connectivity had a major contribution only for occupancy data. The independent effects of resources and microclimate were higher on butterfly abundance than on occupancy. However, a considerable amount of variation in the butterfly patterns was accounted for by the joint effects of predictors and may thus be causally related to two or all three groups of variables. Abundance of the butterfly in the surroundings of the focal grid cell had a significant effect in all analyses, independently of the effects of other predictors. Our results encourage wider applications of partitioning methods in biodiversity studies.

Protected areas alleviate climate change effects on northern bird species of conservation concern.

Global climate change is a major threat to biodiversity, posing increasing pressures on species to adapt in situ or shift their ranges. A protected area network is one of the main instruments to alleviate the negative impacts of climate change. Importantly, protected area networks might be expected to enhance the resilience of regional populations of species of conservation concern, resulting in slower species loss in landscapes with a significant amount of protected habitat compared to unprotected landscapes. Based on national bird atlases compiled in 1974-1989 and 2006-2010, this study examines the recent range shifts in 90 forest, mire, marshland, and Arctic mountain heath bird species of conservation concern in Finland, as well as the changes in their species richness in protected versus unprotected areas. The trends emerging from the atlas data comparisons were also related to the earlier study dealing with predictions of distributional changes for these species for the time slice of 2051-2080, developed using bioclimatic envelope models (BEMs). Our results suggest that the observed changes in bird distributions are in the same direction as the BEM-based predictions, resulting in a decrease in species richness of mire and Arctic mountain heath species and an increase in marshland species. The patterns of changes in species richness between the two time slices are in general parallel in protected and unprotected areas. However, importantly, protected areas maintained a higher level of species richness than unprotected areas. This finding provides support for the significance and resilience provision of protected area networks in preserving species of conservation concern under climate change.

Higher mobility of butterflies than moths connected to habitat suitability and body size in a release experiment.

Mobility is a key factor determining lepidopteran species responses to environmental change. However, direct multispecies comparisons of mobility are rare and empirical comparisons between butterflies and moths have not been previously conducted. Here, we compared mobility between butterflies and diurnal moths and studied species traits affecting butterfly mobility. We experimentally marked and released 2011 butterfly and 2367 moth individuals belonging to 32 and 28 species, respectively, in a 25 m × 25 m release area within an 11-ha, 8-year-old set-aside field. Distance moved and emigration rate from the release habitat were recorded by species. The release experiment produced directly comparable mobility data in 18 butterfly and 9 moth species with almost 500 individuals recaptured. Butterflies were found more mobile than geometroid moths in terms of both distance moved (mean 315 m vs. 63 m, respectively) and emigration rate (mean 54% vs. 17%, respectively). Release habitat suitability had a strong effect on emigration rate and distance moved, because butterflies tended to leave the set-aside, if it was not suitable for breeding. In addition, emigration rate and distance moved increased significantly with increasing body size. When phylogenetic relatedness among species was included in the analyses, the significant effect of body size disappeared, but habitat suitability remained significant for distance moved. The higher mobility of butterflies than geometroid moths can largely be explained by morphological differences, as butterflies are more robust fliers. The important role of release habitat suitability in butterfly mobility was expected, but seems not to have been empirically documented before. The observed positive correlation between butterfly size and mobility is in agreement with our previous findings on butterfly colonization speed in a long-term set-aside experiment and recent meta-analyses on butterfly mobility.

Impacts of land cover data selection and trait parameterisation on dynamic modelling of species' range expansion.

Dynamic models for range expansion provide a promising tool for assessing species' capacity to respond to climate change by shifting their ranges to new areas. However, these models include a number of uncertainties which may affect how successfully they can be applied to climate change oriented conservation planning. We used RangeShifter, a novel dynamic and individual-based modelling platform, to study two potential sources of such uncertainties: the selection of land cover data and the parameterization of key life-history traits. As an example, we modelled the range expansion dynamics of two butterfly species, one habitat specialist (Maniola jurtina) and one generalist (Issoria lathonia). Our results show that projections of total population size, number of occupied grid cells and the mean maximal latitudinal range shift were all clearly dependent on the choice made between using CORINE land cover data vs. using more detailed grassland data from three alternative national databases. Range expansion was also sensitive to the parameterization of the four considered life-history traits (magnitude and probability of long-distance dispersal events, population growth rate and carrying capacity), with carrying capacity and magnitude of long-distance dispersal showing the strongest effect. Our results highlight the sensitivity of dynamic species population models to the selection of existing land cover data and to uncertainty in the model parameters and indicate that these need to be carefully evaluated before the models are applied to conservation planning.

Long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): survey methods, data management, and long-term population trends.

Long-term observational studies conducted at large (regional) spatial scales contribute to better understanding of landscape effects on population and evolutionary dynamics, including the conditions that affect long-term viability of species, but large-scale studies are expensive and logistically challenging to keep running for a long time. Here, we describe the long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia) that has been conducted since 1991 in a large network of 4000 habitat patches (dry meadows) within a study area of 50 by 70 km in the Åland Islands in Finland. We explain how the landscape structure has been described, including definition, delimitation, and mapping of the habitat patches; methods of field survey, including the logistics, cost, and reliability of the survey; and data management using the EarthCape biodiversity platform. We describe the long-term metapopulation dynamics of the Glanville fritillary based on the survey. There has been no long-term change in the overall size of the metapopulation, but the level of spatial synchrony and hence the amplitude of fluctuations in year-to-year metapopulation dynamics have increased over the years, possibly due to increasing frequency of exceptional weather conditions. We discuss the added value of large-scale and long-term population studies, but also emphasize the need to integrate more targeted experimental studies in the context of long-term observational studies. For instance, in the case of the Glanville fritillary project, the long-term study has produced an opportunity to sample individuals for experiments from local populations with a known demographic history. These studies have demonstrated striking differences in dispersal rate and other life-history traits of individuals from newly established local populations (the offspring of colonizers) versus individuals from old, established local populations. The long-term observational study has stimulated the development of metapopulation models and provided an opportunity to test model predictions. This combination of empirical studies and modeling has facilitated the study of key phenomena in spatial dynamics, such as extinction threshold and extinction debt.