Restoring spatiotemporal variability to enhance the capacity for dispersal-limited species to track climate change.

Climate refugia are areas where species can persist through climate change with little to no movement. Among the factors associated with climate refugia are high spatial heterogeneity, such that there is only a short distance between current and future optimal climates, as well as biotic or abiotic environmental factors that buffer against variability in time. However, these types of climate refugia may be declining due to anthropogenic homogenization of environments and degradation of environmental buffers. To quantify the potential for restoration of refugia-like environmental conditions to increase population persistence under climate change, we simulated a population's capacity to track their temperature over space and time given different levels of spatial and temporal variability in temperature. To determine how species traits affected the efficacy of restoring heterogeneity, we explored an array of values for species' dispersal ability, thermal tolerance, and fecundity. We found that species were more likely to persist in environments with higher spatial heterogeneity and lower environmental stochasticity. When simulating a management action that increased the spatial heterogeneity of a previously homogenized environment, species were more likely to persist through climate change, and population sizes were generally higher, but there was little effect with mild temperature change. The benefits of heterogeneity restoration were greatest for species with limited dispersal ability. In contrast, species with longer dispersal but lower fecundity were more likely to benefit from a reduction in environmental stochasticity than an increase in spatial heterogeneity. Our results suggest that restoring environments to refugia-like conditions could promote species' persistence under the influence of climate change in addition to conservation strategies such as assisted migration, corridors, and increased protection.

Unclear relationships between mean survival rate and its environmental variance in vertebrates.

Current environmental changes may increase temporal variability of life history traits of species thus affecting their long-term population growth rate and extinction risk. If there is a general relationship between environmental variances (EVs) and mean annual survival rates of species, that relationship could be used as a guideline for analyses of population growth and extinction risk for populations, where data on EVs are missing. For this purpose, we present a comprehensive compilation of 252 EV estimates from 89 species belonging to five vertebrate taxa (birds, mammals, reptiles, amphibians and fish) covering mean annual survival rates from 0.01 to 0.98. Since variances of survival rates are constrained by their means, particularly for low and high mean survival rates, we assessed whether any observed relationship persisted after applying two types of commonly used variance stabilizing transformations: relativized EVs (observed/mathematical maximum) and logit-scaled EVs. With raw EVs at the arithmetic scale, mean-variance relationships of annual survival rates were hump-shaped with small EVs at low and high mean survival rates and higher (and widely variable) EVs at intermediate mean survival rates. When mean annual survival rates were related to relativized EVs the hump-shaped pattern was less distinct than for raw EVs. When transforming EVs to logit scale the relationship between mean annual survival rates and EVs largely disappeared. The within-species juvenile-adult slopes were mainly positive at low (<0.5) and negative at high (>0.5) mean survival rates for raw and relativized variances while these patterns disappeared when EVs were logit transformed. Uncertainties in how to interpret the results of relativized and logit-scaled EVs, and the observed high variation in EV's for similar mean annual survival rates illustrates that extrapolations of observed EVs and tests of life history drivers of survival-EV relationships need to also acknowledge the large variation in these parameters.

Inferring stability and persistence in the vaginal microbiome: A stochastic model of ecological dynamics.

The interplay of stochastic and ecological processes that govern the establishment and persistence of host-associated microbial communities is not well understood. Here we illustrate the conceptual and practical advantages of fitting stochastic population dynamics models to multi-species bacterial time series data. We show how the stability properties, fluctuation regimes and persistence probabilities of human vaginal microbial communities can be better understood by explicitly accommodating three sources of variability in ecological stochastic models of multi-species abundances: 1) stochastic biotic and abiotic forces, 2) ecological feedback and 3) sampling error. Rooting our modeling tool in stochastic population dynamics modeling theory was key to apply standardized measures of a community's reaction to environmental variation that ultimately depends on the nature and intensity of the intra-specific and inter-specific interaction strengths. Using estimates of model parameters, we developed a Risk Prediction Monitoring (RPM) tool that estimates temporal changes in persistence probabilities for any bacterial group of interest. This method mirrors approaches that are often used in conservation biology in which a measure of extinction risks is periodically updated with any change in a population or community. Additionally, we show how to use estimates of interaction strengths and persistence probabilities to formulate hypotheses regarding the molecular mechanisms and genetic composition that underpin different types of interactions. Instead of seeking a definition of "dysbiosis" we propose to translate concepts of theoretical ecology and conservation biology methods into practical approaches for the management of human-associated bacterial communities.

