Next Article in Journal
Landscape Preferences of Visitors to the Danube Floodplains National Park, Vienna
Next Article in Special Issue
Vertical Boundary Mixing Events during Stratification Govern Heat and Nutrient Dynamics in a Windy Tropical Reservoir Lake with Important Water-Level Fluctuations: A Long-Term (2001–2021) Study
Previous Article in Journal
Seasonal Variations in the Biochemical Compositions of Phytoplankton and Transparent Exopolymer Particles (TEPs) at Jang Bogo Station (Terra Nova Bay, Ross Sea), 2017–2018
Previous Article in Special Issue
Methane Production and Oxidation in Mangrove Soils Assessed by Stable Isotope Mass Balances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Páramo Lakes of Colombia: An Overview of Their Geographical Distribution and Physicochemical Characteristics

by
Angela Zapata
1,*,
Carlos A. Rivera-Rondón
1,
Daivan Valoyes
1,
Claudia L. Muñoz-López
1,
Marino Mejía-Rocha
2 and
Jordi Catalan
3,4
1
Laboratorio de Limnología, Unidad de Ecología y Sistemática-UNESIS, Departamento de Biología, Pontificia Universidad Javeriana, Bogotá 110231, Colombia
2
Instituto de Salud Pública, Pontificia Universidad Javeriana, Bogotá 110231, Colombia
3
CREAF, Universitat Autònoma de Barcelona Campus, Edifici C, E-08193 Cerdanyola del Vallès, Spain
4
CSIC, Universitat Autònoma de Barcelona Campus, E-08193 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Water 2021, 13(16), 2175; https://doi.org/10.3390/w13162175
Submission received: 31 May 2021 / Revised: 13 July 2021 / Accepted: 20 July 2021 / Published: 9 August 2021

Abstract

:
The páramo lakes, Colombia, situated in the neotropical region, are of great value both as natural heritage and as water source to the most populated areas of the Andes, but are threatened by expanding agriculture, livestock, mining, and landscape fragmentation. Nonetheless, a general assessment of the lakes’ distribution and biogeochemical characteristics was lacking. We made a complete inventory of the Colombian páramo lakes and characterized their morphometry and water chemistry based on a survey of 51 lakes in the Eastern Cordillera. There are 3250 lakes distributed across 28 páramo complexes in Colombia, mainly located between 3600 and 4400 m a.s.l. The lakes are usually small (<10 ha) and shallow (<10 m). Most of them are slightly acidic (average pH~6), with high nutrient (total phosphorus ~0.6 mg/L, total nitrogen ~1 mg/L) and total organic carbon (~6 mg/L) concentrations and low oxygen (~3.5 mg/L) at the bottom. Water chemistry varies according to two main independent gradients related to watershed bedrock geology and trophic state. Global change pressures may be challenging the preservation of these unique ecosystems. Increasing the protection of more lake watersheds should reduce these potential impacts by mitigating negative synergies with local pressures.

1. Introduction

The páramo aquatic ecosystems are a unique element of the planet’s natural heritage and are socially crucial because they provide water to the most densely populated areas of the neotropical region of Colombia [1]. The páramo extends across high mountains from Ecuador to Venezuela. In these areas, air temperature is relatively constant throughout the year and relative humidity high. These conditions have favored the formation of peat and a large number of wetlands and lakes. The páramo usually presents deep and organic-rich soils. In Colombia, the Andes branch into three mountain ranges in which there are differentiated páramo complexes (Figure 1) above 3000 m altitude. The three mountain ranges show different amounts of precipitation (1000–4500 mm) [2]. The largest area of páramo ecosystems in Colombia (70%) is in the Eastern Cordillera followed by the Central Cordillera [3]. The Western Cordillera is much lower in altitude and, therefore, holds fewer complexes. Tertiary and Quaternary sedimentary rocks characterize the Eastern Cordillera, while the Central combines not only sedimentary material but also plutonic rocks and a history of high volcanism [4]. Paramo vegetation is dominated by grasses (Festuca and Calamagrostis spp.) and rosette plants of the Asteraceae family (Espeletia spp.), which are missing from the central and southern Andes [5]. Páramo complexes are “sky islands” of the northern Andes and host endemic biodiversity [6]. Despite their extraordinary natural and societal relevance, there has been limited research on Colombian páramo lakes.
Existing studies have primarily focused on descriptions of the aquatic biota and relationships with the environment e.g., [7,8,9,10]. Very few studies have analyzed functional aspects such as primary production and trophic food webs [11,12,13,14]. However, extensive surveys of the assessment of the general limnological characteristics of these lakes are missing. Such studies are required to develop a reference framework for conserving the natural heritage and evaluating global change pressures, including climate shift, atmospheric pollution, and landscape fragmentation [15]. In the present study, we aim to cover the gap by making an inventory of the Colombian páramo lakes and characterizing the water chemistry based on a survey of lakes in the Eastern Cordillera.
The Eastern Cordillera holds the greatest páramo extension in Colombia, although it is also the one with the highest human landscape transformation, between 49% and 74% in some areas during the last two centuries [3,16]. Potato crops and livestock are the two activities that historically had the greatest impact on the páramo ecosystems [17], resulting in a significant landscape fragmentation [3]. Compared to the temperate zone, the tropical high-mountain sustains more population, and, thus, some parts have been submitted to local human pressures for centuries. Páramo ecosystems store a significant amount of carbon in their vegetation and soils [18], offering two key services in the current climate change scenario; carbon sequestration and water regulation. Since they are sensitive to disturbances, an understanding of the natural variability of their aquatic ecosystems could have significant implications for land management decisions.
The protection of páramo areas in Colombia is relatively recent, mostly occurring in the late 1970s. There are still many lakes outside protected areas and, even in them, many watersheds have had some degree of intervention in the last 50 years. Therefore, we considered lakes both in protected and non-protected areas in our survey and planned to discuss whether protection measures are sufficient at present and which are the main threats under ongoing global change [19]. On the other hand, this work provides a baseline of the physical and chemical characteristics of these high mountain lakes. This information, in addition to being a reference point for monitoring in the scenario of global change, will allow the identification of páramo complexes that deserve greater attention in the definition of management policies and ecosystems that have a high value for their conservation.

2. Methods

