Reduced variability in range‐edge butterfly populations over three decades of climate warming

Populations at the high latitude edge of species’ geographical ranges are thought to show larger interannual population fluctuations, with subsequent higher local extinction risk, than those within the ‘core’ climatic range. As climate envelopes shift northward under climate warming, however, we would expect populations to show dampened variability. We test this hypothesis using annual abundance indices from 19 butterfly species across 79 British monitoring sites between 1976 and 2009, a period of climatic warming. We found that populations in the latter (warmer) half of the recording period show reduced interannual population variability. Species with more southerly European distributions showed the greatest dampening in population variability over time. Our results suggest that increases in population variability occur towards climatic range boundaries. British sites, previously existing at the margins of suitable climate space, now appear to fall closer to the core climatic range for many butterfly species.     

Obligate annual and successive facultative diapause establish a bet‐hedging strategy of Rhagoletis cerasi (Diptera: Tephritidae) in seasonally unpredictable environments

To cope with temporal and spatial heterogeneity of habitats, herbivorous insects in the temperate zone usually enter diapause that facilitates synchronization of their life cycle with specific stages of host plants, such as fruit ripening. In the present study, we address those factors regulating dormancy responses as part of a ‘longer strategy’ to persist and thrive in temperate environments, focusing on Rhagoletis cerasi, a univoltine, oligophagous species, which overwinters as pupae and emerges when host fruits are available for oviposition at local scale. To ensure population survival and reproduction at habitats with ecological heterogeneity, R. cerasi has evolved a sophisticated diapause strategy based on a combination of local adaptation and diversified bet‐hedging strategies. Diapause duration is determined both by (i) the adaptive response to local host fruit phenology patterns (annual diapause) and (ii) the plastic responses to unpredictable inter‐annual (temporal) climatic variability that drives a proportion of the populations to extend dormancy by entering a second, successive, facultative cycle of prolonged diapause as part of a bet‐hedging strategy. Besides the dormant periods, post‐diapause development (which varies among populations) exerts ‘fine tune’ adjustments that assure synchronization and may correct possible errors. Adults emerging from pupae with prolonged diapause are larger in body size compared with counterparts emerging during the first year of diapause. However, female fecundity rates are reduced, followed by an extended post‐oviposition period, whereas adult longevity remains unaffected. Overall, it appears that R. cerasi populations are adapted to ecological conditions of local habitats and respond plastically to unpredictable environmental (climatic) conditions.     

The role of demography,intra‐species variation,and species distribution models in species' projections under climate change

Organisms are projected to shift their distribution ranges under climate change. The typical way to assess range shifts is by species distribution models (SDMs), which predict species’ responses to climate based solely on projected climatic suitability. However, life history traits can impact species’ responses to shifting habitat suitability. Additionally, it remains unclear if differences in vital rates across populations within a species can offset or exacerbate the effects of predicted changes in climatic suitability on population viability. In order to obtain a fuller understanding of the response of one species to projected climatic changes, we coupled demographic processes with predicted changes in suitable habitat for the monocarpic thistle Carlina vulgaris across northern Europe. We first developed a life history model with species‐specific average fecundity and survival rates and linked it to a SDM that predicted changes in habitat suitability through time with changes in climatic variables. We then varied the demographic parameters based upon observed vital rates of local populations from a translocation experiment. Despite the fact that the SDM alone predicted C. vulgaris to be a climate ‘winner’ overall, coupling the model with changes in demography and small‐scale habitat suitability resulted in a matrix of stable, declining, and increasing patches. For populations predicted to experience declines or increases in abundance due to changes in habitat suitability, altered fecundity and survival rates can reverse projected population trends.     

How universal are reserve design rules? A test using butterflies and their life history traits

