Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Measuring resilience and assessing vulnerability of terrestrial ecosystems to climate change in South America

  • Luciano J. S. Anjos ,

    Contributed equally to this work with: Luciano J. S. Anjos, Peter Mann de Toledo

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Writing – original draft, Writing – review & editing

    ljsanjos@gmail.com

    Affiliations Campus de Parauapebas, Universidade Federal Rural da Amazônia - UFRA, Parauapebas, Pará, Brazil, Programa de Pós-Graduação em Ciências Ambientais - PPGCA, Instituto de Geociências, Universidade Federal do Pará - UFPA, Belém, Pará, Brazil

  • Peter Mann de Toledo

    Contributed equally to this work with: Luciano J. S. Anjos, Peter Mann de Toledo

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliations Programa de Pós-Graduação em Ciências Ambientais - PPGCA, Instituto de Geociências, Universidade Federal do Pará - UFPA, Belém, Pará, Brazil, Centro de Ciência do Sistema Terrestre - CCST, Instituto Nacional de Pesquisas Espaciais - INPE, São José dos Campos, São Paulo, Brazil

Abstract

Climate change has been identified as the primary threat to the integrity and functioning of ecosystems in this century, although there is still much uncertainty about its effects and the degree of vulnerability for different ecosystems to this threat. Here we propose a new methodological approach capable of measuring and mapping the resilience of terrestrial ecosystems at large scales based on their climatic niche. To do this, we used high spatial resolution remote sensing data and ecological niche modeling techniques to calculate and spatialize the resilience of three stable states of ecosystems in South America: forest, savanna, and grassland. Also, we evaluated the sensitivity of ecosystems to climate stress, the likelihood of exposure to non-analogous climatic conditions, and their respective adaptive capacities in the face of climate change. Our results indicate that forests, the most productive and biodiverse terrestrial ecosystems on the earth, are more vulnerable to climate change than savannas or grasslands. Forests showed less resistance to climate stress and a higher chance of exposure to non-analogous climatic conditions. If this scenario occurs, the forest ecosystems would have less chance of adaptation compared to savannas or grasslands because of their narrow climate niche. Therefore, we can conclude that a possible consolidation of non-analogous climatic conditions would lead to a loss of resilience in the forest ecosystem, significantly increasing the chance of a critical transition event to another stable state with a lower density of vegetation cover (e.g., savanna or grassland).

Introduction

Empirical data show that in 2016 there was a drop in the temperature record [1]. It was also the driest year recorded for South America since observations began in 1900 [2]. Ongoing climate changes have been identified as the main threat to the integrity and functioning of terrestrial ecosystems in the 21st century [3]. In particular, projections for South America indicate that by the end of this century the continent is likely to be subjected to non-analogous climatic conditions, having little overlap with the current climate [4,5]. This represents a high risk to the extraordinary biodiversity on the continent due to the high climatic sensitivity of its ecosystems, which could lead to the compromise of a vast set of goods and services provided to humans [6,7].

Maintaining the structure and functioning of ecosystems over a broad-scale climate gradient range [8] suggests that each stable state of the ecosystem must be strictly adapted and evolved under a specific subset of abiotic conditions, enclosed in the multidimensional space of its climatic niche. Thus, under such a combination of conditions, the ecosystem is expected to have a higher capacity to absorb disturbances and recover the condition closer to the original state in the shortest time after a disturbance. Also, it is important to note that there is a strong correlation between ecosystem resilience and ecosystem structure.

The theory of ecological stability predicts that if changes to environmental conditions occur and a critical threshold of resilience is surpassed, catastrophic transitions between stable ecosystem states can be abruptly precipitated [9]. In this sense, it is crucial to measure and monitor ecosystem resilience, because this attribute plays a crucial role in mediating transitional events between different ecosystems [10]. Recent studies have been devoted to investigating and understanding the mechanisms behind critical transition events in tropical terrestrial ecosystems [1114].

Despite recent theoretical and methodological advances, resilience is a concept that remains qualitative [15], and quantification at large scales can be challenging [16]. To overcome this limitation, we propose the methodological coupling of two significant bodies of theoretical knowledge in ecology: ecological stability theory [17] and ecological niche theory [1820]. Some previous studies have already proposed measuring ecosystem resilience at large scales by calculating the likelihood of finding an ecosystem at a given level of precipitation [8,21,22]. However, none of the authors considered the robust theoretical framework and the recent methodological advances of ecological niche modeling [23,24]. With this approach, it is possible to construct spatially explicit models that measure the resilience of ecosystems through a metric of climate suitability, based on the multidimensional niche preferably occupied by them.

Overall, there are still many uncertainties about what effects climate change will have on the resilience and response of terrestrial ecosystems in structural and functional terms [25,26]. In this sense, it is fundamental to evaluate the intrinsic vulnerability of terrestrial ecosystems through a decomposition into three different axes: sensitivity, exposure, and adaptive capacity of the ecosystem [27]. Here, we propose to objectively estimate the comparative sensitivities of terrestrial ecosystems in South America from an analysis of ecosystem resistance to a gradual increase in climatic stress. Also, we will briefly assess the likelihood of exposure to unfavorable climatic conditions by identifying the propensity of each stable climate state to cope with non-analogous climate changes within or outside its optimal climatic niche, taking into account observed and simulated climate trends. Finally, to infer the adaptive capacity of ecosystems to changes in climate, we will assess the amount of climate variability that the ecosystem can tolerate without losing its structure and function.

Materials and methods

Study area