Pace and parity predict the short-term persistence of small plant populations.

Life history traits are used to predict asymptotic odds of extinction from dynamic conditions. Less is known about how life history traits interact with stochasticity and population structure of finite populations to predict near-term odds of extinction. Through empirically parameterized matrix population models, we study the impact of life history (reproduction, pace), stochasticity (environmental, demographic), and population history (existing, novel) on the transient population dynamics of finite populations of plant species. Among fast and slow pace and either a uniform or increasing reproductive intensity or short or long reproductive lifespan, slow, semelparous species are at the greatest risk of extinction. Long reproductive lifespans buffer existing populations from extinction while the odds of extinction of novel populations decrease when the reproductive effort is uniformly spread across the reproductive lifespan. Our study highlights the importance of population structure, pace, and two distinct aspects of parity for predicting near-term odds of extinction.

Why so many polyploids? Accounting for environmental stochasticity in unreduced gamete formation lowers the perceived barriers to polyploid establishment.

While polyploids are common in nature, existing models suggest that polyploid establishment should be difficult and rare. We explore this apparent paradox by focussing on the role of unreduced gametes, as their union is the main route for the formation of neopolyploids. Production of such gametes is affected by genetic and environmental factors, resulting in variation in the formation rate of unreduced gametes (u). Once formed, neopolyploids face minority cytotype exclusion (MCE) due to a lack of viable mating opportunities. More than a dozen theoretical models have explored factors that could permit neopolyploids to overcome MCE and become established. Until now, however, none have explored variability in u and its consequences for the rate of polyploid establishment. Here, we determine the distribution that best fits the available empirical data on u. We perform a global sensitivity analysis exploring the consequences of using empirical distributions of u to investigate effects on polyploid establishment. We determined that in many cases, u is best fit by a log-normal distribution. We found environmental stochasticity in u dramatically impacts model predictions when compared to a static u. Our results help reconcile previous modelling results suggesting high barriers to the polyploid establishment with the observation that polyploids are common in nature.

Assessing the sensitivity and resistance of communities in a changing environment.

We propose that the ecological resilience of communities to permanent changes of the environment can be based on how variation in the overall abundance of individuals affects the number of species. Community sensitivity is defined as the ratio between the rate of change in the log expected number of species and the rate of change in the log expected number of individuals in the community. High community sensitivity means that small changes in the total abundance strongly impact the number of species. Community resistance is the proportional reduction in expected number of individuals that the community can sustain before expecting to lose one species. A small value of community resistance means that the community can only endure a small reduction in abundance before it is expected to lose one species. Based on long-term studies of four bird communities in European deciduous forests at different latitudes large differences were found in the resilience to environmental perturbations. Estimating the variance components of the species abundance distribution revealed how different processes contributed to the community sensitivity and resistance. Species heterogeneity in the population dynamics was the largest component, but its proportion varied among communities. Species-specific response to environmental fluctuations was the second major component of the variation in abundance. Estimates of community sensitivity and resistance based on data only from a single year were in general larger than those based on estimates from longer time series. Thus, our approach can provide rapid and conservative assessment of the resilience of communities to environmental changes also including only short-term data. This study shows that a general ecological mechanism, caused by increased strength of density dependence due to reduction in resource availability, can provide an intuitive measure of community resilience to environmental variation. Our analyses also illustrate the importance of including specific assumptions about how different processes affect community dynamics. For example, if stochastic fluctuations in the environment affect all species in a similar way, the sensitivity and resistance of the community to environmental changes will be different from communities in which all species show independent responses.

Processes governing species richness in communities exposed to temporal environmental stochasticity: A review and synthesis of modelling approaches.

