Next Article in Journal
Metamitron and Shade Effects on Leaf Physiology and Thinning Efficacy of Malus × domestica Borkh
Next Article in Special Issue
Defoliation Dynamics in Kikuyugrass Pastures Subjected to Intensities of Defoliation
Previous Article in Journal
Adult Desert Locust Swarms, Schistocerca gregaria, Preferentially Roost in the Tallest Plants at Any Given Site in the Sahara Desert
Previous Article in Special Issue
Long-Term Yield and Quality Performance of Perennial Energy Grasses (Agropyron spp.) on Marginal Land
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Páramo Ecosystems in Ecuador’s Southern Region: Conservation State and Restoration

by
Víctor J. García
1,2,*,
Carmen O. Márquez
3,4,
Marco V. Rodríguez
5,
Jonathan J. Orozco
5,
Christian D. Aguilar
6 and
Anita C. Ríos
3
1
Facultad de Ingeniería, Carrera de Ingeniería Civil, Universidad Nacional de Chimborazo, Riobamba, Provincia de Chimborazo 060150, Ecuador
2
Facultad de Ciencias, Universidad de Los Andes, Mérida, Estado Mérida 5101, Venezuela
3
Facultad de Ingeniería, Carrera de Ingeniería Ambiental, Universidad Nacional de Chimborazo, Riobamba, Provincia de Chimborazo 060150, Ecuador
4
Facultad de Ciencias Forestales y Ambientales, Universidad de Los Andes, Mérida, Estado Mérida 5101, Venezuela
5
Dirección de Investigación, Universidad Nacional de Chimborazo Riobamba, Provincia de Chimborazo 060150, Ecuador
6
Mancomunidad Bosque Seco, Bosque Seco, Provincia de Loja 110104, Ecuador
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(12), 1922; https://doi.org/10.3390/agronomy10121922
Submission received: 29 October 2020 / Revised: 27 November 2020 / Accepted: 2 December 2020 / Published: 7 December 2020
(This article belongs to the Special Issue Sustainable Management and Utilization of Permanent Grassland)

Abstract

:
The páramo is home to a significant proportion of global biodiversity and provides essential services for the development of life for millions of people in Ecuador. However, land use/land cover (LULC) changes threaten biodiversity and modify its functioning. The objectives of this study were: (1) to evaluate the conservation status of the herbaceous páramo (HP) ecosystem by analyzing its LULC in Ecuador’s southern region. (2) to identify possible regions where the native páramo ecosystem is being restored. We analyzed Landsat 8 images using Object-Based Image Analysis (OBIA) and a Classifier Decision Tree (CDT) to achieve these objectives. The results show that the native herbaceous páramo (NHP) ecosystem is being transformed into an anthropogenic HP (AHP). The area covered by the NHP ecosystem (296,964 ha) has been reduced by 50% (149,834 ha). Nevertheless, we identified five regions where the NHP is upgrading. These regions are relevant for studying NHP regeneration in Ecuador’s southern region, where soils are mostly andosols. The LU of the páramo, with cycles of exploitation, abandonment, and regeneration in a secondary páramo, is transforming the NHP ecosystem. These exploitation practices, global climate change, and lack of knowledge about the NHP ecosystem’s regeneration and its soils’ recovery threaten to substantially reduce the NHP area, its functionality, and its ecosystem services.

1. Introduction

The páramos are humid tropical ecosystems located from the forest’s upper limit with herbaceous and shrub vegetation [1,2]. The páramos is home to a significant proportion of the Ecuadorian biodiversity and provide essential services for the development of life for millions of people in Ecuador (e.g., supplying water and storing carbon) [2,3,4,5]. These ecosystemic services are possible due to a favorable relationship between biodiversity and ecosystem functioning [6].
When the native páramo ecosystem (PE) is destroyed, biodiversity recovery can take hundreds to thousands of years. Although the number of different species may increase over time, the plant community colonizing the secondary páramo remains distinct from that of the native páramo [7]. Therefore, recognizing that the native páramo is different from the secondary páramo will improve our understanding of the relationship between biodiversity, community assembly time, and ecological functioning. Knowing these differences shed light on understanding human activities’ ecological consequences and global climate change on the native and secondary páramo ecosystems. A previous study assessed the páramos of central Ecuador’s conservation state based on the differentiation between natural and secondary (anthropogenic) PE [8]. Given the global importance of the páramos, the páramo ecosystems’ dynamics must be understood to advance our comprehension of the páramo ecosystems’ response to environmental changes and anthropogenic disturbances such as farming and agriculture.
The Andean corridor hosts the Ecuadorian PEs [9]. Their cold and humid climate has allowed several plant species to develop unique adaptive strategies to cope with reasonably hostile environmental conditions [10]. These particular conditions have led to high endemism and a remarkable diversity of plants, the highest among all alpine landscapes in the world [11]. Another characteristic of the PE is its capacity to store water due to the high content of organic matter (OM) of its soils [12,13]. The high content of OM in its soils is favored by environmental conditions such as a low intensity but long duration of annual rainfall and high temperature during the day and low temperature during the night [14,15]. The alteration of native PEs caused by land use/land cover (LULC) change is the main threat to biodiversity and a predominant modifier of ecosystem functioning. These alterations continue to increase, driven primarily by economic and demographic imperatives associated with intensive exploitation of essential resources needed for human needs [7,16]. Therefore, the restoration of degraded PE is recognized as a global priority.
There are few reports about how the intensive exploitation of the páramo soils in Ecuador’s southern region has affected the páramo ecosystems [3,4,5,17,18]. Very little is known about the regeneration of native PE. Sustainable development of the páramos will fail if we do not know how the exploitation of the páramos affects these ecosystems and what factors determine the regeneration of the native PE.
In previous studies, García et al. [8] reported a methodology based on the OBIA and CART algorithm to analyze the conservation state of the páramos in the central region of Ecuador. These authors found a decision tree with six predictor variables to classify the LULC at the PE. These authors reported that the methodology is simple, fast, and can be implemented in other Ecuador regions as an exploratory tool for evaluating and planning the páramos.
The objectives of this study were: (1) to evaluate the conservation state of the herbaceous paramo ecosystem by analyzing the LULC of the soils under the páramo ecosystems in Ecuador’s southern region; and (2) to identify possible regions where the native herbaceous páramo (NHP) is being restored. To achieve these objectives, we analyzed Landsat 8 images using the object-based methodology and decision tree to classify objects. We generated an updated database of LULC in the páramos ecosystems in Ecuador’s southern region.