2.1. Distribution of Colombian High-Mountain Lakes

The Colombian high-mountain lakes are clustered in isolated páramo complexes distributed across the mountain ranges and covering 2,900,000 km2 (Figure 1). Using available 30-m resolution digital maps [20,21,22,23], we extracted layers of páramo complexes, protected areas, and lakes and intersected them to obtain lake inventories and statistics for each páramo complex. ArcGIS® software by Esri (Redlands, CA, USA) was used to perform all cartographic spatial analyses.

2.2. Eastern Cordillera Sampling and Analyses

The Eastern Cordillera extends between approximately 1°05′ and 8°40′ N and is predominantly tectonic in origin [24]. The lithology predominantly dates back to Cretaceous and Eocene periods, consisting of sedimentary rocks of continental or transitional marine-continental origin, with some Jurassic plutonic rocks and Quaternary modeling at high altitudes [25]. We sampled 51 lakes at the Eastern Cordillera >0.3 ha and above 3000 m a.s.l., selected to cover variation across altitude, geology, protection level, and land use. The lake survey was performed during the drought period (between January and March 2017), when radiation is the highest, and the effect of dilution of ions and nutrients caused by rainfall is reduced. The selected lakes were mainly located in areas with current protection status (86%). There is no detailed information on land use in the lake basins, but field observations indicate that among the lakes with protection status, 35% present some type of current or very recent affectation (dams, introduced trout, presence of cattle, and vestiges of agriculture in the basin).
Lake depth was measured using an acoustic probe and GPS in points evenly distributed in each lake. A minimum of 30 points was measured in the smaller lakes (<1 ha) and at least 100 points in the larger ones. The lake bathymetry was obtained using kriging interpolation. Lake area was calculated from orthophoto maps and total volume by trapezoidal integration of sections: 0.2-m layers for lakes with depth <2 m, 0.5-m for 5 m depth, and 1-m for the deepest lakes.
The transparency of the water was measured using the Secchi disk depth (SD) and a photic index calculated as the ratio between the photic zone depth (i.e., SD × 2.3) and the maximum depth (Zmax). Vertical profiles of temperature and oxygen were performed, and surface water was sampled in the central area of the lake for chemical analyses of alkalinity, Ca2+, Mg2+, Mn2+, Na+, K+, Cl, total phosphorus (TP), SO42−, PO43−, total organic carbon (TOC), total Kjeldahl nitrogen (TKN), NH4+, NO3, NO2, and SiO2. Chemical analyses were performed following standard methods for the examination of water and wastewater [26]. The main environmental gradients were assessed by principal component analysis (PCA). The environmental variables were transformed using Log10 or square root. Temperature, altitude and pH were not transformed because they already had a central distribution. PCA is an ordination technique that arranges environmental variables and samples in a space with a reduced number of dimensions. The PCA allows analysis of environmental gradients, correlations between variables and the distribution of the samples with respect to the set of gradients of all the variables. The set of variables that has a greater weight in one of the dimensions (axes) can be interpreted together as a response to a latent factor.

3. Results

3.1. Distribution of Páramo Lakes

The analysis of digital maps showed that there are 3250 lakes in 28 out of 36 Colombian páramo complexes (Table 1), with 71% in protected areas. The highest number occurs in the Eastern Cordillera (44%), 39.6% in the Central Cordillera, and 13.7% in the Sierra Nevada de Santa Marta. The highest density of lakes occurs in Los Nevados, Sierra Nevada de Santa Marta, Las Hermosas, and Pisba complexes 0.003, 0.0029, 0.0025, and 0.0025 lakes ha−1, respectively).
The lake’s altitude ranges from 3000 to 5100 m.a.s.l., showing a bimodal distribution, with a primary mode at 3700 m and a secondary at 4300 m (Figure 2). The highest mode corresponds mainly to the Cocuy, Sierra Nevada de Santa Marta, and Los Nevados complexes.
Lake area ranges between 0.005 and 217 ha (Figure 2). The distribution is unimodal, with the mode slightly lower than 1 ha. The largest lakes mainly occur in three páramo complexes, Chiles-Cumbal, Sotará, and Cocuy.

3.2. Characteristics of Eastern Cordillera Lakes

The altitudinal distribution of the lakes surveyed in the Eastern Cordillera followed the primary mode of the whole set of páramo lakes (Figure 3), whereas the lake area showed a larger median (~2 ha), with centered distribution and a tail of a few large lakes (Table 2, Figure 3). There is no spatial pattern in the distribution of lakes according to altitude or area (Figure 4). Maximum depth was usually <10 m, only few lakes showed >20 m. Likewise, the mean depth was low (median < 2 m). Consequently, the lake volume was relatively small compared to the lake surface. Morphometry differences between páramo complexes were irrelevant, except in the Sumapaz complex, where the lakes were exceptionally shallow (Figure 4).
Water transparency was low (average SD 3.5 m). However, the photic zone extended to bottom in 65% of the lakes because of the shallowness (Table 2). The lakes located in the northern zone of the survey (Pisba and Guantiva-Russia páramos) showed the highest water transparency (Figure 4).
Surface water temperature was usually <15 °C. The average temperature difference between surface and bottom was 1.2 °C, but in the deepest lakes, it reached up to 4 °C. The oxygen saturation was high at surface (>78%) but highly variable at bottom, ranging from oversaturation to anoxia (Table 2).
Alkalinity was low in most lakes (<200 µeq/L), except in a few with values ~1 meq/L. The dominant cation was Ca2+ but with relatively low values (Table 2). Consequently, lake waters showed low buffering capacity, and pH was commonly well below 7, but with only a few cases of extremely low values (<5). Total organic carbon was generally relatively high (~6 mg/L), although it did not correlate significantly with water acidity. TOC did not show a spatial pattern (Figure 4).
Nitrogen concentration was low. The largely dominant dissolved inorganic form was NH4+ (Table 2), usually accounting for more than 50% of the TKN. TP and PO43− concentrations were high in most of the investigated complexes, indicating general mesotrophic and eutrophic conditions, with a few hypereutrophic exceptions. TP presented higher values in the páramos of the southern area sampled (Cruz Verde-Sumapaz and Chingaza).
The PCA summarized two latent factors of similar contribution to the water characteristics (Figure 5). One of the factors related to bedrock geology, which provides the major water chemicals and buffering capacity (i.e., pH, alkalinity, Ca2+, Mg2+, and SO42−). The second factor encompassed many variables defining the trophic state. Phosphorus variables showed the highest contribution to this factor. Variation in both factors was relatively high in every páramo complex. Lakes in Pisba and Guantiva-La Rusia were less productive and tended to hold softer waters; Tota-Rabanal, the hardest waters; and Sumapaz and Guerrero were he lakes of the highest trophic state.

4. Discussion

4.1. Distribution of Páramo Lakes

The altitudinal distribution of páramo lakes responds to the glacial modeling of the Andes of Colombia. During the last glacial maximum, glaciers reached 4300 m a.s.l. [27,28], whereas the formation of lake basins at lower altitudes can be traced back to older cold periods [29]. These dual erosive glacial phases could be the explanation for the bimodal lake distribution found. The páramo complexes, where glaciers still exist today, are the ones with lakes at the higher altitudes, suggesting a more recent basin formation. These lakes show small and silted basins. Nonetheless, there is no significant general relationship between altitude and lake number or size.
The complexes of Las Hermosas, Sierra Nevada de Santa Marta, and Los Nevados are of greater interest from the conservation perspective because of the higher number of lakes than in other complexes and the large proportion with protection status. In contrast, a small number of lakes are part of protected areas in the complexes of Cruz Verde-Sumapaz, Santurbán, and Pisba.

4.2. Eastern Cordillera Lake Characteristics

Generally, glacier modeling in temperate mountain regions has caused relatively deep basins compared to the lake area, in contrast to ice modeling in subarctic areas [30]. The average depth of lakes in the Pyrenees is around 15–20 m [31]. Lakes in the Andean equatorial region tend to show lower average depth, e.g., around 10 m in the Cajas massive in Southern Ecuador [32]. Existing information indicates that the lake aspect ratio for the Colombian Andes tends to be closer to subarctic lakes than those of the temperate zone [15]. Our survey confirms this tendency. The reason may be related to an erosive process in a more open high-mountain landscape than in other ranges with narrower valleys.
Unlike most alpine high-mountain watersheds, the páramo is covered by dense shrubby vegetation [33]. High air humidity throughout most of the year and average temperature values between 3–13 °C [34,35] promote extensive marsh vegetation around the lakes, which are limited by an extensive belt of terrestrial–aquatic transition zone [35]. Consequently, the organic matter load is potentially high, explaining the consistently high TOC values in the lakes studied [36,37]. This feature also makes the páramo lakes more similar to tundra lakes than temperate alpine lakes.
The Secchi disk values indicated that the transparency of the studied lakes in the Eastern Cordillera is very low compared to that in the temperate zone [30,38] and mountain lakes from Southern Ecuador [37]. However, despite the low water transparency, the light penetration is sufficient for macrophyte and algal growth at the bottom of most of the lakes due to the shallowness. This fact has major implications for the ecosystem as it increases the lake’s primary productivity [39,40].
The investigated lakes showed pH values mainly in the acidic range, with low buffer capacity because of the low cation concentration [41]. TOC and pH were not correlated (r = 0.1, p = 0.5, n = 60), indicating a low contribution of organic acids to acidity, which was primarily related to the inorganic ionic balance. All in all, the high diversity in TOC concentration among lakes merits further investigation about the TOC specific nature and the process determining its levels.
Total nitrogen and phosphorus levels were significantly higher than those reported for high-mountain lakes in other areas of the world [42,43]. Ammonium accumulation at levels above nitrate is particularly remarkable, which may point to some limitations for nitrification or a continuous flux from deep layers and sediments. This issue, and the nitrogen cycle as a whole, considering atmospheric deposition, watershed loadings, and fixation, deserve further research in these lakes [44].
The proportion of PO43− respective to TP fluctuates strongly but is remarkably high (>20%) in many lakes. This feature and the relatively high nitrogen levels points to a light limitation in some of these lakes. Thus, there is a significant relationship between the ratio PO43−/TP and the TOC (r = 0.32, p = 0.02, n = 51). In any case, TP is high, indicating mesotrophic to eutrophic conditions, in contrast to the usual oligotrophic character of high-mountain lakes [10]. On the one hand, the relatively low oxygen values at the bottom may facilitate phosphorus diffusion from recycling in the sediments [45]. On the other hand, the shallowness of the lakes fosters the relevance of the littoral and bottom primary producers and the influence of thycoplankton in open waters [46,47]. Although the characteristics pointed here refer to the páramo lakes of the Eastern Cordillera, there is no reason to expect that they will not apply to the whole set of páramo lakes in Colombia.

4.3. Páramo Lakes in the Context of Global Change

The main factors associated with global change that could affect the páramo lakes are the increase in air temperature, changes in the rainfall regime, changes in land use, and atmospheric transport of pollutants. These processes have been widely documented in remote mountain lakes [48,49,50]. The potential impacts of these processes range from in-lake moderate changes (e.g., water column stability, oxygen demand) to large watershed and external loading modifications (e.g., water level, organic loads, pollutant accumulation), eventually leading to marked changes in communities and food webs (Figure 6).
Studies have projected an increase in air temperature of between 3 and 4 °C in the Andes high mountains [51], resulting in some of these lakes possibly exceeding 20 °C during periods of high radiation. Higher temperatures will accelerate organism metabolism and, consequently, biogeochemical processes. Respiration is more sensitive than photosynthesis to temperature, but higher recycling may result in a higher trophic state of the lakes [52]. Higher biomass and metabolism will increase oxygen demand [53], and warmer waters and stronger stratification may limit the supply [54]; thus, the lakes may become more similar to tropical lakes of lower altitudes [55]. Oxygen values in hypolimnion already tend to be low in the lakes studied, so an increase in hypoxia will enhance phosphorus mobilization. Lakes that currently have the highest phosphorus values could experience significant water quality deterioration in the coming decades. This is the case of lakes near the main urban centers of Colombia. These systems, due to their proximity to the city of Bogotá, presented a history of high affectation during the last two centuries and, therefore, could have a low capacity to mitigate the impact of the increase in temperature. On the other hand, during the dry season all these effects can be enhanced by substantially increasing the metabolism of the lakes. Since our study is based on a sampling carried out during the dry period, the chemical composition of the lakes in future scenarios, with longer periods of drought, could be similar to the more enriched lakes that were observed. During the rainy season, dilution could reduce the concentration of nutrients, but could increase the contribution of organic matter from the basin.
Climate change may also affect rainfall seasonality, with more extended and more frequent periods of drought expected [1]. Our results show that the lakes of the Eastern Cordillera are very shallow ecosystems, so that longer periods of drought will have a great impact. Even though there are a high number of small bodies of water, the shallow depth determines that the water storage capacity will be low. Aquatic plant communities, especially those that inhabit the terrestrial–aquatic transition and submerged macrophytes without rapid regeneration, will be quickly affected by these changes [56,57]. Most of the lakes studied showed a photic layer that reaches the sediment, so the reduction in depth could increase the role of algal productivity in the deeper lakes, while the shallower lakes could suffer a progressive invasion of littoral plants. More importantly, páramo peatlands may experience periods of high release of organic material to the lakes, enhancing lake metabolism and oxygen demand. The shifting conditions may last for long periods, and human land use in the watersheds may enhance the processes. The long-range transport of pollutants to these lakes, another component of global change, is largely unknown [15]. The lake watershed protection should reduce undesirable synergies between global change and local land use and provide better stewardship of any unavoidable shift produced by climate transitions [58].

4.4. Páramo Lake Conservation

The results of the lake chemical analyses did not show a pattern, from which it could be concluded that a low lake trophic state is related to protection status (Figure 7). Despite the imbalance in the number of lakes considered between conserved and non-conserved areas, the high variability and the high concentration of TP and TOC in lakes of some National Parks such as Chingaza and Sumapaz suggests that this result shows more the degree of affectation suffered by the different complexes than the current conservation status. The System of Protected Areas in Colombia is relatively recent. Páramo protection is mainly under the jurisdiction of the national park system (33%), and smaller areas are managed by municipal or private regional reserves (12%). The majority of the páramo ecosystems were declared protected in 1977 as part of the National Natural Parks System of Colombia (decree 2811–1974 of the Republic of Colombia, 1975). Many of these areas previously had agricultural, livestock, and mining activities. Currently, most of them still have some type of agricultural management with low impact, such as extensive livestock. Surveying the past is possible with paleolimnological techniques [59]. The study of lake sediments can evaluate the effectiveness of the watershed conservation status and define reference conditions to manage human land use and climate interaction (e.g., land erosion and siltation of the lakes).
Strategic planning should also consider increasing the protected areas. Páramo ecosystems are a unique natural heritage for biodiversity. Fish introduction is one of the leading food web alterations in high-mountain lakes worldwide [60], including trout in the Andean tropical lakes [61]. Protection should reduce trout spreading and other undesirable fish species. Furthermore, if climate change modifies water availability, hydrological interventions (e.g., building dams, diverting streams) could be a pressure that prevention and planning should anticipate. The results presented here constitute a baseline for the definition of the trophic status of the lakes and for their monitoring. The shallow depth and the high content of organic matter and nutrients of most ecosystems suggest that they are ecosystems that require specific actions even in protected areas. The total elimination of livestock, the control of the introduction of trout and the assessment of the real potential of some lakes to provide water to local communities are actions that should be implemented urgently. The regular measurement of the water transparency, TP, TKN and the oxygen concentration in the hypolimnion would allow a response from these ecosystems to these actions during the progression of climate change.

5. Conclusions

The páramo lakes of Colombia are relatively small and fairly shallow lakes with limnological characteristics (e.g., low buffering capacity, high organic content, high nutrient concentration, and low oxygen) that make them extremely sensitive to local environmental pressures in their watersheds and during ongoing climate change. It becomes an urgent matter to increase the protected area and focus future research on assessing whether the current management of the watersheds is enough to guarantee a good ecological status and respond to challenging shifting global conditions.

Author Contributions