Research into the processes governing species richness has often assumed that the environment is fixed, whereas realistic environments are often characterised by random fluctuations over time. This temporal environmental stochasticity (TES) changes the demographic rates of species populations, with cascading effects on community dynamics and species richness. Theoretical and applied studies have used process-based mathematical models to determine how TES affects species richness, but under a variety of frameworks. Here, we critically review such studies to synthesise their findings and draw general conclusions. We first provide a broad mathematical framework encompassing the different ways in which TES has been modelled. We then review studies that have analysed models with TES under the assumption of negligible interspecific interactions, such that a community is conceptualised as the sum of independent species populations. These analyses have highlighted how TES can reduce species richness by increasing the frequency at which a species becomes rare and therefore prone to extinction. Next, we review studies that have relaxed the assumption of negligible interspecific interactions. To simplify the corresponding models and make them analytically tractable, such studies have used mean-field theory to derive fixed parameters representing the typical strength of interspecific interactions under TES. The resulting analyses have highlighted community-level effects that determine how TES affects species richness, for species that compete for a common limiting resource. With short temporal correlations of environmental conditions, a non-linear averaging effect of interspecific competition strength over time gives an increase in species richness. In contrast, with long temporal correlations of environmental conditions, strong selection favouring the fittest species between changes in environmental conditions results in a decrease in species richness. We compare such results with those from invasion analysis, which examines invasion growth rates (IGRs) instead of species richness directly. Qualitative differences sometimes arise because the IGR is the expected growth rate of a species when it is rare, which does not capture the variation around this mean or the probability of the species becoming rare. Our review elucidates key processes that have been found to mediate the negative and positive effects of TES on species richness, and by doing so highlights key areas for future research.

Evolutionary rescue under demographic and environmental stochasticity.

Populations suffer two types of stochasticity: demographic stochasticity, from sampling error in offspring number, and environmental stochasticity, from temporal variation in the growth rate. By modelling evolution through phenotypic selection following an abrupt environmental change, we investigate how genetic and demographic dynamics, as well as effects on population survival of the genetic variance and of the strength of stabilizing selection, differ under the two types of stochasticity. We show that population survival probability declines sharply with stronger stabilizing selection under demographic stochasticity, but declines more continuously when environmental stochasticity is strengthened. However, the genetic variance that confers the highest population survival probability differs little under demographic and environmental stochasticity. Since the influence of demographic stochasticity is stronger when population size is smaller, a slow initial decline of genetic variance, which allows quicker evolution, is important for population persistence. In contrast, the influence of environmental stochasticity is population-size-independent, so higher initial fitness becomes important for survival under strong environmental stochasticity. The two types of stochasticity interact in a more than multiplicative way in reducing the population survival probability. Our work suggests the importance of explicitly distinguishing and measuring the forms of stochasticity during evolutionary rescue.

Coevolution of Patch Selection in Stochastic Environments.

AbstractSpecies interact in landscapes where environmental conditions vary in time and space. This variability impacts how species select habitat patches. Under equilibrium conditions, evolution of this patch selection can result in ideal free distributions where per capita growth rates are zero in occupied patches and negative in unoccupied patches. These ideal free distributions, however, do not explain why species occupy sink patches, why competitors have overlapping spatial ranges, or why predators avoid highly productive patches. To understand these patterns, we solve for coevolutionarily stable strategies (coESSs) of patch selection for multispecies stochastic Lotka-Volterra models accounting for spatial and temporal heterogeneity. In occupied patches at the coESS, we show that the differences between the local contributions to the mean and the variance of the long-term population growth rate are equalized. Applying this characterization to models of antagonistic interactions reveals that environmental stochasticity can partially exorcize the ghost of competition past, select for new forms of enemy-free and victimless space, and generate hydra effects over evolutionary timescales. Viewing our results through the economic lens of modern portfolio theory highlights why the coESS for patch selection is often a bet-hedging strategy coupling stochastic sink populations. Our results highlight how environmental stochasticity can reverse or amplify evolutionary outcomes as a result of species interactions or spatial heterogeneity.

Partitioning variance in population growth for models with environmental and demographic stochasticity.