South America is located in the Neotropical biogeographic domain and has an area of almost 18 million km2 distributed among 12 countries. The continent has an extensive environmental heterogeneity, encompassing over 100 different ecoregions [28]. Among these are some of the most productive and megadiverse ecosystems on the planet, including forests (e.g., Amazon rainforest and Atlantic forest), savanna ecosystems (e.g., Cerrado, Beni, Gran-savanna), seasonally dry tropical forests (e.g., Caatinga, Chiquitano), and even deserts (e.g., Atacama).

The high variability can partially explain this diversity of ecosystems in latitudinal and altitudinal gradients, which directly influence the patterns of moisture and energy availability throughout the continent. Based on the climatic dataset of this study, we observed that average annual temperatures could vary between negative values at high altitudes (e.g., the Andes Mountain Range) to 30 °C in the low-altitude equatorial tropical zone. The annual range of temperature varies mainly as a function of latitude, being low near the equator, and can reach 18 °C in high latitudes. In some regions of the tropical zone, it can rain more than 6000 mm·year-1 with low seasonality throughout the year (e.g., Amazonian West), while other areas are exceptionally arid, with a marked dry season throughout the year and less than 100 mm·year-1 of precipitation.

Vegetational cover data and determination of stable states

To define the terrestrial ecosystems of South America as stable states, we used the variable called tree cover for the year 2001, from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor. This variable describes the percentage of vegetation cover varying between zero and 100% [29] and can be interpreted as the density or abundance of trees. The transition limits between ecosystem’s stable states were inferred from the frequency distribution of the tree cover [8], which presented a trimodal distribution. Each distribution mode comprises a stable state of a terrestrial ecological system. The grasslands had zero to 5% tree cover; savannas had 5 to 60%, and forests had values above 60%. The original dataset, spanning all of South America, was resampled from a 500-m to a 6-km spatial resolution.

Due to the effects of historical changes in land use, mistaken pixel classifications occurred among the three ecosystem classes. Notably, the forest ecosystem (with tree cover between 60 and 100%) was the most sensitive to spectral and radiometric changes following the intrusion of alternative land uses (e.g., selective logging, monocultures, or pastures). To minimize such classification bias, and thus to match the original distribution of the forest ecosystem as closely as possible, we evaluated the classification accuracy by comparing the tree cover with a census database of high-resolution land cover [30]. Whenever a pixel was classified with a tree cover percentage of less than 60% but classified as forest based on the consensus dataset, the pixel was reclassified to forest.

Presence-absence data for ecological niche modeling

After correction of land use classification inaccuracies, the raster database was converted to a presence-absence binary, where each map pixel was converted to a point in vector format. In total, n = 37763 samples were generated for South America, of which 53% were classified as savanna; 38% as forest; and 9% as grassland. These ecosystems are mutually exclusive and geographic substitutes, so when a particular type of ecosystem was recorded as present (1) for a given area, the others were automatically recorded as absent (0) in that same area. These presence-absence data were used as an input in the ecological niche modeling procedure described below. Each point was assigned a geo-referenced signature with longitude and latitude information.

Climatic variables and selection of bioclimatic predictors

For this work, the climatic conditions of South America were defined along two principal axes: availability of moisture (1) and energy (2). At the scale used for this study, these are considered the primary predictors of climate response patterns in terrestrial ecosystems [31,32].

To describe the pattern of moisture availability, we used current precipitation data for South America from the climatic dataset CHPclim (v.1.0), produced by the Climate Hazards Group’s Precipitation Climatology. It is a combination of satellite data, physiographic indicators, and standard in situ climatological data, with data available from 1980 to 2009. The final product is a monthly global climatic precipitation measure with a spatial resolution of 0.05° (~6 km) [33]. This database was copied from the address 'http://chg.geog.ucsb.edu/data/CHPclim/' in ‘tiff’ format and then cut to the region of interest. CHPclim has a critical comparative advantage over other precipitation climatic bases available for South America (e.g., WorldClim): It has a better fit for regions with low densities of meteorological stations and high variability of precipitation, such as the Amazon, and responds satisfactorily to complex terrain, such as in the Andes Mountain Range [33].

To describe energy availability across the continent, we used the WorldClim temperature database [34], which consists of a climatological average between the years 1960 and 1990, with data from stations around the globe. The spatial resolution of the dataset is 0.041° (~ 5 km) and was copied from 'http://www.worldclim.org/version1' in ‘tiff’ format. Although it has the same limitation of a low density of meteorological stations for the Neotropical region, the variable in question does not present significant variability near the equator, so that the quality of the observed data is not compromised. Finally, we selected four bioclimatic predictors related to energy availability and moisture, which originate from the ecophysiological point of view of ecosystems: (1) annual cumulative precipitation (ACP), (2) precipitation seasonality coefficient (PSC), (3) average annual temperature (AAT), and (4) annual range of temperature (ART).

Modeling the climatic niche of ecosystems

Ecological niche modeling finds strong support in ecological niche theory [35], which predicts that each organism has a multidimensional environmental space with optimal conditions for its survival, growth, and reproduction [3639]. The various modeling methods quantify the relationships between occurrences (presence/absence) and environmental predictors to delimit and adjust the multidimensional environmental niche. The result is the quantification and spatialization of the suitability of the environment for the modeled ecosystem.

To model the distribution of ecosystems, we used the biomod2 package implemented in R software [40]. Distribution models were calibrated using presence-absence data from each ecosystem and climate predictors for the South American continent. We have adopted the ensemble strategy that emphasizes the most consensus in predictions among different modeling methods [23,41], thus minimizing the effect of uncertainties on model prediction [42]. The models were run using 10 different methods: Bioclim (SRE), Classification Tree Analysis (CTA) [43], Maxent [44,45], Random Forest (RF) [46], Generalized Linear Models (GLM) [47], Generalized Additive Models (GAM) [48], Generalized Boosted Regression Models (GBM) [49], Function Discriminant Analysis (FDA) [50], Artificial Neural Networks (ANN) [50], and Multiple Additive Regression Splines (MARS) [51].

