Journal Information
Vol. 22. Issue 3.
Pages 240-249 (July - September 2024)
Share
Share
Download PDF
More article options
Visits
551
Vol. 22. Issue 3.
Pages 240-249 (July - September 2024)
Essays and Perspectives
Full text access
Climate influence on future suitability of high-altitude wetlands in two natural protected areas from the Central Andes of Argentina
Visits
551
Bárbara Ventoa,
Corresponding author
bvento@mendoza-conicet.gov.ar

Corresponding author.
, Juan Riverab, Marcela Ontiveroa
a Instituto Argentino de Investigaciones en Zonas Áridas (CCT-Mendoza/CONICET), Av. Ruiz Leal s/n, Parque General San Martín, Mendoza 5500, Argentina
b Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales (CCT-Mendoza/CONICET), Av. Ruiz Leal s/n, Parque General San Martín, Mendoza 5500, Argentina
Highlights

  • High-altitude wetlands suitability will be reduced under future climate conditions.

  • Projected changes in temperature and precipitation will affect future distribution of wetlands.

  • Greater reductions in wetlands suitability are projected for elevations lower than 4,000 m.

  • Increases in wetlands suitability are projected for elevations higher than 4,200 m.

This item has received
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (7)
Show moreShow less
Tables (2)
Table 1. Biolimatic variables used in the modeling analysis.
Table 2. Climate variable contribution and optimal range values for the five uncorrelated and most relevant variables.
Show moreShow less
Abstract

Climate conditions have a strong influence on the distribution of many natural communities; thus, the influence of global climate change may alter ecosystems. High-altitude wetlands in the Central Andes of Argentina (CAA) provide relevant ecosystem benefits and promote human activities. However, a degradation of Andean wetlands has been observed and documented in the last decades. Projecting the impacts of climate change on future distribution of wetlands is an essential subject in ecological studies. In this work, the present and future suitability for wetland systems in a pristine section of the CAA using spatial distribution modeling under low and high-emission scenarios are explored. The studied wetlands are strongly driven by bioclimatic variables such as mean annual temperature, precipitation, and its seasonality. Projections show that most of the currently occupied areas will modify under future climate conditions. Changes in temperature and precipitation patterns will decrease the potential suitability in low elevation areas for the next decades for the species inhabiting there, especially south of 29 °S. Additionally, future warmer climatic conditions, greater temperature variability, and reduction in precipitation would probably affect the snow cover and the available water supply which are key limiting factors for the distribution of Andean wetlands. This research is a contribution to understanding possible effects of climate change on high-altitude ecosystems. Mitigation measures for conservation of wetlands in CAA are immediately required to compensate for the impact of climate change under future environmental conditions.

Keywords:
Andes
Bioclimatic variables
Climate change
Future projections
Suitability
Wetlands
Full Text
Introduction

High-altitude wetlands are one of the most important productive mountain ecosystems worldwide, being hotspots of biodiversity and relevant for the regional water cycle, and provide fundamental ecosystem benefits such as freshwater for many human settlements, livestock, and refuge for plants and animals inhabiting there. They are spread in mountain areas around the world and possess an especial ecological relevance due to their species richness and endemism (Mitsch et al., 2013; Rubio et al., 2017).

In the Central Andes of Argentina (CAA), wetlands are distributed at high altitudes and are characterized by extreme low temperatures, large thermal amplitude, snow cover during several months of the year, frozen soils (permafrost), and frequent strong winds (Martínez-Carretero and Ontivero, 2016; Mitsch et al., 2013). Andean wetlands are supplied by a continuous water flow that offers optimal conditions for associated plant communities (Carilla et al., 2018; Izquierdo et al., 2022; Ruthsatz et al., 2020). They provide numerous local benefits such as flooding control, water supply and regulation, aquifer recharge, nutrient recycling, erosion control, terrestrial carbon storage, fixation, and soil formation (Mitsch et al., 2013). Moreover, they are important for regional climate regulation with a strong influence in the water cycle and rain production (Bullock and Acreman, 2003). The presence and ecological importance of wetlands is greater in places where water resources are scarce as the case of arid and semi-arid ecosystems in southern South America (Martínez-Carretero and Ontivero, 2016). The geographic position of high-altitude wetlands makes them highly vulnerable to changes in temperature and precipitation patterns (Cuesta et al., 2023; Ontivero, 2015).

Climate change is the major threat to mountain ecosystems at a global scale (Huggel et al., 2019; Pauli and Halloy, 2019) and is expected to become the main driver of biodiversity and ecosystem degradation in the next decades (Malhi et al., 2020). Currently, the hydrological cycle in the Andean region is seriously affected by the increase of temperature causing glacier retreat and alterations in the precipitation patterns, thus placing the Andean wetlands at a high vulnerability level (Bradley et al., 2006; Hugonnet et al., 2021; Rivera et al., 2021). High-altitude wetlands are strongly sensitive to global warming because they are mainly handled by low temperatures throughout the year. The degradation of these fragile ecosystems has been documented at a global scale (Dangles et al., 2017; Durigan et al., 2022). Thus, they are seriously threatened all over the world and a call to avoid a declining availability of fresh water has been made (Finlayson et al., 2019). A modification and loss of these areas is projected for the next decades with a consequent negative impact in natural habitats and a high risk of disappearance of many endemic species (Cuesta et al., 2023; Sala et al., 2000). Studies affirm that the mean temperature in the CCA will increase significantly under future climate conditions, particularly under high emissions scenarios (Rivera et al., 2020). Therefore, it is expected that climate changes may produce serious structural and functional modifications of wetland dynamics with the limitation in the provision of ecosystem services (Cuesta et al., 2023; Salimi et al., 2021). In the CAA, the increase in the aridization processes of high-altitude wetlands in natural protected areas was highlighted (Ontivero, 2015; Martínez-Carretero and Ontivero 2016). Species divergent responses and diversity losses in Andean biomes have been projected under future climate scenarios (Tovar et al., 2022). However, the quantification of impacts and potential distribution of high-altitude wetlands under future climate conditions are still unclear. For this reason, awareness of their present suitable distribution areas in the CAA related to climate variables is essential for understanding the future perspectives of these delicate mountain ecosystems (Carilla et al., 2018; Dangles et al., 2017). Knowing the current distribution area occupied by high-altitude wetlands in the CAA is relevant for comprehending the structure and functioning of these ecosystems under global climate change. Here, the principal hypothesis deals with a change of suitable areas for high-altitude wetlands in CAA under future climate conditions. Therefore, in this research, the present and future distribution of high-altitude wetlands in two natural protected areas from the CAA were explored. Projected changes under future climate conditions using two greenhouse gasses emissions scenarios were analyzed, and the effects of prospective impacts on their structure and functioning were discussed.