How demographic factors lead to variation or change in growth rates can be investigated using life table response experiments (LTRE) based on structured population models. Traditionally, LTREs focused on decomposing the asymptotic growth rate, but more recently decompositions of annual 'realized' growth rates using 'transient' LTREs have gained in popularity. Transient LTREs have been used particularly to understand how variation in vital rates translate into variation in growth for populations under long-term study. For these, complete population models may be constructed to investigate how temporal variation in environmental drivers affect vital rates. Such investigations have usually come down to estimating covariate coefficients for the effects of environmental variables on vital rates, but formal ways of assessing how they lead to variation in growth rates have been lacking. We extend transient LTREs to further partition the contributions from vital rates into contributions from temporally varying factors that affect them. The decomposition allows one to compare the resultant effect on the growth rate of different environmental factors, as well as density dependence, which may each act via multiple vital rates. We also show how realized growth rates can be decomposed into separate components from environmental and demographic stochasticity. The latter is typically omitted in LTRE analyses. We illustrate these extensions with an integrated population model (IPM) for data from a 26 years study on northern wheatears (Oenanthe oenanthe), a migratory passerine bird breeding in an agricultural landscape. For this population, consisting of around 50-120 breeding pairs per year, we partition variation in realized growth rates into environmental contributions from temperature, rainfall, population density and unexplained random variation via multiple vital rates, and from demographic stochasticity. The case study suggests that variation in first year survival via the unexplained random component, and adult survival via temperature are two main factors behind environmental variation in growth rates. More than half of the variation in growth rates is suggested to come from demographic stochasticity, demonstrating the importance of this factor for populations of moderate size.

Density regulation amplifies environmentally induced population fluctuations.

Background: Density-dependent regulation is ubiquitous in population dynamics, and its potential interaction with environmental stochasticity complicates the characterization of the random component of population dynamics. Yet, this issue has not received attention commensurate with its relevance for descriptive and predictive modeling of population dynamics. Here we use a Bayesian modeling approach to investigate the contribution of density regulation to population variability in stochastic environments. Methods: We analytically derive a formula linking the stationary variance of population abundance/density under Gompertz regulation in a stochastic environment with constant variance to the environmental variance and the strength of density feedback, to investigate whether and how density regulation affects the stationary variance. We examine through simulations whether the relationship between stationary variance and density regulation inferred analytically under the Gompertz model carries over to the Ricker model, widely used in population dynamics modeling. Results: The analytical decomposition of the stationary variance under stochastic Gompertz dynamics implies higher variability for strongly regulated populations. Simulation results demonstrate that the pattern of increasing population variability with increasing density feedback found under the Gompertz model holds for the Ricker model as well, and is expected to be a general phenomenon with stochastic population models. We also analytically established and empirically validated that the square of the autoregressive parameter of the Gompertz model in AR(1) form represents the proportion of stationary variance due to density dependence. Discussion: Our results suggest that neither environmental stochasticity nor density regulation can alone explain the patterns of population variability in stochastic environments, as these two components of temporal variation interact, with a tendency for density regulation to amplify the magnitude of environmentally induced population fluctuations. This finding has far-reaching implications for population viability. It implies that intense intra-specific resource competition increases the risk of environment-driven population collapse at high density, making opportune harvesting a sensible practice for improving the resistance of managed populations such as fish stocks to environmental perturbations. The separation of density-dependent and density-independent processes will help improve population dynamics modeling, while providing a basis for evaluating the relative importance of these two categories of processes that remains a topic of long-standing controversy among ecologists.

Generation Time in Stage-Structured Populations under Fluctuating Environments.

AbstractA common measure of generation time is the average distance between two recruitment events along a genetic lineage. In populations with stage structure that live in a constant environment, this generation time can be computed from the elasticities of stable population growth to fecundities, and it is equivalent to another common measure of generation time: the average parental age of reproductive-value-weighted offspring. Here, we show three things. First, when the environment fluctuates, the average distance between two recruitment events along a genetic lineage is computed from the elasticities of the stochastic growth rate to fecundities. Second, under environmental stochasticity, this measure of generation time remains equivalent to the average parental age of reproductive-value-weighted offspring. Third, the generation time of a population in a fluctuating environment may deviate from the generation time the population would have in the average environment.

Applying stochastic and Bayesian integral projection modeling to amphibian population viability analysis.

Integral projection models (IPMs) can estimate the population dynamics of species for which both discrete life stages and continuous variables influence demographic rates. Stochastic IPMs for imperiled species, in turn, can facilitate population viability analyses (PVAs) to guide conservation decision-making. Biphasic amphibians are globally distributed, often highly imperiled, and ecologically well suited to the IPM approach. Herein, we present a stochastic size- and stage-structured IPM for a biphasic amphibian, the U.S. federally threatened California tiger salamander (CTS) (Ambystoma californiense). This Bayesian model reveals that CTS population dynamics show greatest elasticity to changes in juvenile and metamorph growth and that populations are likely to experience rapid growth at low density. We integrated this IPM with climatic drivers of CTS demography to develop a PVA and examined CTS extinction risk under the primary threats of habitat loss and climate change. The PVA indicated that long-term viability is possible with surprisingly high (20%-50%) terrestrial mortality but simultaneously identified likely minimum terrestrial buffer requirements of 600-1000 m while accounting for numerous parameter uncertainties through the Bayesian framework. These analyses underscore the value of stochastic and Bayesian IPMs for understanding both climate-dependent taxa and those with cryptic life histories (e.g., biphasic amphibians) in service of ecological discovery and biodiversity conservation. In addition to providing guidance for CTS recovery, the contributed IPM and PVA supply a framework for applying these tools to investigations of ecologically similar species.

Seasonality, density dependence, and spatial population synchrony.

Studies of spatial population synchrony constitute a central approach for understanding the drivers of ecological dynamics. Recently, identifying the ecological impacts of climate change has emerged as a new important focus in population synchrony studies. However, while it is well known that climatic seasonality and sequential density dependence influences local population dynamics, the role of season-specific density dependence in shaping large-scale population synchrony has not received attention. Here, we present a widely applicable analytical protocol that allows us to account for both season and geographic context-specific density dependence to better elucidate the relative roles of deterministic and stochastic sources of population synchrony, including the renowned Moran effect. We exemplify our protocol by analyzing time series of seasonal (spring and fall) abundance estimates of cyclic rodent populations, revealing that season-specific density dependence is a major component of population synchrony. By accounting for deterministic sources of synchrony (in particular season-specific density dependence), we are able to identify stochastic components. These stochastic components include mild winter weather events, which are expected to increase in frequency under climate warming in boreal and Arctic ecosystems. Interestingly, these weather effects act both directly and delayed on the vole populations, thus enhancing the Moran effect. Our study demonstrates how different drivers of population synchrony, presently altered by climate warming, can be disentangled based on seasonally sampled population time-series data and adequate population models.

Laplace's demon in biology: Models of evolutionary prediction.

Our ability to predict natural phenomena can be limited by incomplete information. This issue is exemplified by "Laplace's demon," an imaginary creature proposed in the 18th century, who knew everything about everything, and thus could predict the full nature of the universe forward or backward in time. Quantum mechanics, among other things, has cast doubt on the possibility of Laplace's demon in the full sense, but the idea still serves as a useful metaphor for thinking about the extent to which prediction is limited by incomplete information on deterministic processes versus random factors. Here, we use simple analytical models and computer simulations to illustrate how data limits can be captured in a Bayesian framework, and how they influence our ability to predict evolution. We show how uncertainty in measurements of natural selection, or low predictability of external environmental factors affecting selection, can greatly reduce predictive power, often swamping the influence of intrinsic randomness caused by genetic drift. Thus, more accurate knowledge concerning the causes and action of natural selection is key to improving prediction. Fortunately, our analyses and simulations show quantitatively that reasonable improvements in data quantity and quality can meaningfully increase predictability.

Effects of predator modulation and vector preference on pathogen transmission in plant populations.

Biological control programs frequently rely on predators to control vector-borne pathogens by consumptive effects on vector abundance in agroecosystems. Meanwhile, the spread of vectored disease depends on the vector preference for host status (healthy or infected hosts). Yet, it is unclear how vector preferences alter the controlled effectivity of predators in pathogen transmission. Therefore, we here addressed the plant-vector-pathogen models assessing how pathogen transmission in plant was affected by variable predation rates and vector preferences for host status. Specifically, we discussed effects of predators on vector abundance and pathogen transmission under both a non-spatial model and a spatially structured metapopulation model. We showed that predators can decrease the vector abundance and inhibit pathogen prevalence, whereas vector preference contributes profoundly to the controlled effectivity of predators on the spread of vector-borne pathogens. Moreover, predation can increase oscillation amplitude of the pathogen prevalence in both plant and vector; suggesting that the inclusion of predator can amplify the effects of environmental stochasticity on pathogen dynamics. In conclusion, our results support the prediction of theoretical disease models showing predator can be a natural enemy for pathogen control, and also extend that predatory interactions interacting with vector preferences play the singularly joint effects on the spread of vector-borne pathogens.