We used butterfly species lists available for a set of 125 Czech Republic National Nature Reserves and Monuments, the highest small‐sized conservation category in the country encompassing practically all biotope types existing in central Europe, to test the validity of generally agreed “reserve design rules” using multivariate ordination analyses. Further, we used ordination analysis of butterfly life history traits to seek for biological mechanisms responsible for butterfly community responses to essentially geometric reserves characteristics. Reserve area, relative perimeter, within reserve habitat heterogeneity, and surrounding landscape compositional and configurational heterogeneity all affected the composition of butterfly assemblages after controlling for effects of geographical position and prevailing biotope type. Species inclining towards large reserves displayed low mobility and high local population density, probably because they require large habitat areas to maintain self‐sustaining populations; such species tend to have restricted distribution in the country and threatened status. Reserves with relatively long boundaries hosted species with high mobility, broad trophic range and long adult period; faunas of such reserves contain high proportions of widespread generalists. Species with narrow trophic ranges inclined towards reserves containing diverse habitats, probably due to requirements for high floristic diversity. Species with short adult flight, low generations number and overwintering in early stages inclined towards reserves situated amidst diverse landscapes, perhaps because such species require finely‐grained mosaics for metapopulation dynamics. Commonly agreed reserve design rules thus hold for Central European butterflies, but different design characteristics are important for individual species, depending on their life histories.     

Components of breeding performance in two competing species: habitat heterogeneity, individual quality and density-dependence

Density‐dependent breeding performance due to habitat heterogeneity has been shown to regulate populations of territorial species, since the progressive occupation of low quality territories as breeding density increases may cause a decline in the mean per capita fecundity of a population while variation in fecundity increases. Although the preemptive use of sites may relegate low quality individuals to sites of progressively lower suitability, few studies on density dependence have tried to separate the effects of territory quality from individual quality, and none have simultaneously considered the effects of heterospecific competitors. Using two long‐term monitored populations, we assessed the relative contribution of habitat heterogeneity and bird quality (in terms of age) on the productivity of sympatric golden Aquila chrysaetos and Bonelli's eagles Hieraaetus fasciatus under different scenarios of intra‐ and inter‐specific competition. Productivity (number of offspring fledged) varied among territories and average annual productivity was negatively related to its variability in both species and populations, thus giving some support to the habitat heterogeneity hypothesis. However, the effect of habitat heterogeneity on productivity became non‐significant when parental age and local density estimators were included in multivariate analyses. Therefore, temporal changes in bird quality (age) combined with intra‐ and interspecific competition explained variability in territory productivity rather than habitat heterogeneity among territories per se. The recruitment of subadult breeders, a surrogate of mortality in eagles, strongly varied among territories. Habitat heterogeneity in productivity may thus arise not because sites differ in suitability for reproduction but because of differences in factors affecting survival. Territories associated with high mortality risks have a higher probability of being occupied by young birds, whose lower quality, interacting with the density competitors, leads to a reduction of productivity. Site‐dependent variability in adult survival and interspecific competition may be extensive, but so far largely overlooked, factors to be seriously considered for the site‐dependent population regulation framework.     

Population resilience to an extreme drought is influenced by habitat area and fragmentation in the local landscape

Most studies on the biological impact of climate change have focussed on incremental climate warming, rather than extreme events. Yet responses of species’ populations to climatic extremes may be one of the primary drivers of ecological change. We assess the resilience of individual populations in terms of their sensitivity to‐ and ability to recover from‐ environmental perturbation. We demonstrate the method using a model species, the ringlet butterfly Aphantopus hyperantus, and analyse the effects of an extreme drought event using data from 79 British sites over 10 yr. We find that populations crashed most severely in drier regions but, additionally, the landscape structure around sites influenced population responses. Larger and more connected patches of woodland habitat reduced population sensitivity to the drought event and also facilitated faster recovery. Having enough, sufficiently connected habitat appears essential for species’ populations to be resilient to the increased climatic variability predicted under future scenarios.     

Climate and habitat configuration limit range expansion and patterns of dispersal in a non‐native lizard