For each method, ten replicates with 75% and 25% partitions were run for training and testing, respectively. Quality of the models produced by the different methods was evaluated using True Skill Statistics (TSS) and Receiver Operating Characteristic (ROC) metrics (S1 Table). We calculated the contribution of each bioclimatic predictor to build and explain the patterns of response variable which varies between zero (lower importance) to 1 (higher importance) (S1 Table). The models selected to compose the ensemble were those with the best TSS, which measures combined sensitivity and specificity of the model [52]. For the threshold effect, only the models with TSS values equal to or greater than 0.7 were considered for inclusion in the ensemble. The consensus distribution model was then obtained as the arithmetic average taken for the selected models [53].

Measuring and mapping the resilience of terrestrial ecosystems

Here we present an alternative interpretation for the climate suitability map produced from ecological niche modeling. For the first time, a methodological convergence between predictions of ecological stability and ecological niche theory is explicitly suggested. The phenotypic adaptation of an ecosystem to a given set of climatic conditions has a direct and robust relationship with the resilience of the ecosystem itself. In this sense, the ecosystem is expected to lose resilience as it moves away from the optimal climatic conditions of its niche. Therefore, the product of niche modeling can be interpreted as a continuous and objective measure of resilience, ranging from zero (minimum resilience) to 1000 (maximum resilience) for terrestrial ecosystems in South America.

Sensitivity of terrestrial ecosystems to climate stress

We assume that climate stress increases in ecosystems as they move away from their optimal niche. Since the climate suitability gradient, predicted by niche modeling, can be interpreted as an approximation of the fundamental ecosystem niche, we estimate climate stress from a simple inversion of the gradient as predicted by the models. Thus, a band with low climatic suitability is assigned a high level of climatic stress and vice versa.

Using this climate stress vector (x) as the predictor of the relative frequency of observations for each ecosystem (y), we fit an exponential model according to the equation below: (1) where a and b are the estimated coefficients of the exponential model. The objective was to use the coefficient b, which defines the curve fit of the exponential model to the observed data, as an indicator of the ecosystem resistance to climatic stress. In other words, if a small increment of climatic stress significant changes the ecosystem, we can assume that the ecosystem has low resistance to climatic stress, and consequently a higher value of the coefficient b. Conversely, resistant ecosystems will have a lower value of the coefficient b.

Exposure of ecosystems to non-analogous climatic conditions and their respective adaptive capacities

In a two-dimensional climatic space representing the availability of moisture and energy present in South America, we evaluated the propensity of ecosystems for exposure to non-analogous climatic conditions and their respective adaptive capacities based on the climate space occupied by each ecosystem and its respective observed resilience gradient. We have taken as a reference the observed current climate trends [2,54,55] and simulations of future climate scenarios [5,56,57] to assess the likely progression of climate changes in South America. In general, both observed and simulated data point to an increase in climatic seasonality associated with increased aridity and temperatures for South America. We, therefore, delimited a polygon for the bidimensional climate space to indicate the likely orientation climate changes on a continental scale. Observed ecosystem resilience data was then compared with the polygon to determine which ecosystems would be most abruptly exposed to unfavorable climatic conditions. In the same sense, it was also possible to observe the response width of ecosystem resilience and, in this way, infer the capacity to adapt to a new climatic reality.

Results

Mapping ecosystem resilience

Here we present a continuous metric of spatially explicit ecosystem resilience, based on the climatic niche to which each ecosystem is adapted (Fig 1). Grasslands have high resilience at high altitudes and arid environments, such as in the Andes, the Atacama Desert, and Patagonia. Savannas are widely distributed throughout the continent and exhibit high resilience along the continental diagonal polygon as well as to the north of South America. Forests have good resilience near the equator, covering almost entirely the Amazon basin to the foothills of the Andes. Forests also occupy a narrow strip following the Atlantic coast (Atlantic Forest) and the eastern region lying near latitude 50° (Araucaria forest). All of these areas show a lower level of resilience than do the forests of the Amazon basin. In the southwestern portion of the continent, the model predicted the presence of a vertical band of temperate forest (Valdivia).

thumbnail
Fig 1. Ecosystem resilience based on the climatic niche projected under geographic space.

A Grasslands; B Savannas; C Forest. High values (~1000) indicate a high recovery capacity after a disturbance, whereas low values (~zero) mean slower recovery capacity after a disturbance.

https://doi.org/10.1371/journal.pone.0194654.g001

There is a spatially structured pattern of loss of resilience toward the transition zones between the three ecosystems. Models based on the climatic niche of ecosystems were highly sensitive in detecting that ecotone zones tend to have a lower resilience value than do areas in an optimal climate niche. This pattern reinforces the idea of geographic ecosystem replacement caused by transitional events through the loss of resilience through geographic space and climatic gradient.

Sensitivity of terrestrial ecosystems to climate stress gradient

In general, the three ecosystems presented similar statistical behaviors, with a decrease in their relative frequencies accompanying a gradual increase in climatic stress (Fig 2). However, there are significant differences concerning the forms of the response curves observed for the three ecosystems. For savannas and grasslands, the distribution of observations was more uniform, while forests account for almost 80% of the observations in the range of the lowest climatic stress (close to the climatic optimum). Because of this, savannas and grasslands have a smoother response curve as climate stress increases than do forests. This pattern was also reflected in the value of the exponential model parameter, with forests (b = -0.024) having a value almost double those for grasslands (b = -0.010) and savannas (b = -0.007), indicating that savannas and grasslands have better resistance to climatic stress than do forests.