Impacts of demographic and environmental stochasticity on population dynamics with cooperative effects.

Different types of stochasticity play essential roles in shaping complex population dynamics. This paper presents a novel approach to model demographic and environmental stochasticity in a single-species model with cooperative components that are measured by component Allee effects. Our work provides rigorous mathematical proof on stochastic persistence and extinction, ergodicity (i.e., the existence of a unique stationary distribution) and the existence of a nontrivial periodic solution to study the impacts of demographic and environmental stochasticity on population dynamics. The theoretical and numerical results suggest that stochasticity may affect the population system in a variety of ways, specifically: (i) In the weak Allee effects case (e.g., strong cooperative efforts), the demographic stochasticity from the attack rate contributes to the expansion of the population size, while the demographic stochasticity from the handling rate and the environmental stochasticity have the opposite role, and may even lead to population extinction; (ii) In the strong Allee effects case (cooperative efforts not strong enough), both demographic and environmental stochasticity play a similar role in the survival of population, and are related to the initial population level: if the initial population level is large enough, demographic stochasticity and environmental stochasticity may be detrimental to the survival of population, otherwise if the initial population level is small enough, demographic stochasticity and environmental stochasticity may bring survival opportunities for the population that deterministically would extinct indefinitely; (iii) In the extinction case, demographic and environmental stochasticity cannot change the trend of population extinction, but they can delay or promote population extinction.

Plasticity's role in adaptive evolution depends on environmental change components.

To forecast extinction risks of natural populations under climate change and direct human impacts, an integrative understanding of both phenotypic plasticity and adaptive evolution is essential. To date, the evidence for whether, when, and how much plasticity facilitates adaptive responses in changing environments is contradictory. We argue that explicitly considering three key environmental change components - rate of change, variance, and temporal autocorrelation - affords a unifying framework of the impact of plasticity on adaptive evolution. These environmental components each distinctively effect evolutionary and ecological processes underpinning population viability. Using this framework, we develop expectations regarding the interplay between plasticity and adaptive evolution in natural populations. This framework has the potential to improve predictions of population viability in a changing world.

Temporally auto-correlated predator attacks structure ecological communities.

For species primarily regulated by a common predator, the P* rule of Holt & Lawton (Holt & Lawton, 1993. Am. Nat.142, 623-645. (doi:10.1086/285561)) predicts that the prey species that supports the highest mean predator density (P*) excludes the other prey species. This prediction is re-examined in the presence of temporal fluctuations in the vital rates of the interacting species including predator attack rates. When the fluctuations in predator attack rates are temporally uncorrelated, the P* rule still holds even when the other vital rates are temporally auto-correlated. However, when temporal auto-correlations in attack rates are positive but not too strong, the prey species can coexist due to the emergence of a positive covariance between predator density and prey vulnerability. This coexistence mechanism is similar to the storage effect for species regulated by a common resource. Negative or strongly positive auto-correlations in attack rates generate a negative covariance between predator density and prey vulnerability and a stochastic priority effect can emerge: with non-zero probability either prey species is excluded. These results highlight how temporally auto-correlated species' interaction rates impact the structure and dynamics of ecological communities.

A Trade-Off between Robustness to Environmental Fluctuations and Speed of Evolution.

AbstractUnderstanding how a species' life history affects its capacity to cope with environmental changes is important in the context of rapid climate changes. Reinterpreting previous results from a well-developed theoretical framework, we show that a trade-off exists between a species' ability to genetically adapt to long-term gradual environmental changes and its ability to demographically resist short-term environmental perturbations, causing variation in its vital rates. Surprisingly, this important insight has not been made formally explicit before. Choosing archetypal life histories along the fast-slow pace-of-life continuum and modeling their eco-evolutionary dynamics, we further show that long-lived species have larger demographic robustness to interannual fluctuations but limited trait evolutionary responses in gradually changing environments. In contrast, short-lived species had larger evolvability but reduced demographic robustness. This trade-off bears heavily on extinction probabilities of populations tracking fast trait changes in stochastic environments. Faster trait evolution in short-lived species came at the expense of their higher sensitivity to short-term fluctuations, causing higher extinction rates than for long-lived species. Long-lived species persisted better on short timescales but built maladaptation and an extinction debt over time. This work shows how modeling species' eco-evolutionary dynamics can help to assess species vulnerability to environmental changes.