2. Materials and Methods

2.1. Area of Study

The study area is in Ecuador’s southern region. It extends from the north of Cañar province, between the coordinates UTM 18 S, Datum WGS84 (X = 742,040; Y = 973,6413), and UTM 17 South (X = 684,253; Y = 9,463,053) to the south of Loja province (Figure 1). The study area was segmented into 62 regions of 20 × 20 km (40,000 ha) each region was processed individually.
The study area has a diverse topography where the elevations vary from 2600 to 4480 m above sea level (m.a.s.l.), which gives rise to various ecological interactions and plant formations [19,20]. Cañar and Azuay’s region has a cold and humid climate influenced by the continental air masses of the Amazon basin and the west’s dry and cold air masses. The annual average precipitation is 1210 ± 101 mm per year [17]. In the provinces of Loja and Zamora, the average annual temperature varies from 7 to 25 °C, with a rainfall of 500 to 800 mm per year. The rain predominates in the afternoon, and part of it is of low intensity, known as a drizzle. The Andean mountain range determines the precipitation regime of the entire study area. Precipitation is low but constant throughout the year, so paramos’ ecosystems play the role of water regulators by continually providing the water produced by rain. This water flows down the slopes and riparian areas [21,22].
The geology of the study area is dominated by the Tarqui Formation, occupying 51% of the area (Agglomerates, tuff agglomerates, and tufas of intermediate to acid composition), while 10% is occupied by the Zamora series (Sandstones, siltstones, and shales), 8% by the Saraguro Formation (Pyroclasts: tufas and thick agglomerates with lava blocks and lavas: andesites with alternating porphyritic texture), 4% metamorphic rocks, 5% glacial deposits, and the remaining 22% is dominated by a variety of geological units, each with less than 3% representation in the paramo ecosystems in the area of study [18].
Most of the soils under páramo located in the south are classified as Andosols [3]. Regardless of their parent material, páramo’s soils have very dark epipedons consisting of Ah and A horizons with little differentiation, light overlying colored subsoil (C horizon) rich in minerals (dominated by clay) [17]. Among the most important soil properties are high organic C content (170–200 g/kg), very low bulk density (400 kg/m3), and absence of allophane. The soils have developed strong hydric properties, with water retention at −1500 kPa of 1.4 g/g [18].

2.2. The Herbaceous Páramo

The HP is characterized by species of the genera Calamagrostis, Agrostis, Festuca, Cortaderia, Stipa, and an abundant diversity of lifeforms [20]. In south-central Ecuador, in the Cajas páramo, the plant community varies as humidity conditions drop, and associations are created between Calamagrostis sp. and Viola humboldtii. In the southern part of Ecuador, the paramo ecosystem descends to 2600 m.a.s.l. and is composed of Agrostis breviculmis, Calamagrostis spp., Festuca asplundii, and Stipa ichu; in areas with steep slopes greater than 70% [2]. In lower altitude areas, the HP ecosystem is usually in association with shrub species.
The Ministry of Environment of Ecuador (MEE) defines the HP according to vegetation types, precipitation, and soils. The HP covers 73% of the study area, where the HP ecosystem (61%) and the HP ecosystem with evergreen shrubs (12%). A semi-deciduous shrub ecosystem covers 15% of the area. Different ecosystems cover the remaining 12% of the area [20]. For the sake of comparison, Figure 5a shows the HP páramo ecosystem distribution reported by the MEE [20], which served as a baseline. The methodology implemented by the MEE [23] includes six environmental variables: bioclimate, ecological floors, flood regime, relief, phenology, and biogeography. The bioclimatic model was built from the variables bioclimate, ambrotypes, thermotypes, and ecological floors. The flooding and fluvial origins are combined with the general relief map, macro, and mesorelief to give the geoforms model. The ecosystems’ map was built from the bioclimatic model, geoforms model, phenology, biogeography, and land cover (natural vegetation) map. Landsat and Rapideye images acquired in 2007 by the MEE were used to generate the land cover map (natural vegetation). These images were pre-processed and corrected. A supervised classification was performed; the resulting environmental units were verified in the field. The vegetation map used the classes defined by the ecosystem type developed by the ecosystem classification system for continental Ecuador [20]. Finally, the ecosystem map was subjected to a process of (1) spatial debugging to eliminate spurious or unrepresentative patterns, considering a minimum mapping unit of 25 ha., and (2) verification of thematic consistency with floristic data and expert criteria.

2.3. Workflow

The methodology comprised four main stages: (1) image pre-processing, implementation of OBIA, (2) attribute extraction, (3) validation and performance evaluation of the decision tree, and (4) land cover and land use classification in the study area, map generation, database, and analysis (Figure 2).
Image pre-processing included downloading Landsat 8 images, panoramic sharpness, and brightness, and reflectance correction. The implementation of OBIA included image segmentation and object generation. Attribute extraction included extracting the maximum, mean, and minimum reflectance values, spectral indices, and ancillary data for each object. Decision tree validation involved sampling 610 random objects categorized by the decision tree and validating them using the Google Earth platform. Classification and analysis included a Python code that allowed us to classify each object according to the category and generate a map with the different LULC categories.