thumbnail
Fig 2. Model of ecosystem responses to a climate stress gradient based on the prediction of relative frequencies (%).

The parameters of the exponential model indicate that there are differences regarding resistance to climatic stress among the three ecosystems. Grassland n = 598; Savanna n = 889; Forest n = 891.

https://doi.org/10.1371/journal.pone.0194654.g002

Propensity for exposure to unfavorable climatic conditions

The moisture availability gradient showed a strong interaction with climate resilience. Forests present high resilience to rainfall above 2000 mm·year-1 with a relatively low seasonality (<100%) during the year (Fig 3). Savannas are adapted to a precipitation range between 500 and 1800 mm·year-1, so that resilience tends to decrease with increasing seasonality of precipitation and annual rainfall volumes. In contrast, grasslands show high resilience in the gradient range where a more arid climate prevails, with low rainfall volumes and high rainfall seasonality throughout the year. In this sense, if the observed trends are assessed alongside the simulated future climate scenarios, both savannas and grasslands would be favored, while forests would inevitably lose resilience due to moisture reduction.

thumbnail
Fig 3. Observed patterns between climatic variation and ecosystem resilience.

In the two-dimensional gradient of moisture availability (A) and energy availability gradient (B), if ecosystems are exposed to non-analogous climatic conditions (gray polygon), it is possible to assess what the impacts would be and whether they are likely to adapt to climate change. Forest observations: n = 14,466 samples—38.30%; Savana: n = 19,923 samples- 52.75%; Grassland: n = 3,374 samples—8.93%.

https://doi.org/10.1371/journal.pone.0194654.g003

For the energy availability gradient, forests are adapted at high levels of resilience to a narrow range of high average annual temperatures (between 20–28 °C), but tolerate only low variability throughout the year (between 0–5 °C·year-1). Outside of this range, there was a drastic loss in ecosystem resilience. On the other hand, savannas have high resilience at average annual temperatures above 10 °C, also supporting temperatures above 25 °C, while their resilience tends to decrease at temperatures below 10 °C. Savannas also persist in places with considerable variation in temperature throughout the year, with ranges up to 20 °C. Grasslands are adapted for temperatures in a wide range, from negative to 20 °C average annual temperatures, and with a large annual thermal amplitude (5 to 20 °C·year-1). The observed and simulated average temperature trends indicate abrupt increases, as delimited by the gray polygon, which would favor forest ecosystems, provided there is no increase in the annual temperature range. If there is an increase in annual temperature variability, the scenario would be more favorable to savannas at the expense of forests, which would probably lose resilience.

Adaptive capacity of terrestrial ecosystems

Forests and savannas have strong niche-specific adaptive responses to moisture availability persisting over a specific range of the moisture gradient (Fig 3). This suggests a low adaptive capacity to new conditions outside the optimum climatic niche for forest and savanna ecosystems. On the other hand, grasslands show high levels of resilience to a wide range of moisture availability. The pattern is further enhanced by the seasonal axis of precipitation, with grasslands maintaining high levels of resilience for a wide range of moisture availability. The most significant adaptive constraint regarding grassland resilience is regions with rainfall volumes over 2000 mm·year-1.

Savannas and grasslands would have higher adaptive capacities to changing temperature gradients because they show high resilience at a wide range of temperature combinations. The forest ecosystem, however, has a relatively narrow climatic niche and demonstrates high resilience only for high temperatures with low annual variability. This suggests that savannas and grasslands have a higher adaptive capacity to changes in the energy availability than do forests, which have a narrow climatic niche.

Discussion

In the context of robust transformations induced by climate change, the resilience estimated from the climatic niche emerges as an alternative to the traditional metrics of biological diversity (such as richness or species composition). This new metric can be interpreted as a measure of ecosystem health [58], and may be useful in the design and implementation of conservation strategies at different spatial scales. For example, ecosystem resilience could be monitored in real time to predict catastrophic transitional events between stable states as a function of observed or simulated climate changes [13]. This method is already available, since remote sensing data are widely available at high spatial, spectral, and temporal resolutions [59].

Another potential application would be in actions to support the restoration of degraded ecosystems, considering that the success of an ecological restoration process is directly linked to the ecosystem resilience of the area to be restored. In other words, recovering a forest in an area with low resilience (for this type of ecosystem) would be more costly regarding resources and time, if not infeasible, which could jeopardize successful project implementation.

Our results indicate that the terrestrial ecosystems of South America are sensitive to the expected gradual increases in climatic stress, since the relative frequency of ecosystems decrease as conditions become marginal for the climate niche of the ecosystem in question. However, there are differences in the form of responses between different vegetation classes. When the exponential model parameters were compared between ecosystems, forests presented the least resistance to climate stress, indicating a higher sensitivity intrinsic to changes in climate. These results suggest that forests are more sensitive to climate variability [11,60], mainly due to decreasing availability of moisture [6163]. This result highlights the fragility of these ecosystems to climate change since vulnerability is positively correlated with sensitivity [27].

The response curve for the stable state of forests also indicates that an incremental change in the climatic stress gradient could trigger large-scale transformations, particularly if a rupture threshold is reached [64,65]. This factor may be critical in the short term, with the accelerated levels of transformation currently being experienced, so it is critical to the identify thresholds for rapid forest decline. Even so, it may take decades for forests to restore the services they provide [58]. In contrast, a larger amount of climate stress would be required for savannas or grasslands to become more susceptible to large-scale transformations. It is possible to explain their higher tolerance to climatic stress due to better adaptation to water deficit and the fire effect on biota [66,67].

However, it is important to note that when forests are subjected to some degree of climatic stress, particularly to moisture deficit, their capacity to retain CO2 in biomass is compromised [6871]. The reduction of moisture availability throughout the year affects the floristic composition of forest species [72] and is also associated with high tree mortality rates [73,74]. Hydric stress, combined with factors such as fire, can lead to species loss and biodiversity erosion [75,76]. Furthermore, changes in land use at large scales, e.g., deforestation or selective logging, have affected biophysical climatic parameters such as air temperature [77], as well as the precipitation regime [78,79]. Since temperature is an essential predictor of the maintenance of structure and function of terrestrial ecosystems, these bioclimatic temperature variations are expected to directly and negatively affect forest resilience.

Our results indicate that forests in South America are more likely to face unfavorable climatic conditions in the near future than are savannas and grasslands, which have a lower percentage of vegetation cover. Forest ecosystems display greater vulnerability due to their higher probability of exposure to non-analogous climates [27]. Empirical evidence shows that forests in border areas of the Amazon southwest have presented slower biomass recovery after extreme drought events [80,81], suggesting that they have recently declined in resilience.

On the other hand, these new climatic conditions are likely to favor other ecosystem types that display increasing resilience. In this way, the dynamics of resilience within an ecological system with multistability are closely related to transition events between stable states and depend directly on the orientation of change, thus favoring a particular stable state. For example, at the end of the Upper Miocene, there was a reduction in global temperature, and grassland ecosystems dominated most of South America [82]. However, in the Eocene, with higher temperatures and precipitation, forests thrived to the detriment of grasslands and savannas [83].

Tropical forests are known to have a narrow thermal tolerance [84,85], evidence of their higher vulnerability to climate change, representing a low adaptive capacity compared to other ecosystems that can maintain high resilience values under a wide range of climatic conditions. Moreover, in theory, for an adaptation event to occur, the adaptive strategy should be within the limits of the diversity of responses, i.e., phenotypic plasticity, displayed in the current climate. Outside the boundaries of the climatic niche, there would be slow adaptation, and the ecosystem would succumb to a loss of structure and function until it reached a point of rupture or another stable state.

Our results indicate that, in this sense, savannas demonstrate an advantage over forests. Savannas are successful in a wide range of moisture and temperature availability and tolerate strong seasonal variations. Savannas can, therefore, maintain a steady tree cover state despite climatic variation, which confers a tremendous adaptive potential to the ecosystem if it is exposed to non-analogous climatic conditions.

There are no predicted energy restrictions in the near future, since all simulated projections point to an increase in temperature, regardless of the emissions scenario [86]. An increase in atmospheric CO2 concentration could favor forests by increasing the availability of CO2 as a resource [87,88]. However, the most critical factor for forest persistence is the widespread need for moisture [62,89,90]. In contrast to the hypothesis that forests are resilient to climate change [87], our results suggest that, due to the close adaptive relationship of ecophysiological and evolutionary patterns between terrestrial ecosystems and climate, the high values of resiliency to changing climatic conditions can only be maintained at the conditions of the optimum climatic niche occupied by the ecosystem. This dependence would thus be reflected in ecosystem responses following a disturbance event [91].

Conclusions

Our models, with strong support in niche theory and ecological stability, have shown a significant sensitivity and may be useful in several practical applications within the science of conservation, including in assessing the vulnerability of ecosystems to climate change. We found some worrying evidence that South American forest ecosystems are intrinsically more vulnerable to climate change than other ecosystems. Ongoing climate change can accelerate the loss of ecosystem resilience by promoting erosion of forest biodiversity and leading to another stable state with a lower density of vegetation cover.

Supporting information

S1 Table. Importance metrics of the predictors to build the ecological niche models and evaluation metrics by different methods for each ecosystem.

Average values were calculated utilizing 10 replicates, as well as the evaluation metrics.

https://doi.org/10.1371/journal.pone.0194654.s001

(PDF)

Acknowledgments

We are grateful to Ima Vieira, Aline Lima, Roberta Cerqueira, Tereza Giannini, Alexandre Aleixo, Everaldo Souza and Edson da Rocha for their valuable contributions and suggestions throughout the development of this work.

References

  1. 1. Mann ME, Miller SK, Rahmstorf S, Steinman BA, Tingley M. Record temperature streak bears anthropogenic fingerprint. Geophys Res Lett. 2017; 7936–7944.
  2. 2. Erfanian A, Wang G, Fomenko L. Unprecedented drought over tropical South America in 2016: significantly under-predicted by tropical SST. Sci Rep. Springer US; 2017;7: 5811. pmid:28724896
  3. 3. Bellard C, Bertelsmeier C, Leadley P, Thuiller W, Courchamp F. Impacts of climate change on the future of biodiversity. Ecol Lett. 2012;15: 365–377. pmid:22257223
  4. 4. Williams JW, Jackson ST, Kutzbach JE. Projected distributions of novel and disappearing climates by 2100 AD. 2007;
  5. 5. Garcia RA, Cabeza M, Rahbek C, Araujo MB. Multiple Dimensions of Climate Change and Their Implications for Biodiversity. Science (80-). 2014;344: 1247579–1247579. pmid:24786084
  6. 6. Cardinale BJ, Duffy JE, Gonzalez A, Hooper DU, Perrings C, Venail P, et al. Biodiversity loss and its impact on humanity. Nature. 2012;489: 326–326.
  7. 7. Pecl GT, Araujo MB, Bell JD, Blanchard J, Bonebrake TC, Chen I-C, et al. Biodiversity redistribution under climate change: impacts on ecosystems and human well-being. Science (80-). 2017; in press. pmid:28360268
  8. 8. Hirota M, Holmgren M, Van Nes EH, Scheffer M. Global Resilience of Tropical Forest and Savanna to Critical Transitions. Science (80-). 2011;334: 232–235. pmid:21998390
  9. 9. Scheffer M, Carpenter SR. Catastrophic regime shifts in ecosystems: Linking theory to observation. Trends Ecol Evol. 2003;18: 648–656.
  10. 10. Gunderson L. Ecological resilience in theory and application. Annu Rev Ecol Syst. 2000;31: 425–439.
  11. 11. Holmgren M, Hirota M, van Nes EH, Scheffer M. Effects of interannual climate variability on tropical tree cover. Nat Clim Chang. Nature Publishing Group; 2013;3: 755–758.
  12. 12. Dantas V de L, Hirota M, Oliveira RS, Pausas JG, Dantas DL, Dantas V de L, et al. Disturbance maintains alternative biome states. Ecol Lett. 2016;19: 12–19. pmid:26493189
  13. 13. Verbesselt J, Umlauf N, Hirota M, Holmgren M, Van Nes EH, Herold M, et al. Remotely sensed resilience of tropical forests. Nat Clim Chang. 2016;1: 1–5.
  14. 14. Flores BM, Holmgren M, Xu C, van Nes EH, Jakovac CC, Mesquita RCG, et al. Floodplains as an Achilles’ heel of Amazonian forest resilience. Proc Natl Acad Sci. 2017; 201617988. pmid:28396440
  15. 15. Scheffer M, Carpenter SR, Dakos V, Van Nes EH. Generic Indicators of Ecological Resilience: Inferring the Chance of a Critical Transition. Annu Rev Ecol Evol Syst. 2015;
  16. 16. Reyer CPO, Brouwers N, Rammig A, Brook BW, Epila J, Grant RF, et al. Forest resilience and tipping points at different spatio-temporal scales: Approaches and challenges. J Ecol. 2015;103: 5–15.
  17. 17. Holling CS. Resilience and Stability of Ecological Systems. Annu Rev Ecol Syst. 1973;4: 1–23.
  18. 18. Elton CS. Animal Ecology [Internet]. Animal ecology. Charles Elton; 1927.
  19. 19. Grinnell J. The Niche-Relationships of the California Thrasher. Auk. 1917;34: 427–433.
  20. 20. Hutchinson GE. Concluding remarks. Cold Spring Harb Symp Quant Biol. 1957;22: 415–427.
  21. 21. Zemp DC, Schleussner CF, Barbosa HMJ, Rammig A. Deforestation effects on Amazon forest resilience. Geophys Res Lett. 2017;44: 6182–6190.
  22. 22. Zemp DC, Schleussner C-F, Barbosa HMJ, Hirota M, Montade V, Sampaio G, et al. Self-amplified Amazon forest loss due to vegetation-atmosphere feedbacks. Nat Commun. Nature Publishing Group; 2017;8: 14681. pmid:28287104
  23. 23. Franklin J. Mapping species distributions: spatial inference and prediction. 2010;
  24. 24. Peterson AT, Soberón J, Pearson RG, Anderson RP, Martínez-Meyer E, Nakamura M, et al. Ecological niches and geographic distributions [Internet]. Choice Reviews Online. 2011.
  25. 25. Urban MC. Accelerating extinction risk from climate change. Foreign Policy. 2015;348.
  26. 26. Urban MC, Bocedi G, Hendry AP, Mihoub J-B, Peer G, Singer A, et al. Improving the forecast for biodiversity under climate change. Science (80-). 2016;353: aad8466–aad8466. pmid:27609898
  27. 27. Dawson TP, Jackson ST, House JI, Prentice IC, Mace GM. Beyond Predictions: Biodiversity Conservation in a Changing Climate. 2011;332.
  28. 28. Olson DM, Dinerstein E, Wikramanayake ED, Burgess ND, Powell GVN, Underwood EC, et al. Terrestrial Ecoregions of the World: A New Map of Life on Earth. Bioscience. 2001;51: 933.
  29. 29. Townshend JR., Hansen MC, Carroll M, DiMiceli C, Sohlberg R, Huang C. Vegetation Continuous Fields MOD44B, 2010 Percent Tree Cover. Collect 5, Univ Maryland, Coll Park Maryl. 2011; 1–12. http://glcf.umiacs.umd.edu/data/vcf/
  30. 30. Tuanmu MN, Jetz W. A global 1-km consensus land-cover product for biodiversity and ecosystem modelling. Glob Ecol Biogeogr. 2014;23: 1031–1045.
  31. 31. Lehmann CER, Anderson TM, Sankaran M, Higgins SI, Archibald S, Hoffmann WA, et al. Savanna Vegetation-Fire-Climate Relationships Differ Among Continents. Science (80-). 2014;343: 548–553. pmid:24482480
  32. 32. Oliveras I, Malhi Y. Many shades of green: The dynamic tropical forest-savanna transitions. Philos Trans Ser B, Biol Sci. 2016;
  33. 33. Funk C, Verdin A, Michaelsen J, Peterson P, Pedreros D, Husak G. A global satellite-assisted precipitation climatology. Earth Syst Sci Data. 2015;7: 1–15.
  34. 34. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25: 1965–1978.
  35. 35. Colwell RK, Rangel TF. Hutchinson’s duality: the once and future niche. Proc Natl Acad Sci U S A. 2009;106 Suppl: 19651–8. pmid:19805163
  36. 36. Guisan A, Zimmerman NE. Predictive habitat distribution models in ecology. Ecol Modell. 2000;135: 147–186.
  37. 37. Guisan A, Thuiller W, Antoine G, Thuiller W. Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005;8: 993–1009.
  38. 38. Hirzel AH, Le Lay G. Habitat suitability modelling and niche theory. J Appl Ecol. 2008;45: 1372–1381.
  39. 39. Soberón J, Nakamura M. Niches and distributional areas: Concepts, methods, and assumptions. Proc Natl Acad Sci. 2009;106: 19644–19650. pmid:19805041
  40. 40. Thuiller W, Lafourcade B, Engler R, Araújo MB, Araujo MB. BIOMOD—a platform for ensemble forecasting of species distributions. Ecography (Cop). 2009;32: 369–373.
  41. 41. Araújo MB, New M. Ensemble forecasting of species distributions. Trends Ecol Evol. 2007;22: 42–47. pmid:17011070
  42. 42. Diniz-Filho JAF, Bini LM, Rangel TF, Loyola RD, Hof C, Nogue D, et al. Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change. Ecography (Cop). 2009;32: 897–906.
  43. 43. Scull P, Franklin J, Chadwick OA. The application of classification tree analysis to soil type prediction in a desert landscape. Ecol Modell. 2005;181: 1–15.
  44. 44. Phillips SJ, Anderson RP, Schapire RE, Phillips ENCANDPR, J S, Anderson RP, et al. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190: 231–259.
  45. 45. Phillips SJ, Anderson RP, Dudík M, Schapire RE, Blair ME. Opening the black box: an open-source release of Maxent. Ecography (Cop). 2017;40: 887–893.
  46. 46. Breiman L. Random forests. Mach Learn. 2001;45: 5–32.
  47. 47. Guisan A, Edwards TC, Hastie T. Generalized linear and generalized additive models in studies of species distributions: setting the scene. Ecol Modell. 2002;157.
  48. 48. Hastie T, Tibshirani R. Generalized Additive Models. Statistics (Ber). 1986;10: 409–435.
  49. 49. Friedman JH. Greedy Function Approximation: a Gradient Boosting Machine. Statistics (Ber). 2001;29: 1189–1232.
  50. 50. Manel S, Dias JM, Ormerod SJ. Comparing discriminant analysis, neural networks and logistic regression for predicting species distributions: A case study with a Himalayan river bird. Ecol Modell. 1999;120: 337–347.
  51. 51. Friedman JH. Multivariate Adaptive Regression Splines. Ann Stat. 1991;19: 1–67.
  52. 52. Allouche O, Tsoar A, Kadmon R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS). J Appl Ecol. 2006;43: 1223–1232.
  53. 53. Diniz-Filho JAF, Ferro VG, Santos T, Nabout JC, Dobrovolski R, De Marco P. The three phases of the ensemble forecasting of niche models: geographic range and shifts in climatically suitable areas of Utetheisa ornatrix. Rev Bras Entomol. 2010;54: 339–349.
  54. 54. Fu R, Yin L, Li W, Arias PA, Dickinson RE, Huang L, et al. Increased dry-season length over southern Amazonia in recent decades and its implication for future climate projection. Proc Natl Acad Sci. 2013;110: 18110–18115. pmid:24145443
  55. 55. Fu R. Global warming-accelerated drying in the tropics. Proc Natl Acad Sci. 2015;112: 201503231. pmid:25770219
  56. 56. Malhi Y, Aragao LEOC, Galbraith D, Huntingford C, Fisher R, Zelazowski P, et al. Exploring the likelihood and mechanism of a climate-change-induced dieback of the Amazon rainforest. Proc Natl Acad Sci. 2009;106: 20610–20615. pmid:19218454
  57. 57. Boisier JP, Ciais P, Ducharne A, Guimberteau M. Projected strengthening of Amazonian dry season by constrained climate model simulations. Nat Clim Chang. 2015;5: 656–660.
  58. 58. Trumbore S, Brando P, Hartmann H. Forest health and global change. Science (80-). 2015;349: 814–818. pmid:26293952
  59. 59. Turner W, Spector S, Gardiner N, Fladeland M, Sterling E, Steininger M. Remote sensing for biodiversity science and conservation. 2003;
  60. 60. Seddon AW, Macias-Fauria M, Long PR, Benz D, Willis KJ. Sensitivity of global terrestrial ecosystems to climate variability. Nature. Nature Publishing Group; 2016;531: 229–232. pmid:26886790
  61. 61. Hilker T, Lyapustin AI, Tucker CJ, Hall FG, Myneni RB, Wang Y, et al. Vegetation dynamics and rainfall sensitivity of the Amazon. Proc Natl Acad Sci. 2014;111: 16041–16046. pmid:25349419
  62. 62. Malhi Y, Roberts JT, Betts RA, Killeen TJ, Li W, Nobre CA. Climate change, deforestation, and the fate of the Amazon. Science. 2008;319: 169–72. pmid:18048654
  63. 63. Seidl R, Thom D, Kautz M, Martin-Benito D, Peltoniemi M, Vacchiano G, et al. Forest disturbances under climate change. Nat Clim Chang. Nature Publishing Group; 2017;7: 395–402. pmid:28861124
  64. 64. Nobre CA, Borma LDS. “Tipping points” for the Amazon forest. Curr Opin Environ Sustain. 2009;1: 28–36.
  65. 65. Salazar LF, Nobre CA. Climate change and thresholds of biome shifts in Amazonia. Geophys Res Lett. 2010;37: 1–5.
  66. 66. Bond WJ, Keeley JE. Fire as a global “herbivore”: the ecology and evolution of flammable ecosystems. Trends Ecol Evol. 2005;20: 387–94. pmid:16701401
  67. 67. Jolly WM, Cochrane MA, Freeborn PH, Holden ZA, Brown TJ, Williamson GJ, et al. Climate-induced variations in global wildfire danger from 1979 to 2013. Nat Commun. Nature Publishing Group; 2015;6: 1–11. pmid:26172867
  68. 68. Baccini A, Walker W, Carvahlo L, Farina M, Sulla-Menashe D, Houghton R. Tropical forests are a net carbon source based on new measurements of gain and loss. Rev. 2017;5962: 1–11.
  69. 69. Feldpausch TR, Phillips OL, Brienen RJWW, Gloor E, Lloyd J, Malhi Y, et al. Amazon forest response to repeated droughts. Global Biogeochem Cycles. 2016;30: 964–982.
  70. 70. Brienen RJW, Phillips OL, Feldpausch TR, Gloor E, Baker TR, Lloyd J, et al. Long-term decline of the Amazon carbon sink. Nature. Nature Publishing Group; 2015;519: 344–8. pmid:25788097
  71. 71. Phillips OL, Aragão LEOC, Lewis SL, Fisher JB, Lloyd J, López-González G, et al. Drought Sensitivity of the Amazon Rainforest. Science (80-). 2009;323: 1344–1347. pmid:19265020
  72. 72. Esquivel-Muelbert A, Galbraith D, Dexter KG, Baker TR, Lewis SL, Meir P, et al. Biogeographic distributions of neotropical trees reflect their directly measured drought tolerances. Sci Rep. Springer US; 2017;7: 8334. pmid:28827613
  73. 73. Greenwood S, Ruiz-Benito P, Martínez-Vilalta J, Lloret F, Kitzberger T, Allen CD, et al. Tree mortality across biomes is promoted by drought intensity, lower wood density and higher specific leaf area. Ecol Lett. 2017; pmid:28220612
  74. 74. Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M, et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For Ecol Manage. 2010;259: 660–684.
  75. 75. Brando PM, Balch JK, Nepstad DC, Morton DC, Putz FE, Coe MT, et al. Abrupt increases in Amazonian tree mortality due to drought-fire interactions. Proc Natl Acad Sci U S A. 2014;111: 6347–52. pmid:24733937
  76. 76. Barlow J, Peres CA. Fire-mediated dieback and compositional cascade in an Amazonian forest. Philos Trans R Soc Lond B Biol Sci. 2008;363: 1787–1794. pmid:18267911
  77. 77. Alkama R, Cescatti A. Biophysical climate impacts of recent changes in global forest cover. Science (80-). 2016;351: 600–604.
  78. 78. Khanna J, Medvigy D, Fueglistaler S, Walko R. Regional dry-season climate changes due to three decades of Amazonian deforestation. 2017;
  79. 79. Chambers JQ, Artaxo P. Biosphere–atmosphere interactions: Deforestation size influences rainfall. Nat Publ Gr. Nature Publishing Group; 2017;
  80. 80. Morton DC, Le Page Y, DeFries R, Collatz GJ, Hurtt GC. Understorey fire frequency and the fate of burned forests in southern Amazonia. Philos Trans R Soc B Biol Sci. 2013;368: 20120163–20120163. pmid:23610169
  81. 81. Saatchi S, Asefi-Najafabady S, Malhi Y, Aragao LEOC, Anderson LO, Myneni RB, et al. Persistent effects of a severe drought on Amazonian forest canopy. Proc Natl Acad Sci. 2013;110: 565–570. pmid:23267086
  82. 82. Herbert TD, Lawrence KT, Tzanova A, Peterson LC, Caballero-Gill R, Kelly CS. Late Miocene global cooling and the rise of modern ecosystems. Nat Geosci. 2016;9.
  83. 83. Fine PVA, Ree RH. Evidence for a Time-Integrated Species-Area Effect on the Latitudinal Gradient in Tree Diversity. Am Nat. 2006;168: 796–804. pmid:17109321
  84. 84. Perez TM, Stroud JT, Feeley KJ. Thermal trouble in the tropics. Science (80-). 2016;351: 1392–1393. pmid:27013713
  85. 85. Janzen DH. Why Mountain Passes are Higher in the Tropics. The American Naturalist. 1967. p. 233.
  86. 86. IPCC. Climate Change 2014 Synthesis Report. Contrib Work Groups I, II III to Fifth Assess Rep Intergov Panel Clim Chang. 2014; 1–151.
  87. 87. Huntingford C, Zelazowski P, Galbraith D, Mercado LM, Sitch S, Fisher R, et al. Simulated resilience of tropical rainforests to CO2-induced climate change. Nat Geosci. Nature Publishing Group; 2013;6: 268–273.
  88. 88. Lloyd J, Farquhar GD. Effects of rising temperatures and [CO2] on the physiology of tropical forest trees. Philos Trans R Soc B Biol Sci. 2008;363: 1811–1817. pmid:18267901
  89. 89. Betts R, Sanderson M, Woodward S. Effects of large-scale Amazon forest degradation on climate and air quality through fluxes of carbon dioxide, water, energy, mineral dust and isoprene. Philos Trans R Soc B Biol Sci. 2008;363: 1873–1880. pmid:18267906
  90. 90. Davidson EA, de Araújo AC, Artaxo P, Balch JK, Brown IF, C Bustamante MM, et al. The Amazon basin in transition. Nature. Nature Publishing Group; 2012;481: 321–328. pmid:22258611
  91. 91. Poorter L, Bongers F, Aide TM, Zambrano AMA, Balvanera P, Becknell JM, et al. Biomass resilience of Neotropical secondary forests. Nature. Nature Publishing Group; 2016; 1–15. pmid:26840632