The analysis of climate projections and the response of biological processes help to understand future impacts on natural systems (Parmesan et al., 2000). Moreover, the knowledge about the potential and future distribution of Andean wetlands will contribute to the implementation of conservation practices, placing information available for the decision makers (Crego et al., 2014; Ontivero et al., 2022). Thus, the effectiveness of management practices demands a better knowledge of how these ecosystems will change and adapt under future variations in regional temperature and precipitation patterns.

MethodsStudy area

High-altitude wetlands in the CAA are mainly supplied by groundwater along the Puna and the Altoandina phytogeographic provinces (Martínez-Carretero and Ontivero, 2016). The study area comprises two adjacent natural protected areas in the CAA (Fig. 1). The San Guillermo protected area (28º 27' S, 69º 32' W) is located in the northwest of San Juan province, at an altitude of 1,700–5,000 m. On the other hand, the Laguna Brava protected area (28º 28' S, 69º 10' W) is located in the southwest of La Rioja province and characterized by a system of shallow lagoons of saline and hypersaline water at altitudes of around 2,500–4,500 m (Wurstten et al., 2014). The protected areas encompass nearly 18,000 km2 and have become exposed, in many sections, to anthropogenic disturbances such as mining, livestock, and tourism (Didier et al., 2009; Silvani et al., 2017). Due to their regional importance and fragility, they were incorporated into conservation programs for protecting wetlands and their associated diversity. The protected area of San Guillermo was declared as Reserve of the Biosphere by UNESCO in 1980 while Laguna Brava was incorporated into the Ramsar Convention in 2003. These areas protect and preserve a relevant biodiversity associated with the local wetlands and plant communities, promoting the presence and reproduction of flamingos and camelids.

Fig. 1.

Map of the study area showing the location of the occurrences, the protected areas, and the regional topography.

(0.63MB).

The regional climate conditions are strongly influenced by the pronounced altitudinal gradient and the complex topography. The area is characterized by a wide thermal amplitude, with mean monthly temperatures ranging from 9 °C in summer to −7 °C in winter, strong cross-barrier winds, and most of the scarce precipitation (approximately 100 mm per year) falling as rain during summer months, with occasional snowfall during the winter season (Wurstten et al., 2014). The vegetation cover forms an extensive meadow, distributed in an altitudinal range from 3,800 to 4500 m with dominance of microthermal floristic components (Ontivero et al., 2022). It is composed of an association of grass steppes and shrubby species mostly in the area of the slopes with wetland ecosystems distributed in the lowlands.

Occurrence data

A database of georeferenced occurrences of the high-altitude wetlands was constructed from remote sensing using Sentinel 2B images from 2019 at a scale 1:50.000. The occurrences were selected on the images following the wetlands classification based on the local geomorphology and hydrology proposed by Ontivero et al. (2022). High-altitude wetlands are groundwater dependent systems in a close association with vegetation communities sharing high climate and environment affinities (Méndez, 2007). Recently, distribution models have been used to analyze wetlands as entities (Raney and Leopold, 2018). Therefore, high-altitude wetlands placed in the studied area were analyzed as entities, recording their coordinates for conducting the present study. Additionally, many of the identified wetlands were visited during different field trips to the area and their conservation status was discussed by Ontivero (2015). A total of 310 data points were recorded for both natural protected areas in zones with reduced mining activities and grazing (Fig. 1).

Climate variable selection

An initial set of 19 bioclimatic variables in spatial resolution of 30 arc s (∼1 km2) from the CHELSA v2.1 database (Karger et al., 2020) were used to characterize the distribution of high-altitude wetlands in the study area. This database has higher accuracy in mountain regions than other precipitation and temperature datasets, given that it better represents orographic precipitation effects (Bobrowski et al., 2021). The bioclimatic variables (Table 1) were derived from monthly temperature (maximum and minimum) and precipitation data that represent mean annual values, seasonality, and extreme conditions (O’Donnell and Ignizio, 2012). To avoid collinearity and select the most relevant bioclimatic parameters, the total of variables was reduced using the Spearman’s correlation, with a threshold of |r| = 0.80. Correlated variables were removed based on their highest contribution.

Table 1.

Biolimatic variables used in the modeling analysis.

Bioclimatic Index  Acronym  Unit 
Mean Annual Temperature  MAT  °C 
Mean Diurnal Range  AMDR  °C 
Isothermality  ISO  °C 
Temperature Seasonality  TS  (%) 
Max Temperature of the Warmest Month  MTWM  °C 
Min Temperature of the Coldest Month  mTCM  °C 
Annual Temperature Range  ATR  °C 
Mean Temperature of the Wettest Quarter  MTWEQ  °C 
Mean Temperature of the Driest Quarter  MTDQ  °C 
Mean Temperature of the Warmest Quarter  MTWQ  °C 
Mean Temperature of the Coldest Quarter  MTCQ  °C 
Annual Precipitation  AP  mm 
Precipitation of the Wettest Month  PPWEM  mm 
Precipitation of the Driest Month  PPDM  mm 
Precipitation Seasonality  PPS  Dimensionless 
Precipitation of the Wettest Quarter  PPWEQ  mm 
Precipitation of the Driest Quarter  PPDQ  mm 
Precipitation of the Warmest Quarter  PPWQ  mm 
Precipitation of the Coldest Quarter:  PPCQ  mm 
Potential distribution modeling

The potential distribution model was conducted using Wallace EcoMod 2.1 (Kass et al., 2018) in R (R Core Team, 2021) and the algorithm implemented in the maximum entropy (MaxEnt) method (Phillips, 2021). As the program works with georeferenced data points; the 310 sites corresponding to the occurrence data were added into the program for the analysis. A spatial filtering (spThin R package) at a minimum distance of 1 km among occurrence data was applied to reduce clustering caused by sampling bias (Aiello-Lammens et al., 2015). The background extent was set to minimum convex polygon. The buffer distance was 0.1 degrees and 10,000 background points were sampled. Block spatial partitioning was used and occurrence data were partitioned in four groups based on spatial clustering of occurrence points (Brown, 2014). To avoid model overfitting, regularization multiplier values ranging from 1 to 5 (with 0.5 steps) were set. A combination of four feature classes: linear (L), quadratic (Q), product (P), and hinge (H) were tested. The model evaluation was run through the ENMeval R package (Kass et al., 2021) and the selection was made considering the lowest Corrected Akaike Information Criterion (AICc) that corrects for sample size (Symonds and Moussalli, 2011; Warren and Seifert, 2011). The tenth percentile training presence threshold (10%TP) to provide a more comprehensive evaluation and quantification of the projected changes was applied. The presentation of the results in continuous grids were converted to binary grids (presence/absence) to discriminate between suitable and unsuitable areas. Although the threshold selection is still controversial, the 10%TP results are one of the most used and less restrictive, offering a good performance (Escalante et al., 2013).

Climate projections

To evaluate the current and future distributions of the wetlands the selected set of bioclimatic variables available in the CHELSA v2.1 dataset were used. Bioclimatic layers were constructed using five general circulation models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) (Eyring et al., 2016), in particular those selected for the Intersectoral Impact Model Intercomparison Project version 3b (ISIMIP3b, Lange, 2019).

The models from which data were taken are: GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, and UKESM1-0-LL, for a spatial resolution of 30-arc-second. Climatic data of the selected bioclimatic variables for a reference period (1981–2010) and for the mid-term time horizon (2041–2070) using the SSP1-2.6 and SSP5-8.5 scenarios as a forcing pathway. The SSP1-2.6 scenario will produce a multi-model mean of significantly less than 2 °C warming by 2100, and therefore can support analyses of this policy goal (O’Neill et al., 2016). The SSP5-8.5 represents a high fossil-fuel development throughout the 21st century, with a maximum projected global temperature change of almost 5 °C for 2100 (Meinshausen et al., 2020).

A multi-model ensemble considering the average of the five GCMs was obtained in order to reduce the uncertainty related to the model biases. This allowed us to obtain a single climatology for the reference period that was compared to the future climate in order to obtain differences among the bioclimatic indicators and, thus, the potential distribution of the high-altitude wetlands in the CAA.

ResultsClimate variables and model selection

The set of uncorrelated and best predictor variables after the distribution model analysis were: MAT, AMDR, ATR, AP, and PPS (Table 2). The model with the highest performance (lowest AICc = 4852.98) identified the optimal beta-multiplier and set of the most relevant bioclimatic variables, which were used to project habitat suitability of wetland ecosystems in the CAA both the reference and the future periods. According to the Maxent models, the lowest AICc scores and highest area under curve (AUC) for the training dataset were assigned to the model with a beta-multiplayer of LQ4. The calculated value of the AUC was 0.72 which indicates a good model performance. The climatic response curves obtained from the analysis indicated the optimal values of the five bioclimatic variables, which model the current climate conditions for the studied wetlands (Table 2).

Table 2.

Climate variable contribution and optimal range values for the five uncorrelated and most relevant variables.

Bioclimatic Variable  Contribution (%)  Optimal Range 
MAT (C°)  56.6  −8– (−1) 
AP (mm)  16.1  <150 
AMDR (C°)  14.1  10–15 
ATR (C°)  10  <24 
PPS (mm)  3.2  50–60 

MAT: mean annual temperature; AP: annual precipitation; AMDR: annual mean diurnal range; ATR: annual temperature range; PPS: precipitation seasonality.

Suitability of wetlands in the CAA for the reference period (1981–2010)

Considering the selected bioclimatic indices, a suitability map showing the potential distribution of high-altitude wetlands was obtained. Areas with higher suitability values are mainly located in an altitudinal range from 3,000 to 4,500 m, while those areas with low suitability (0–0.2) are located in the southeastern part of the San Guillermo protected area which corresponds to low altitudinal range low suitability areas are also observed in the western edge of this protected area, a factor that is linked to the presence of high-elevation areas with altitudes above 4,500 m. On the other hand, Laguna Brava protected area shows a high suitability for wetlands in most of its area, especially in the central part of the reserve, i.e., with altitudes ranging from 3,500 to 4,500 m, where the highest suitability values were obtained (Fig. 2).

Fig. 2.

Suitability map for the distribution of high-altitude wetlands in the CAA for the reference period (1981–2010).

(0.35MB).
Future changes in Andean wetlands

The resulting suitability map for high-altitude wetlands in the CAA using the projected bioclimatic indices for the SSP1-2.6 scenario shows a general reduction in the suitable areas for the period 2041-2070 compared to the reference period, particularly over the San Guillermo area (Fig. 3A). Much of the central and southernmost parts of the San Guillermo area has a projected reduction in suitability compared to the present, mainly for the areas located at elevation lower than 4,000 m. The projected reductions over the Laguna Brava protected area are located in the western portion of the region (Fig. 3B), where a large number of wetlands are placed (Fig. 1). The projected change in suitability shows higher percentages (i.e. more suitable future conditions for high-altitude wetlands) in those zones constrained to elevation values higher than 4,000 m (Figs. 1 and 3B).

Fig. 3.

(A) Future suitability for high-altitude wetlands in the CAA for the period 2041–2070, based on the SSP1-2.6 emissions scenario. (B) Projected suitability change for the distribution of wetlands in the study area for the same period and emissions scenario.

(0.6MB).

Based on the SSP5-8.5 scenario, results show a similar spatial pattern in terms of positive and negative projected changes in suitable areas for the period 2041–2070 (Fig. 4A–B). The most suitable areas are restricted to the northern part of the study area, mainly along the Laguna Brava protected area and the northeastern portion of San Guillermo protected area. Areas with low suitability (0–0.2) are projected mostly for the western and southeastern regions of San Guillermo protected area (Fig. 4A), which also correspond to low suitability areas based on the reference period (Fig. 2). The projected change in suitability shows higher percentages in the highest zones of the study area (Fig. 4B). This becomes more evident for the eastern section of Laguna Brava protected area, and the northern portion of San Guillermo protected area. The regions with elevation higher than 4,000 exhibit high percentage values, which means a projected positive change in suitability. Most of the San Guillermo area shows a negative change in suitability, particularly for those areas with altitude values lower than 4,000 m. The area with the highest density of sites in the western portion of Laguna Brava, also shows a projected decrease in suitability.

Fig. 4.

(A) Future suitability for high-altitude wetlands in the CAA for the period 2041–2070, based on the SSP5-8.5 emissions scenario. (B) Projected suitability change for the distribution of wetlands in the study area for the same period and emissions scenarios.

(0.59MB).

Projected climate change using a threshold confirms the placement of future optimal regions in the northern part of the studied area and deterioration towards the southern part. The projected changes in the distribution of high-altitude wetlands display an elevation dependence (Fig. 5A–C). The suitability for the studied wetlands was discriminated for the reference as well as for both future periods by the altitudinal bands (Fig. 6). The projected decrease in suitability can be linked to elevations below the 4,200 m a.s.l. with abrupt changes projected for the altitudinal range of 3,000 to 3,400 m a.s.l., and 3,600 to 3,800 m a.s.l. It is observed a projected increase in suitability in those areas located at elevations higher than 4,200 m a.s.l. The most notorious difference is shown between the projected future scenarios and the reference period for the regions located above 4,400 m a.s.l. It must be noted that the figure shows an optimum present suitability range for wetlands placed between 3,600–4,000 m a.s.l., sharply declining to elevations higher than 4,400 m a.s.l. (Fig. 6).

Fig. 5.

Presence of high-altitude wetlands in CAA using the tenth percentile training presence threshold (10% TP). (A) Reference period (1981–2010). (B) Period 2041–2070, based on the SSP1-2.6 emissions scenario. (C) Period 2041–2070, based on the SSP5-8.5 emissions scenario.

(0.43MB).
Fig. 6.

Average suitability for the high-altitude wetlands considering the selected sites, as a function of the elevation of the study area. Green bars show the values for the reference period while yellow and red bars show the values for the mid-term time horizon in both emissions scenarios.

(0.18MB).
DiscussionClimate influence on wetlands in the CAA

High-altitude wetlands in CAA are a combination of evergreen meadows in a hydromorphic soil supply by groundwater and running on a permanent water flow (Martínez-Carretero and Ontivero, 2016; Ontivero, 2015). Therefore, the distribution and presence of these ecosystems are in association with the regional available water resources and vegetations. Plant communities associated with wetlands in CAA are frequently distributed in altitudinal zonations with vegetation records rare or scarce at elevations between 4,200–4,300 m a.s.l. and practically absent for elevations over that range (Méndez, 2004; Ontivero, 2015).

Recent studies showed that changes in temperature and precipitation have disturbed the habitat of the Andean vegetation, limiting its distribution to small favorable refuges (Ruthsatz et al., 2020) with the consequent loss of native species (Schuchardt et al., 2023). In the studied protected areas, an increase of dried wetlands and habitat modifications with loss plant coverage and diversity, plus soil deterioration, have been observed during the last years. The status classification and analysis proposed by Ontivero (2015) for six wetlands in the CAA indicated that at least four were in a regular and intermediate conservation status with a high aridization percentage. In the San Guillermo natural protected area, high suitability areas are placed at an elevation range from 3,000 to 4,000 m. Low suitability ranges were in concordance with the lowlands (<3,000 m) placed in the southeastern part of the area, where the water supply for the development of wetlands is scant. Additionally, areas with low suitability are spread in the western part of the CAA at elevations above 4,500 m, a region characterized by rock glaciers and permafrost, with a geological setup that precludes the abundance of wetlands and vegetation. On the other hand, Laguna Brava natural protected area has a more equal distribution of highlands and suitability ranges (Fig. 2).

The present distribution of high-altitude wetlands in the CAA is strongly influenced by bioclimatic variables linked to temperature and precipitation (Fig. 2; Table 2). Seasonal precipitation is regularly stored as snowpack, playing a major role in determining plant composition. Precipitation either as rainfall or snow exerts a large influence in shaping wetlands in the Andes (Cooper et al., 2019; Feeley et al., 2020). Temperature combined with the spatial and temporal pattern of seasonal precipitation are relevant in determining the growing season in mountain ecosystems where plants spend a long dormancy due to frozen days and snow cover (Cavieres and Sierra-Almeida, 2018; Wipf and Rixen, 2010). In high-altitude ecosystems, snowmelt defines the onset and length of the growing warm season as well as the available water and nutrients that influence soil moisture and nutrient status for plants (Keller et al., 2005; Körner and Hiltbrunner, 2021). The regional temperature interacts with many components of the hydrological cycle, especially in determining the duration of the snowpack, a factor that affects the cover, diversity, and productivity of high-altitude plant communities (Wipf and Rixen, 2010). The MAT and ATR are also linked to water supply from glaciers, mainly in periods when snow has been completely melted (Crespo et al., 2020). Thermal daily and seasonal range have a relevant role in many physiological processes as flowering and seed production for high-altitude plants (Hülber et al., 2010). This relationship at daily and seasonal scale highlights the role of the temperature in the suitability of wetlands probably due to the modulation of water availability. In this sense, the warming recorded during the last decades in the Alps has caused plants with low thermal requirements to retreat towards more elevated areas (Körner and Hiltbrunner, 2021). In South America, several plant communities have experienced a thermophilization, positively correlated with warming rates (Feeley et al., 2020). Warming effect has caused a long-term reduction in glaciers along the CAA, a situation that has been exacerbated during the last decade (Dussaillant et al., 2019), along with the occurrence of extreme drought conditions that affected regional water availability (Rivera et al., 2021). Glacier retreat has a direct impact on wetlands (Bradley et al., 2006), and the sustained warming is causing a decline of permafrost (Rangecroft et al., 2016), a cryosphere component that plays a significant role as a water source in the study area (Esper-Angillieri, 2017).

Perspectives of Andean wetlands in the forthcoming decades under climate change scenarios

An increase in warming is projected for the next decades with serious effects for the highest altitudes of the CAA (Rivera et al., 2020), which implies an increase of regional temperatures from 1 °C to 4 °C by the end of the 21st century (Almazroui et al., 2021). Moreover, a reduction in the precipitation, runoff, and snow cover are projected for the rest of the century, even when the most optimistic emission scenarios are considered (Rivera et al., 2020; Almazroui et al., 2021). Analyzing the projected changes according to the selected bioclimatic indices for the two greenhouse gasses emissions scenarios, a negative change in the suitability for high-altitude wetlands is projected, particularly towards the south of 29 °S and the lower elevation areas (Figs. 3A and 4A). Projections for wetlands in the southernmost Andes of Patagonia indicated that suitable areas will decline under future climate conditions (Crego et al., 2014). A warmer climate will produce changes in life conditions for high mountain ecosystems at very small scales, mainly as a consequence of the reduction in snow duration, with species moving towards more optimal thermal habitats (Körner and Hiltbrunner, 2021). In this sense, species with a cold period required for seed germination may be endangered in the future (Cavieres and Sierra-Almeida, 2018).

Suitable areas in San Guillermo and Laguna Brava protected areas will be restricted at high-altitudes in the CAA (Figs. 3A and 4A) with a probable reduction of vegetation cover and composition. Recently, changes in vegetation cover and species richness for the CAAs in the north of Argentina were recorded with the migration of species towards wider thermal niches (Cuesta et al., 2023). Environmental changes, linked to regional climate variability, were detected in wetlands from the northwestern part of Argentina and southeastern Bolivia. These studies recorded a severe aridization trend in the Andean wetlands and a decrease in the plant productivity several decades ago (Carilla et al., 2013; Morales et al., 2015). Vegetation communities will probably reshape highland ecosystems in the next decades with different distribution patterns, abundances, and specific composition (Feeley et al., 2020).

The reduction of wetlands is the result of soil dryness and aridization with the loss of organic matter and the disappearance of the plant communities that formed these ecosystems. Therefore, the deterioration of wetlands in response to temperature encourages the loss of biodiversity. Additionally, an early snowmelt under future climate would imply snow-free mountain areas during prolonged periods that could promote the invasion by species potentially more efficient in flowering than natives (Keller et al., 2005; Martínez-Carretero and Ontivero, 2016; Schuchardt et al., 2023). According to the projected models, the most suitable areas with a greater presence of wetlands in the future will be restricted to altitudes of more than 4,200 m a.s.l. (Figs. 5A–C and 6). These areas will keep adequate climate conditions for high-altitude wetlands and probably for vegetation migration.

However, there are permanently frozen soils and pronounced slopes with intense cryogenic processes and soil mobilization in the CAA at altitudes above 4,600 m a.s.l. that makes the rooting of plants difficult (Méndez, 2007).

Even though the future changes in suitability indicate better climate conditions for the north part of the studied area in the higher altitudes, these ecosystems may experience a reduction with a high risk of disappearing in the next decades (Figs. 3A, 4A, and 5A–C). In the CAA, most of the wetlands face limited provincial or national environmental regulations in addition to human impacts such as mining, tourism or livestock. There are a few studies which attempt to define local parameters for evaluating their present and future conditions. Monitoring wetlands is an essential task that should be considered for preservation of these natural areas (Carilla et al., 2013; Cuesta et al., 2023; Ontivero, 2015).

The conservation of these sensitive high-mountain ecosystems will ensure that water resources originating in the Andes will be sustained in the future under the increased global warming conditions (Finlayson et al., 2019).

Conclusions

Climate change is a global concern related to human activities. Drastic consequences on the Earth’s ecosystems will rapidly increase as temperature continues rising.

For the Andean wetlands, placed in two natural protected areas from the CAA, changes in their distribution areas are projected for the period 2041–2070 under a low and a high-emissions scenario relative to the 1981–2010 period. Undoubtedly, the present distribution of high-altitude wetlands in the CAA will be modified under future climate change. The projected increase in temperature, together with a decrease in precipitation will have a direct effect on the suitable conditions, influencing the current distribution of Andean wetlands. Most of the negative changes in the distribution of wetlands will occur at altitudes below 4,200 m. The projected changes, considering the most relevant bioclimatic indices, will affect the future distribution of high-altitude wetlands, probably related to the reduction in water supply at lower elevations. These changes probably affect the plant communities, diversity, and their distribution in Andean wetlands following the reduction of the snow cover and cryosphere components such as glaciers and permafrost. High optimum areas will be placed at altitudes of more than 4,200 m. a.s.l. However, the vegetation and the associated life forms in areas of more than 4,500 m a.s.l. are currently restricted to the topography. It means a narrow range of possibilities for the existence of wetlands and green lands in the long term and a high risk of the disappearance of many endemic and native plants. Understanding their present conditions and any potentially threatening such as global climate change are essential to take action. Management practices for preserving these ecosystems are necessary in the short term. Our study represents the first step to explore the present and future distribution suitability of high-altitude wetland systems in CAA as regional entities under changing climate conditions. The intrinsic value of high Andean wetlands, the provision of ecosystem benefits, and the climate threats have led to increased awareness about the necessity of actions to protect and preserve these natural environments.

Conflict of interest

The authors declare not any conflict of interest.

Financial support

This research was supported by the National Agency for Scientific and Technological Promotion (ANPCyT) and the University Juan Agustín Maza (Grant No. PICTO-UUMM-2019-00004).