Conceptualization, C.A.R.-R., A.Z. and J.C.; methodology, A.Z., C.A.R.-R.; D.V., C.L.M.-L. and M.M.-R.; data analysis: A.Z., C.A.R.-R.; D.V., C.L.M.-L., M.M.-R. and J.C.; writing—original draft preparation, A.Z., C.A.R.-R., D.V., C.L.M.-L. and M.M.-R.; writing—review and editing, A.Z., C.A.R.-R. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Ministerio de Ciencia, Tecnología e Innovación of Colombia (Colciencias 50000-050-2016) and Pontificia Universidad Javeriana (IDPPTA 6789 and ID PPTA 7066). C. Muñoz-López was supported by a PhD scholarship MinCiencias Colciencias-Colfuturo (No. 779 2017-Regiones Boyacá).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Thanks to Parques Nacionales Naturales de Colombia (Resolución No. 141-241016) and to the staff of the PNN Chingaza, PNN Sumapaz, SFF Guanentá Alto Río Fonce, SFF lguaque, and PNN Pisba, for their invaluable help during the fieldwork. Thanks to the project fieldwork team: Y. Caleño, D. Castillo, M. Guerrero, J. Jimenez, and L. Garzón for their invaluable help exploring the mountains lakes.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Herzog, S.; Martínez, R.; Jørgensen, P.; Tiessen, H. (Eds.) Climate Change and Biodiversity in the Tropical Andes; Inter-American Institute for Global Change Research (IAI); Scientific Committee on Problems of the Environment (SCOPE): Paris, France, 2011; p. 348. [Google Scholar]
  2. Van der Hammen, T.; Otero, J. Los páramos: Archipiélagos terrestres en el norte d elos Andes. In Atlas de Páramos de Colombia; Morales, M., Otero, J., Van der Hammen, T., Torres, A., Cadena, C., Pedraza, C., Rodríguez, N., Franco, C., Betancourt, J.C., Olaya, E., et al., Eds.; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia, 2007; p. 208. [Google Scholar]
  3. Cabrera, M.; Ramirez, W. (Eds.) Restauración Ecológica de los Páramos de Colombia: Transformación y Herramientas para su Conservación; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia, 2014; p. 296. [Google Scholar]
  4. Gómez, J.; Montes, N.E. Mapa Geológico de Colombia 2020, Escala 1:1,000,000; Servicio Geológico Colombiano, 2 Hojas: Bogotá, Colombia, 2020; p. 2. [Google Scholar]
  5. Rangel, J.; Pinto-Zárate, J. Colombian Páramo Vegetation Database (CPVD)–The database on high Andean páramo vegetation in Colombia. Biodivers. Ecol. 2012, 4, 275–286. [Google Scholar] [CrossRef]
  6. Peyre, G.; Lenoir, J.; Karger, D.N.; Gomez, M.; González, A.; Broennimann, O.; Guisan, A. The fate of páramo plant assemblages in the sky islands of the northern Andes. J. Veg. Sci. 2020, 31, 967–980. [Google Scholar] [CrossRef]
  7. Aranguren-Riaño, N.; Guisande, C.; Ospina, R. Factors controlling crustacean zooplankton species richness in Neotropical lakes. J. Plankton Res. 2011, 33, 1295–1303. [Google Scholar] [CrossRef]
  8. Donato, J.; González, L.E.; Rodríguez, C.L. Ecología de dos Sistemas Acuáticos de Páramo; Academia Colombiana de Ciencias Exactas, Físicas y Naturales: Bogotá, Colombia, 1996. [Google Scholar]
  9. Meneses-Ortegón, L.; Herrera, Y. Bacterioplancton de tres humedales altoandinos de la cordillera Oriental de Colombia. Rev. Biota Colomb. 2015, 16, 1–10. [Google Scholar]
  10. Ramos, C.; Cárdenas, N.; Herrera, Y. Caracterización de la comunidad de Macrofitas acuáticas en lagunas del Páramo de La Rusia (Boyacá-Colombia). Cienc. Desarro. 2013, 4, 73–82. [Google Scholar]
  11. Donato, J.; Jimenez, H.; Reynolds, C. Limnology in El Dorado: Some surprising aspects of the regulation of phytoplankton productive capacity in a high-altitude Andean lake (Laguna de Guatavita, Colombia). Rev. Biol. Trop. 2012, 60, 967–980. [Google Scholar] [CrossRef] [PubMed]
  12. Pinilla-A, G.A.; Donato-R, J.C.; Rivera-R, C.A. Photosynthetic efficiency of phytoplankton in a tropical mountain lake. Caldasia 2006, 28, 57–66. [Google Scholar]
  13. Zapata, A.; Rivera Rondón, C.A.; Donato, J. Dynamics of photosynthetic pigments in an Andean lake in Colombia. Lakes Reserv. Res. Manag. 2006, 11, 29–38. [Google Scholar] [CrossRef]
  14. Zapata, A. Variaciones Diarias y Mensuales de la Productividad Primaria en un Lago Andino (Laguna de Guatavita. Cundinamarca). Master’s Thesis, Pontificia Universidad Javeriana, Bogotá, Colombia, 2001. [Google Scholar]
  15. Catalan, J.; Donato Rondón, J.C. Perspectives for an integrated understanding of tropical and temperate high-mountain lakes. J. Limnol. 2016, 75, 19. [Google Scholar] [CrossRef] [Green Version]
  16. Spehn, E.M.; Liberman, M.; Korner, C. Land Use Change and Mountain Biodiversity; CRC Press: Boca Raton, FL, USA, 2006. [Google Scholar]
  17. Sánchez-Cuervo, A.M.; Aide, T.M.; Clark, M.L.; Etter, A. Land Cover Change in Colombia: Surprising Forest Recovery Trends between 2001 and 2010. PLoS ONE 2012, 7, e43943. [Google Scholar] [CrossRef] [Green Version]
  18. Tonneijck, F.H.; Jansen, B.; Nierop, K.G.J.; Verstraten, J.M.; Sevink, J.; De Lange, L. Towards understanding of carbon stocks and stabilization in volcanic ash soils in natural Andean ecosystems of northern Ecuador. Eur. J. Soil Sci. 2010, 61, 392–405. [Google Scholar] [CrossRef]
  19. Steffen, W.; Richardson, K.; Rockström, J.; Cornell, S.E.; Fetzer, I.; Bennett, E.M.; Biggs, R.; Carpenter, S.R.; de Vries, W.; de Wit, C.A.; et al. Planetary boundaries: Guiding human development on a changing planet. Science 2015, 347, 1259855. [Google Scholar] [CrossRef] [Green Version]
  20. Instituto Geográfico Agustín Codazzi. Datos Abiertos Cartografía y Geografía. Available online: https://geoportal.igac.gov.co/contenido/datos-abiertos-cartografia-y-geografia (accessed on 11 November 2019).
  21. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. Una Nueva Cartografía de los Páramos Colombianos: Diversidad, Territorio e Historia; Instituto Alexander von Humboldt: Bogotá, Colombia, 2013. [Google Scholar]
  22. Registro Único de Áreas Protegidas (RUNAP). Runap en Cifras. Available online: https://runap.parquesnacionales.gov.co/cifras (accessed on 22 July 2021).
  23. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. Base de Datos Geográfica de la Propuesta de Límite de los Complejos de Páramo a Escala 1:25.000. Available online: http://geonetwork.humboldt.org.co/geonetwork/srv/spa/catalog.search#/metadata/36139b13-b15e-445e-b44d-dd7a5dbe8185 (accessed on 22 July 2021).
  24. Caballero, V.; Parra, M.; Bohorquez, A.R.M. Levantamiento de la cordillera oriental de Colombia durante el eoceno tardío-oligoceno temprano: Proveniencia sedimentaria en el sinclinal de nuevo mundo, cuenca Valle medio del Magdalena. Bol. Geol. 2010, 32, 45. [Google Scholar]
  25. Gómez, J.; Nivia, A.; Montes, N.E.; Jiménez, D.M.; Tejada, M.L.; Sepúlveda, M.J.; Osorio, J.A.; Gaona, T.; Diederix, H.; Uribe, H.; et al. Mapa Geológico de Colombia. Escala 1:2’800.000, 2nd ed.; INGEOMINAS: Bogotá, Colombia, 2007.
  26. APHA; AWWA; WEF. Standard Methods for the Examination of Water and Wastewater, 21st ed.; American Public Health Association: Washington, DC, USA, 2005; p. 1325. [Google Scholar]
  27. Thouret, J.-C.; Van der Hammen, T.; Salomons, B.; Juvigné, E. Paleoenvironmental Changes and Glacial Stades of the Last 50,000 Years in the Cordillera Central, Colombia. Quat. Res. 1996, 46, 1–18. [Google Scholar] [CrossRef]
  28. Helmens, K. The Quaternary glacial record of the Colombian Andes. Dev. Quat. Sci. 2004, 2, 115–134. [Google Scholar] [CrossRef]
  29. Angel, I.; Guzman, O.; Carcaillet, J. Pleistocene Glaciations in the Northern Tropical Andes, South America (Venezuela, Colombia and Ecuador). Cuad. Investig. Geogr. 2017, 43, 20. [Google Scholar] [CrossRef] [Green Version]
  30. Catalan, J.; Curtis, C.J.; Kernan, M. Remote European mountain lake ecosystems: Regionalisation and ecological status. Freshw. Biol. 2009, 54, 2419–2432. [Google Scholar] [CrossRef]
  31. Rivera-Rondón, C.A.; Catalan, J. Diatom diversity in the lakes of the Pyrenees: An iconographic reference. Limnetica 2017, 36, 127–396. [Google Scholar] [CrossRef]
  32. Mosquera, P.V.; Hampel, H.; Vázquez, R.F.; Alonso, M.; Catalan, J. Abundance and morphometry changes across the high-mountain lake-size gradient in the tropical Andes of Southern Ecuador. Water Resour. Res. 2017, 53, 7269–7280. [Google Scholar] [CrossRef] [Green Version]
  33. Buytaert, W.; Célleri, R.; De Bièvre, B.; Cisneros, F.; Wyseure, G.; Deckers, J.; Hofstede, R. Human impact on the hydrology of the Andean páramos. Earth-Sci. Rev. 2006, 79, 53–72. [Google Scholar] [CrossRef]
  34. Castro, C.; Agualimpia, Y.; Sánchez, F. Modelo climático de los páramos de la cordillera Oriental colombiana aplicado a regímenes de temperatura del suelo. Perspect. Geogr. 2016, 21, 33–62. [Google Scholar] [CrossRef] [Green Version]
  35. Schmidt-Mumm, U.; Vargas Ríos, O. Comunidades vegetales de las transiciones terrestre-acuáticas del páramo de Chingaza, Colombia. Rev. Biol. Trop. 2012, 60, 35–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Curtis, C.J.; Juggins, S.; Clarke, G.; Battarbee, R.W.; Kernan, M.; Catalan, J.; Thompson, R.; Posch, M. Regional influence of acid deposition and climate change in European mountain lakes assessed using diatom transfer functions. Freshw. Biol. 2009, 54, 2555–2572. [Google Scholar] [CrossRef]
  37. Van Colen, W.R.; Mosquera, P.; Vanderstukken, M.; Goiris, K.; Carrasco, M.-C.; Decaestecker, E.; Alonso, M.; León-Tamariz, F.; Muylaert, K. Limnology and trophic status of glacial lakes in the tropical Andes (Cajas National Park, Ecuador). Freshw. Biol. 2017, 62, 458–473. [Google Scholar] [CrossRef] [Green Version]
  38. Rivera-Rondón, C.A.; Catalan, J. The ratio between chrysophycean cysts and diatoms in temperate, mountain lakes: Some recommendations for its use in paleolimnology. J. Paleolimnol. 2017, 57, 273–285. [Google Scholar] [CrossRef]
  39. Rooney, N.; Kalff, J.; Habel, C. The Role of Submerged Macrophyte Beds in Phosphorus and Sediment Accumulation in Lake Memphremagog, Quebec, Canada. Limnol. Oceanogr. 2003, 48, 1927–1937. [Google Scholar] [CrossRef] [Green Version]
  40. Sayer, C.D.; Burgess, A.M.Y.; Kari, K.; Davidson, T.A.; Peglar, S.; Yang, H.; Rose, N. Long-term dynamics of submerged macrophytes and algae in a small and shallow, eutrophic lake: Implications for the stability of macrophyte-dominance. Freshw. Biol. 2010, 55, 565–583. [Google Scholar] [CrossRef]
  41. Psenner, R.; Catalan, J. Chemical composition of lakes in crystalline basins: A combination of atmospheric deposition, geologic background, biological activity and human action. In Limnology Now: A Paradigm of Planetary Problems, 1st ed.; Margalef, R., Ed.; Elsevier Science: Amsterdam, The Netherlands, 1994; pp. 255–314. [Google Scholar]
  42. Rivera-Rondón, C.A.; Catalan, J. Diatoms as indicators of the multivariate environment of mountain lakes. Sci. Total Environ. 2020, 703, 135517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Cantonati, M.; Lange-Bertalot, H. Diatom monitors of close-to-pristine, very-low alkalinity habitats: Three new eunotia species from springs in nature parks of the south-eastern alps. J. Limnol. 2011, 70, 209–221. [Google Scholar] [CrossRef]
  44. Canham, C.D.; Pace, M.L.; Weathers, K.C.; McNeil, E.W.; Bedford, B.L.; Murphy, L.; Quinn, S. Nitrogen deposition and lake nitrogen concentrations: A regional analysis of terrestrial controls and aquatic linkages. Ecosphere 2012, 3, art66. [Google Scholar] [CrossRef]
  45. Davidson, T.A.; Jeppesen, E. The role of palaeolimnology in assessing eutrophication and its impact on lakes. J. Paleolimnol. 2013, 49, 391–410. [Google Scholar] [CrossRef]
  46. Vadeboncoeur, Y.; Peterson, G.; Vander Zanden, M.J.; Kalff, J. Benthic algal production across lake size gradients: Interactions among morphometry, nutrients, and light. Ecology 2008, 89, 2542–2552. [Google Scholar] [CrossRef]
  47. Donato, J. Phytoplankton of Andean Lakes in Northern South America (Colombia); Witkowski, A., Ed.; A.R.G. Gantner Verlag K. G.: Ruggell, Liechtenstein, 2010; Volume 11, p. 419. [Google Scholar]
  48. Catalan, J.; Camarero, L.; Felip, M.; Pla, S.; Ventura, M.; Buchaca, T.; Bartumeus, F.; de Mendoza, G.; Miró, A.; Casamayor, E.; et al. High mountain lakes: Extreme habitats and witnesses of environmental changes. Limnetica 2006, 25, 551–584. [Google Scholar]
  49. Catalan, J.; Bartrons, M.; Camarero, L.; Grimalt, J. Mountain Waters as Witnesses of Global Pollution. In Living with Water: Targeting Quality in a Dynamic World; Pechan, P., de Vries, G.E., Eds.; Springer: New York, NY, USA, 2013; pp. 31–67. [Google Scholar] [CrossRef]
  50. Williams, J.J.; Lynch, J.A.; Saros, J.E.; Labou, S.G. Critical loads of atmospheric N deposition for phytoplankton nutrient limitation shifts in western U.S. mountain lakes. Ecosphere 2017, 8, e01955. [Google Scholar] [CrossRef] [Green Version]
  51. Bradley, R.S.; Vuille, M.; Diaz, H.F.; Vergara, W. Threats to water supplies in the tropical andes. Science 2006, 312, 1755–1756. [Google Scholar] [CrossRef] [PubMed]
  52. Striebel, M.; Schabhüttl, S.; Hodapp, D.; Hingsamer, P.; Hillebrand, H. Phytoplankton responses to temperature increases are constrained by abiotic conditions and community composition. Oecologia 2016, 182, 815–827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Burns, N.; Gibbs, M.; Hickman, M. Measurement of oxygen production and demand in lake waters. N. Z. J. Mar. Freshw. Res. 1996, 30, 127–133. [Google Scholar] [CrossRef]
  54. Fang, X.; Stefan, H.G. Simulations of climate effects on water temperature, dissolved oxygen, and ice and snow covers in lakes of the contiguous U.S. under past and future climate scenarios. Limnol. Oceanogr. 2009, 54, 2359–2370. [Google Scholar] [CrossRef]
  55. Lewis, W.M.J. Tropical lakes: How latitude makes a difference. In Perspectives in Tropical Limnology; Schiemer, F., Boland, K.T., Eds.; SPB Academic Publishing: Amsterdam, The Netherlands, 1996; pp. 43–64. [Google Scholar]
  56. Beklioglu, M.; Altinayar, G.; Tan, C. Water level control over submerged macrophyte development in five Mediterranean Turkey. Arch. Hydrobiol. 2006, 166, 535–556. [Google Scholar] [CrossRef]
  57. Zhao, D.; Jiang, H.; Cai, Y.; An, S. Artificial regulation of water level and its effect on aquatic macrophyte distribution in Taihu Lake. PLoS ONE 2012, 7, e44836. [Google Scholar] [CrossRef]
  58. Catalan, J.; Ninot, J.M.; Aniz, M.M. The High Mountain Conservation in a Changing World. In High Mountain Conservation in a Changing World; Catalan, J., Ninot, J.M., Aniz, M.M., Eds.; Springer: Cham, Switzerland, 2017; pp. 3–36. [Google Scholar] [CrossRef] [Green Version]
  59. Bennion, H.; Battarbee, R.W.; Sayer, C.D.; Simpson, G.L.; Davidson, T.A. Defining reference conditions and restoration targets for lake ecosystems using palaeolimnology: A synthesis. J. Paleolimnol. 2011, 45, 533–544. [Google Scholar] [CrossRef]
  60. Ventura, M.; Tiberti, R.; Buchaca, T.; Buñay, D.; Sabás, I.; Miró, A. Why should we preserve fishless high-mountain lakes? In High Mountain Conservation in a Changing World; Catalan, J., Ninot, J.M., Aniz, M.M., Eds.; Springer: Cham, Switzerland, 2017; pp. 181–206. [Google Scholar]
  61. Herrera-Martínez, Y.; Paggi, J.; Garcia, C. Cascading effect of exotic fish fry on plankton community in a tropical Andean high mountain lake: A mesocosm experiment. J. Limnol. 2017, 76, 397–408. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Distribution of páramo complexes in Colombia. Adapted from Instituto de Investigación de Recursos Biológicos Alexander von Humboldt (2013, 2021). The box shows the distribution of the lakes sampled for the physical and chemical study.
Figure 1. Distribution of páramo complexes in Colombia. Adapted from Instituto de Investigación de Recursos Biológicos Alexander von Humboldt (2013, 2021). The box shows the distribution of the lakes sampled for the physical and chemical study.
Water 13 02175 g001
Figure 2. Area and altitude distributions of the Colombian páramo lakes.
Figure 2. Area and altitude distributions of the Colombian páramo lakes.
Water 13 02175 g002
Figure 3. Morphometry of the 51 lakes studied in the Eastern Cordillera of the Colombian Andes.
Figure 3. Morphometry of the 51 lakes studied in the Eastern Cordillera of the Colombian Andes.
Water 13 02175 g003
Figure 4. Spatial variation of geographical, geomorphological, physical and chemical variables of the 51 lakes sampled.
Figure 4. Spatial variation of geographical, geomorphological, physical and chemical variables of the 51 lakes sampled.
Water 13 02175 g004
Figure 5. Biplot of the principal component analyses of selected environmental variables. The upper figure presents the relationship between the environmental variables. The direction of the arrows indicates the direction of the gradient and the position of samples with the highest values of each variable. The closer the arrows, the more correlated are the environmental variables. The second figure shows the position of each lake with respect to the environmental variables. The colors of the lakes correspond to the different páramo complexes studied.
Figure 5. Biplot of the principal component analyses of selected environmental variables. The upper figure presents the relationship between the environmental variables. The direction of the arrows indicates the direction of the gradient and the position of samples with the highest values of each variable. The closer the arrows, the more correlated are the environmental variables. The second figure shows the position of each lake with respect to the environmental variables. The colors of the lakes correspond to the different páramo complexes studied.
Water 13 02175 g005
Figure 6. Conceptual sketch of the potential effects of global change on páramo lakes (Arrows indicate drivers (left) enhancement of the processes (right).
Figure 6. Conceptual sketch of the potential effects of global change on páramo lakes (Arrows indicate drivers (left) enhancement of the processes (right).
Water 13 02175 g006
Figure 7. Boxplot of the second principal components (PCA 2) showed in Figure 5, grouping lakes by non-protected (n = 7) and protected areas (n = 44). Dots represent outliers.
Figure 7. Boxplot of the second principal components (PCA 2) showed in Figure 5, grouping lakes by non-protected (n = 7) and protected areas (n = 44). Dots represent outliers.
Water 13 02175 g007
Table 1. Páramo complexes of Colombia and lake distribution. C: Cordillera.
Table 1. Páramo complexes of Colombia and lake distribution. C: Cordillera.
IDPáramo ComplexSectorTotal Páramo Area (ha)Protected Fraction (%)Protection under National Park System (%)Number of LakesNumber of Lakes in Protected Areas Area Covered by Lakes (ha)
1AlmorzaderoEastern C.156,5520.4 128 301
2Altiplano CundiboyacenseEastern C.46574
3BelmiraWestern C.10,622100 110.3
4Cerro PlateadoWestern C.17,07017
5Chiles-CumbalNariño-Putumayo63,22315 3 247
6Chilí-BarragánCentral C.80,70829 7811267
7ChingazaEastern C.109,9567344124103245
8CitaráWestern C.11,23356 112
9Cruz Verde-SumapazEastern C.333,42052423242201095
10Doña Juana-ChimayoyNariño-Putumayo69,2635035281065
11El DuendeWestern C.445434
12Farallones de CaliWestern C.4545100100
13Frontino-UrraoWestern C.13,921858227
14Guanacas-Puracé-CoconucosCentral C.137,67723238139285
15Guantiva-La RusiaEastern C.119,7504539288139
16GuerreroEastern C.42,32575 715165
17Iguaque-MerchánEastern C.28,3114618151217
18Santurbán-BerlínEastern C.142,60835 20792290
19La Cocha-PatascoyNariño-Putumayo145,5393545028179
20Las HermosasCentral C.208,01164555264291437
21Los NevadosCentral C.146,0274939438395373
22Los PicachosEastern C.23,7258856112
23MirafloresEastern C.19,92885
24Nevado del Huila-MorasCentral C.150,538787714996465
25ParamilloWestern C.6744100100226
26PerijáEastern C.29,72781 1 0.1
27PisbaEastern C.106,2432727265140253
28Rabanal y río BogotáEastern C.24,65066 337
29Sierra Nevada de Santa MartaSierra Nevada151,02195954454451600
30Sierra Nevada del CocuyEastern C.271,03365657847576
31SonsónCentral C.870740
32SotaráCentral C.80,92938311512146
33TamáEastern C.16,33975532 3
34TatamáWestern C.10,930100100
35Tota-Bijagual-MamapachaEastern C.151,49847 12073208
36YariguíesEastern C.42529996
Total 2,906,13750.833.3325023038281
Table 2. Summary of the morphometrical and chemical characteristics of 51 páramo lakes in the Eastern Cordillera of the Colombian Andes.
Table 2. Summary of the morphometrical and chemical characteristics of 51 páramo lakes in the Eastern Cordillera of the Colombian Andes.
VariableMinimumMaximumMedianMeanCoefficient of Variation
Altitude (m a.s.l.)30143830361235800.05
Area (ha)0.3104.82.56.32.45
Maximum length (m)7717812453430.82
Volume (×103 m3)1.013,30338.6438.64.28
Zmax (m)0.444.36.38.20.97
Zmean (m)0.117.01.82.81.18
Secchi disk (m)0.38.43.23.50.64
Photic index0.292.351.631.390.49
Conductivity (µS/cm)2.072.05.010.71.26
Surface temperature (°C)9.418.812.112.50.15
Vertical temperature gradient (°C)04.211.210.94
O2 at surface (mg/L)58.26.526.720.10
O2 at bottom (mg/L)06.94.083.490.73
O2 saturation at surface (%)78117.199.999.00.10
O2 saturation at bottom (%)0110.656.650.50.74
pH4.78.06.16.10.16
Alkalinity (µeq/L CaCO3) 6848301001.56
SO42− (mg/L)0.00344.10.0031.3224.67
Ca2+ (mg/L)0.177.820.621.381.38
Mg2+ (mg/L)0.041.110.160.260.97
Na+ (mg/L)0.371.730.660.730.33
Mn2+ (mg/L) 0.000.090.010.020.88
K+ (mg/L)0.030.860.120.140.94
Total Fe (mg/L)0.042.200.210.321.09
Cl (mg/L)0.033.730.180.342.03
Total organic carbon (mg/L)2.614.06.26.30.40
Total Kjeldahl nitrogen (mg/L)0.701.400.981.010.17
NH4+ (mg/L)0.280.840.420.510.36
NO3 (mg/L)0.0010.1790.0070.0161.83
Total phosphorus (mg/L)0.022.050.400.550.79
PO43− (mg/L)0.0031.1250.0140.1641.77
Si (mg/L)0.133.751.201.370.62
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zapata, A.; Rivera-Rondón, C.A.; Valoyes, D.; Muñoz-López, C.L.; Mejía-Rocha, M.; Catalan, J. Páramo Lakes of Colombia: An Overview of Their Geographical Distribution and Physicochemical Characteristics. Water 2021, 13, 2175. https://doi.org/10.3390/w13162175

AMA Style

Zapata A, Rivera-Rondón CA, Valoyes D, Muñoz-López CL, Mejía-Rocha M, Catalan J. Páramo Lakes of Colombia: An Overview of Their Geographical Distribution and Physicochemical Characteristics. Water. 2021; 13(16):2175. https://doi.org/10.3390/w13162175

Chicago/Turabian Style

Zapata, Angela, Carlos A. Rivera-Rondón, Daivan Valoyes, Claudia L. Muñoz-López, Marino Mejía-Rocha, and Jordi Catalan. 2021. "Páramo Lakes of Colombia: An Overview of Their Geographical Distribution and Physicochemical Characteristics" Water 13, no. 16: 2175. https://doi.org/10.3390/w13162175

APA Style

Zapata, A., Rivera-Rondón, C. A., Valoyes, D., Muñoz-López, C. L., Mejía-Rocha, M., & Catalan, J. (2021). Páramo Lakes of Colombia: An Overview of Their Geographical Distribution and Physicochemical Characteristics. Water, 13(16), 2175. https://doi.org/10.3390/w13162175

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