2.4. Landsat-8 Images

The Landsat 8 images were downloaded from the official NASA website, “earthexplorer.usgs.gov.” Cloudiness of less than 30% was considered acceptable to make the visible field as complete as possible. The date of the downloaded images was 20 November 2016, and the grid number on the map (path = 10/row = 62 for the South Central zone and path = 10/row = 63 for the South zone).
The study area contemplated the use of two satellite images with the following labels: (1) LC80100612016325LGN01 and LC80100622016325LGN01; (2) time of year: summer; (3) route/row: 010/062 and 010/063; (4) azimuth of the sun (degrees): 128.09132614 and 125.72777466; (5) sun elevation (degrees): 61.31854934 and 62.02116942; (6) scene center time: 15:27:26 and 15:27:50; (7) cloud cover: 27.75% and 5.44%; (8) sensors: OLI_TIRS; and, (9) datum—projection: WGS 84—UTM ZONE 17S.
The geometry accuracy of Landsat 8 images was verified using the topographic charts and base cartography of rivers and roads, at a scale of 1:50,000, georeferenced in the UTM Datum WGS84 projection of the Ecuadorian Military Geographical Institute [24].

2.5. Categorization of Land Use and Land Cover

Land use and the cover were categorized into five classes: NHP, anthropogenic herbaceous páramo (AHP), forest (FRS), water bodies (WTR), and snow (SNW). Figure 3 shows the classes of LULC based on visual observation of the landscape and imagery: Water—watercourses and water bodies, including rivers, small lakes, and reservoirs. The water class includes the lake system and adjacent areas with high susceptibility to flooding. Forest—forest, and shrubs. The forests in the province of Cañar y Azuay include forests and shrubs located in low windy and humid areas. They are dominated by the genus Polylepis with heights ranging from 5 to 7 m. The forests in the provinces of Loja y Zamora are mainly evergreen shrubs that reach heights up to 3 m and found from 2800 to 3000 m.a.s.l.; the NHP is dominated by Calamagrostis spp. In regions well to the south, the NHP is associated with shrub species, extending from 2800 to 3900 m.a.s.l. The AHP is the páramo that has been transformed by humans. It includes farmland, pasture, bare soil, built/infrastructure areas, roads, burned areas, and regenerated grasslands with other species’ association and exhibits alterations in its phenology. This class is in areas with altitudes lower than 3800 m.a.s.l. and slopes lower than 11%. They are buffer zones and human settlements that allow the transition from Andean forests to crops, plantations, or infrastructure. The snow class corresponds to the areas with an altitude greater than 3827 m.a.s.l., with small deposits of ice crystals giving the highest values of reflectance.

2.6. Object-Based Image Analysis