Invasive species are one of the main causes of biodiversity loss worldwide. As introduced, populations increase in abundance and geographical range, so does the potential for negative impacts on native communities. As such, there is a need to better understand the processes driving range expansion as species become established in recipient landscapes. Through an investigation into capacity for population growth and range expansion of introduced populations of a non‐native lizard (Podarcis muralis), we aimed to demonstrate how multi‐scale factors influence spatial spread, population growth, and invasion potential in introduced species. We collated location records of P. muralis presence in England, UK through data collected from field surveys and a citizen science campaign. We used these data as input for presence‐background models to predict areas of climate suitability at a national‐scale (5 km resolution), and fine‐scale habitat suitability at the local scale (2 m resolution). We then integrated local models into an individual‐based modeling platform to simulate population dynamics and forecast range expansion for 10 populations in heterogeneous landscapes. National‐scale models indicated climate suitability has restricted the species to the southern parts of the UK, primarily by a latitudinal cline in overwintering conditions. Patterns of population growth and range expansion were related to differences in local landscape configuration and heterogeneity. Growth curves suggest populations could be in the early stages of exponential growth. However, annual rates of range expansion are predicted to be low (5–16 m). We conclude that extensive nationwide range expansion through secondary introduction is likely to be restricted by currently unsuitable climate beyond southern regions of the UK. However, exponential growth of local populations in habitats providing transport pathways is likely to increase opportunities for regional expansion. The broad habitat niche of P. muralis, coupled with configuration of habitat patches in the landscape, allows populations to increase locally with minimal dispersal.     

Contrasting growth forecasts across the geographical range of Scots pine due to altitudinal and latitudinal differences in climatic sensitivity

Ongoing changes in global climate are altering ecological conditions for many species. The consequences of such changes are typically most evident at the edge of a species’ geographical distribution, where differences in growth or population dynamics may result in range expansions or contractions. Understanding population responses to different climatic drivers along wide latitudinal and altitudinal gradients is necessary in order to gain a better understanding of plant responses to ongoing increases in global temperature and drought severity. We selected Scots pine (Pinus sylvestris L.) as a model species to explore growth responses to climatic variability (seasonal temperature and precipitation) over the last century through dendrochronological methods. We developed linear models based on age, climate and previous growth to forecast growth trends up to year 2100 using climatic predictions. Populations were located at the treeline across a latitudinal gradient covering the northern, central and southernmost populations and across an altitudinal gradient at the southern edge of the distribution (treeline, medium and lower elevations). Radial growth was maximal at medium altitude and treeline of the southernmost populations. Temperature was the main factor controlling growth variability along the gradients, although the timing and strength of climatic variables affecting growth shifted with latitude and altitude. Predictive models forecast a general increase in Scots pine growth at treeline across the latitudinal distribution, with southern populations increasing growth up to year 2050, when it stabilizes. The highest responsiveness appeared at central latitude, and moderate growth increase is projected at the northern limit. Contrastingly, the model forecasted growth declines at lowland‐southern populations, suggesting an upslope range displacement over the coming decades. Our results give insight into the geographical responses of tree species to climate change and demonstrate the importance of incorporating biogeographical variability into predictive models for an accurate prediction of species dynamics as climate changes.     

Butterfly species richness in the north-west Mediterranean Basin: the role of natural and human-induced factors

Aim We developed a model enabling us to evaluate the contribution of both natural and human‐related factors to butterfly species richness in Catalonia, a Mediterranean area that harbours one of the most diverse butterfly faunas in Europe. Location The study was carried out in Catalonia (north‐east Iberian Peninsula), a region of 31,930 km² lying between the Pyrenees, the Ebro depression and the Mediterranean sea. Methods Data from the Catalan Butterfly Monitoring Scheme were used to assess butterfly species richness from 55 transects spread all over the region. Three groups of environmental variables likely to affect the presence of butterfly species were calculated, above all from geographic information system data: (1) climatology and topography, (2) vegetation structure and (3) human disturbance. Because climatic and topographic variables are expected to be strongly correlated, we first performed a principal component analysis (PCA) to create a summarizing factor that would account for most of the variance within this set of variables. Subsequently, a backward stepwise multiple regression was performed in order to assess the effects of environmental factors on butterfly species richness. Results A total of 131 species were detected in the monitoring transects, representing 75.7% of the butterfly fauna known from Catalonia. Mean species richness per transect and per year was 41.4, although values varied greatly among sites (range: 14–76.8). The final regression model explained more than 80% of the total variance, which indicated a strong association between butterfly species richness and the studied environmental factors. The model revealed the very important contribution of climatic and topographic variables, which were combined into a single factor in the PCA. In contrast to what has been found in other, more northerly countries, species richness was negatively correlated with temperature and positively correlated with rainfall, except for extreme cold and wet conditions. This may be a consequence of the predictably adverse effects of the Mediterranean summer drought on herbivorous insects, and the fact that the limits of distribution of many butterflies correlate well with climatic variables. Human disturbance (defined as the proportion of urban and agricultural landscape cover in buffer areas of 5 km around the transect sites) was the second most important predictor for species richness. We found that a significant decrease in species numbers was associated with an increase in human pressure, a finding that indicates that not only building development, but also modern‐day agricultural practices are detrimental to the conservation of Mediterranean butterflies. Surprisingly, vegetation variables had an almost negligible effect on butterfly species richness. Main conclusions Our findings strongly indicate that the current motors of global change will have a negative effect on Mediterranean butterfly assemblages. First, changes in land‐use are transforming and fragmenting the landscape into an inhospitable and less permeable matrix for butterflies. Secondly, the negative correlation between species richness and temperature will lead to a predictable loss of diversity over the coming years, as predicted in the most plausible scenarios of climate change. Considering the high butterfly richness characterizing the Mediterranean Basin, this future trend will pose a serious threat to biodiversity.     

Genetic structure of a recent climate change‐driven range extension

The life‐history strategies of some species make them strong candidates for rapid exploitation of novel habitat under new climate regimes. Some early‐responding species may be considered invasive, and negatively impact on ‘naïve’ ecosystems. The barrens‐forming sea urchin Centrostephanus rodgersii is one such species, having a high dispersal capability and a high‐latitude range margin limited only by a developmental temperature threshold. Within this species’ range in eastern Australian waters, sea temperatures have increased at greater than double the global average rate. The coinciding poleward range extension of C. rodgersii has caused major ecological changes, threatening reef biodiversity and fisheries productivity. We investigated microsatellite diversity and population structure associated with range expansion by this species. Generalized linear model analyses revealed no reduction in genetic diversity in the newly colonized region. A ‘seascape genetics’ analysis of genetic distances found no spatial genetic structure associated with the range extension. The distinctive genetic characteristic of the extension zone populations was reduced population‐specific F_ST, consistent with very rapid population expansion. Demographic and genetic simulations support our inference of high connectivity between pre‐ and post‐extension zones. Thus, the range shift appears to be a poleward extension of the highly‐connected rangewide population of C. rodgersii. This is consistent with advection of larvae by the intensified warm water East Australian current, which has also increased Tasmanian Sea temperatures above the species’ lower developmental threshold. Thus, ocean circulation changes have improved the climatic suitability of novel habitat for C. rodgersii and provided the supply of recruits necessary for colonization.     

Turnover and trends in butterfly communities on two British tidal islands: stochastic influences and deterministic factors

Aim Community trends were investigated for two small islands and two local mainland butterfly communities within the UK over a period of 20–30 years. Location Hilbre Island off the Wirral Peninsula at 53.33° N, 3.10° W; Lindisfarne, an island off the Northumberland coast at 56.41° N, 1.48° W; Leighton Moss at 54.08° N, 2.26° W; Wyre Forest at 52.23° N, 2.14° W, UK. Methods Butterfly species data were collected on Hilbre and two mainland sites (Leighton Moss and Wyre Forest) from 1983 to 2006, and on Lindisfarne from 1977 to 2006, as part of the National Habitat Survey, the UK Butterfly Monitoring Scheme and ‘Butterflies for the New Millennium Atlas’ recording. Matrices of associations (Sokal and Michener’s matching coefficient S_SM; resemblance coefficient) were computed between years and subject to non‐metric multidimensional scaling (NMDS) and Mantel tests. The pattern of extinctions and colonizations at sites were examined, their heterogeneity tested by applying a Friedman test to fractional incidences for the same years. Regression analysis (multiple regression and logit regression) was used to relate butterfly numbers and incidences to climate variables, time and previous records. Results Significant community trends based on population counts and species’ incidences were found for all four sites. There was a significant climatic signal for Hilbre; although this was not apparent for the remaining sites, significant associations occurred between records for a number of species and climatic variables at all sites. Substantial turnover of species on the islands was inversely related to numbers of records for species but not to their conspicuousness to recorders. Main conclusions We argue that time trends are widespread in butterfly communities, even for relatively short periods; they are largely generated by stochastic influences rather than by more substantive factors such as climate change. Potential biases in surveying and recording history are shown to be unlikely. A clear climate signal was found only for the small Hilbre Island, for which there was also evidence for the significant influence of colonization capability of individual source species. We conclude that for many species, small islands will be sinks or pseudosinks and their ‘populations’ vulnerable to small changes in source–sink dynamics.     

Landscape structure and the genetic effects of a population collapse

Both landscape structure and population size fluctuations influence population genetics. While independent effects of these factors on genetic patterns and processes are well studied, a key challenge is to understand their interaction, as populations are simultaneously exposed to habitat fragmentation and climatic changes that increase variability in population size. In a population network of an alpine butterfly, abundance declined 60–100% in 2003 because of low over-winter survival. Across the network, mean microsatellite genetic diversity did not change. However, patch connectivity and local severity of the collapse interacted to determine allelic richness change within populations, indicating that patch connectivity can mediate genetic response to a demographic collapse. The collapse strongly affected spatial genetic structure, leading to a breakdown of isolation-by-distance and loss of landscape genetic pattern. Our study reveals important interactions between landscape structure and temporal demographic variability on the genetic diversity and genetic differentiation of populations. Projected future changes to both landscape and climate may lead to loss of genetic variability from the studied populations, and selection acting on adaptive variation will likely occur within the context of an increasing influence of genetic drift.     

Topographic,latitudinal and climatic distribution of Pinus coulteri: geographic range limits are not at the edge of the climate envelope

With changing climate, many species are projected to move poleward or to higher elevations to track suitable climates. The prediction that species will move poleward assumes that geographically marginal populations are at the edge of the species' climatic range. We studied Pinus coulteri from the center to the northern (poleward) edge of its range, and examined three scenarios regarding the relationship between the geographic and climatic margins of a species' range. We used herbarium and iNaturalist.org records to identify P. coulteri sites, generated a species distribution model based on temperature, precipitation, climatic water deficit, and actual evapotranspiration, and projected suitability under future climate scenarios. In fourteen populations from the central to northern portions of the range, we conducted field studies and recorded elevation, slope and aspect (to estimate solar insolation) to examine relationships between local and regional distributions. We found that northern populations of P. coulteri do not occupy the cold or wet edge of the species' climatic range; mid‐latitude, high elevation populations occupy the cold margin. Aspect and insolation of P. coulteri populations changed significantly across latitudes and elevations. Unexpectedly, northern, low‐elevation stands occupy north‐facing aspects and receive low insolation, while central, high‐elevation stands grow on more south‐facing aspects that receive higher insolation. Modeled future climate suitability is projected to be highest in the central, high elevation portion of the species range, and in low‐lying coastal regions under some scenarios, with declining suitability in northern areas under most future scenarios. For P. coulteri, the lack of high elevation habitat combined with a major dispersal barrier may limit northward movement in response to a warming climate. Our analyses demonstrate the importance of distinguishing geographically vs. climatically marginal populations, and the importance of quantitative analysis of the realized climate space to understand species range limits.     

Problems with using large‐scale oceanic climate indices to compare climatic sensitivities across populations and species

To understand which populations and species are most sensitive to climate change, studies correlate time series of climate variables with those of traits important for population dynamics, and subsequently compare which aspects of a species’ ecology or life‐history best explain variation in climate sensitivity. Often large‐scale oceanic climate indices (LOCIs) are used as a proxy for local climatic drivers, with many studies reporting geographic gradients in climate sensitivity to LOCIs (e.g. suggesting that species living further from the equator are relatively climate sensitive). However, the relationship between LOCIs and local weather variables also varies geographically, raising the possibility that apparent intra‐ and inter‐specific differences in climate sensitivity to LOCIs could also reflect geographic variation in how well LOCIs function as a proxy for local climatic drivers. This hypothesis is rarely tested due to lack of knowledge about the specific local climatic drivers. Here we show, using reproductive and climate data from 16 long‐term population studies of 7 Australian fairy‐wren species (Malurus genus), that the use of LOCIs can result in 1) strong overestimation of the amount of inter‐specific variation in climate sensitivity and 2) spurious patterns, particularly geographic gradients. Consequently a paradox emerges: LOCIs often explain much of the temporal variation in traits important for population dynamics, but the common usage of LOCIs may prevent meaningful intra‐ and inter‐specific comparisons of climate sensitivities over large spatial scales. Our results thus may offer an alternative interpretation of the widely reported geographic gradients in sensitivity to LOCIs. Future progress will likely require better knowledge about the identity and temporal features of local environmental drivers of population dynamics.     

Less favourable climates constrain demographic strategies in plants

Correlative species distribution models are based on the observed relationship between species’ occurrence and macroclimate or other environmental variables. In climates predicted less favourable populations are expected to decline, and in favourable climates they are expected to persist. However, little comparative empirical support exists for a relationship between predicted climate suitability and population performance. We found that the performance of 93 populations of 34 plant species worldwide – as measured by in situ population growth rate, its temporal variation and extinction risk – was not correlated with climate suitability. However, correlations of demographic processes underpinning population performance with climate suitability indicated both resistance and vulnerability pathways of population responses to climate: in less suitable climates, plants experienced greater retrogression (resistance pathway) and greater variability in some demographic rates (vulnerability pathway). While a range of demographic strategies occur within species’ climatic niches, demographic strategies are more constrained in climates predicted to be less suitable.     

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.     

Historical and event‐based bioclimatic suitability predicts regional forest vulnerability to compound effects of severe drought and bark beetle infestation

Vulnerability to climate change, and particularly to climate extreme events, is expected to vary across species ranges. Thus, we need tools to standardize the variability in regional climatic legacy and extreme climate across populations and species. Extreme climate events (e.g., droughts) can erode populations close to the limits of species' climatic tolerance. Populations in climatic‐core locations may also become vulnerable because they have developed a greater demand for resources (i.e., water) that cannot be enough satisfied during the periods of scarcity. These mechanisms can become exacerbated in tree populations when combined with antagonistic biotic interactions, such as insect infestation. We used climatic suitability indices derived from Species Distribution Models (SDMs) to standardize the climatic conditions experienced across Pinus edulis populations in southwestern North America, during a historical period (1972–2000) and during an extreme event (2001–2007), when the compound effect of hot drought and bark beetle infestation caused widespread die‐off and mortality. Pinus edulis climatic suitability diminished dramatically during the die‐off period, with remarkable variation between years. P. edulis die‐off occurred mainly not just in sites that experienced lower climatic suitability during the drought but also where climatic suitability was higher during the historical period. The combined effect of historically high climatic suitability and a marked decrease in the climatic suitability during the drought best explained the range‐wide mortality. Lagged effects of climatic suitability loss in previous years and co‐occurrence of Juniperus monosperma also explained P. edulis die‐off in particular years. Overall, the study shows that past climatic legacy, likely determining acclimation, together with competitive interactions plays a major role in responses to extreme drought. It also provides a new approach to standardize the magnitude of climatic variability across populations using SDMs, improving our capacity to predict population's or species' vulnerability to climatic change.     

Molecular substitution rate increases with latitude in butterflies: evidence for a trans‐glacial latitudinal layering of populations?

A well‐documented consequence of repeated global ice ages is the negative relationship between latitude and intraspecific genetic diversity. However, little is known about additional effects of such major climatic events on population genetic structure. Here we studied the phylogeographic structure of five lycaenid butterfly species with varied ecological adaptations, sampled across a latitudinal gradient in the Holarctic region. We found a positive correlation between latitude and substitution rate of mitochondrial DNA sequences in all species investigated. We propose that this result is the signal of increased genetic drift and founder effects during post‐glacial recolonization of northern populations. Given that phylogenetic branch length is the result of a cumulative process over evolutionary time, we hypothesize that a latitudinal layering of populations has generally been maintained during repeated cycles of glaciation, possibly due to a neutral spatial effect and/or local adaptive advantage. This trans‐glacial latitudinal layering could be viewed as a particular case of the more general phenomenon of intraspecific structuring that is created and maintained in a fluctuating environmental gradient.     

A unified framework for diversity gradients: the adaptive trait continuum

Aim Adaptive trait continua are axes of covariation observed in multivariate trait data for a given taxonomic group. These continua quantify and summarize life‐history variation at the inter‐specific level in multi‐specific assemblages. Here we examine whether trait continua can provide a useful framework to link life‐history variation with demographic and evolutionary processes in species richness gradients. Taking an altitudinal species richness gradient for Mediterranean butterflies as a study case, we examined a suite of traits (larval diet breadth, adult phenology, dispersal capacity and wing length) and species‐specific habitat measures (temperature and aridity breadth). We tested whether traits and species‐specific habitat measures tend to co‐vary, whether they are phylogenetically conserved, and whether they are able to explain species distributions and spatial genetic variation in a large number of butterfly assemblages. Location Catalonia, Spain. Methods We formulated predictions associated with species richness gradients and adaptive trait continua. We applied principal components analyses (PCAs), structural equation modelling and phylogenetic generalized least squares models. Results We found that traits and species‐specific habitat measures covaried along a main PCA axis, ranging from multivoltine trophic generalists with high dispersal capacity to univoltine (i.e. one generation per year), trophic specialist species with low dispersal capacity. This trait continuum was closely associated with the observed distributions along the altitudinal gradient and predicted inter‐specific differences in patterns of spatial genetic variability (F_ST and genetic distances), population responses to the impacts of global change and local turnover dynamics. Main conclusions The adaptive trait continuum of Mediterranean butterflies provides an integrative and mechanistic framework to: (1) analyse geographical gradients in species richness, (2) explain inter‐specific differences in population abundances, spatial distributions and demographic trends, (3) explain inter‐specific differences in patterns of genetic variation (F_ST and genetic distances), and (4) study specialist–generalist life‐history transitions frequently involved in butterfly diversification processes.     

BioMove – an integrated platform simulating the dynamic response of species to environmental change

BioMove simulates plant species' geographic range shifts in response to climate, habitat structure and disturbance, at annual time steps. This spatially explicit approach integrates species' bioclimatic suitability and population‐level demographic rates with simulation of landscape‐level processes (dispersal, disturbance, species' response to dynamic dominant vegetation structure). Species population dynamics are simulated through matrix modelling that includes scaling demographic rates by climatic suitability. Dispersal functions simulate population spread. User‐specified plant functional types (PFTs) provide vegetation structure that determines resource competition and disturbance. PFTs respond annually through dispersal, inter‐PFT competition and demographic shifts. BioMove provides a rich framework for dynamic range simulations.