Next Article in Journal
A Treatment Reliability-Based Method for Supporting Infrastructure Asset Management of Wastewater Treatment Plants
Previous Article in Journal
Development of a New Dolomite-Based Adsorbent with Phosphorus and the Adsorption Characteristics of Arsenic (III) in an Aqueous Solution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Regime and Water Balance of Two Tropical High-Mountain Lakes in the Nevado de Toluca Volcano, Mexico

by
Anatoliy Filonov
1,
María del Refugio Barba-López
2,
Iryna Tereshchenko
1,
Lydia B. Ladah
3,*,
Diego A. Pantoja
1,
Javier Alcocer
4 and
Federico Ángel Velázquez-Muñoz
1
1
Department of Physics, University of Guadalajara, Blvd. Marcelino García Barragán 1421, Guadalajara CP 44430, Jalisco, Mexico
2
Posgrado en Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, Av. Ciudad Universitaria 3000, Coyoacán, Ciudad de México CP 04510, Mexico
3
Department of Biological Oceanography, CICESE, Centro de Investigación Científica y Educación Superior de Ensenada, Carretera Ensenada-Tijuana 3918, Zona Playitas, Ensenada CP 22860, Baja California, Mexico
4
Grupo de Investigación en Limnología Tropical, FES Iztacala, Universidad Nacional Autónoma de México, Av. de los Barrios No.1, Los Reyes Iztacala, Tlanepantla CP 54090, Estado de Mexico, Mexico
*
Author to whom correspondence should be addressed.
Water 2022, 14(7), 1104; https://doi.org/10.3390/w14071104
Submission received: 15 February 2022 / Revised: 18 March 2022 / Accepted: 25 March 2022 / Published: 30 March 2022
(This article belongs to the Section Hydrology)

Abstract

:
High-mountain lakes are unique ecosystems with very few examples at tropical latitudes for experimentation. A two-year, high-frequency meteorological and water-column dataset from the crater of the Nevado de Toluca volcano, in Mexico, at an altitude of about 4200 m above sea level, allowed for the study of temporal changes in the thermal structure, water level, and water balance in the Lakes El Sol and La Luna, separated by about 500 m by a lava dome. Annual fluctuations in the water level of the lakes and calculations of the annual variability of the water balance showed that the lakes accumulated almost four times less water from rainfall than expected. Furthermore, the temperature measured at a depth of 15 cm in the bottom sediments of Lake El Sol revealed an unexpected warming during the cold season. Estimated heat fluxes through the lake bottom were less than 0.3 W m−2 during the winter and less than 0.1 W m−2 during the rest of the year. Although the variability of the hydrometeorological regime of high-mountain lakes remains relatively poorly understood, our results significantly improve the understanding of these complex processes of stratification and mixing in these unique lake ecosystems.

1. Introduction

High-mountain lakes (HML) occur above the tree line [1]. The altitude at which the tree line is observed depends on the latitude, and it occurs at higher altitudes at lower latitudes [2,3]. Mountain lakes share many features, regardless of their latitude, including catchments with steep topographic gradients, cold temperatures, and high-incident solar and ultraviolet radiation (UVR) [4,5].
Warming above the planetary average associated with climate change is higher in mountain areas, and warming rate increases with elevation. High-mountain areas experience more rapid changes in temperature than environments at lower elevations [5,6,7], and for that reason, HML are useful sentinels of global change effects, particularly those located in remote areas with limited direct human influence [8,9]. Moreover, variation in lake water level can be directly associated with the precipitation/evaporation relationship present in the lake basin [10,11,12].
Most studies on HML have focused on temperate mountain ranges with limited information from tropical systems in the literature. Tropical HML occur predominantly in the Andes, and to a lesser extent in East Africa [13,14]. The main distinction between temperate and tropical HML arises from the seasonal patterns of heat exchange and the external loadings [15], with notable water-column mixing and stratification differences. Tropical lakes are generally more sensitive than temperate lakes to environmental alterations. For example, nutrient-load increases can lead to rapid eutrophication, reduction of water quality, and loss of deep-water dissolved oxygen concentration, all of which affect aquatic biota. The differences between the limnology of temperate and tropical lakes might require different conservation strategies and are important to consider in management plans [16].
Due to the lower latitude and their high altitude, tropical HML are exposed to intense solar and UV radiation [17]. The water temperature distribution of any lake is directly controlled by solar radiation and is one of their most basic intrinsic characteristics. The upper two-meter water layer absorbs over one-half of incoming radiation, while light energy decreases exponentially with increasing water depth [18,19]. In tropical areas, the heat balance in a lake is driven principally by daily variations.
Although the impacts of local climatic conditions and complex topographical conditions on the physical behavior of lakes has been studied in general [20,21,22,23,24], there are only a few studies [25] exploring the effects of meteorological parameters on tropical HML. There are only two tropical HML in Mexico: El Sol and La Luna, both within the crater of the Nevado de Toluca volcano. The geochemistry and biology of these lakes have been studied with much less attention paid to their physics [26,27,28,29,30,31]. Recent research [25,32] shows alternating patterns of heating and cooling for ~30 days during the warmest months related with the solar irradiance reaching the lake surface. However, despite being within the same crater and very close to each other (about 500 m), these lakes show significantly different thermal responses. Lake El Sol has a discontinuous warm polymictic regime where, between each heating and cooling episode, it develops a stable mixed layer and a deeper layer separated by a weak thermal gradient. Lake La Luna also showed the same heating and cooling cycle, but it mixes and stratifies daily and not over multiple days, showing a continuous warm polymictic regime [32].
Barba-López et al. [25] used long-term time series of meteorological and water-column parameters to study the dynamics of these two HML, showing that the wind regime and solar radiation are different when measured on the outer slope of the volcano versus inside the crater, where Lake El Sol is located. At the center of the volcanic crater, the wind speed weakens, the wind direction changes, and the solar radiation is reduced by almost 20% [25]. Soto et al. [33] showed that topography affects climatic conditions by using climate series in the Nevado de Toluca volcano area. These findings highlight the relevance of studying the spatial–temporal variation of temperature, precipitation, evaporation, and wind inside the volcano to understand how these factors affect the thermal and hydrodynamic behavior of the lakes [25].
The aim of the current study was to reveal how prevailing climate variability affects the hydrophysical behavior of the lakes, including temporal changes in the thermal structure, water balance fluctuations, and surface currents, from more accurate in situ data and a numerical simulation for the last two parameters. We use a continuous two-year dataset of water-column temperature as the basis for analyzing the thermal structure dynamics from the surface to the bottom of both lakes and the water balance in Lake El Sol. A further goal was to explore the thermal interactions between the sediments and the bottom water layer.
We show that in Lake El Sol, similar to shallow lakes of middle latitudes [34], there is an exchange of heat between the bottom sediments and the water bottom layers of water. In winter, heat is transferred from bottom sediments to the water layers, and in summer, heat is accumulated in sediments. Herein, temperature fluctuations in the bottom water layers of a shallow crater lake have been measured for the first time, providing critical new information on bottom stratification and its effect on lake mixing.

2. Study Area

Lakes El Sol and La Luna are in Central Mexico (19°06″ N, 99°45″ W, 4200 m above sea level, see Figure 1). Average monthly mean temperatures ranged between 2.8 °C in February and 5.8 °C in April, with an annual mean temperature of 4.2 °C. Total annual precipitation is 1243.5 mm, ranging from 17.2 mm in December to 270 mm in July [26]. The waters of these two lakes are transparent and cold. Their highest temperatures (~11 °C) occur in the warm, rainy season, and are lowest in the cold, dry winter (~4 °C). Both lakes are well-oxygenated, with concentrations close to saturation (~7 mg/L) throughout the entire water column [10,31]. They are well-mixed lakes with a warm polymictic regime [25,32]. Both lakes have an acidic pH (4.9–5.6) and reduced buffering capacity. Their conductivity is low (18–24 µS/cm) with low dissolved and suspended organic matter [26,28].
The volcano inner slopes are shaped like a horseshoe around the two lakes in the center, where there is a circular lava dome rising almost 115 m above the lakes, separating them. The western slope of the crater is about 345 m higher than its lower eastern part, through which wind currents can enter the crater.
The hydrological parameters for Lakes El Sol and La Luna, according to Filonov et al. [35], are described in Table 1.

3. Data and Methods

3.1. Field Measurements

3.1.1. Floating Weather Station

Previous measurements [25,35] suggested an issue with estimating the wind direction and speed within the volcano crater associated with the presence of the lava dome in the center of the crater preventing their accurate estimations. In this study, we used a novel approach to avoid this issue by installing a floating weather station at the center of Lake El Sol.
The floating weather station was mounted on a metal pipe about 6 m long. The anemometer was at the upper end of the pipe, 2 m above the water level. The other end of the pipe was submerged in the water to a depth of 4 m with a 30-kg weight attached to the lower end (Figure 2). In the middle part, three buoys were attached to the pipe, keeping the station afloat. They were fixed by three ties preventing the station from horizontal oscillations and thereby avoiding errors in wind-direction measurement. Fifteen thermographs (HOBO © Water TempProV2, accuracy ±0.2 °C) were attached to the station. Two of them were above the water level (0.1 m and 2 m), and the other thirteen thermographs below the lake level, distributed every 0.5 m from the surface to 3 m depth, and then every meter to the bottom. Thermographs located above the water level were in two heat-reflecting white plastic cylinders nested inside each other, avoiding inaccuracies due to logger resin heating during the day.
Wind and air temperature measurements were carried out with a time rate of 15 min from 19 June to 25 October 2019; however, wind speed records were unavailable for the last two months of the record.

3.1.2. National Meteorological Service (NMS) Weather Station

To analyze the meteorological regime during the period of our experiment, we use long-term data (2017–2019) from the weather station 15062-Nevado de Toluca, provided by the (NMS 00015062, 19°7.171′ N; 99°44.851′ O, with altitude 4139 m a.s.l.). This weather station is located on the outer slope of the Nevado de Toluca volcano and recorded temperature, wind speed and direction, relative humidity, atmospheric pressure, solar radiation, and rainfall every 10 min. Wind data was not used due to poor performance of the speed sensor.

3.1.3. Temperature on the Internal Slope of the Nevado de Toluca Volcano

One HOBO thermistor (HOBO © Water TempProV2) was placed on the western internal slope of the volcano. The device recorded the temperature every 15 min from June to November 2018 (the measurement point is marked in Figure 1a, number 6).

3.1.4. Lake Level and Bottom Sediment Temperature Measurements

Because the heat flow in shallow lakes reaches the bottom, we assessed its role in the variability of bottom stratification [34,37]. To measure heat fluxes through the bottom of the lake (sediment water), an RBR TDR-2050 temperature and depth recorder was buried into the bottom sediment at a depth of about 15 cm. The temperature channel of the TDR-2050 was calibrated to an accuracy of ±0.002 °C (ITS-90) over a range of −5 to +35 °C. The depth channel was calibrated to an accuracy of 0.05% of full scale. We used an instrument with a depth-measurement range from 0 to 50 m, which ensured the accuracy of our measurements of 2.5 mm. These measurements were taken continuously at a sampling interval of every 2 min for 15 months (from June 2018 to August 2019).

3.1.5. Water Temperature Measurements

A mooring with a chain of HOBO thermistors (HOBO © Water TempProV2) was placed in the deepest part of each lake. The measurement period for both moorings was more than two years, from 10 May 2017, to 8 September 2019. The chain in Lake El Sol had 13 thermographs (attached to the floating weather station, see above), with 12 thermographs in Lake La Luna. All thermographs were distributed vertically from the surface to 3 m depth with seven devices every 0.5 m, and then to the bottom every meter and recorded temperature every 15 min.

3.1.6. Hydrodynamic Modeling

To estimate the circulation in Lake El Sol, we used the hydrodynamic model Delft3D, which is an open-source numerical model developed by WL/Delft Hydraulics and the Delft University of Technology [38]. The Delft3D model solves the Navier-Stokes equations for an incompressible fluid under the shallow water and the Boussinesq assumptions. The model includes the depth-averaged horizontal momentum equations
u t + u · u + g η x f v + u u C 2 d + η F x ρ d + η ν 2 u = 0
v t + u · v + g η y + f u + v u C 2 d + η F y ρ d + η ν 2 v = 0
and the depth-averaged continuity equation
η t + d + η u x + d + η v y = Q d + η
and the vertical momentum equation, which reduces to the hydrostatic pressure relationship via the Boussinesq approximation: p z = ρ g , where u , v are the depth-averaged velocity in the x and y directions, C is the Chézy coefficient, d is the water depth, η is the free surface elevation above the reference plane (at z = 0 ), u is two-dimensional current vector, with · the Euclidean norm, Q are sinks or sources of water, f is the Coriolis force, F x , y is the Reynolds stress, g is the gravity, ν is the horizontal eddy viscosity, p is the pressure, and ρ is the water density.
The RFGRID module of Delft3D was used to generate the grid file. A mesh grid of 61 × 60 cells and 1 layer were used in the horizontal and vertical direction, respectively. Bathymetry measurements obtained during the field campaigns were used to interpolate the depth file by using the QUICKIN module of Delft3D. The time step allowed was ∆t = 0.6 seg. The model was run for 34 days, with the first 4 days used as warm-up in the model. The other physical parameters had typical values: gravity, 9.81 m/s2, water density 1000 kg/m3, air density, 1 kg/m3, uniform Chézy roughness, 65 m1/2/s, background horizontal and vertical eddy viscosity and diffusivity of 1 and 10 m2/s, respectively.
The Murakami model for heat flux was applied with the default values of 75% sky cloudiness and Dalton number of 0.0013. The Secchi depth was set to 2 m.
Delft3D is a model that has been tested and used globally for the study of lake systems [39,40,41]. In addition, it has been used specifically in the hydrodynamic analysis of Mexican crater lakes [42,43].

3.2. Data Analysis

3.2.1. Meteorological Variables

We analyzed the meteorological variables at the floating weather station in Lake El Sol and the nearby weather station of the National Meteorological Service, separated by a horizontal distance of 1.7 km. To show how wind speed and direction changed during the day at the floating weather station, as well as at the NMS weather station, time series with hourly data were presented in the form of two-dimensional matrices, in which the analyzed values measured during the day were plotted along the vertical axis, and the values, by days of measurement, were located along the horizontal axis (these matrices are graphically shown in Figure 3a,b and Figure 6a,b).
Cross-spectral analysis was used to calculate coherence matrices and phase differences between the time series of all measured values. This process was carried out with the objective of recognizing the relationship between the meteorological characteristics and the temperature regime of the lakes [44].
A comparison was made between the air-temperature fluctuations on the western slope of the volcano and the temperature fluctuations in the surface water layer in Lake El Sol (near the floating weather station). The distance between the thermographs was about 200 m.

3.2.2. Water Temperature

The initial series were smoothed with a cosine filter, and the hourly time series were analyzed by using auto- and cross-spectral analysis. The spectral functions were calculated by using the Fourier transform with subsequent smoothing of the periodograms in frequency. The spectra were calculated to explore the periods and amplitudes of the dominant oscillations, as well as the coherence and phase difference between them [11,25,44,45].

3.2.3. Heat Flux between the Water-Bottom Sediments Interface

Based on water temperature measurements at the bottom and temperature in bottom sediments for the 2018–2019 time period, heat fluxes near the water-bottom interface were calculated by the gradient method [34,37] according to the formula Q = λ   T z 1 , where Q is the heat flow near the water-bottom boundary in W/m2; λ = W/m °C is the coefficient of molecular thermal conductivity of water in the 0–10 m layer 0.5813; ∂T/(∂z)-temperature gradient [46].

3.2.4. Rainfall, Evaporation and Level Fluctuations

To analyze the relationship of precipitation, evaporation from the surface of the water and the area of runoff with fluctuations in the level of Lake El-Sol, data from the NMS weather station and time series of the water level measured by the TDR-2050 instrument were used. To estimate the evaporation rate, we used the simplified version of the Penman [47] method proposed by Valiantzas [48].
For the calculations, we used a catchment area of Lake El Sol at 2.17 km2 and Lake La Luna at 2.0 km2, which allowed us to evaluate the response of Lake El Sol level response in one case of heavy precipitation.
The horizontal distance between the weather station and the mirror of Lake El Sol in a straight line is 1.7 km, and it is 47 m below the lake, that is, they are located at almost the same height [25].

4. Results and Discussion

4.1. Long-Term Series. NMS Weather Station (2017–2019)

Figure 3c shows the annual precipitation amount, and the total daily (1) and average monthly precipitation (2). The red polygon in (d) and (e) shows the period of intense rainfall, and with the effect on the level of Lake El Sol discussed in Section 4.5.1.
The influx of solar radiation on the slope of the volcano reaches a maximum value of 1000 W/m2 at about 15:00, and from 18:00 to 9:00 in the morning it is zero. Barba-López et al. [25] showed that due to the high walls, 20% less solar energy penetrates into the crater of the volcano than on its outer slope, and this must be taken into account in thermodynamic calculations. The supply of solar energy to these crater lakes is at maximum in March and April and at minimum from June to August, and low also during the rainy season due to cloud cover. Intra-annual fluctuations in air temperature in the volcano area depend mainly on the influx of solar radiation and are consistent with its variability. Fluctuations in air temperature have their maximum values in April and May in the early afternoon hours of the day (2 to 3 p.m.). The minimum temperatures were observed during the rainy season and during the winter months (Figure 3a,b).
Barba-López et al. [25] presented an analysis of an eight-year series of observations of wind and other meteorological elements at the NMS meteorological station located on the southern outer slope of the Nevado de Toluca volcano. The three-year field measurements analyzed herein, presented in Figure 3a,b, confirm previously obtained results [25].
The values of annual precipitation for the area surrounding the volcano (1162 mm, 1359 mm and 1198 mm) (Figure 3c) are within the range of normal precipitation values mentioned in Soto et al. [33] for the same study area. Precipitation was the main source of water entering the lakes. Due to their specific orography, they have a very small runoff area (Figure 3d,e). According to observations at the NMS station, an average of 1.2771 m3 of precipitation, and 0.9708 m3 evaporation, occurs per year.

4.2. Data Obtained at the Floating Weather Station

4.2.1. Water Temperature

Temperature fluctuations between 6 °C and 14 °C occurred for Lake El Sol, while on Lake La Luna values ranged between 6 °C and 13 °C (Figure 4). These values are similar to those reported by Barba-López et al. [25] and the maximum values found in Alcocer et al. [32] for these same lakes.

4.2.2. Relationship between Air Temperature Fluctuations on the Western Slope of the Volcano and Surface Temperature Fluctuations of Lake El Sol

The cross-spectral analysis of the annual time series (Figure 5) was dominated by diurnal oscillations and their subharmonics (Figure 5b) with a 0.920 agreement, with a confidence interval of 95%. The phase shift of the fluctuations was 3.5 h, with the fluctuations in water temperature leading the fluctuations in air on land (which lagged water fluctuations). That is, unexpectedly, the water warmed and cooled faster than the inner walls of the volcano.

4.2.3. Wind Speed and Direction

In Figure 6, two-dimensional plots (hours of the day/days of measurements) of the hourly variability of speed (a) and wind direction (b) measured at the floating weather station in July–October 2018 are shown. Due to the failure of the wind speed sensor, data was only available during the first two months. Nevertheless, for the first time, we obtained reliable data of the wind speed over the water surface of the lake, which was at maximum from 11:00 to 19:00 and did not exceed 3 m/s, during which time the wind had an almost stable direction of 210°, i.e., blowing over the lake from the north to its south side. The rest of the time, the wind speed was less than 1 m/s and at times, zero. The diurnal wind direction had a bimodal structure with a direction of 90–360° at night and in the morning until 10:00, also shown by the histograms of wind speed and direction measured from June–October 2018 at the station.
The values recorded by our station neared 3 m/s, somewhat lower than the average value reported by Alcocer et al. [32] (4.2 m/s).

4.2.4. Fluctuations in Air Temperature above Water and at the Near-Surface Layer of the Lake

Air-temperature fluctuations from September to May (Figure 7a) had a large amplitude (as well as a greater value by 5–7 °C) at the surface than at a height of 2 m above the water. The minimum temperatures during these months dropped to 2–3 °C at night.
From October to March, there was a balance between the temperature of the water and the air near the surface of the lake. A significant drop in air temperature at night to −5–8 °C was compensated by the influx of heat to the surface from the deep layers of the lake, which did not allow the water to freeze and form an ice cover. Winter vertical circulation was observed here, typical for lakes of middle and high latitudes [49]. In the near-surface layers of the lake during the winter months, the temperature was stable and near 4 °C, i.e., close to the temperature of the maximum density of fresh water. Seasonal fluctuations in water temperature in the near-surface layer of the lake occurred synchronously and did not create significant vertical gradients (Figure 7b).
The air temperature presented a diurnal oscillation, which is typical in tropical areas [50,51]. Besides a daily variation, a seasonal variation with greater amplitude in the months of September to May (2 m above the surface of the lake) and smaller amplitudes to 0.1 m above the surface of the lake for these same months also occurred in air temperature. For both time series of air temperature, the maximum and minimum values (−9 °C, 22 °C) are within the range of temperatures observed in other mountain areas [50,52].

4.3. Modeling Results of Currents and Level Fluctuations in Lake El Sol

In Filonov et al. [35], field measurements of currents carried out in Lake El Sol by using an ADP Sontek 1.0 MHz current meter with a vertical resolution of 0.5 m and a time rate of 1 min were presented. A monthly time series of currents made it possible to explore their structure and circulation in the lake, showing that wind currents exist only in the upper 2-m layer and do not penetrate deeper. The measured currents were weak, did not exceed 10–15 mm/s, and did not demonstrate the strict daily periodicity inherent in temperature fluctuations in the surface layers of the lake. To expand on these data, we estimated the circulation in Lake El Sol with the hydrodynamic model Delft3D.
The use of the floating weather station allowed for the reliable measurement of the direction and speed of the wind directly on the lakes without the confounding influence of the lava dome to then be used for hydrodynamic modeling (Figure 8). Numerical modeling showed that the lake contains two closed circulation patterns with opposite directions of rotation. The more intense one is in the southern, deep-water part of the lake, with a clockwise rotation with current amplitudes up to 8–10 mm/s. The second circulation pattern almost completely covers the central and northeastern parts of the lake. It is weaker, only 2–4 mm/s, and has a counterclockwise rotation. Fluctuations in current speed and direction correlate well with changes in wind speed and direction measured at the floating weather station (Figure 6) and penetrate the near-surface layer to no more than 2 m [35]. Thus, because the lakes under study are protected by high crater walls, the currents in them are very weak and mainly thermal processes drive vertical mixing.
Level fluctuations also correlated well with fluctuations in wind and currents in the lake. They have a maximum height of about ±25 mm, but their role can be significant in the mixing of lake water.
Modeling of currents and water levels in Lake La Luna were not conducted as we did not have data on the wind and fluctuations in the level of this lake.

4.4. Lake-Level Fluctuations and Processes Occurring at Its Bottom

By using the TDR-2050 RBR temperature and level meter immersed in the bottom sludge, it was possible to obtain an accurate time series of more than 13 months (Figure 9), due to the high resolution of the temperature and level sensor (accuracy of ±0.002 °C and accuracy of 0.05% full scale for pressure), which is two orders of magnitude higher compared to temperature measured with the Hobo V2. During the winter months, the temperature near the bottom of the lake drops to a minimum value of 4.7 °C, i.e., it becomes close to the temperature of the maximum density of water. The diurnal variation is well-manifested at the very bottom of the lake, where it has a larger range from November to May (0.2–0.3 °C) than the rest of the year (Figure 9a,b).

4.4.1. Heat Fluxes between the Water-Bottom Sediments Interface

In the autumn–winter period, Lake El Sol was in a state of weak thermal stratification, followed by homothermal conditions. The bottom sediments accumulated heat in summer and released it back into the water column in winter. This is the reason for the well-manifested annual temperature variation of the bottom layer of water and the upper layer of bottom sediments of the lake. In summer and autumn, the heat flow is directed from the water to the bottom sediments; in winter, on the contrary, the heat flow is directed from the bottom sediments to the water layers (Figure 9b).
Calculations showed that heat flowed through the bottom of the lake during the entire measurement period. In spring, summer and autumn, the heat flux from the bottom layers of water to the bottom sediments did not exceed 0.1 W/m2; that is, it was 10−4 times less than the flux of solar radiation penetrating the lake surface [25], where it has a maximum value of about 1000 W/m2. This is due to the dissipation and absorption of solar energy in the water column. On the contrary, in the winter months, the heat flux was directed from bottom sediments into the water column. During these months, it reached almost 0.3 W/m2.
Our results are distinct from the exchange pattern described by O’Niell et al. [53] for Lake Biwa (Japan). The difference is because the heat transfer patterns are often complex and result from a wide range of characteristics. The steady-state temperatures of volcanic lakes are largely determined by the magnitude of the volcanic heat influx relative to the surface area of the lake [54].

4.4.2. Variability of Level Fluctuations and Stratification in Lake El Sol

The annual variation of the level of Lake El Sol is quite smooth, sometimes disturbed by small sharp jumps due to precipitation. Two such time intervals are shown in Figure 9a (see rectangles). The figure shows that the annual minimum temperature in Lake El Sol at the water-bottom sediment boundary in the deep-water part of the lake reached 4.8 °C in early January 2018, where the temperature at the bottom increased slightly to 5.7 °C and dropped again in early February to 4.5 °C. In summer, the temperature at the bottom was close to 10 °C.

4.4.3. Spectra of the Level and Temperature Fluctuations

The frequency spectra of temperature and level fluctuations are shown in Figure 10. Oscillations with periods of 8 and 4 days dominate the spectra, as well as the daily period of 24 h and its overtone at 12 h. Moreover, the semidiurnal peak is higher than the diurnal one in the level spectrum due to the two maxima in atmospheric pressure near 10 and 22 h, and two minima near 4 and 16 h. The presence of these peaks in pressure variability is especially noticeable in tropical latitudes [55] and is generally associated with daily fluctuations in air temperature.
Spectral analysis has been used to study the oscillations in some alpine lakes [56], in Mexican lakes [11,12,43,57], and particularly in the Lakes El Sol and La Luna by Barba-López et al. [25]. Our results are in accordance with Barba-López et al. [25]. For another closed body of water, Arnon et al. [58] detected frequencies of 8, 12, and 24 h, as well as 2 and 3 days. The differences in the dominant frequencies of the parameters analyzed by spectral analysis may be due to the morphology of the aquatic bodies in question [59,60] and variables such as precipitation and runoff area [61].

4.5. Fluctuation of the Water Balance of the Lake Depending on Precipitation and Evaporation

The high slopes of the Nevado de Toluca crater form a drainage basin that fills with atmospheric precipitation and reaches equilibrium through the evaporation and filtration of water. The inflow of groundwater into crater lakes is not possible due to the high altitude of the volcano in comparison with the surrounding space. Thus, to calculate the water balance of a closed lake, it is necessary to know the area of its catchment, the mode of precipitation and evaporation, as well as the volume of water filtration. Of the listed values, we know only the amount of precipitation and the catchment area. This makes it possible to calculate the volume of lost water due to its filtration through the bottom.
Because the crater of the volcano has very steep southern and western slopes that shade the lakes, the penetration of solar energy in the morning and evening hours is limited with about 20% less solar radiation than during the day [25]. We took this factor into account in our calculations and obtained an annual evaporation rate of 940.1 mm/year. The evaporation calculation algorithm is detailed in Appendix A.

4.5.1. Lake El Sol Level Response in One Case of Heavy Precipitation

It is known that a thunderstorm cloud has, on average, a horizontal size of 8–15 km and rises up to 12–15 km [62,63]. That is, the height of the clouds is 2–2.5 times greater than the height of the volcano, and the dimensions of its crater are small (about 2 km) to affect the local microclimate and the patron of precipitation. Therefore, we assumed that the precipitation measurements at the meteorological station on the slope of the volcano will, on average, be the same as in its crater.
To assess the response of the lake level to intense short-term precipitation, we took such a case that occurred at the end of 2018, when the lake level almost reached its minimum annual value. During the day on 29 November, more than 50 mm of precipitation fell; a rise in the lake level by 11 cm was recorded (Figure 11). The initial data used for the calculations and their result are shown in Table 2 and Appendix A (Table A2). Note that considering evaporation from the water surface changes the calculated level rise by an insignificant amount, by only 4%.
Measurements and calculations showed that the real rise in the level of Lake El Sol as a result of heavy rains on 29 November 2018 (11 cm, Table 3) should be about five times higher than that measured by the device. When rainwater enters the lake, it is filtered through its bottom and side walls and seeps through loose volcanic rocks down the slope.
Hydrogeologists often use radioactive or fluorescent tracers and measure their presence in the water and in a drilled section across the slope of the volcano as an accurate method of confirming these calculations. Unfortunately, we were unable to make such measurements and calculations of level fluctuations in Lake La Luna. For this reason, we do not calculate the water budget in Lake La Luna.

4.6. Annual Water Budget of Lake El Sol

The data obtained as a result of our field measurements allow for rough estimates of the water balance of the lakes located in the crater of the Nevado de Toluca volcano (Table 3).
The calculation gave results close to those obtained in the previous paragraph. Over the course of the annual inflow of precipitation into the closed lake, its level in 2018 should have increased by 4.5 m, but the measured level change was only 1.2 m (Table 3).
Williams and Pelletier [61] indicate that a determining factor in these fluctuations may be the size of the lake when the analysis corresponds to long time series (tens of years), although most lake levels are studied under the assumption that they respond quickly and linearly to precipitation and runoff (dependent parameters of precipitation runoff, snowmelt, and groundwater discharge into streams and lake bottoms), in such a way that they are directly attributable to climate variability and climate change. These and many other studies show links between water level and climate variability or climate change [64,65]. However, other studies show that while water levels and climate cycles may be correlated, it is difficult to either isolate the effect of individual forcing or attribute water-level trends to that forcing [66,67].

5. Conclusions

The analysis for the years 2017–2019 of meteorological variables, such as solar radiation and air temperature outside the Nevado de Toluca volcano crater are similar to the values reported for the years 2000–2007 by Soto et al. [33] and Barba-López et al. [25].
However, the magnitude and direction of the wind directly on the surface of Lake El Sol turned out to be much smaller than on the exterior of the crater due to the presence of the walls of the volcano and their shielding of the lake. This process is well-described in works such as [20,22,68,69,70]. The values recorded herein are somewhat lower than the average value reported in the study by Alcocer et al. [32], potentially due to the difference in the measurement period and the location of the weather stations.
In our study, the air temperature presented a diurnal oscillation, which is typical in tropical areas [50,51]. The maximum and minimum values are within the range of temperatures observed in other mountain areas [50,52]. The difference in range oscillations is due to the real relationship between temperature and altitude, which varies both temporally and spatially depending on climatic conditions (moisture, wind, and radiation) and topography [21,71].
The use of spectral analysis allowed us to observe temperature and level oscillations in Lake El Sol. The morphology of the lake basin produces spectral peaks that differ from other lakes.
The numerical model allowed us to observe for the first time the wind-driven free surface level and circulation patterns. The analysis of such physical processes is of key importance as one of the first steps in the evaluation of the chemical and biological processes present in this body of water.
By using in situ measurements, unexpected changes in the temperature near the bottom of the lake were detected, with heat exchange between the bottom water layer and the sediments (in the winter the sediments give heat to the water column).
Finally, our study found that the level of the lake increased less than the theoretical calculations would suggest. The analysis of the level of the lake corresponds to a time series of two years and the difference between the measured levels and the theoretical calculations is observed for the first time at this site.
Our results highlight the difficulty of planning and maintaining equipment in high-mountain areas, in addition to the complicated relationships between meteorological drivers and the thermal regime and water balances in these unique mountain lake ecosystems.

Author Contributions

A.F.: conceptualization, investigation, writing original draft, methodology, software, validation, visualization, project administration. M.d.R.B.-L.: writing original draft, formal analysis, investigation, data processing and modeling. L.B.L.: investigation, writing—review & editing. I.T.: data analysis, investigation, review & editing. D.A.P.: field work, data processing and modeling. J.A.: financial support, review & editing. F.Á.V.-M.: investigation, data processing and modeling. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Universidad Nacional Autónoma de México DGAPA/PAPIIT through projects IN105009 and ES209301, and Programa de Investigación en Cambio Climático (PINCC 2020–2021).

Data Availability Statement

Data are available from the authors upon reasonable request.

Acknowledgments

The authors express their gratitude to each of the expedition participants who made it possible to complete this work. Omar Mireles-Loera provided assistance in the field work. We are also very grateful to the graduate students at the University of Guadalajara, Arturo Orozco-Estrada and Andrea Manrique-Cantillo who, in cold weather, rain or snow, performed the difficult work of measuring at the top of the volcano. Special thanks to graduate student Carlos Villarreal-Olavarrieta for his work as an equipment diver. We also thank the Comisión Estatal de Parques Naturales y de la Fauna (CEPANAF, Secretaría de Ecología, Gobierno del Estado de México) for providing the permit for scientific research at the Área de Protección de Flora y Fauna Nevado de Toluca.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A. Calculation of Evaporation from the Surface of the Lake

Various authors have developed methods and equations for calculating evaporation through adjustments, correlations of climatic variables, as well as through theories of energy balance, water balance, aerodynamic methods and others [72]. One of the most used methods is the one proposed by Penman [47], where evaporation is calculated from the effects of wind speed and radiation solar. To do this, he suggested using the following formula:
P E = Δ Δ + γ + R n λ + γ Δ + γ 6.43 f u D λ
where PE is the potential evaporation (mm/ day); Rn is the net radiation on the surface of the water body (MJ/m2day), ∆ is the water vapor saturation curve (kPa/°C); γ corresponds to the psychrometric constant (kPa/°C); λ represents the latent heat of vaporization (MJ/kg) f u is the function of the wind, and D is the vapor pressure deficit. Valiantzas [48] simplified Equation (A1) by calculating evaporation from more readily available measured data. This new approximation is equivalent in accuracy to Equation (A1) and allows calculating the daily evaporation rate for each month of the year.
The simplified Penman formula proposed by Valiantzas [48] is represented:
P E S 0.047 R S T + 9.5 2.4 R s R A 2 + 0.09 T + 20 1 R H 100
where R s is the solar radiation (MJ/m2 day, T is the monthly average temperature (°C) at 2m above the surface, RH is the relative humidity and (MJ/m2 day) is the extraterrestrial radiation. In this case, the extraterrestrial radiation values can be traced using the formula:
R A = 24 60 π G S C d r [ ω s sin Φ sin δ + cos Φ cos δ sin ω s ]
where RA is the extraterrestrial radiation (MJ/m2 day), G S C is the solar constant G S C = 0.0820 MJ/m2 min, d r is the mean distance from the Earth to the Sun, ω s corresponds to the sunset angle of sun (radians), Φ is latitude (radians) and δ solar declination (radians). To calculate the evaporation in the crater Lakes El Sol and La Luna, the values given in Table A1 were used.
Table A1. Parameters for estimating the evaporation rate in Lakes El Sol and La Luna.
Table A1. Parameters for estimating the evaporation rate in Lakes El Sol and La Luna.
YearMonthT2 m Average (°C)RS Average (MJ/m2 day)RA (MJ/m2 day)HR Average (%)
2018July5.214.426.885.8
August4.716.730.698.3
September5.414.834.785.3
October5.011.837.979.4
November3.014.939.361.9
December2.617.439.565.2
2019January2.117.939.362.3
February3.219.638.352.6
March4.323.335.864.4
April4.622.431.862.5
May6.422.627.779.1
June5.813.925.673.4
Where average T2 m corresponds to the average monthly temperature at 2 m above the surface of the lake, RS average is the monthly solar radiation calculated from data recorded by the automatic meteorological station Nevado de Toluca (NMS) with the monthly extraterrestrial solar radiation values calculated from Equation (A3), and average RH corresponds to the monthly average values of relative humidity recorded by a meteorological station located inside the crater of the volcano. Substituting the values from Table 2 in Equation (A2); we obtain the evaporation values in Table A2.
Table A2. Monthly and annual evaporation in Lakes El Sol and La Luna.
Table A2. Monthly and annual evaporation in Lakes El Sol and La Luna.
Monthly
Evaporation
(mm/month)
July 201868.8
August 201870.7
September 201877.4
October 201872.5
November 201887.9
December 201895.5
January 201996.7
February 2019112.2
March 2019107.2
April 2019107.7
May 201997.0
June 201976.3
Annual
evaporation
1070.0
(mm/year)
CRediT authorship contribution statement.

References

  1. Löffler, H. Contribution to the Limnology of High Montain Lakes in Central America. Int. Rev. Hydrobiol. Hydrogr. 1972, 57, 397–408. [Google Scholar] [CrossRef]
  2. Hutchinson, G.E. A Treatise on Limnology. In Geography, Physics and Chemistry; John Wiley and Sons: New York, NY, USA, 1957; Volume 1, p. 1015. [Google Scholar]
  3. Sommaruga, R. The role of solar UV radiation in the ecology of alpine lakes. J. Photochem. Photobiol. B Biol. 2001, 62, 35–42. [Google Scholar] [CrossRef]
  4. Hobbs, W.O.; Telford, R.J.; Birks, H.J.B.; Saros, J.E.; Hazewinkel, R.R.O.; Perren, B.B.; Saulnier-Talbot, É.; Wolfe, A.P. Quantifying recent ecological changes in remote lakes of North America and greenland using sediment diatom assemblages. PLoS ONE 2010, 5, e10026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Moser, K.A.; Baron, J.S.; Brahney, J.; Oleksy, I.A.; Saros, J.E.; Hundey, E.J.; Sadro, S.A.; Kopáček, J.; Sommaruga, R.; Kainz, M.J.; et al. Mountain lakes: Eyes on global environmental change. Glob. Planet. Change 2019, 178, 77–95. [Google Scholar] [CrossRef] [Green Version]
  6. Pepin, N.; Bradley, R.S.; Diaz, H.F.; Baraer, M.; Caceres, E.B.; Forsythe, N.; Fowler, H.; Greenwood, G.; Hashmi, M.Z.; Liu, X.D.; et al. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Change 2015, 5, 424–430. [Google Scholar]
  7. Schirpke, U.; Ebner, M.; Pritsch, H.; Fontana, V.; Kurmayer, R. Quantifying ecosystem services of high mountain lakes across different socio-ecological contexts. Sustainability 2021, 13, 6051. [Google Scholar] [CrossRef]
  8. Battarbee, R.W.; Grytnes, J.A.; Thompson, R.; Appleby, P.G.; Catalan, J.; Korhola, A.; Lami, A. Comparing palaeolimnological and instrumental evidence of climate change for remote mountain lakes over the last 200 years. J. Paleolimnol. 2002, 28, 161–179. [Google Scholar] [CrossRef]
  9. Catalan, J.; Pla-Rabés, S.; Wolfe, A.P.; Smol, J.P.; Rühland, K.M.; Anderson, N.J.; Kopáček, J.; Stuchlík, E.; Schmidt, R.; Koinig, K.A.; et al. Global change revealed by palaeolimnological records from remote lakes: A review. J. Paleolimnol. 2013, 49, 513–535. [Google Scholar] [CrossRef]
  10. Alcocer, J. Las aguas celestiales. Nevado de Toluca. In Las Aguas Celestiales. Nevado de Toluca; Luna, P., Montero, A., Junco, R., Eds.; Instituto Nacional de Antropología e Historia: Ciudad de México, Mexico, 2009; Volume 16, pp. 10–13. [Google Scholar]
  11. Filonov, A. On the dynamical response of Lake Chapala, Mexico to lake breeze forcing. Hydrobiologia 2002, 141–157. [Google Scholar] [CrossRef]
  12. Filonov, A.; Tereshchenko, I.; Alcocer, J. Dynamic response to mountain breeze circulation in Alchichica, a crater lake in Mexico. Geophys. Res. Lett. 2006, 33, 2–5. [Google Scholar] [CrossRef]
  13. Löffler, H. The limnology of tropical high-mountain lakes Internationale. Vereinigung Theor. Angew. Limnol. Verh. 1964, 15, 176–193. [Google Scholar]
  14. Payne, A.I. The ecology of tropical lakes and rivers. Ecol. Trop. 1991, 75, 582. [Google Scholar] [CrossRef]
  15. Catalan, J.; Donato Rondón, J.C. Perspectives for an integrated understanding of tropical and temperate high-mountain lakes. J. Limnol. 2016, 75, 215–234. [Google Scholar] [CrossRef] [Green Version]
  16. Alcocer, J.; Oseguera, L.A.; Ibarra-Morales, D.; Escobar, E.; García-Cid, L. Responses of benthic macroinvertebrate communities of two tropical, high-mountain lakes to climate change and deacidification. Diversity 2021, 13, 243. [Google Scholar] [CrossRef]
  17. Kinzie, R.A.; Banaszak, A.T.; Lesser, M.P. Effects of ultraviolet radiation on primary productivity in a high altitude tropical lake. Hydrobiologia 1998, 385, 23–32. [Google Scholar] [CrossRef]
  18. Wetzel, R. Limnology. Lake and River Ecosystems; Academic Press: Cambridge, MA, USA, 2001. [Google Scholar]
  19. Wang, J.; Huang, L.; Ju, J.; Daut, G.; Wang, Y.; Ma, Q.; Zhu, L.; Haberzettl, T.; Baade, J.; Mäusbacher, R. Spatial and temporal variations in water temperature in a high-altitude deep dimictic mountain lake (Nam Co), central Tibetan Plateau. J. Great Lakes Res. 2019, 45, 212–223. [Google Scholar] [CrossRef]
  20. Baron, J.S.; Caine, N. Temporal coherence of two alpine lake basins of the Colorado Front Range, U.S.A. Freshw. Biol. 2000, 43, 463–476. [Google Scholar] [CrossRef]
  21. Becker, A.; Bugmann, H. Predicting Global Change Impacts on Mountain Hydrology and Ecology: Integrated Catchment Hydrology/Altitudinal Gradient Studies: A Workshop report. IGBP Rep. 1997, 43, 7–9. [Google Scholar]
  22. Bitan, A. Lake kinneret (sea of galilee) and its exceptional wind system. Bound.-Layer Meteorol. 1981, 21, 477–487. [Google Scholar] [CrossRef]
  23. Bueche, T.; Vetter, M. Simulating water temperatures and stratification of a pre-alpine lake with a hydrodynamic model: Calibration and sensitivity analysis of climatic input parameters. Hydrol. Process. 2014, 28, 1450–1464. [Google Scholar] [CrossRef]
  24. Solarski, M.; Szumny, M. Conditions of spatiotemporal variability of the thickness of the ice cover on lakes in the Tatra Mountains. J. Mt. Sci. 2020, 17, 2369–2386. [Google Scholar] [CrossRef]
  25. Barba-López, M.D.R.; Filonov, A.; Alcocer, J.; Ladah, L.; Tereshchenko, I. Complex interactions between meteorological parameters and the thermal regime of two tropical, high altitude crater lakes (Nevado de Toluca, Mexico). Sci. Total Environ. 2021, 771, 145310. [Google Scholar] [CrossRef] [PubMed]
  26. Alcocer, J.; Oseguera, L.A.; Escobar, E.; Peralta, L.; Lugo, A. Phytoplankton biomass and water chemistry in two high-mountain Tropical Lakes in central Mexico. Arctic. Antarct. Alp. Res. 2004, 36, 342–346. [Google Scholar] [CrossRef] [Green Version]
  27. Armienta, M.A.; De La Cruz-Reyna, S.; Macías, J.L. Chemical characteristics of the crater lakes of Popocatetetl, El Chichon, and Nevado de Toluca volcanoes, Mexico. J. Volcanol. Geotherm. Res. 2000, 97, 105–125. [Google Scholar] [CrossRef]
  28. Caballero-Miranda, M. The diatom flora of two acid lakes in central Mexico. Diatom Res. 1996, 11, 227–240. [Google Scholar] [CrossRef]
  29. Cuna, E.; Caballero, M.; Zawisza, E.; Ruiz, C. Historia Ambiental De Un Lago Alpino En El Centro De México (1230–2010). Tip 2015, 18, 97–106. [Google Scholar] [CrossRef] [Green Version]
  30. García-Palomo, A.; Macías, J.L.; Arce, J.L.; Capra, L.; Garduño, V.H.; Espíndola, J.M. Geology of the Nevado de Toluca Volcano and Surrounding Areas, Central Mexico. 2002. Available online: https://www.researchgate.net/publication/265539610_Geology_of_Nevado_de_Toluca_Volcano_and_surrounding_areas_central_Mexico (accessed on 14 February 2022).
  31. Oseguera, L.A.; Alcocer, J.; Escobar, E. Macroinvertebrados bentónicos de dos lagos tropicales de alta montaña en el volcán Nevado de Toluca, en la región central de México. Hidrobiologica 2016, 26, 419–432. [Google Scholar]
  32. Alcocer, J.; Roberson, J.; Oseguera, L.A.; Lewis, W.M., Jr. Rhythmic episodes of heating and cooling control thermal stratification of two tropical high mountain lakes. Aquat. Sci. 2020, 82, 58. [Google Scholar] [CrossRef]
  33. Soto, V.H.; Yoshikawa, K.; Schörghofer, N. Climatic variation in the high mountains of central Mexico: Temperature and precipitation indices at Nevado de Toluca volcano. Atmósfera 2020, 33, 301–318. [Google Scholar] [CrossRef]
  34. Gavrilenko, G.G.; Zdorovennova, G.E.; Zdorovennov, R.E.; Palshin, N.I.; Mitrokhov, A.V.; Terzhevik, A.Y. Heat Flux At the Water-Sediment Interface in a Shallow Lake. Proc. Karelian Res. Cent. Russ. Acad. Sci. 2015, 3, 3. [Google Scholar] [CrossRef] [Green Version]
  35. Filonov, A.; Barba-López, M.D.R.; Ladah, L.; Tereshchenko, I.; Palacios-Hernández, E.; Alcocer, J.; Ávalos-Cueva, D. The dynamics and temperature regime of the crater lakes in the Nevado de Toluca volcano, Mexico. Limnologica 2022, 92, 125938. [Google Scholar] [CrossRef]
  36. Guía de México: Información, Turismo y Viajes a México. Available online: https://www.turismomexico.es/ (accessed on 14 February 2022).
  37. Golosov, S.; Kirillin, G. A parameterized model of heat storage by lake sediments. Environ. Model. Softw. 2010, 25, 793–801. [Google Scholar] [CrossRef]
  38. Deltares. Simulation of Multi-Dimensional Hydrodynamic Flows and Transport Phenomena, Including Sediments; User Manual Delft3D-FLOW; Version 3.15; Deltares: Delft, The Netherlands, 2020. [Google Scholar]
  39. Baracchini, T.; Hummel, S.; Verlaan, M.; Cimatoribus, A.; Wüest, A.; Bouffard, D. An automated calibration framework and open source tools for 3D lake hydrodynamic models. Environ. Model. Softw. 2020, 134, 104787. [Google Scholar] [CrossRef]
  40. Kranenburg, W.; Tiessen, M.; Veenstra, J.; de Graaff, R.; Uittenbogaard, R.; Bouffard, D.; Sakindi, G.; Umutoni, A.; Van de Walle, J.; Thiery, W.; et al. 3D-modelling of Lake Kivu: Horizontal and vertical flow and temperature structure under spatially variable atmospheric forcing. J. Great Lakes Res. 2020, 46, 947–960. [Google Scholar] [CrossRef]
  41. Razmi, A.M.; Barry, D.A.; Bakhtyar, R.; Le Dantec, N.; Dastgheib, A.; Lemmin, U.; Wüest, A. Current variability in a wide and open lacustrine embayment in Lake Geneva (Switzerland). J. Great Lakes Res. 2013, 39, 455–465. [Google Scholar] [CrossRef]
  42. Cardoso-Mohedano, J.G.; Sanchez-Cabeza, J.A.; Ruiz-Fernández, A.C.; Pérez-Bernal, L.H.; Lima-Rego, J.; Giralt, S. Fast deep water warming of a subtropical crater lake. Sci. Total Environ. 2019, 691, 1353–1361. [Google Scholar] [CrossRef]
  43. Gasca-Ortiz, T.; Pantoja, D.A.; Filonov, A.; Domínguez-Mota, F.; Alcocer, J. Numerical and observational analysis of the hydro-dynamical variability in a small lake: The case of Lake Zirahuén, México. Water 2020, 12, 1658. [Google Scholar] [CrossRef]
  44. Emery, W.J.; Thomson, R.E. Data Analysis Methods in Physical Oceanography; Elsevier: Amsterdam, The Netherlands, 1997. [Google Scholar]
  45. Jenkins, D.M.; Watts, D.G. Spectral Analysis and Its Applications; Card Nr. 67-13840; Holden-Day, Inc.: San Francisco, CA, USA, 1969. [Google Scholar]
  46. Thermal Conductivity—Online Converter. Available online: https://www.engineeringtoolbox.com/thermal-conductivity-calculator-d_857.html (accessed on 14 February 2022).
  47. Penman, H.L. Natural evaporation from open water, bare soil and grass. Proc. R. Soc. Lond. Math. Phys. Sci. 1948, 193, 120–145. [Google Scholar]
  48. Valiantzas, J.D. Simplified versions for the Penman evaporation equation using routine weather data. J. Hydrol. 2006, 331, 690–702. [Google Scholar] [CrossRef]
  49. Rodgers, G.K. The Thermal Bar in Lake Ontario, Spring 1965 and Winter 1965–66. Gt. Lakes Res. Div. Univ. Michigan Ann. Arbor. 1966, 15, 369–374. [Google Scholar]
  50. Labaj, A.L.; Michelutti, N.; Smol, J.P. Annual stratification patterns in tropical mountain lakes reflect altered thermal regimes in response to climate change. Fundam. Appl. Limnol. 2018, 191, 267–275. [Google Scholar] [CrossRef]
  51. Michelutti, N.; Labaj, A.L.; Grooms, C.; Smol, J.P. Equatorial mountain lakes show extended periods of thermal stratification with recent climate change. J. Limnol. 2016, 75, 403–408. [Google Scholar] [CrossRef]
  52. Livingstone, D.M.; Lotter, A.F.; Walker, I.R. The decrease in summer surface water temperature with altitude in Swiss Alpine lakes: A comparison with air temperature lapse rates. Arctic. Antarct. Alp. Res. 1999, 31, 341–352. [Google Scholar] [CrossRef]
  53. O’Neill, K.; Ashton, G.D. Bottom Heat Transfer to Water Bodies in Winter. 1981. Available online: https://erdc-library.erdc.dren.mil/jspui/handle/11681/12042 (accessed on 14 February 2022).
  54. Pasternack, G.B.; Varekamp, J.C. Volcanic lake systematics I. Physical constraints. Bull. Volcanol. 1997, 58, 528–538. [Google Scholar] [CrossRef]
  55. Berberan-Santos, M.N.; Bodunov, E.N.; Pogliani, L. On the barometric formula inside the Earth. Am. J. Phys. 1997, 65, 404–412. [Google Scholar] [CrossRef]
  56. van Haren, H.; Dijkstra, H.A. Convection under internal waves in an alpine lake. Environ. Fluid Mech. 2021, 21, 305–316. [Google Scholar] [CrossRef]
  57. Serrano, D. Dynamic response to valley breeze circulation in Santa María del Oro, a volcanic lake in Mexico. Geophys. Res. Lett. 2002, 29, 13–16. [Google Scholar] [CrossRef]
  58. Arnon, A.; Brenner, S.; Selker, J.S.; Gertman, I.; Lensky, N.G. Seasonal dynamics of internal waves governed by stratification stability and wind: Analysis of high-resolution observations from the Dead Sea. Limnol. Oceanogr. 2019, 64, 1864–1882. [Google Scholar] [CrossRef]
  59. Gloor, M.; Wüest, A.; Münnich, M. Benthic boundary mixing and resuspension induced by internal seiches. Hydrobiologia 1994, 284, 59–68. [Google Scholar] [CrossRef]
  60. Rueda, F.J.; Schladow, S.G. Surface seiches in lakes of complex geometry. Limnol. Oceanogr. 2002, 47, 906–910. [Google Scholar] [CrossRef]
  61. Williams, Z.C.; Pelletier, J.D. Self-affinity and surface-area-dependent fluctuations of lake-level time series. Water Resour. Res. 2015, 51, 7258–7269. [Google Scholar] [CrossRef]
  62. Oliver, J.E. Encyclopedia of World Climatology; Springer: New York, NY, USA, 2005; ISBN 978-1-4020-3264-6. [Google Scholar]
  63. Yau, M.K.; Rogers, R.R. A Short Course in Cloud Physics, 3rd ed.; Butterworth-Heinemann: Oxford, UK, 1989; ISBN 9780750632157. [Google Scholar]
  64. Mercier, F.; Cazenave, A.; Maheu, C. Interannual lake level fluctuations (1993–1999) in Africa from Topex/Poseidon: Connections with ocean-atmosphere interactions over the Indian Ocean. Glob. Planet. Change 2002, 32, 141–163. [Google Scholar] [CrossRef]
  65. Wilcox, D.A.; Thompson, T.A.; Booth, R.K.; Nicholas, J.R. Lake-level variability and water availability in the Great Lakes. US Geol. Surv. Circ. 2007, 1311, 25. [Google Scholar] [CrossRef]
  66. Pasquini, A.I.; Lecomte, K.L.; Depetris, P.J. Climate change and recent water level variability in Patagonian proglacial lakes, Argentina. Glob. Planet. Change 2008, 63. [Google Scholar] [CrossRef]
  67. Sellinger, C.E.; Stow, C.A.; Lamon, E.C.; Qian, S.S. Recent water level declines in the Lake Michigan-Huron system. Environ. Sci. Technol. 2008, 42, 367–373. [Google Scholar] [CrossRef] [PubMed]
  68. Kossmann, M.; Sturman, A.P.; Zawar-Reza, P.; McGowan, H.A.; Oliphant, A.J.; Owens, I.F.; Spronken-Smith, R.A. Analysis of the wind field and heat budget in an alpine lake basin during summertime fair weather conditions. Meteorol. Atmos. Phys. 2002, 81, 27–52. [Google Scholar] [CrossRef]
  69. McGowan, H.A.; Owens, I.F.; Sturman, A.P. Thermal and dynamic characteristics of alpine lake breezes, Lake Tekapo, New Zealand. Bound.-Layer Meteorol. 1995, 76, 3–24. [Google Scholar] [CrossRef]
  70. Santo, M.A.; Toffolon, M.; Zanier, G.; Giovannini, L.; Armenio, V. Large eddy simulation (LES) of wind-driven circulation in a peri-alpine lake: Detection of turbulent structures and implications of a complex surrounding orography. J. Geophys. Res. Ocean. 2017, 122, 4704–4722. [Google Scholar] [CrossRef] [Green Version]
  71. Barry, R.G. Mountain Weather and Climate, 3rd ed.; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  72. Barco, J.; Cuartas, A.; Mesa, O.; Poveda, G.; Vélez, J.I.; Hoyos, C.; Mejía, J.F.; Botero, B. Estimación de la evaporación en colombia. Av. EN Recur. Hidraul. 2000, 7, 43–51. [Google Scholar]
Figure 1. (a) Geographical location of the Nevado de Toluca volcano and Lakes El Sol and La Luna. Photograph of the volcano crater is available in [36]. The numbers indicate: 1—Lake El Sol, 2—Lake La Luna, 3—Lava dome, 4—the position of the floating weather station, 5—the position of the meteorological station installed on the shore in 2007–2010, 6—location of the Hobo V-2 thermograph on the western slope of the crater. (2). (b) A topographical chart of the volcano. The red dots indicate the location of the floating weather station in Lake El Sol (1), and the location of the SMN weather station on the outer slope. Volcano diagram adapted from [29].
Figure 1. (a) Geographical location of the Nevado de Toluca volcano and Lakes El Sol and La Luna. Photograph of the volcano crater is available in [36]. The numbers indicate: 1—Lake El Sol, 2—Lake La Luna, 3—Lava dome, 4—the position of the floating weather station, 5—the position of the meteorological station installed on the shore in 2007–2010, 6—location of the Hobo V-2 thermograph on the western slope of the crater. (2). (b) A topographical chart of the volcano. The red dots indicate the location of the floating weather station in Lake El Sol (1), and the location of the SMN weather station on the outer slope. Volcano diagram adapted from [29].
Water 14 01104 g001
Figure 2. Structural diagram of a floating weather station. In the right corner is a photograph of the part of the station above the water level.
Figure 2. Structural diagram of a floating weather station. In the right corner is a photograph of the part of the station above the water level.
Water 14 01104 g002
Figure 3. Hourly variability of solar radiation (a), air temperature (b) and precipitation (c) measured during 2017–2019 at the NMS station. In (d,e), the red curves show the runoff areas in the crater lakes of the Nevado de Toluca volcano.
Figure 3. Hourly variability of solar radiation (a), air temperature (b) and precipitation (c) measured during 2017–2019 at the NMS station. In (d,e), the red curves show the runoff areas in the crater lakes of the Nevado de Toluca volcano.
Water 14 01104 g003
Figure 4. Hourly variability of temperature (a) at thirteen depths of Lake El Sol (a) and at twelve depths of Lake La Luna (b) for the same time period (instrumental depths are indicated in the legend).
Figure 4. Hourly variability of temperature (a) at thirteen depths of Lake El Sol (a) and at twelve depths of Lake La Luna (b) for the same time period (instrumental depths are indicated in the legend).
Water 14 01104 g004
Figure 5. (a) Variability of water temperature in the surface layer of Lake El Sol in the area of the floating meteorological station (1) and air temperature on the western slope of the volcano crater (2), from June-November 2018. (b) The corresponding spectra of these time series. The vertical line in Figure 5b shows the 95% confidence interval.
Figure 5. (a) Variability of water temperature in the surface layer of Lake El Sol in the area of the floating meteorological station (1) and air temperature on the western slope of the volcano crater (2), from June-November 2018. (b) The corresponding spectra of these time series. The vertical line in Figure 5b shows the 95% confidence interval.
Water 14 01104 g005
Figure 6. Hourly variability of speed (a) and wind direction (b) measured at the floating weather station located in the center of Lake El Sol. (c,d) Representations of the histograms of (a,b). The red vertical line shows the day for which the numerical modeling of fluctuations of currents and the level of Lake El Sol was carried out.
Figure 6. Hourly variability of speed (a) and wind direction (b) measured at the floating weather station located in the center of Lake El Sol. (c,d) Representations of the histograms of (a,b). The red vertical line shows the day for which the numerical modeling of fluctuations of currents and the level of Lake El Sol was carried out.
Water 14 01104 g006
Figure 7. Hourly variability of air and water temperature in the area of the floating weather station. (a) At a height of 0.1 m above the level of the lake (1) and at a height of 2 m (2). (b) Temperature fluctuations directly below the water surface (3) and at a depth of 2 m (4).
Figure 7. Hourly variability of air and water temperature in the area of the floating weather station. (a) At a height of 0.1 m above the level of the lake (1) and at a height of 2 m (2). (b) Temperature fluctuations directly below the water surface (3) and at a depth of 2 m (4).
Water 14 01104 g007
Figure 8. Modeling results of currents and level fluctuation in Lake El Sol from observations of wind and air temperature on 5 September 2018. In the upper left, a bathymetric map of the lake shows the floating meteorological station marked with a white circle (in the rest of the figures, the position of the weather station is marked with a black dot). The remaining rectangles show simulated currents and level fluctuations every two hours. In the lower left, black and red pattern arrows indicate the corresponding scales of current velocities and measured wind.
Figure 8. Modeling results of currents and level fluctuation in Lake El Sol from observations of wind and air temperature on 5 September 2018. In the upper left, a bathymetric map of the lake shows the floating meteorological station marked with a white circle (in the rest of the figures, the position of the weather station is marked with a black dot). The remaining rectangles show simulated currents and level fluctuations every two hours. In the lower left, black and red pattern arrows indicate the corresponding scales of current velocities and measured wind.
Water 14 01104 g008
Figure 9. (a) Annual series of temperature fluctuations (1) of bottom sediments and (2) the level of Lake El Sol in the area of the floating meteorological station. Moments of a sharp change in the measured characteristics are marked by black and red rectangles. (b) Vertical gradient of water temperature between the bottom and 1 m above the bottom (9 and 8 m depth) (1). The red curve (2) is the calculated heat flux across the water-bottom sediment interface (positive values are heat flux from bottom sediments to water layers).
Figure 9. (a) Annual series of temperature fluctuations (1) of bottom sediments and (2) the level of Lake El Sol in the area of the floating meteorological station. Moments of a sharp change in the measured characteristics are marked by black and red rectangles. (b) Vertical gradient of water temperature between the bottom and 1 m above the bottom (9 and 8 m depth) (1). The red curve (2) is the calculated heat flux across the water-bottom sediment interface (positive values are heat flux from bottom sediments to water layers).
Water 14 01104 g009
Figure 10. Spectra of lake level fluctuations (1) and temperature fluctuations (2) in bottom sediments in 2018–2019 in the area of the floating weather station (see Figure 2). The vertical line in Figure 10 shows the 95% confidence interval.
Figure 10. Spectra of lake level fluctuations (1) and temperature fluctuations (2) in bottom sediments in 2018–2019 in the area of the floating weather station (see Figure 2). The vertical line in Figure 10 shows the 95% confidence interval.
Water 14 01104 g010
Figure 11. Daily precipitation (NMS data) and level values (TDR 2050 data) measured on 27–30 November 2018 at the time of heavy precipitation. This segment of observations corresponds to the red rectangle inn Figure 9a.
Figure 11. Daily precipitation (NMS data) and level values (TDR 2050 data) measured on 27–30 November 2018 at the time of heavy precipitation. This segment of observations corresponds to the red rectangle inn Figure 9a.
Water 14 01104 g011
Table 1. Hydrological parameters for Lake El Sol and Lake La Luna. Modified from Filonov et al. [35].
Table 1. Hydrological parameters for Lake El Sol and Lake La Luna. Modified from Filonov et al. [35].
At the End of the Rainy Season (October)At the End of the the Dry Season (May)
El Sol LakeLa Luna LakeEl Sol LakeLa Luna Lake
Surface area2.020 × 105 m20.203 × 105 m22.003 × 105 m20.202 × 105 m2
Maximum depth13.8 m10.2 m10.9 m7.6 m
Average depth6.4 m5.7 m4.0 m3.5 m
Volume13.022 × 105 m31.172 × 105 m38.102 × 105 m30.725 × 105 m3
Table 2. Parameters used to calculate the response of Lake El Sol to heavy rainfall on 29 November 2018.
Table 2. Parameters used to calculate the response of Lake El Sol to heavy rainfall on 29 November 2018.
Dry season area of the lake El Sol [35]200,330 m2
Runoff area2,170,000 m2
Evaporation rate (November 2018)3.33 mm/day
Rain case50.8 mm
Real level rise0.11 m
Accumulated in runoff area110,236 m3
Expected level rise with evaporation0.550 m
Expected level rise without evaporation0.586 m
Table 3. Annual average values (with CNA measures).
Table 3. Annual average values (with CNA measures).
Lake El Sol average area for the rainy and dry seasons [35]201,165 m2
Runoff area2,170,000 m2
Rain rate1.2771 m
(NMS measures)
1.359 m
(calculations)
Evaporation rate0.9708 m
(NMS measures)
0.941 m
(calculations)
Accumulated volume907,060 m3
Expected level rise with evaporation4.509 m
Real level rise1.20 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Filonov, A.; Barba-López, M.d.R.; Tereshchenko, I.; Ladah, L.B.; Pantoja, D.A.; Alcocer, J.; Velázquez-Muñoz, F.Á. Thermal Regime and Water Balance of Two Tropical High-Mountain Lakes in the Nevado de Toluca Volcano, Mexico. Water 2022, 14, 1104. https://doi.org/10.3390/w14071104

AMA Style

Filonov A, Barba-López MdR, Tereshchenko I, Ladah LB, Pantoja DA, Alcocer J, Velázquez-Muñoz FÁ. Thermal Regime and Water Balance of Two Tropical High-Mountain Lakes in the Nevado de Toluca Volcano, Mexico. Water. 2022; 14(7):1104. https://doi.org/10.3390/w14071104

Chicago/Turabian Style

Filonov, Anatoliy, María del Refugio Barba-López, Iryna Tereshchenko, Lydia B. Ladah, Diego A. Pantoja, Javier Alcocer, and Federico Ángel Velázquez-Muñoz. 2022. "Thermal Regime and Water Balance of Two Tropical High-Mountain Lakes in the Nevado de Toluca Volcano, Mexico" Water 14, no. 7: 1104. https://doi.org/10.3390/w14071104

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