References
[Aiello-Lammens et al., 2015]
M.E. Aiello-Lammens, R.A. Boria, A. Radosavljevic, B. Vilela, R.P. Anderson.
spThin: an R package for spatial thinning of species occurrence records for use in ecological niche models.
Ecography, 38 (2015), pp. 541-545
[Almazroui et al., 2021]
M. Almazroui, M. Ashfaq, M.N. Islam, I.U. Rashid, K. Shahzad, M.A. Abid, E. O’Brien, M. Ismail, M.S. Reboita, A.A. Sörensson, P.A. Arias, L. Alves, M.K. Tippett, S. Saeed, R. Haarsma, F.J. Doblas-Reyes, F. Saeed, F. Kucharski, I. Nadeem, Y. Silva-Vidal, J.A. Rivera, M.A. Ehsan, D. Martínez-Castro, A.G. Muñoz, M.A. Ali, E. Coppola, M.B. Sylla.
Assessment of CMIP6 performance and projected temperature and precipitation changes Over South America.
Earth Syst. Environ., 5 (2021), pp. 155-183
[Bobrowski et al., 2021]
M. Bobrowski, J. Weidinger, U. Schickhoff.
Is new always better? Frontiers in global climate datasets for modeling treeline species in the Himalayas.
Atmosphere, 12 (2021), pp. 543
[Bradley et al., 2006]
R.S. Bradley, M. Vuille, H.F. Díaz, W. Vergara.
Threats to water supplies in the tropical Andes.
Science, 312 (2006), pp. 1755-1756
[Brown, 2014]
J.L. Brown.
SDMtoolbox: a python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses.
Methods Ecol. Evol., 5 (2014), pp. 694-700
[Bullock and Acreman, 2003]
A. Bullock, M. Acreman.
The role of wetlands in the hydrological cycle.
Hydrol. Earth Syst. Sci., 7 (2003), pp. 358-389
[Carilla et al., 2013]
J. Carilla, H.R. Grau, L. Paolini, M. Mariano.
Lake fluctuations, plant productivity, and long-term variability in high-elevation tropical Andean ecosystems.
Arctic Antarctic Alpine Res., 45 (2013), pp. 179-189
[Carilla et al., 2018]
J. Carilla, S. Halloy, S. Cuello, A. Grau, A. Malizia, F. Cuesta.
Vegetation trends over eleven years on mountain summits in NW Argentina.
Ecol. Evol., 8 (2018), pp. 11554-11567
[Cavieres and Sierra-Almeida, 2018]
L.A. Cavieres, A. Sierra-Almeida.
Assessing the importance of cold-stratification for seed germination in alpine plant species of the High-Andes of central Chile.
Perspect. Plant Ecol. Evol. Syst., 30 (2018), pp. 125-131
[Cooper et al., 2019]
D.J. Cooper, J. Sueltenfuss, E. Oyague, K. Yager, D. Slayback, E.M.C. Caballero, J. Argollo, B.G. Mark.
Drivers of peatland water table dynamics in the central Andes, Bolivia and Peru.
Hydrol. Processes, 33 (2019), pp. 1913-1925
[Crego et al., 2014]
R.D. Crego, C.K. Nielsen, K.A. Didier.
Climate change and conservation implications for wet meadows in dry Patagonia.
Environ. Conserv., 41 (2014), pp. 122-131
[Crespo et al., 2020]
S.A. Crespo, F. Fernandoy, L. Cara, S. Klarian, C. Lavergne.
First snow, glacier and groundwater contribution quantification in the upper Mendoza River basin using stable water isotopes.
Isot. Environ. Health Stud., 56 (2020), pp. 566-585
[Cuesta et al., 2023]
F. Cuesta, J. Carilla, L.D. LLambí, P. Muriel, M.V. Lencinas, R.I. Meneses, K.J. Feeley, H. Pauli, N. Aguirre, S. Beck, A. Bernardi, S. Cuello, S.A. Duchicela, P. Eguiguren, L.E. Gamez, S. Halloy, L. Hudson, R. Jaramillo, P.L. Peri, L.A. Ramírez, P. Rosero-Añazco, N. Thompson, K. Jagger, C. Tovar.
Compositional shifts of alpine plant communities across the high Andes.
Global Ecol. Biogeogr., 32 (2023), pp. 1591-1606
[Dangles et al., 2017]
O. Dangles, A. Rabatel, M. Kraemer, G. Zeballos, A. Soruco, D. Jacobsen, F. Anthelme.
Ecosystem sentinels for climate change? Evidence of wetland cover changes over the last 30 years in the tropical Andes.
[Didier et al., 2009]
K. Didier, M. Glennon, A. Novaro, E. Sanderson, S. Strindberg, S. Walker, S. Di Martino.
The landscape species approach: spatially-explicit conservation planning applied in the Adirondacks, USA, and San Guillermo-Laguna Brava, Argentina, Landscapes.
Oryx, 43 (2009), pp. 476-487
[Durigan et al., 2022]
G. Durigan, C.B. Munhoz, M.J.B. Zakia, R.S. Oliveira, N. Pilon, R.S.T. do Valle, B.M.T. Walter, E.A. Honda, A. Pott.
Cerrado wetlands: multiple ecosystems deserving legal protection as a unique and irreplaceable treasure.
Perspect. Ecol. Conserv., 20 (2022), pp. 185-196
[Dussaillant et al., 2019]
I. Dussaillant, E. Berthier, F. Brun, M. Masiokas, R. Hugonnet, V. Favier, A. Rabatel, P. Pitte, L. Ruiz.
Two decades of glacier mass loss along the Andes.
Nat. Geosci., 12 (2019), pp. 802-808
[Escalante et al., 2013]
T. Escalante, G. Rodríguez-Tapia, M. Linaje, P. Illoldi-Rangel, R. González-López.
Identification of areas of endemism from species distribution models: threshold selection and Nearctic mammals.
TIP Revista Especializada en Ciencias Químico-Biológicas, 16 (2013), pp. 5-17
[Esper-Angillieri, 2017]
M.Y. Esper-Angillieri.
Permafrost distribution map of San Juan Dry Andes (Argentina) based on rock glacier sites.
J. South Amer. Earth Sci., 73 (2017), pp. 42-49
[Eyring et al., 2016]
V. Eyring, S. Bony, G.A. Meehl, C.A. Senior, B. Stevens, R.J. Stouffer, K.E. Taylor.
Overview of the coupled model intercomparison project phase 6 (CMIP6) experimental design and organization.
Geosci. Model Dev., 9 (2016), pp. 1937-1958
[Feeley et al., 2020]
K.J. Feeley, C. Bravo-Avila, B. Fadrique, T.M. Perez, D. Zuleta.
Climate-driven changes in the composition of New World plant communities.
Nat. Clim. Change, 10 (2020), pp. 965-970
[Finlayson et al., 2019]
C.M. Finlayson, G.T. Davies, W.R. Moomaw, G.L. Chmura, S.M. Natali, J.E. Perry, N. Roulet, A.E. Sutton-Grier.
The second warning to humanity–providing a context for wetland management and policy.
Wetlands, 39 (2019), pp. 1-5
[Huggel et al., 2019]
C. Huggel, V. Muccione, M. Carey, R. James, C. Jurt, R. Mechle.
Loss and damage in the mountain cryosphere.
Reg. Environ. Change, 19 (2019), pp. 1387-1399
[Hugonnet et al., 2021]
R. Hugonnet, R. McNabb, E. Berthier, B. Menounos, C. Nuth, L. Girod, D. Farinotti, M. Huss, I. Dussaillant, F. Brun, A. Kääb.
Accelerated global glacier mass loss in the early twenty-first century.
Nature, 592 (2021), pp. 726-731
[Hülber et al., 2010]
K. Hülber, M. Winkler, G. Grabherr.
Intraseasonal climate and habitat‐specific variability controls the flowering phenology of high alpine plant species.
Funct. Ecol., 24 (2010), pp. 245-252
[Izquierdo et al., 2022]
A.E. Izquierdo, C. Blundo, J. Carilla, J. Foguet, C.J. Navarro, E. Casagranda, M.F. Chiappero, M.V. Vaieretti.
Floristic types of high‐Andean wetlands from northwest Argentina and their remote‐sensed characterization at a regional scale.
Appl. Veg. Sci., 25 (2022),
[Karger et al., 2020]
D.N. Karger, D. Schmatz, D. Dettling, N.E. Zimmermann.
High resolution monthly precipitation and temperature time series for the period 2006-2100.
[Kass et al., 2018]
J.M. Kass, B. Vilela, M.E. Aiello-Lammens, R. Muscarella, C. Me­row, R.P. Anderson.
Wallace: a flexible platform for reproducible modeling of species niches and distributions built for community expansion.
Methods Ecol. Evol., 9 (2018), pp. 1151-1156
[Kass et al., 2021]
J. Kass, R. Muscarella, P.J. Galante, C.L. Bohl, G.E. Pinilla-Buitrago, R.A. Boria, M. Soley-Guardia, R.P. Anderson.
ENMeval: Automated Tuning and Evaluations of Ecological Niche Models.
[Keller et al., 2005]
F. Keller, S. Goyette, M. Beniston.
Sensitivity analysis of snow cover to climate change scenarios and their impact on plant habitats in alpine terrain.
Clim. Change, 72 (2005), pp. 299-319
[Körner and Hiltbrunner, 2021]
C. Körner, E. Hiltbrunner.
Why is the alpine flora comparatively robust against climatic warming?.
Diversity, 13 (2021), pp. 383
[Lange, 2019]
S. Lange.
Trend-preserving bias adjustment and statistical downscaling with ISIMIP3BASD (v1.0).
Geosci. Model Dev., 12 (2019), pp. 3055-3070
[Malhi et al., 2020]
Y. Malhi, J. Franklin, N. Seddon, M. Solan, M.G. Turner, C.B. Field, N. Knowlton.
Climate change and ecosystems: threats, opportunities and solutions.
Philos. Trans. R. Soc. B, 375 (2020),
[Martínez-Carretero and Ontivero, 2016]
E. Martínez-Carretero, M. Ontivero.
Ecosistema altoandino de importancia biológica, ecológica y socio-económica, pp. 187-197
[Meinshausen et al., 2020]
M. Meinshausen, Z.R.J. Nicholls, J. Lewis, M.J. Gidden, E. Vogel, M. Freund, U. Beyerle, C. Gessner, A. Nauels, N. Bauer, J.G. Canadell, J.S. Daniel, A. John, P.B. Krummel, G. Luderer, N. Meinshausen, S.A. Montzka, P.J. Rayner, S. Reimann, S.J. Smith, M. van den Berg, G.J.M. Velders, M.K. Vollmer, R.H.J. Wang.
The shared socio-economic pathway (SSP) greenhouse gas concentrations and their extensions to 2500.
Geosci. Model Dev., 13 (2020), pp. 3571-3605
[Méndez, 2004]
E. Méndez.
La vegetación de los Altos Andes I: Pisos de vegetación del flanco oriental del Cordón del Plata (Mendoza, Argentina).
Boletín de la Sociedad Argentina de Botánica, 39 (2004), pp. 227-253
[Méndez, 2007]
E. Méndez.
La vegetación de los Altos Andes II: Las Vegas del flanco oriental del Cordón del Plata (Mendoza, Argentina).
Boletín de la Sociedad Argentina de Botánica, 42 (2007), pp. 273-294
[Mitsch et al., 2013]
W.J. Mitsch, B. Bernal, A.M. Nahlik, Ü. Mander, L. Zhang, C.J. Anderson, E. Jørgensen, H. Brix.
Wetlands, carbon, and climate change.
Landscape Ecol., 28 (2013), pp. 583-597
[Morales et al., 2015]
M.S. Morales, J. Carilla, H.R. Grau, R. Villalba.
Multi-century lake area changes in the Southern Altiplano: a tree-ring-based reconstruction.
Clim. Past, 11 (2015), pp. 1139-1152
[O’Donnell and Ignizio, 2012]
M.S. O’Donnell, D.A. Ignizio.
Bioclimatic predictors for supporting ecological applications in the conterminous United States.
US Geological Survey Data Series, 691 (2012), pp. 4-9
[O’Neill et al., 2016]
B.C. O’Neill, C. Tebaldi, D.P. van Vuuren, V. Eyring, P. Friedlingstein, G. Hurtt, R. Knutti, E. Kriegler, J.F. Lamarque, J. Lowe, G.A. Meehl, R. Moss, K. Riahi, B.M. Sanderson.
The Scenario Model intercomparison project (ScenarioMIP) for CMIP6.
Geosci. Model Dev., 9 (2016), pp. 3461-3482
[Ontivero, 2015]
M. Ontivero.
Caracterización biofísica de vegas de los altos andes centrales de Argentina. Tesis doctoral.
Universidad Nacional de Córdoba, (2015), pp. 1-171
[Ontivero et al., 2022]
M. Ontivero, E. Martínez-Carretero, L. Perucca.
Clasificación de humedales de montaña (Vegas) en los Andes Centrales de Argentina.
Boletín de la Sociedad Argentina de Botánica, 57 (2022), pp. 41-50
[Parmesan et al., 2000]
C. Parmesan, T.L. Root, M.R. Willig.
Impacts of extreme weather and climate on terrestrial biota.
Bull. Am. Meteorol. Soc., 81 (2000), pp. 443-450
[Pauli and Halloy, 2019]
H. Pauli, S.R. Halloy.
High mountain ecosystems under climate change.
Oxford Research Encyclopedias of Climate Science, Oxford University Press, United Kingdom, (2019), http://dx.doi.org/10.1093/acrefore/9780190228620.013.764
[Phillips, 2021]
S.J. Phillips.
Maxnet: Fitting ‘Maxent’ Species Distribution Models with ‘glmnet’. CRAN. R package version 0.4.1.
[R Core Team, 2021]
R Core Team.
R: A language and environment for statistical computing.
[Raney and Leopold, 2018]
P.A. Raney, D.J. Leopold.
Fantastic wetlands and where to find them: modeling rich fen distribution in New York State with Maxent.
Wetlands, 38 (2018), pp. 81-93
[Rangecroft et al., 2016]
S. Rangecroft, A.J. Suggitt, K. Anderson, J. Magrath, A.P. Castel, P. Pacheco.
Future climate warming and changes to mountain permafrost in the Bolivian Andes.
Clim. Change, 137 (2016), pp. 231-243
[Rivera et al., 2020]
J.A. Rivera, E. Naranjo-Tamayo, M. Viale.
Water resources change in Central-Western Argentina under the Paris Agreement warming targets.
[Rivera et al., 2021]
J.A. Rivera, S. Otta, C. Lauro, N. Zazulie.
A decade of hydrological drought in Central-Western Argentina.
[Rubio et al., 2017]
M.C. Rubio, C. Rubio, M.A. Salomón, E. Abraham.
Conservación de servicios ecosistémicos en humedales altoandinos: la participación social en la creación de un área natural protegida.
Ecología Austral, 27 (2017), pp. 177-192
[Ruthsatz et al., 2020]
B. Ruthsatz, K. Schittek, B. Backes.
The vegetation of cushion peatlands in the Argentine Andes and changes in their floristic composition across a latitudinal gradient from 39 S to 22 S.
Phytocoenologia, 50 (2020), pp. 249-278
[Sala et al., 2000]
O.E. Sala, F. Stuart-Chapin, J.J. Armesto, E. Berlow, J. Bloomfield, R. Dirzo, E. Huber-Sanwald, L.F. Huenneke, R.B. Jackson, A. Kinzig, R. Leemans, D.M. Lodge, H.A. Mooney, M. Oesterheld, N. Leroy-Poff, M.T. Sykes, B.H. Walker, M. Walker, D.H. Wall.
Global biodiversity scenarios for the year 2100.
Science, 287 (2000), pp. 1770-1774
[Salimi et al., 2021]
S. Salimi, S.A.A.A.N. Almuktar, M. Scholz.
Impact of climate change on wetland ecosystems: a critical review of experimental wetlands.
J. Environ. Manag., 286 (2021),
[Schuchardt et al., 2023]
M.A. Schuchardt, B.J. Berauer, J. Giejsztowt, A.V. Heßberg, Y. Niu, M. Bahn, A. Jentsch.
Drought erodes mountain plant community resistance to novel species under a warming climate.
Arctic Antarctic Alpine Res., 55 (2023),
[Silvani et al., 2017]
V.A. Silvani, R.P. Colombo, M.V. Scorza, L. Fernández-Bidondo, C.P. Rothen, A. Scotti, S. Fracchia, A. Godeas.
Arbuscular mycorrhizal fungal diversity in high-altitude hypersaline Andean wetlands studied by 454-sequencing and morphological approaches.
Symbiosis, 72 (2017), pp. 143-152
[Symonds and Moussalli, 2011]
M.R. Symonds, A. Moussalli.
A brief guide to model selection, multimodel inference and model averaging in behavioral ecology using Akaike’s information criterion.
Behav. Ecol. Sociobiol., 65 (2011), pp. 13-21
[Tovar et al., 2022]
C. Tovar, A.F. Carril, A.G. Gutiérrez, A. Ahrends, L. Fita, P. Zaninelli, P. Flombaum, A.M. Abarzúa, D. Alarcón, V. Aschero, S. Báez, A. Barros, J. Carilla, M.E. Ferrero, S.G. Flantua, P. Gonzáles, C.G. Menéndez, O.A. Pérez-Escobar, A. Pauchard, R.C. Ruscica, T. Särkinen, A.A. Sörensson, A. Srur, R. Villalba, P.M. Hollingsworth.
Understanding climate change impacts on biome and plant distributions in the Andes: challenges and opportunities.
J. Biogeogr., 49 (2022), pp. 1420-1442
[Warren and Seifert, 2011]
D.L. Warren, S.N. Seifert.
Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria.
Ecol. Appl., 21 (2011), pp. 335-342
[Wipf and Rixen, 2010]
S. Wipf, C. Rixen.
A review of snow manipulation experiments in Arctic and alpine tundra ecosystems.
Polar Res., 29 (2010), pp. 95-109
[Wurstten et al., 2014]
A. Wurstten, A.J. Novaro, R.S. Walker.
Habitat use and preference by guanacos, vicuñas, and livestock in an altitudinal gradient in Northwest Argentina.
Eur. J. Wildlife Res., 60 (2014), pp. 35-43
Copyright © 2024. Associação Brasileira de Ciência Ecológica e Conservação
Perspectives in Ecology and Conservation
Article options
Tools