The images of each scene were processed in ArcGIS software 10.3, Redlands, California 92373 8100, USA (http://www.esri.com/software/arcgis). We followed the workflow suggested by Urbanski [25] to convert Landsat-8 OLI images into a land cover map where the fundamental units are geographical objects. In total, 1,876,324 geographical objects with common attributes were obtained for the 62 regions under study.
A map of LULC in the study area was generated using the CDT shown in Figure 4. This CDT served as an exploratory tool for data-based discovery and prediction and gaining new knowledge about LULC under the paramo ecosystems in Ecuador’s central region [8]. The CDT requires only six predictor variables (reflectance, VARIg, altitude, EVI2, NBR2, and slope). In this study, we used this CDT to predict the LULC in the páramo ecosystems in Ecuador’s southern region.

2.7. Predictor Variables

The predictor variables were recovered from the Landsat 8 bands and the DEM. A total of 6 predictor variables were used, including the necessary spectral information on reflectance, three spectral vegetation indices (SVI), and two indicators derived from DEM (Table 1). Band algebra was used through the ENVI Math Band tool to generate the detailed indices according to their corresponding Landsat 8 bands. The SVI were used as predictor variables because they are the primary source of information for monitoring and evaluating land cover [8,12,26].

2.8. Decision Tree

Although other classifiers could achieve better classification results, the CDT represents a clear set of rules with values that define thresholds used in similar contexts without training [30]. Cao et al. [31] reported that decision trees could provide stable performance and reliable machine learning and data mining research results. Thus, decision trees have been used to identify spectral bands with the highest class discrimination capabilities and low misclassification rates [31]. Decision trees have also been used to identify greenhouse crops from remote sensing data [30]. However, in general, CDT rules are restricted by the root variable. In this study, only rules that are indicated by “Reflectance” were extracted. The CDT in Figure 4 exhibited excellent performance in the learning and validation phase for the five categories used in the characterization of LULC in the páramo ecosystem in Ecuador’s central region [8]. In the CDT’s performance, the predictor variables’ relevance was 100% for the reflectance, 63% for the VARIg index, 56% for the altitude, 46% for the EVI2 index, 34% for the NBR2 index, and 23% for the slope.

2.9. Validation of the Classification Decision Tree

The use of the CDT in Figure 4 was validated in Ecuador’s southern region using 610 random objects. Among the 610 objects, 50 belong to the category WTR, 100 to FRS, 225 to NHP, 225 to AHP, and 10 to SNW. The allocation of the number of samples per category was proportional to the total number of objects in each category. We choose 610 objects for the sake of practicality and statistical reasons. In reality, it could be any number larger than 385, the number of samples statistically representative of a population of 1,876,324 objects, with a confidence level of 95% and a confidence interval of 5%. The allocation of samples per category was proportional to the total number of objects in each category.
The classification results by the CDT were verified with the Google Earth platform’s information. Table 2 shows the table of confusion obtained in the validation of the CDT in the study area. The CDT performance was evaluated by a binary crosstabulation derived from the confusion matrix (Table 3).
The CDT’s overall accuracy was above 94%. The misclassification rate was below 6%. The CDT’s specificity (how well the CDT can recognize negative samples) was above 94%. The probability of an informed classification for each class was above 86%. How much each class is marked as a predictor or possible cause of another was above 70%. Matthew’s correlation coefficient between observed and predicted classification was above 84%. Therefore, the percentage of correct categorizations of the 610 known ground objects was 91% (Table 2). This value is acceptable if we consider the heterogeneity of the study area. In addition, there is little information on land use and land cover under paramo ecosystems in the southern region of Ecuador. Therefore, the CDT can be used as an exploratory tool to classify and differentiate LULC objects under the páramo ecosystem in Ecuador’s southern region. The CDT allowed the mapping of geographical objects’ input space into a predefined class with good performance.

3. Results

Figure 5b shows the current distribution of LULC in the study area. The MEE report is the official tool for the characterization and definition of ecosystems at the national level. The Ecosystem Classification System of Ecuador has been established from continental to fine scales (landscape, local) [20].
The MEE reported in 2013 [20] the existence of 296,964 ha of NHP. There is no report for the AHP surface area. In our study, we found that the surface area covered by NHP, AHP, FRS, WTR, and SNW was 149,834 ha (39%), 146,829 ha (38%), 61,404 ha (16%), 27,916 ha (7%) and 726 ha (0.2%), respectively.
Table 4 shows the distribution of LULC for each region in the study area. The regions with the largest area of NHP, AHP, FRS, WTR, and SNW are region 47, with 27,377 ha, region 40 (12,431 ha), region 29 (5152 ha), region 22 (3210 ha), and region 46 (313 ha), respectively (Table 4). Each region has an area of 40,000 ha. Therefore, 74% of the surface of region 47 is covered by NHP. Region 40 exhibits the largest anthropogenic activity area, with 35% of its surface area of AHP. Forests cover 33% of the surface area of Region 29. Water covers 19% of Region 22, and 3% of Region 46 is covered by snow.
There are 14 regions (4, 9, 11, 13, 24, 25, 27, 31, 38, 42, 51, 57,58, and 59) where the MEE map and our analysis agree that there is no area covered by the NHP ecosystem. Meanwhile, in regions 1, 3, 6, 10, and 16, the MEE map reports no NHP, and our results report the existence of 42, 837, 133, 517, and 539 ha of NHP, respectively. This result suggests that the CDT identifies that the HP ecosystem is upgrading in these regions. Similarly, in regions 12 and 26, the MEE reports 1021 and 108 ha of NHP, and our results show that in these regions, there is no area covered by NHP. On the other hand, the CDT identified 23 regions (18, 22, 23, 28, 29, 32, 33, 34, 35, 36, 40, 41, 43, 44, 46, 48, 49, 50, 53, 54, 55, 60, and 61) where the surface area covered by the NHP ecosystem reported by the MEE has decreased by more than 2000 ha. This result suggests intensive land use that has led to NHP degradation; for example, regions 29 and 18 have lost 15,035 and 14,616 ha, respectively. In regions 2, 3, 5, 10, 16, and 17, the CDT identifies an increase of more than 500 ha in the surface area covered by the NHP ecosystem concerning the values reported by the MEE. Particularly in region 17, the CDT finds an increase of 3042 ha under the NHP. This result suggests that in these regions, the NHP is regenerating.
In summary, the sum of all positive differences between the NHP surface’s value achieved with our analysis and the MEE value suggests a total gain of 7607 ha. The sum of all negative differences between the total area given by OBIA and the reported by MEE in each region suggests a total loss of 154,618 ha. In contrast, regarding the MEE report, the results suggest a total net loss of 147,129 ha (296,963–149,833) of NHP. However, the MEE only reported the existence of the NHP ecosystem.

4. Discussion

This study found that the NHP ecosystem is being transformed into AHP. The results show that the NHP ecosystem (296,964 ha) has been reduced by 50% (149,834 ha). However, the rate of changes in the extent of NHP and AHP cannot be clearly defined. There is a strong possibility that the MEE classified all paramos (NHP and AHP) together as one unit and reported it as NHP. Similarly, we do not rule out the possibility that the 50% loss in NHP occurred over a period that started before the baseline MEE map. Nevertheless, this result agrees with what several authors have suspected [2,32]. Among the factors that may be contributing to the transformation of NHP are (1) intensive land use or expansion of anthropogenic activity, (2) global climate change, and (3) limited knowledge about NHP regeneration, among others.
Intensive land use has modified a significant part of the páramo ecosystems in the study area. In 23 of the regions studied, the NHP ecosystem area decreased by more than 2000 ha. This decrease is because the paramos’ inhabitants, to exploit the soils of the páramo, introduce replacement systems of the NHP ecosystem (Figure 6).
In the páramos, there is a wide variety of production systems with tremendous social and environmental impacts on soil, water quality, and social-environmental dynamics [34]. Human-induced transformations, such as fire and clearing activities, have decreased the area covered by páramo ecosystems and increased their fragmentation [35].
The transformation of the natural paramo due to LULC changes affects soil carbon retention. When soil is cultivated, the A horizon is degraded, and the layer of organic remains disappears due to erosion [36]. When the soil is unproductive, it is abandoned and left to “recover,” leading to a secondary herbaceous páramo (Figure 6) with unpredictable results. However, very little is known about the recovery of the NHP and the páramo’s soils.
Although there may be a wide range of soil rehabilitation recommendations, it is essential to conduct a complete ecosystem analysis before site rehabilitation. The analysis must integrate the description of the developed soils and the slow pedogenic processes [16]. Thus, the soils’ study under páramo is indispensable to understand the páramo ecosystem’s functioning and optimize its ecological rehabilitation efficiently.
The páramos are among the most vulnerable ecosystems to global climate change. It is expected that the páramo temperature change will be lower than in other lower altitude regions. However, a 10% reduction in annual rainfall is expected in the páramo by 2040. The change in temperature can cause an altitudinal displacement of a vast number of species and, with them, of the bioclimatic band. This temperature change may reduce the surface covered initially by páramo ecosystems by 15 to 25% [34]. The growing transformation of the NHP, the increase in temperature, and the decrease in humidity reduce taxonomic and functional diversity, and the functioning and supply of ecosystem services that the páramo provides.
Soils in the southern Ecuadorian páramos are burned after three to four years of exploitation and left to recover. The burning of vegetation to grow grass is a widespread management practice in the region [4]. When burned and abandoned, the secondary HP ecosystem (AHP) regenerates through secondary succession, with limited functionality and reduced taxonomic diversity (Figure 6) [2,33,35]. There is little information on how soil properties under paramo ecosystems impact herbaceous and shrub species’ natural regeneration. However, this secondary HP may or may not lead to an NHP. The difficulty in the NHP regeneration is favored by the alteration of the pedogenesis of its andosol soil [12,13,18]. In some cases, the exploitation has been so intensive that only bare soil remains. In this condition, the ecosystem demands an active restoration with direct intervention on the degraded páramo ecosystem’s structure and characteristics to overcome stressors that prevent regeneration and ensure recovery [2,7,16,18,32,34]. However, few know how to or even if it is possible to regenerate the NHP, the resilience of NHP, and the resilience of andosols to continue supporting the NHP ecosystem. There are few studies on the páramos’ soil regarding the impact resulting from exploiting the páramos or their restoration after exploitation [2,4,18,32].
In our study, the CDT identified five regions where the NHP ecosystem is upgrading. In particular, region 17 shows an increase in the surface covered by HPN in 3000 ha. The CDT used spectral indices of vegetation as biophysical indicators of the LULC. This outcome suggests that in region 17, the NHP is regenerating. Few studies on the recovery of paramo ecosystems after the human intervention have ceased, limiting the implementation of restoration actions [5]. In páramo ecosystems, the geomorphological position regarding LU and management practices conditions the availability of carbon in the soil. It affects pedogenetic processes that compromise the regeneration of NHP. However, pedogenetic processes are prolonged and even blocked by specific environmental parameters [16]. The processes involved in the evolution of intensively exploited soils in the páramos have not yet been addressed. The nature of pedogenetic processes in the páramo remains unexplored. All of the NHP ecosystem’s compartments must be considered, particularly the soil, to know whether to intervene after LU ceases and restore the ecological communities’ functioning in the long term.

5. Conclusions

The intensive LU of the páramo with cycles of exploitation, abandonment, and regeneration in a secondary páramo transforms the NHP ecosystem. These exploitation practices, global climate change, and the lack of knowledge about the regeneration of the ecosystem and recovery of its soils threaten to substantially reduce the ecosystem’s area, functionality, and provision of the NHP ecosystem services. Analysis of the LULC database in the páramo ecosystem of Ecuador’s southern region shows that 39% (149,834/386,709) of the total study area remains as a native HP (NHP). However, 38% (146,829/386,709) have been transformed into anthropogenic HP (AHP), and 23% ((61,404 + 27,916 + 726)/386,709) remains FRS, WTR, and SNW. Overall, 386,709 ha is the sum of all classes’ area determined in our study (see Table 4). Regarding the MEE report, 50% (149,834 ha) of the surface area of the páramo ecosystem (296,964 ha) remain as a native HP. Among the factors that may be contributing to the transformation of NHP are (1) intensive land use or expansion of anthropogenic activity, (2) global climate change, and (3) limited knowledge about NHP regeneration, among others.
However, we found evidence of regeneration of the NHP, indicating its resilience. Additional studies to integrate new technological tools based on artificial intelligence algorithms, evolutionary processes, functional diversity, community dynamics, and ecosystem services of the NHP and AHP are necessary for exploring NHP’s response to global climate change and all páramo communities’ adaptability to such change. The future of NHP depends on the study of the factors that favor NHP regeneration. Thus, the páramos’ conservation depends on the new knowledge we can acquire by studying areas where NHP regeneration has been identified. This study will be primarily useful to scientists and decision-makers interested in páramo ecosystems in Ecuador’s southern region.

Author Contributions

V.J.G. and C.O.M. conceived and designed the experiments; M.V.R. and J.J.O. performed the experiment; V.J.G. and C.O.M. analyzed the data; A.C.R. and J.J.O. and C.D.A. contributed reagents/materials/analysis tools; V.J.G., C.O.M. and M.V.R. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors express their gratitude to the Vice-rectorate of Postgraduate Studies and Research of the National University of Chimborazo (Unach) through the group “Interdisciplinary Studies” and the project “Soil organic carbon evaluation and sequestration in Ecuadorian páramos ecosistem.”

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Acevedo, C.; Alarcon, L.; Miranda-Esquivel, D. Páramos Neotropicales como unidades biogeográficas. Biol. Trop. 2020, 68, 503–516. [Google Scholar]
  2. Hofstede, R.; Calles, J.; López, V. Los Páramos Andinos Qué Sabemos? Estado de Conocimiento Sobre el Impacto del Cambio Climático en el Ecosistema Páramo. 2014. Available online: https://portals.iucn.org/library/node/44760 (accessed on 17 July 2018).
  3. Poulenard, J.; Podwojewski, P.; Janeau, J.-L.; Collinet, J. Runoff and soil erosion under rainfall simulation of Andisols from the Ecuadorian Páramo: Effect of tillage and burning. Catena 2001, 45, 185–207. [Google Scholar] [CrossRef]
  4. López, S.; López-Sandoval, M.F.; Gerique, A.; Salazar, J. Landscape change in Southern Ecuador: An indicator-based and multi-temporal evaluation of land use and land cover in a mixed-use protected area. Ecol. Indic. 2020, 115, 106357. [Google Scholar] [CrossRef]
  5. Calderón-Loor, M.; Cuesta, F.; Pinto, E.; Gosling, W.D. Carbon sequestration rates indicate ecosystem recovery following human disturbance in the equatorial Andes. PLoS ONE 2020, 15, e0230612. [Google Scholar] [CrossRef]
  6. Binder, S.; Isbell, F.; Polasky, S.; Catford, J.A.; Tilman, D. Grassland biodiversity can pay. Proc. Natl. Acad. Sci. USA 2018, 115, 3876–3881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Nerlekar, A.N.; Veldman, J.W. High plant diversity and slow assembly of old-growth grasslands. Proc. Natl. Acad. Sci. USA 2020, 117, 18550–18556. [Google Scholar] [CrossRef] [PubMed]
  8. García, V.J.; Márquez, C.O.; Isenhart, T.M.; Rodríguez, M.; Crespo, S.D.; Cifuentes, A.G. Evaluating the conservation state of the páramo ecosystem: An object-based image analysis and CART algorithm approach for central Ecuador. Heliyon 2019, 5, e02701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Joslin, A.J.; Jepson, W.E. Territory and authority of water fund payments for ecosystem services in Ecuador’s Andes. Geoforum 2018, 91, 10–20. [Google Scholar] [CrossRef]
  10. Barta, B.; Mouillet, C.; Espinosa, R.; Andino, P.; Jacobsen, D.; Christoffersen, K.S. Glacial-fed and páramo lake ecosystems in the tropical high Andes. Hydrobiologia 2017, 813, 19–32. [Google Scholar] [CrossRef]
  11. Madriñán, S.; Cortés, A.J.; Richardson, J.E. Páramo is the world’s fastest evolving and coolest biodiversity hotspot. Front. Genet. 2013, 4, 192. [Google Scholar] [CrossRef] [Green Version]
  12. Ayala-Izurieta, J.E.; Márquez, C.O.; García, V.J.; Recalde, C.; Llerena, M.V.R.; Damián-Carrión, D.A. Land Cover Classification in an Ecuadorian Mountain Geosystem Using a Random Forest Classifier, Spectral Vegetation Indices, and Ancillary Geographic Data. Geoscience 2017, 7, 34. [Google Scholar] [CrossRef] [Green Version]
  13. Prat, C. The Soils of Ecuador. In World Soils Book Series; Springer Science and Business Media LLC.: Berlin, Germany, 2018; pp. 11–49. [Google Scholar]
  14. Hribljan, J.A.; Suárez, E.R.; Heckman, K.A.; Lilleskov, E.A.; Chimner, R.A. Peatland carbon stocks and accumulation rates in the Ecuadorian páramo. Wetl. Ecol. Manag. 2016, 24, 113–127. [Google Scholar] [CrossRef]
  15. Wiesmeier, M.; Hübner, R.; Barthold, F.K.; Spörlein, P.; Geuß, U.; Hangen, E.; Reischl, A.; Schilling, B.; Von Lützow, M.; Kögel-Knabner, I. Amount, distribution and driving factors of soil organic carbon and nitrogen in cropland and grassland soils of southeast Germany (Bavaria). Agric. Ecosyst. Environ. 2013, 176, 39–52. [Google Scholar] [CrossRef]
  16. Chenot, J.; Jaunatre, R.; Buisson, E.; Bureau, F.; Dutoit, T. Impact of quarry exploitation and disuse on pedogenesis. Catena 2018, 160, 354–365. [Google Scholar] [CrossRef]
  17. Carrillo-Rojas, G.; Silva, B.; Rollenbeck, R.; Célleri, R.; Bendix, J. The breathing of the Andean highlands: Net ecosystem exchange and evapotranspiration over the páramo of southern Ecuador. Agric. For. Meteorol. 2019, 265, 30–47. [Google Scholar] [CrossRef]
  18. Buytaert, W.; Sevink, J.; De Leeuw, B.; Deckers, J. Clay mineralogy of the soils in the south Ecuadorian páramo region. Geoderma 2005, 127, 114–129. [Google Scholar] [CrossRef]
  19. Rivera, J.I.; Bonilla, C.A. Predicting soil aggregate stability using readily available soil properties and machine learning techniques. Catena 2020, 187, 104408. [Google Scholar] [CrossRef]
  20. Ministerio del Ambiente del Ecuador. Sistema de Clasificación de los Ecosistemas del Ecuador Continental; Subsecretaría de Patrinmonio Natural: Quito, Ecuador, 2013.
  21. Saldarriaga, J.C.; Hoyo, D.Á.; Correa, M.Á. Evaluación De Procesos Biológicos Unitarios En La Remoción Simultánea De Nutrientes Para Minimizar La Eutrofización. DialnetUniriojaEs 2011, 15, 129–140. [Google Scholar]
  22. Severino, G.; De Bartolo, S.; Toraldo, G.; Srinivasan, G.; Viswanathan, H. Water Resources Research. Water Resour. Res. 2012, 53, 5998–6017. [Google Scholar]
  23. Ministerio del Ambiente del Ecuador. Metodología Para la Representación Cartográfica de los Ecossitemas del Ecuador Continental; Subsecretaría de Patrimonio Natural: Quito, Ecuador, 2012.
  24. Instituto Geográfico Militar. Instituto Geográfico Militar. Available online: http://www.igm.gob.ec/cms/files/cartabase/enie/ENIEV_C2.htm (accessed on 15 December 2013).
  25. Urbanski, J. Integration of GEOBIA with GIS for semi-automatic land cover mapping from Landsat 8 imagery. In Proceedings of the 5th GEOBIA conference, Thessaloniki, Poland, 21–24 May 2014. [Google Scholar]
  26. Yengoh, G.T.; Dent, D.; Olsson, L.; Tengberg, A.; Iii, C.J.T. The Use of the Normalized Difference Vegetation Index (NDVI) to Assess Land Degradation at Multiple Scales: Current Status, Future Trends, and Practical Considerations; Springer: New York, NY, USA, 2015. [Google Scholar]
  27. Gitelson, A.A.; Kaufman, Y.J.; Stark, R.; Rundquist, D. Novel algorithms for remote estimation of vegetation fraction. Remote. Sens. Environ. 2002, 80, 76–87. [Google Scholar] [CrossRef] [Green Version]
  28. Jiang, Z.; Huete, A.; Didan, K.; Miura, T. Development of a two-band enhanced vegetation index without a blue band. Remote Sens. Environ. 2008, 112, 3833–3845. [Google Scholar] [CrossRef]
  29. Storey, E.A.; Stow, D.A.; O’Leary, J.F. Assessing postfire recovery of chamise chaparral using multi-temporal spectral vegetation index trajectories derived from Landsat imagery. Remote Sens. Environ. 2016, 183, 53–64. [Google Scholar] [CrossRef] [Green Version]
  30. Nemmaoui, A.; Aguilar, M.A.; Aguilar, F.J.; Novelli, A.; Lorca, A.M.G. Greenhouse Crop Identification from Multi-Temporal Multi-Sensor Satellite Imagery Using Object-Based Approach: A Case Study from Almería (Spain). Remote. Sens. 2018, 10, 1751. [Google Scholar] [CrossRef] [Green Version]
  31. Cao, J.; Leng, W.; Liu, K.; Liu, L.; He, Z.; Zhu, Y. Object-Based Mangrove Species Classification Using Unmanned Aerial Vehicle Hyperspectral Images and Digital Surface Models. Remote Sens. 2018, 10, 89. [Google Scholar] [CrossRef] [Green Version]
  32. Dahik, C.Q.; Marín, F.; Arias, R.; Crespo, P.J.; Weber, M.; Palomeque, X. Comparison of Natural Regeneration in Natural Grassland and Pine Plantations across an Elevational Gradient in the Páramo Ecosystem of Southern Ecuador. Forests 2019, 10, 745. [Google Scholar] [CrossRef] [Green Version]
  33. Llambí, L.D.; Teresa, M.; Independiente, B. Building a Strategy for the Integrated Monitoring of High Mountain Ecosystems in Colombia. Biodivers Práctica 2019, 4, 150–172. [Google Scholar]
  34. Sarmiento, C.; Osejo, A.; Ungar, P. Páramos habitados: Desafíos para la gobernanza ambiental de la alta montaña en Colombia. Biodivers Práctica 2017, 2, 122–145. [Google Scholar]
  35. Wiesmeier, M.; Urbanski, L.; Hobley, E. Páramos Neotropicales como unidades biogeográficas. Glob. Chang. Biol. 2019, 54, 651–661. [Google Scholar]
  36. Gutiérrez-Salazar, P.; Medrano-Vizcaíno, P. The effects of climate change on decomposition processes in andean paramo ecosystem–synthesis, a systematic review. Appl. Ecol. Environ. Res. 2019, 17, 4957–4970. [Google Scholar] [CrossRef]
Figure 1. Ecuador’s southern region.
Figure 1. Ecuador’s southern region.
Agronomy 10 01922 g001
Figure 2. Workflow and methodology of land use and land cover classification.
Figure 2. Workflow and methodology of land use and land cover classification.
Agronomy 10 01922 g002
Figure 3. Typical characteristics of the LULC class identified from field testing and recognized in Lansat 8 Imagery. (a) Native herbaceous páramo (NHP); (b) anthropogenic herbaceous páramo (AHP); (c) forest (FRS); (d) water (WRT); and (e) snow (SNW).
Figure 3. Typical characteristics of the LULC class identified from field testing and recognized in Lansat 8 Imagery. (a) Native herbaceous páramo (NHP); (b) anthropogenic herbaceous páramo (AHP); (c) forest (FRS); (d) water (WRT); and (e) snow (SNW).
Agronomy 10 01922 g003
Figure 4. The classifier decision tree.
Figure 4. The classifier decision tree.
Agronomy 10 01922 g004
Figure 5. (a) The native herbaceous páramo reported in Ecuador’s ecosystem map [20]. (b) The native herbaceous páramo and anthropogenic herbaceous páramo predicted by the classification decision tree in the study area.
Figure 5. (a) The native herbaceous páramo reported in Ecuador’s ecosystem map [20]. (b) The native herbaceous páramo and anthropogenic herbaceous páramo predicted by the classification decision tree in the study area.
Agronomy 10 01922 g005
Figure 6. Possible evolution paths of taxonomic–functional diversity and ecosystem functioning services in different land-use scenarios [33].
Figure 6. Possible evolution paths of taxonomic–functional diversity and ecosystem functioning services in different land-use scenarios [33].
Agronomy 10 01922 g006
Table 1. Features used for classification [8].
Table 1. Features used for classification [8].
Object FeaturesDescription
Reflectance basic spectral information (Landsat 8 bands B2, B3, B4, and B8).The average of the ground REFLECTANCE values at the bottom of the atmosphere for all pixels within an object.
Spectral vegetation indexes derived from Landsat 8 B2, B3, B4, B5, B6, and B7 bands.EVI2, VARIg, y NBR2.
Topographic indexes derived from DEM.ALTITUDE y SLOPE.
Spectral vegetation indexes derived from Landsat 8 bands B2, B3, B4, B5, B6, and B7.
IndexFormulationDescription
VARIg: Visible Atmospherically Resistant
Vegetation Index Green.
( G R ) / ( G + R B ) Lineally sensitive to vegetation fraction, it
exhibits a good correlation with nitrogen contents [27].
EVI2: Enhanced Vegetation Index 2 2.5 × ( NIR R ) ( NIR + 2.4 × R + 1 ) Provide greater sensitivity in regions with high biomass while minimizing the soil’s influence and atmosphere [28].
NBR2: Normalized Burn Ratio 2 ( SWIR 1 SWIR 2 ) / ( SWIR 1 + SWIR 2 ) NBR2 is useful for the postfire recovery assessment [29].
Table 2. The resulting confusion matrix from the validation process in the study area.
Table 2. The resulting confusion matrix from the validation process in the study area.
Clase Predict
Current Class
Water (50)Forest (100)Native Herbaceous páramo (225)Anthropogenic Herbaceous páramo (225)Snow (10)
Water (35)350000
Forest (98)288080
Native herbaceous páramo (254)512221151
Anthropogenic herbaceous páramo (215)8042021
Snow (8)00008
Table 3. Performance of the classification decision tree classifying LULC in Ecuador’s southern region.
Table 3. Performance of the classification decision tree classifying LULC in Ecuador’s southern region.
WaterForestNative Herbaceous páramoAnthropogenic Herbaceous páramoSnow
Sensitivity (User accuracy)1.000.900.870.941.00
Specificity0.970.980.990.941.00
Presicion (Producer accuracy)0.700.880.980.900.80
Accuracy (Overall accuracy)0.980.960.940.941.00
Misclassification rate0.020.040.060.060.00
Informedness0.970.870.860.881.00
Markedness0.700.860.900.860.80
Matthew’s correlation0.840.870.840.880.89
Table 4. The distribution of land use/land cover for the complete study area and each region.
Table 4. The distribution of land use/land cover for the complete study area and each region.
RegionCoordinates
UTM—Zone 17
Southern Hemisphere
Surface Area (ha)
Native Herbaceous páramo
(MEE 2012)
Native Herbaceous páramoAnthropogenic Herbaceous páramoForestWaterSnow
XY
16870309455278042201213550
2667030947527807029007022050
3687030947527808372962320411570
46670309495278000000
56870309495278941363430431967930
6707030949527801337258922410
7687030951527853453167167740
870,0309515278346360307210
96870309535278000000
1070703095352780517171918295790
116870309555278000000
127070309555278102100000
136670309575278000000
1468703095752781507123699180
157070309575278168121288175160
1666703095952780540413152680
1768703095952782543295621716318350
18707030959527816,142152710,505369511490
19727030959527855934250413210
2066703096152783758571665394840
21687030961527813513416559690
22707030961527894486908421474832100
23727030961527892251366508041519160
247470309615278700800
256470309635278000000
26667030963527810800000
2768703096352781100000
2870703096352788010273585210195640
29727030963527817,62925946117515216900
30747030963527815572137553863060
3164703096552780229035330
32667030965527896242678723413338130
33687030965527819,22013,606642610293376
34707030965527880824121362317541310
357270309655278370339625599653200
36747030965527811,07325744753349920920
377670309655278596194343671810
386470309675278000000
396670309675278175814852053248
40687030967527826,46418,37412,4311590279341
41707030967527817,64313,526363933908840
427270309675278000000
437470309675278496748793712554220
44767030967527849431372773441380
45667030969527815051656682125070
46687030969527812,57710,03670859842486313
47707030969527829,11427,37747212760188547
48727030969527881064167235812562120
497470309695278392776414471061830
5076703096952784038373255357690
516670309715278000000
526870309715278575623313490
53707030971527812,80376315862295121941
5472703097152789896720453725424423
55747030971527815,1065627735047804890
5676703097152781709916234621600
576670309735278000000
586870309735278000000
597070309735278000000
60727030973527876775279230323410027
61747030973527812,73283024222190350
62767030973527882864222839130
Total296,964149,834146,82961,40427,916726
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

García, V.J.; Márquez, C.O.; Rodríguez, M.V.; Orozco, J.J.; Aguilar, C.D.; Ríos, A.C. Páramo Ecosystems in Ecuador’s Southern Region: Conservation State and Restoration. Agronomy 2020, 10, 1922. https://doi.org/10.3390/agronomy10121922

AMA Style

García VJ, Márquez CO, Rodríguez MV, Orozco JJ, Aguilar CD, Ríos AC. Páramo Ecosystems in Ecuador’s Southern Region: Conservation State and Restoration. Agronomy. 2020; 10(12):1922. https://doi.org/10.3390/agronomy10121922

Chicago/Turabian Style

García, Víctor J., Carmen O. Márquez, Marco V. Rodríguez, Jonathan J. Orozco, Christian D. Aguilar, and Anita C. Ríos. 2020. "Páramo Ecosystems in Ecuador’s Southern Region: Conservation State and Restoration" Agronomy 10, no. 12: 1922. https://doi.org/10.3390/agronomy10121922

APA Style

García, V. J., Márquez, C. O., Rodríguez, M. V., Orozco, J. J., Aguilar, C. D., & Ríos, A. C. (2020). Páramo Ecosystems in Ecuador’s Southern Region: Conservation State and Restoration. Agronomy, 10(12), 1922. https://doi.org/10.3390/agronomy10121922

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop