Range shifting on a fragmented landscape

Where and why are species' range shifts hampered by unsuitable landscapes?

There is widespread concern that species will fail to track climate change if habitat is too scarce or insufficiently connected. Targeted restoration has been advocated to help species adapt, and a "conductance" metric has been proposed, based on simulation studies, to predict effective habitat configurations. However, until now there is very little empirical evidence on how the configuration of habitat is affecting expansion at species' cool range margins. We analysed the colonisation events that have occurred in continuously monitored trap locations for 54 species of southerly distributed moths in Britain between 1985 and 2011. We tested whether the time until colonisation was affected by attributes of each species, and of intervening landcover and climate between the trap and the baseline distribution (1965-1985). For woodland species, the time until colonisation of new locations was predicted by the "conductance" of woodland habitat, and this relationship was general, regardless of species' exact dispersal distances and habitat needs. This shows that contemporary range shifts are being influenced by habitat configuration as well as simple habitat extent. For species associated with farmland or suburban habitats, colonisation was significantly slower through landscapes with a high variance in elevation and/or temperature. Therefore, it is not safe to assume that such relatively tolerant species face no geographical barriers to range expansion. We thus elucidate how species' attributes interact with landscape characteristics to create highly heterogeneous patterns of shifting at cool range margins. Conductance, and other predictors of range shifts, can provide a foundation for developing coherent conservation strategies to manage range shifts for entire communities.

The collective effect of finite-sized inhomogeneities on the spatial spread of populations in two dimensions

The dynamics of a population expanding into unoccupied habitat has been primarily studied for situations in which growth and dispersal parameters are uniform in space or vary in one dimension. Here, we study the influence of finite-sized individual inhomogeneities and their collective effect on front speed if randomly placed in a two-dimensional habitat. We use an individual-based model to investigate the front dynamics for a region in which dispersal or growth of individuals is reduced to zero (obstacles) or increased above the background (hotspots), respectively. In a regime where front dynamics is determined by a local front speed only, a principle of least time can be employed to predict front speed and shape. The resulting analytical solutions motivate an event-based algorithm illustrating the effects of several obstacles or hotspots. We finally apply the principle of least time to large heterogeneous environments by solving the Eikonal equation numerically. Obstacles lead to a slow-down that is dominated by the number density and width of obstacles, but not by their precise shape. Hotspots result in a speed-up, which we characterize as function of hotspot strength and density. Our findings emphasize the importance of taking the dimensionality of the environment into account.

Changing environmental gradients over forty years alter ecomorphological variation in Guadalupe Bass Micropterus treculii throughout a river basin

Understanding the degree of intraspecific variation within and among populations is a key aspect of predicting the capacity of a species to respond to anthropogenic disturbances. However, intraspecific variation is usually assessed at either limited temporal, but broad spatial scales or vice versa, which can make assessing changes in response to long-term disturbances challenging. We evaluated the relationship between the longitudinal gradient of changing flow regimes and land use/land cover patterns since 1980 and morphological variation of Guadalupe Bass Micropterus treculii throughout the Colorado River Basin of central Texas. The Colorado River Basin in Texas has experienced major alterations to the hydrologic regime due to changing land- and water-use patterns. Historical collections of Guadalupe Bass prior to rapid human-induced change present the unique opportunity to study the response of populations to varying environmental conditions through space and time. Morphological differentiation of Guadalupe Bass associated with temporal changes in flow regimes and land use/land cover patterns suggests that they are exhibiting intraspecific trait variability, with contemporary individuals showing increased body depth, in response to environmental alteration through time (specifically related to an increase in herbaceous land cover, maximum flows, and the number of low pulses and high pulses). Additionally, individuals from tributaries with increased hydrologic alteration associated with urbanization or agricultural withdrawals tended to have a greater distance between the anal and caudal fin. These results reveal trait variation that may help to buffer populations under conditions of increased urbanization and sprawl, human population growth, and climate risk, all of which impose novel selective pressures, especially on endemic species like Guadalupe Bass. Our results contribute an understanding of the adaptability and capacity of an endemic population to respond to expected future changes based on demographic or climatic projection.

Expansion Under Climate Change: The Genetic Consequences

Range expansion and range shifts are crucial population responses to climate change. Genetic consequences are not well understood but are clearly coupled to ecological dynamics that, in turn, are driven by shifting climate conditions. We model a population with a deterministic reaction-diffusion model coupled to a heterogeneous environment that develops in time due to climate change. We decompose the resulting travelling wave solution into neutral genetic components to analyse the spatio-temporal dynamics of its genetic structure. Our analysis shows that range expansions and range shifts under slow climate change preserve genetic diversity. This is because slow climate change creates range boundaries that promote spatial mixing of genetic components. Mathematically, the mixing leads to so-called pushed travelling wave solutions. This mixing phenomenon is not seen in spatially homogeneous environments, where range expansion reduces genetic diversity through gene surfing arising from pulled travelling wave solutions. However, the preservation of diversity is diminished when climate change occurs too quickly. Using diversity indices, we show that fast expansions and range shifts erode genetic diversity more than slow range expansions and range shifts. Our study provides analytical insight into the dynamics of travelling wave solutions in heterogeneous environments.

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.

Pushing the pace of tree species migration

Plants and animals have responded to past climate changes by migrating with habitable environments, sometimes shifting the boundaries of their geographic ranges by tens of kilometers per year or more. Species migrating in response to present climate conditions, however, must contend with landscapes fragmented by anthropogenic disturbance. We consider this problem in the context of wind-dispersed tree species. Mechanisms of long-distance seed dispersal make these species capable of rapid migration rates. Models of species-front migration suggest that even tree species with the capacity for long-distance dispersal will be unable to keep pace with future spatial changes in temperature gradients, exclusive of habitat fragmentation effects. Here we present a numerical model that captures the salient dynamics of migration by long-distance dispersal for a generic tree species. We then use the model to explore the possible effects of assisted colonization within a fragmented landscape under a simulated tree-planting scheme. Our results suggest that an assisted-colonization program could accelerate species-front migration rates enough to match the speed of climate change, but such a program would involve an environmental-sustainability intervention at a massive scale.

Between migration load and evolutionary rescue: dispersal, adaptation and the response of spatially structured populations to environmental change

The evolutionary potential of populations is mainly determined by population size and available genetic variance. However, the adaptability of spatially structured populations may also be affected by dispersal: positively by spreading beneficial mutations across sub-populations, but negatively by moving locally adapted alleles between demes. We develop an individual-based, two-patch, allelic model to investigate the balance between these opposing effects on a population's evolutionary response to rapid climate change. Individual fitness is controlled by two polygenic traits coding for local adaptation either to the environment or to climate. Under conditions of selection that favour the evolution of a generalist phenotype (i.e. weak divergent selection between patches) dispersal has an overall positive effect on the persistence of the population. However, when selection favours locally adapted specialists, the beneficial effects of dispersal outweigh the associated increase in maladaptation for a narrow range of parameter space only (intermediate selection strength and low linkage among loci), where the spread of beneficial climate alleles is not strongly hampered by selection against non-specialists. Given that local selection across heterogeneous and fragmented landscapes is common, the complex effect of dispersal that we describe will play an important role in determining the evolutionary dynamics of many species under rapidly changing climate.

Disequilibrium vegetation dynamics under future climate change

PREMISE OF THE STUDY

Near-future climate changes are likely to elicit major vegetation changes. Disequilibrium dynamics, which occur when vegetation comes out of equilibrium with climate, are potentially a key facet of these. Understanding these dynamics is crucial for making accurate predictions, informing conservation planning, and understanding likely changes in ecosystem function on time scales relevant to society. However, many predictive studies have instead focused on equilibrium end-points with little consideration of the transient trajectories.

METHODS

We review what we should expect in terms of disequilibrium vegetation dynamics over the next 50-200 yr, covering a broad range of research fields including paleoecology, macroecology, landscape ecology, vegetation science, plant ecology, invasion biology, global change biology, and ecosystem ecology.

KEY RESULTS

The expected climate changes are likely to induce marked vegetation disequilibrium with climate at both leading and trailing edges, with leading-edge disequilibrium dynamics due to lags in migration at continental to landscape scales, in local population build-up and succession, in local evolutionary responses, and in ecosystem development, and trailing-edge disequilibrium dynamics involving delayed local extinctions and slow losses of ecosystem structural components. Interactions with habitat loss and invasive pests and pathogens are likely to further contribute to disequilibrium dynamics. Predictive modeling and climate-change experiments are increasingly representing disequilibrium dynamics, but with scope for improvement.

CONCLUSIONS

The likely pervasiveness and complexity of vegetation disequilibrium is a major challenge for forecasting ecological dynamics and, combined with the high ecological importance of vegetation, also constitutes a major challenge for future nature conservation.

Risky movement increases the rate of range expansion

The movement rules used by an individual determine both its survival and dispersal success. Here, we develop a simple model that links inter-patch movement behaviour with population dynamics in order to explore how individual dispersal behaviour influences not only its dispersal and survival, but also the population's rate of range expansion. Whereas dispersers are most likely to survive when they follow nearly straight lines and rapidly orient movement towards a non-natal patch, the most rapid rates of range expansion are obtained for trajectories in which individuals delay biasing their movement towards a non-natal patch. This result is robust to the spatial structure of the landscape. Importantly, in a set of evolutionary simulations, we also demonstrate that the movement strategy that evolves at an expanding front is much closer to that maximizing the rate of range expansion than that which maximizes the survival of dispersers. Our results suggest that if one of our conservation goals is the facilitation of range-shifting, then current indices of connectivity need to be complemented by the development and utilization of new indices providing a measure of the ease with which a species spreads across a landscape.

Comparing organic farming and land sparing: optimizing yield and butterfly populations at a landscape scale

Organic farming aims to be wildlife-friendly, but it may not benefit wildlife overall if much greater areas are needed to produce a given quantity of food. We measured the density and species richness of butterflies on organic farms, conventional farms and grassland nature reserves in 16 landscapes. Organic farms supported a higher density of butterflies than conventional farms, but a lower density than reserves. Using our data, we predict the optimal land-use strategy to maintain yield whilst maximizing butterfly abundance under different scenarios. Farming conventionally and sparing land as nature reserves is better for butterflies when the organic yield per hectare falls below 87% of conventional yield. However, if the spared land is simply extra field margins, organic farming is optimal whenever organic yields are over 35% of conventional yields. The optimal balance of land sparing and wildlife-friendly farming to maintain production and biodiversity will differ between landscapes.

How range shifts induced by climate change affect neutral evolution

We investigate neutral evolution during range shifts in a strategic model of a metapopulation occupying a climate gradient. Using heritable, neutral markers, we track the spatio-temporal fate of lineages. Owing to iterated founder effects ('mutation surfing'), survival of lineages derived from the leading range limit is enhanced. At trailing limits, where habitat suitability decreases, survival is reduced (mutations 'wipe out'). These processes alter (i) the spatial spread of mutations, (ii) origins of persisting mutations and (iii) the generation of diversity. We show that large changes in neutral evolution can be a direct consequence of range shifting.

Evolved dispersal strategies at range margins

Dispersal is a key component of a species's ecology and will be under different selection pressures in different parts of the range. For example, a long-distance dispersal strategy suitable for continuous habitat at the range core might not be favoured at the margin, where the habitat is sparse. Using a spatially explicit, individual-based, evolutionary simulation model, the dispersal strategies of an organism that has only one dispersal event in its lifetime, such as a plant or sessile animal, are considered. Within the model, removing habitat, increasing habitat turnover, increasing the cost of dispersal, reducing habitat quality or altering vital rates imposes range limits. In most cases, there is a clear change in the dispersal strategies across the range, although increasing death rate towards the margin has little impact on evolved dispersal strategy across the range. Habitat turnover, reduced birth rate and reduced habitat quality all increase evolved dispersal distances at the margin, while increased cost of dispersal and reduced habitat density lead to lower evolved dispersal distances at the margins. As climate change shifts suitable habitat poleward, species ranges will also start to shift, and it will be the dispersal capabilities of marginal populations, rather than core populations, that will influence the rate of range shifting.

Landscape structure and boundary effects determine the fate of mutations occurring during range expansions

The interplay between the spatial dynamics of range expansion and evolutionary processes is receiving considerable attention. Recent theory has demonstrated that mutations occurring towards the front of a spatially expanding population can sometimes 'surf' to high frequency and spatial extent. Here, we extend this work to consider how the fate of a novel mutation is influenced by where and when it occurs. Specifically, we are interested in establishing how the origin of a mutation relative to a habitat edge influences its dynamics, and in understanding how this is mediated by the behaviour of individuals at those boundaries. Using a coupled-map lattice model, we demonstrate that the survival probability, abundance and spatial extent of surviving mutants can depend on their origin. An edge effect is often observed and can be quite different both qualitatively and quantitatively depending on the behavioural rules assumed. Mutations, especially those that are deleterious, that arise at a habitat edge with reflective boundary conditions can be many more times likely to survive for substantial periods of time than those that arise away from the edge. Conversely, with absorbing boundary conditions, their survival is greater when they arise well away from the edge. Our results clearly illustrate that landscape structure, habitat edges and boundary conditions have a considerable influence on the likely fate of mutations that occur during a period of range expansion.

Landscape as a model: the importance of geometry

In all models, but especially in those used to predict uncertain processes (e.g., climate change and nonnative species establishment), it is important to identify and remove any sources of bias that may confound results. This is critical in models designed to help support decisionmaking. The geometry used to represent virtual landscapes in spatially explicit models is a potential source of bias. The majority of spatial models use regular square geometry, although regular hexagonal landscapes have also been used. However, there are other ways in which space can be represented in spatially explicit models. For the first time, we explicitly compare the range of alternative geometries available to the modeller, and present a mechanism by which uncertainty in the representation of landscapes can be incorporated. We test how geometry can affect cell-to-cell movement across homogeneous virtual landscapes and compare regular geometries with a suite of irregular mosaics. We show that regular geometries have the potential to systematically bias the direction and distance of movement, whereas even individual instances of landscapes with irregular geometry do not. We also examine how geometry can affect the gross representation of real-world landscapes, and again show that individual instances of regular geometries will always create qualitative and quantitative errors. These can be reduced by the use of multiple randomized instances, though this still creates scale-dependent biases. In contrast, virtual landscapes formed using irregular geometries can represent complex real-world landscapes without error. We found that the potential for bias caused by regular geometries can be effectively eliminated by subdividing virtual landscapes using irregular geometry. The use of irregular geometry appears to offer spatial modellers other potential advantages, which are as yet underdeveloped. We recommend their use in all spatially explicit models, but especially for predictive models that are used in decisionmaking.

Many different areas of science try to simulate and predict (model) how processes act across virtual landscapes. Sometimes these models are abstract, but often they are based on real-world landscapes and are used to make real-world planning or management decisions. We considered two separate issues: how movement occurs across landscapes and how uncertainty in spatial data can be represented in the model. Most studies represent the landscape using regular geometries (e.g., squares and hexagons), but we generated landscapes of irregular shapes. We tested and compared how the shapes that make up a landscape affected cell-to-cell movement across it. All of the virtual landscapes formed with regular geometries had the potential to bias the direction and distance of movement. Those formed with irregular geometry did not. We have also shown that describing whole real-world landscapes with regular geometries will lead to errors and bias, whereas virtual landscapes formed with irregular geometries are free from both. We recommend the use of multiple versions of virtual landscapes formed using irregular geometries for all spatially explicit models as a way of minimizing this source of bias and error; this is especially relevant in predictive models (e.g., climate change) that are difficult to test and are designed to help make decisions.