Boletín de la Sociedad Geológica Mexicana Volumen 73, núm. 1, A171220, 2021 http://dx.doi.org/10.18268/BSGM2021v73n1a171220

Environmental variability during the last three millennia in the rain shadows of central Mexico

Variabilidad ambiental durante los últimos tres milenios en las sombras de lluvia de México central

Gustavo Olivares-Casillas¹, Alex Correa-Metrio^2,*, Edyta Zawisza³, Marta Wojewódka-Przybył³, Maarten Blaauw⁴, Francisco Romero²

¹ Posgrado en Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, Ciudad Universitaria, 04510, Coyoacán,CDMX, Mexico.

² Instituto de Geología, Universidad Nacional Autónoma de México, Ciudad Universitaria, 04510, Coyoacán, CDMX, Mexico.

³ Institute of Geological Sciences, Polish Academy of Sciences, Research Centre in Warzaw, 00818, Twarda 51/55, Poland.

⁴ School of Natural and Built Environment, Queen’s University Belfast, Belfast, BT7 1NN, Belfast, United Kingdom.

* Corresponding author: (A. Correa-Metrio) This email address is being protected from spambots. You need JavaScript enabled to view it.

How to cite this article:

Olivares-Casillas, G., Correa-Metrio, A., Zawisza, E., Wojewódka-Przybył, M., Blaauw, M., Romero, F., 2021, Environmental variability during the last three millennia in the rain shadows of central Mexico: Boletín de la Sociedad Geológica Mexicana, 73 (1), A171220. http://dx.doi.org/10.18268/BSGM2021v73n1a171220

ABSTRACT

The last three millennia have been characterized by global temperature oscillations of around one Celsius degree, and high frequency variability on precipitation. Two main temperature anomalies have been reported worldwide, the Medieval Warm Period (MWP) and the Little Ice Age (LIA), characterized by higher and lower than average temperatures, respectively. Precipitation variability has been mostly associated with El Niño anomalies in the Equatorial Pacific. These global variability modes have been modulated by regional factors such as sea surface temperatures and their interaction with continental landmasses. Understanding regional responses to these anomalies would shed light on ecosystem response to environmental variability, a paramount tool for conservation purposes on the light of modern climate change. Here we present a 3,000-year sedimentary record from Lake Metztitlán, located in a Biosphere Reserve under the rain shadow of the Sierra Madre Oriental. Cladoceran and geochemical analyses were used to reconstruct lacustrine dynamics through the time period encompassed by the record. Our record points to highly dynamic lacustrine systems, coupled with global and regional climatic variability. In Metztitlán, the MWP was associated with low lake levels and a high torrentiality of the precipitation reflected in high-frequency peaks of detrital material. The LIA was associated with an enlarged water body, probably as a result of lower evapotranspiration. Overall, global climatic variability resulted in high variability of regional precipitation and detrital input in the Metztitlán region, in turn associated with changes in lake morphometry and depth. Our record highlights the vulnerability of the area to changes in sea surface temperature of the Gulf of Mexico, and to changes in the frequency of El Niño events. Although the effects of global climate change in the region are inescapable, our results emphasize the importance of controlling anthropogenic activities as an additional source of pressure on the regional ecosystems.

Keywords: Medieval Warm Period, Little Ice Age, Cladocera, late Holocene, rain shadow, Barranca de Metztitlán Biosphere Reserve.

RESUMEN

Los últimos tres milenios se han caracterizado por oscilaciones en la temperatura global alrededor de un grado Celsius, y una variabilidad de alta frecuencia de la precipitación. Dos principales anomalías de temperatura han sido reportadas mundialmente, el Periodo Cálido Medieval (PCM) y la Pequeña Edad de Hielo (PEH), caracterizadas por temperaturas más cálidas y más frías que el promedio, respectivamente. La variabilidad de la precipitación ha estado estrechamente relacionada con las anomalías de El Niño en el Pacífico Ecuatorial. A nivel regional, estos modos de variabilidad han sido modulados por factores como la temperatura superficial de los océanos y su interacción con las masas continentales. El entendimiento de las respuestas regionales a estas anomalías ayuda a comprender la respuesta de los ecosistemas a la variabilidad global, un elemento crítico para la conservación en vista del cambio climático moderno. En este trabajo se presenta un registro sedimentario de 3000 años del lago Metztitlán, ubicado en la Reserva de la Biosfera en la Sierra Madre Oriental. Se empleo un análisis geoquímico y cladócero para reconstruir dinámicas lacustre en el período de tiempo que abarca el informe. Nuestro registro apunta a sistemas dinámicos lacustres, asociados con la variabilidad climática global y regional. En Metztitlán, el PCM estuvo asociado con bajos niveles lacustres y una alta torrencialidad de la lluvia, reflejada en picos de alta frecuencia en los indicadores de aporte detrítico. La PEH estuvo asociada con un cuerpo lacustre más extenso, probablemente resultado de una menor evapotranspiración. En general, los efectos de la variabilidad climática global se reflejaron en una alta variabilidad de la precipitación y entrada de detríticos (detrital input) en la región de Metztitlán, a su vez asociada con cambios en la morfometría y la profundidad del lago. En nuestro registro resalta la vulnerabilidad del área a cambios en la temperatura superficial de mar en el Golfo de México, y a cambios en la frecuencia de los eventos El Niño. Aunque los efectos del cambio climático en el área son inevitables, nuestros resultados enfatizan la importancia de controlar las actividades antropogénicas como una fuente adicional de presión sobre los ecosistemas regionales.

Palabras clave: Periodo cálido medieval, , pequeña edad de hielo, cladóceros, Holoceno tardío, sombra de lluvia, Reserva de la Biósfera Barranca de Metztitlán.

1. Introduction

Environmental and ecosystem dynamics are strongly influenced by the interaction of global and regional climatic variability with local processes at diverse spatial and temporal scales (Delcourt and Delcourt, 1991; Bradley, 2015). Through the Holocene, climatic variability has been expressed mainly through changes in regional precipitation regimes, but also through variability in regional temperatures (Mayewski et al., 2004; Bradley, 2015). Throughout the Holocene, changes in precipitation have been mostly associated with precessional changes in insolation, generating millennial-scale patterns of progressive increases or decreases of precipitation (Haug et al., 2001; Bradley, 2015).

Within the last two millennia, the Medieval Warm Period and the Little Ice Age represent the two main changes in global temperatures (Mann et al., 2009) that are reflected by sea surface temperatures and therefore by regional hydrological cycles (Richey et al., 2007). The sub-decadal changes in the temporal distribution of El Niño anomalies in the equatorial Pacific have been an additional source of climatic variability that manifests in regional hydrologic anomalies across the globe (Cane, 2005). The effects of these variability modes on regional environments depend on mechanisms that couple regional and global mechanisms of energy redistribution (Bradley, 2015). At a local scale, ecosystems represent envelopes that encompass biological communities and their interactions with constantly changing environmental conditions (Delcourt and Delcourt, 1991; Bush, 2003).

In recent decades, numerous paleoenvironmental studies have focused on understanding the successional development of lacustrine environments and their relationship with global and regional climatic variability (e.g. Cuna et al., 2014; Rodríguez-Ramírez et al., 2015; Correa-Metrio et al., 2016; Vázquez-Molina et al., 2016; Escobar et al., 2020; Franco‐Gaviria et al., 2020). A general property of lake development is the sedimentation

process that, in the long term, transforms aquatic environments into terrestrial realms through the gradual infilling of the lacustrine basin (Wetzel, 2001). This process is the result of the interaction between the lithosphere and the atmosphere, modulated by climatic and biological factors.

The eutrophication of freshwater ecosystems is closely associated with these factors and has been a major topic of interest, as the availability of water for ecosystem functioning and human consumption is highly dependent on this process. Eutrophication implies a change on lake trophic status of lacustrine environments as a result of the input of exogenic materials producing an increase in autotrophic productivity (Wetzel and Likens, 2000; Wetzel, 2001). This process can be driven by natural causes and varies depending on lake morphometry, regional climates, and characteristics of the catchment basin such as size and vegetation cover (Carpenter et al., 1998; Wetzel, 2001).

However, during the late Holocene, this process has been accelerated by anthropogenic input of nutrients derived from intensive and extensive human activities (e.g. Carpenter et al., 1998; Franco‐Gaviria et al., 2020). Sediment accumulations in lacustrine environments contain biological, chemical, and physical evidence of lake development through time (Cohen, 2003). The study of this sedimentary evidence offers an opportunity to understand processes and factors involved in lake development through time scales beyond the historic observations.

Numerous lakes are distributed across Central Mexico, offering the opportunity to study lake formation and its development in a highly complex landscape, marked by an intricate history of human occupation and natural environmental variability (Caballero‐Rodríguez et al., 2017). Here we present the environmental history of Lake Metztitlán, located in the Biosphere Reserve of the same name. Understanding the temporal environmental dynamics of the region is crucial for its conservations as it has been classified as one of the most vulnerable biosphere reserves in the light of ongoing global climate change (Esperon-Rodriguez et al., 2019).

We retrieved a sediment core and analyzed it for geochemical composition and cladoceran remains. Our aim was to understand regional and local environmental variability through time, and their coupling with global environmental variability through the late Holocene. Overall, we intend to investigate the effects of climatic variability of the last three millennia on regions that are today characterized by an arid climate. Also, we wanted to evaluate the effects of regional climatic variability on the development of local lacustrine conditions.

2. Study area

Lake Metztitlán is located within the Barranca de Metztitlán Biosphere Reserve (20°42’ - 20°39’ N, 98°53’ - 98°49’ W, at an elevation of 1,253 m asl), State of Hidalgo, on the western flank of the Sierra Madre Oriental, Mexico (Figure 1). The regional geology is composed of Mesozoic marine sedimentary rocks, particularly limestones and turbidites overlaid by clastic and basaltic materials of Neogene igneous origin, as well as alluvial Quaternary deposits (Carrasco-Velázquez et al., 2008).

The regional relief is composed of parallel ranges and plateaus delimited by faults, constituting a highly dissected landscape with evidence of multiple slumps and detrital flows (García Arizaga et al., 1996). Lake Metztitlán is a large (~2,949 ha) and shallow (^zmean = 11 m) water body originated by a massive slump of one of the flanks of the valley of the Metztitlán River that resulted in a natural dam (Suter, 2004).

Two regional, E-W striking, south-dipping normal faults feature pronounced scarps at the northeastern margin of the lake which seem associated with the damming slump (Suter, 2004). According to Suter (2004), the damming of the river occurred between 500 and 1,100 years BP. Currently, the slumped material covers an area of approximately 2.5 km² (Figure 1).

Figure 1. Study area. A) Continental context: 1. Lake Metztitlán; 2. Lake Aljojuca (Bhattacharya and Byrne, 2016); 3. Lake Pallcacocha (Moy et al., 2002); 4. Pigmy Basin (Richey et al., 2007).; 5. Cariaco Basin (Haug et al., 2001). B) Regional topographic context of Lake Metztitlán, showing the coring site (red star), and the damming slump (yellow polygon). C) Regional climatology of the region from Metztitlán meteorological station, showing monthly evapotranspiration, precipitation, and mean temperature (data from Servicio Meteorológico Nacional, 2020).

The lake is an eutrophic freshwater body of low transparency, with a polymictic mixing regime that hinders the formation of either a thermocline or an oxycline (Sigala et al., 2017). The lake is fed by Metztitlán River, whereas the main drainage of water occurs by infiltration through the debris of the slumping with a resurgence of water by the northeast side of the lake, around 250 m below the level of the lake water (Suter, 2004).

The regional climate is warm and dry with pronounced torrential summer rains. Generally dry conditions are associated with the rain-shadow effect of the Sierra Madre Oriental. On the other hand, the seasonal precipitation regime derives from the northerly displacement of the Intertropical Convergence Zone during the boreal summer (CONANP, 2003). Additionally, during the hurricane season of the Atlantic, moisture from tropical storms ascends along the eastern flank of the Sierra Madre Oriental, bringing additional precipitation. Around the lake, moisture deficit persists throughout the entire year with an annual precipitation of 406 mm, and a mean evapotranspiration of 1,788 mm (Figure 1; Servicio Meteorológico Nacional, 2020).

The lake maintains water throughout the year because of permanent water input from the higher areas of the catchment basin. Mean annual temperature is 20.7 °C, with monthly mean minimum and maximum of 16.5 and 23.5 °C in January and May, respectively (Figure 1; Servicio Meteorológico Nacional, 2020). The steep altitudinal gradient and the highly dissected landscape have created the conditions for a complex mosaic of arboreal vegetation and shrublands. The mean vegetation types are tropical deciduous forests (1100-1500 m asl), lower montane shrublands (1500-1800 m asl), xerophytic shrublands (1000-1800 m asl), and coniferous forests (1800-2500 m asl) (CONANP, 2003; Rzedowski, 2006).

3. Methods

In May 2017, a 5.32-m-long core was retrieved from Lake Metztitlán using a modified Livingstone piston corer (Colinvaux et al., 1999). The core was stored at the Instituto de Geología at a temperature of 2 °C for preserving the sedimentary evidence. The core was longitudinally sectioned for lithological description that included description of sediment color, qualitative description of sediment texture, and sediment reaction to HCl. Seven samples were extracted at different depths along the core for chronological control of the sedimentary sequence. Organic matter from each sample was dated by analyzing radiocarbon through accelerator mass spectrometry (AMS). Radiocarbon ages were calibrated using the IntCal20 curve (Reimer et al., 2020), and a Bayesian age-depth model was constructed using bacon (Blaauw and Christen, 2011). All ages were expressed in calibrated years before 1950 CE (cal BP, hereafter).

Every two cm along the core, a one-cm³ sample was extracted for geochemical analyses. Samples were freeze-dried and ground to silt texture using a mortar and pestle. Elemental concentrations were determined using X-Ray fluorescence with a Thermo Scientific Niton XL3t, making three repeated measurements per sample. Elemental concentrations of Ti and Zr were transformed to ratios using Ca, accounting for dilution effects and variability associated with elements undetected by the XRF (Löwemark et al., 2011). Elemental concentrations of Fe and Mn were transformed to ratios using Ti, aiming to detect possible post-depositional alterations of these elements, as they are sensitive to changes in oxygen availability in the mud-water interface (Tribovillard et al., 2006; Algeo and Liu, 2020). Also, the ratio Rb/Ti was calculated, aiming to an indicator of dominant sediment grain size, whereas the Sr/Rb ratio was calculated as an indicator of dominance of chemical weathering (Kylander et al., 2011).

The sedimentary sequence was sampled every 10 cm for analysis of cladocerans. From each selected depth, a one-cm³ sample was extracted and treated according to standard protocols (Frey, 1986; Szeroczyńska and Sarmaja-Korjonen, 2007). Whereas 10% KOH was used for deflocculating the samples, excess of carbonates in samples was removed by washing them with 10% HCl. The remaining material was sieved through a 35- µm mesh, taken to a 5 ml volume and tainted with safranine to better observe the morphology of cladoceran remains. For each sample, remains were counted until 200 individuals were reached by assembling headshields, shells, postabdomens, and postabdominal claws. Counts of cladocerans were transformed to concentrations (individuals/cm³). Taxa were identified at species level using illustrated taxonomic guides (Elías-Gutiérrez et al., 2008; Wojewódka et al., 2020a, 2020b).

The geochemical dataset was submitted to a principal component analysis (PCA), aiming to summarize the variability of the dataset and determine the relationships among elemental concentrations and ratios (Jolliffe, 1986). This technique redistributes the variability of the dataset through new orthogonal components that are linear combinations of the original attributes. Attributes are represented by vectors in the newly ordinated space, with vector magnitude representing variability of the original attribute and angles among vector representing association among attributes. Sample scores along PCA axes were used to infer environmental changes through time. All analyses were performed using R (R Core Team, 2020), especially the vegan package (Oksanen et al., 2019).

4. Results

4.1. LITHOLOGICAL PROFILE AND AGE-DEPTH MODEL

The core was mainly composed of a massive layer of clays overlying two thinner layers of sand and lime (Figure 2). The sediments were characterized by a moderate reaction to HCl, indicating a relatively high content of carbonates in the sediment. According to the texture of sediments, the sequence was divided in three units: U1, from 532 to 441 cm below lake floor (blf, hereafter), characterized by dark silt; U2, from 441 to 398 cm blf, composed by coarse dark sands; and U3, from 398 cm blf to the top of the record, composed of homogeneous bands of light clays.

Figure 2. Stratigraphy and age-depth model of the sediment core. The gray rectangle highlights the sediment that represents a probable instant deposit after the damming of the lake.

According to the age-depth model, from 398 cm blf to the top , the sequence comprises a record of the last ~3000 years through the upper 398 cm. Calibrated radiocarbon dates on units 1 and 2 showed reversals and were therefore excluded from the analyses (Figure 2 and Table 1).

Given the contrasting nature of the sediments and the high instability of the hillsides (Suter, 2004), as well as the torrentiality of the regional climate, it is likely that units 1 and 2 were product of an instantaneous event. Thus, we based our temporal reconstruction exclusively on the sediments of U3. The landslide that dammed the river and created the lake was a mass remotion process that extended for more than 2 km across the valley (Figure 1; Suter, 2004).

Table 1. Radiocarbon dates for the sediments of Lake Metztitlán.

This landslide was probably followed by instability of the hillsides and of the removed landmass, creating subsequent mass movements inside the lake, mixing detritus flows with older sediments from the bottom of the lake. Furthermore, the coarser sediments that characterized U1 and U2 suggest these sediments were deposited by processes of more energy than lacustrine progressive sedimentation.

Whereas one of the radiocarbon dates within U1 was statistically undifferentiable from the basal age, the other was younger (Figure 2), probably as a result of illuviation of organic matter through the coarser materials of U1 and U2. Given that with our results we cannot offer a reliable explanation for units U1 and U2, from this point on, we will base our results and discussion only on U3.

4.2. GEOCHEMISTRY

Repeated measurement of Ca, Fe, K, Mn, Pb, Rb, Sr, Ti, Zn, and Zr were statistically consistent and showed significant variability through time. Concentrations of Ca, Sr, K, and Rb were highly variable through the record, showing a general increasing trend towards present (Figure 3). Ti and Zr concentrations were also highly variable, showing an abrupt decrease from 236 to 209 cm blf (from ~1620 to 1350 cal BP).

On the other hand, Fe and Mn concentrations showed little variability through the record, with the latter showing a substantial increase from 46 cm blf to the top of the record (last ~250 years). Pb and Zn concentrations were highly variable from the bottom of the record to 283 cm blf (from ~3000 to 2000 cal BP), whereas they showed relative stability from 283 cm blf to the top of the record (last ~2000 years).

Figure 3. Elemental concentrations and elemental ratios of the sedimentary record of Lake Metztitlán.

PC1 and PC2 explained 41.0 and 20.8% of the variance of the geochemical dataset, respectively, and both were statistically significant under the broken-stick model (Figure 4). Althrough PC3 and PC4 were also statistically significant, we decided to exclude them from the interpretation as they were associated to very low amounts of variance (Figure 4).

Concentrations of Fe, Pb, Ti, Zr, and Zn were positively correlated and independent from the Fe/Ti ratio, Rb concentrations were positively associated with K, Mn, and Fe/Ti ratio. Ca and Sr concentration showed a positively association between them, and were in turn negatively associated with Fe, Pb, Ti, Zr, and Zn.

The negative association between Ca with Ti/Ca and Zr/Ca ratios (Figure 4), implies diluting effects by high concentrations of these elements in their original concentrations (Löwemark et al., 2011).

Figure 4. Correlation biplot of the principal component analysis of the geochemical record of Lake Metztitlán. The inset shows the variance associated with each component (grey bars), and the null variance associated with the broken-stick model.

Ti/Ca and Zr/Ca ratios were high from ~2800 and 1620 cal BP, and from ~1350 and 900 cal BP, whereas they showed low values between ~3000 and 2800 cal BP, ~1620 and 1350 cal BP, and through the last ~900 years. The Rb/Ti ratio was in antiphase with Ti/Ca and Zr/Ca. Fe/Ti and Mn/Ti ratios showed low variability throughout the record, yet showing a progressive increase through the last ~900 years.

4.3. CLADOCERA REMAINS

Nine Cladocera taxa where found in 54 samples analyzed through the entire record. These taxa were classified in three and six planktonic and littoral species, respectively. Whereas the plankton group was composed of Daphnia sp. (Müller, 1785), Liederobosmina sp. (Brtek, 1997), and Bosmina sp. (Baird, 1845), the littoral was composed of Alona sp. (Baird, 1843), A. ossiani (Sinev, 1998), Ovalona glabra (Sars, 1901), Chydorus cf. sphaericus (Müller, 1776), Leydigia cf. louisi louisi (Jenkin, 1934), and L. cf. striata (Birabén, 1939). Two of Leydigia species, L. cf. louisi louisi and L. cf. striata are benthonic bottom dwellers. Liederobosmina sp. and Bosmina sp. were the most abundant taxa through the record, reaching concentrations higher than 2,000 individuals per cm³, with an abrupt increase from 150 cm blf to the top of the record (last ~800 years) (Figure 5). A. ossiani and Daphnia sp. were present only from 381 to 347 and from 145 to 120 cm blf (~2750 to 2500 and ~750 to 650 cal BP, respectively). Samples from 236 to 209 cm blf (~1550 to 1350 cal BP) were barren for Cladocera.

Figure 5.  Cladocera concentrations in the sedimentary record of Lake Metztitlán. Sample scores of the principal component analysis along axes 1 and 2 are shown.

5. Discussion

5.1. DETRITAL AND AUTHIGENIC INPUT

Concentrations of K, Rb, Ti, and Zr in sedimentary deposits are usually associated with the input of detrital materials, mainly associated to erosive processes (Rothwell and Rack, 2006; Kylander et al., 2011). Rb, Ti, and Zr are conservative elements that, after deposited, maintain a signal of the provenance, as most of their abundance corresponds to detritus from parental rocks (Rothwell and Rack, 2006; Boës et al., 2011; Algeo and Liu, 2020). In our record, whereas Ti and Zr were positively associated, Rb was decoupled from these elements, even when concentrations were standardized to Rb/Ti (Figure 4). This observation suggests that concentrations of these elements in our record could have been associated with the energy of the transporting agents causing a size sorting of the transported materials. Whereas Ti and Zr are usually associated with relative coarse and heavy grains, Rb mainly associates with finer and lighter sediment fractions (Rothwell and Rack, 2006; Calvert and Pedersen, 2007; Kylander et al., 2011). Given its high capacity of sorption in clay minerals, Rb is an element of low mobility, but it also replaces K in feldspars of variable grain size (Rothwell and Rack, 2006; Calvert and Pedersen, 2007; Kylander et al., 2011). Thus, the positive correlation between K and Rb in the sediments of Lake Metztitlán (Figure 4) suggests these elements were associated with the input of fine material into the sedimentary deposit.

Ca and Sr concentrations in lacustrine sediments are the result of both clastic input from the erosion of parental materials and authigenic precipitation of carbonates (Rothwell and Rack, 2006; Algeo and Liu, 2020). Carbonate precipitation occurs when dissolved minerals in water reach saturation concentrations, and physical, chemical and biological processes facilitate their deposition in the lake floor (Cohen, 2003). These conditions often associate with low lake levels. In our record, negative correlations among detrital elements (Ti and Zr) and Ca and Sr (Figure 4), indicate that these latter elements are mainly result of authigenic precipitation of carbonates. Thus, whereas Ti/Ca and Zr/Ca ratios reflect the detrital input, incorporating the Ca diluting effect, the authigenic flux is represented by Ca, reflecting increased regional evapotranspiration, and probably low lake levels.

The evident separation of elements, associated with authigenic and detrital flows in the PCA (Figure 4), suggests that PC1 directly reflects erosive processes. Overall, erosion is controlled by weathering and mechanisms of detrital transport, the latter mainly associated with the erosive action of precipitation and wind (Sageman and Lyons, 2003). Rainfall and wind represent the main erosive agents in humid and arid areas, respectively (Rose, 1985; Ahnert, 1998). Although, through the late Holocene, both erosive agents might have played a role in Metztitlán, aeolian erosion in the area was probably diminished by the regional geomorphology, whereas the torrential precipitation regime has surely been associated with large amounts of pluvial laminar and linear erosion (Ahnert, 1998). Thus, detrital input into our record is most likely associated with regional precipitation, and positive scores of PC1 are a result of phases of a wetter and/or a more torrential regional climate. Given the sparse vegetation associated with the regional climate (Rzedowski, 2006), wetter periods would be associated with more regional erosion. Transport and deposition of a high load of sediments would result in a faster infilling of the lacustrine basin, impeding the deepening of the water body. PC2, on the other hand, was negatively correlated with Fe, Fe/Ti, Mn, Rb, Rb/Ti (Figure 4), suggesting it was also an indicator of detrital input, although it was probably associated with finer materials.

5.2. REDOX CONDITIONS

Fe and Mn can be affected by processes of reduction and oxidation that take place after deposition (Boës et al., 2011). Under conditions of oxygen availability in the water column, Mn and Fe precipitate in the sediments as oxyhydroxides of magnesium and ferric iron. On the other hand, under anoxic conditions, these elements redissolve and are liberated back to the water column (Algeo and Liu, 2020). Given their detrital origin, an unaltered signal of these elements would imply positive associations with detrital elements that are not prone to post-depositional processes. In our record, both Fe and Mn are decoupled from Ti and Zr concentrations, implying a large amount of variance caused by diagenetic processes. Thus, the Fe/Ti and Mn/Ti ratios could be used as indicators of post-depositional processes, with low values pointing to redissolution of these elements caused by anoxic conditions (Algeo and Liu, 2020). Additionally, in the record of Lake Metztitlán, peaks of Mn/Ti and Fe/Ti ratios were in phase with high concentrations of Ca, suggesting enrichment and fixation of Mn and Fe through authigenic precipitation of carbonates (Cohen, 2003). Overall, changes in Mn/Ti and Fe/Ti ratios can be interpreted as changes in the oxygen content of the mud-water interface, or as changes in the precipitation of carbonates. These processes are closely related as carbonate precipitation usually takes place under low lake levels, which in turn would improve oxygenation of water by wind mixing.

5.3. LAKE SURFACE, LITTORAL DEVELOPMENT, AND LAKE PRODUCTIVITY

Cladocerans are a highly specialized biological group, systematically distributed across lacustrine environments (Whiteside and Swindoll, 1988; Korhola and Rautio, 2001). Overall, the lake environment can be divided in the littoral and limnetic or open water zones (Wetzel and Likens, 2000), and the distribution of Cladocera species across these areas is driven by their ecological affinities (Korhola and Rautio, 2001). Thus, the turnover of taxa in our record offers elements to reconstruct environmental dynamics of the lake through the time period spanned by our record. Members of the Chydoridae family are mainly distributed across shallow littoral areas with preferences for particular microhabitats such as rocks, mud, and aquatic vegetation (Korhola and Rautio, 2001; Adamczuk, 2014). Contrastingly, open waters are colonized by planktonic species that mostly belong to the families Bosminidae and Daphniidae (Whiteside and Swindoll, 1988; Korhola and Rautio, 2001). Thus, changes of relative abundant of littoral vs. planktonic taxa are informative of changes in factors that modify the habitat of cladocerans, such as lake morphology (Hofmann, 1987). Given the morphology of Lake Metztitlán, abundance increases of planktonic taxa are probably not associated with a deepening of the lake, but more likely an enlargement of the open water area of the lake.

The higher Bosminidae biomass often reflects increase of nutrients (Bays and Crisman, 1983) while the turnover of bosminid taxa is usually associated as an indicator of changes in the trophic status of lakes (Korhola and Rautio, 2001). However, species of Liederobosmina appearing in the region may inhabit, both, environment poor (Liederobosmina tubicen) and rich in nutrients (Liederobosmina hagmanni); and some species (Liederobosmina huaronensis) did not show prominent relation with trophic state, whereas Bosmina sp. is more commonly found under eutrophic conditions (Beaver et al., 2018). On the other hand, independently of trophic status, species group of Bosminidae can coexist in open water environments (Hofmann, 1987; Szeroczyńska, 2002). The aforementioned lack of species turnover among bosminid taxa suggests abundances of Liederobosmina sp. and Bosmina sp. through the record were more associated with the expansion of the open water zone than with changes in the trophic status of the lake. Aquatic vegetation increases the diversity of cladocerans as they offer shelter for taxa adapted to littoral environments (Adamczuk, 2014). Littoral macrophytes can regulate the trophic state of lakes as they retain nutrients, favoring clear waters (Sand-Jensen and Borum, 1991). The Cladocera genus Alona, and O. glabra are commonly associated with aquatic vegetation (Korhola and Rautio, 2001). On the other hand, C. cf. sphaericus is a dualistic species, as it can be found under conditions of low nutrient concentration and also associated with macrophytes (Korhola, 1999). Thus, in our record, littoral species (Alona sp., A. ossiani, C. cf. sphaericus, and O. glabra) are indicating the development of mature littoral areas, probably resulting from advances of the water surface across the southern alluvial plain of Metztitlán valley. During times of expansion of the lake, the alluvial plain was probably covered by a shallow water layer that allowed the establishment of mats of aquatic vegetation that promoted the growth of populations of Alona sp., A. ossiani, C. cf. sphaericus, and O. glabra. Lastly, the benthic species found in the record of Lake Metztitlán, L. cf. striata and L. cf. louisi louisi, have been reported for turbid environments with high nutrient input (Kotov et al., 2003; Kotov and Elías-Gutiérrez, 2004; Wojewódka et al., 2016; Sigala et al., 2017). In our record, these taxa were probably associated with a highly productive environment, which was in turn promoted by nutrients brought by the river together with the detrital input.

5.4. THE LAST 3,000 YEARS AT LAKE METZTITLÁN

According to a visual inspection of the variability of geochemical ratios, we identified five main stages in the evolution of the lake. These stages were characterized by relatively homogeneous conditions and were apparently dominated by changes in the regional moisture availability. The high variability of Cladocera assemblages through time allowed for the identification of an additional stage towards the end of the record. This discretization of the record was supported by PCA scores, which show consistency between the geochemical and biological signals (Figure 6).

5.4.1. ZONE I (398-256 CM BLF; ~3000-1800 CAL BP)

Through this time interval, higher abundances of L. cf. louisi louisi and L. cf. striata suggest a sha-llow environment with rather weakly developed submerged vegetation and probably low water transparency (Figure 6). This zone overlays the instantaneous deposit of coarse material that characterizes the bottom of the core (Figure 2). Such instantaneous event filled the bottom of the lake bringing a high load of nutrients derived from the fresh material recently eroded. Decreasing Rb/Ti ratios (Figure 3) suggest that progressively coarser materials were deposited into the lake, coinciding with the succession of cladoceran assemblages from littoral to benthic and planktonic species. Indeed, the progressive decrease of the Rb/Ti ratio indicates progressively increasing of grain size in the sediments, whereas the slight increase in lake level probably caused anoxic conditions in the mud-water interface, reflecting in the decrease of the Fe/Ti ratio (Figure 6). Increasing PC1 scores reflect the progressive increase of Ti/Ca and Zr/Ca ratios (Figures 3 and 6), suggesting relatively wet conditions.

Relatively wet conditions have been reported for the region (Bhattacharya et al., 2015), which would have allowed the infilling of the lake after the landslide damming. From ~2550 to 1800 cal BP, as benthonic Cladocera taxa diminished, pelagic and littoral taxa increased (Figure 6). Increasing concentrations of Liederobosmina sp. and Bosmina sp. suggest a slight increase in lake extension at that time.

Figure 6. The paleoenvironmental record of Lake Metztitlán in the regional and continental context. A, B, and C, concentrations of planktonic, benthic and littoral Cladocera taxa, respectively, in the sediments of Lake Metztitlán; samples barren of cladocerans are represented by gray diamonds. D. Red intensity of the sediments of Lake Pallcacocha, Ecuadorian Andes, as a proxy of El Niño frequency (Moy et al., 2002). E. Ti concentration in the sediments of Cariaco Basin (Haug et al., 2001). F. Ti/Ca ratios (red) and PC1 sample scores (black) of the sedimentary record of Lake Metztitlán. G. Sea surface temperatures of the Gulf of Mexico reconstructed from the sediments of Pigmy Basin (Richey et al., 2007). H. d18O in the sedimentary record of Lake Aljojuca, Serdán Oriental Basin (Bhattacharya and Byrne, 2016). I. Rb/Ti ratios in the sediments of Lake Metztitlán.

Overall, this stage shows the initial development and progressive infilling of the lacustrine basin, which could have been associated with the damming of the river by the massive landslide in the context of relatively high moisture availability in the region (Bhattacharya et al., 2015).

Nevertheless, as we do not have elements to associate U1 and U2 to the landslide, it is also possible that the lake was much older than the time spanned by our sedimentary record. The relatively wet conditions could be associated with a low frequency of El Niño events in the equatorial Pacific (Moy et al., 2002). Today, El Niño events are a major contributor to precipitation variability in Mexico (Magaña et al., 2003), bringing pronounced droughts to the Metztitlán region (CONANP, 2003). According to the record of Lake Pallcacocha (Moy et al., 2002), El Niño frequency was apparently low through most of zone 1 of our record, probably producing a low frequency of droughts that facilitated the progressive infilling of the lake. Superimposed on the trend towards wetter conditions, peaks of Ti/Ca and Zr/Ca suggest a high variability in the input of detrital material, that could have been associated with either short lived intermissions of even wetter conditions or torrential events.

High variability of Cladocera assemblages further support our interpretation of high variability of lake level and extension. These erosive peaks were probably associated with the high frequency variability mode of the Intertropical Convergence Zone revealed by the Cariaco basin Ti record (Haug et al., 2001).

5.4.2. ZONE II (256-209 CM BLF, ~1800-1350 CAL BP)

The bottom of this zone (~1800-1620 cal BP) showed an abrupt decrease of Cladocera remains, accompanied by low detrital input, reflected in low scores of PC1 (Figure 6). The peak of terrigenic input identified between ~1750 and 1620 cal BP probably corresponds to torrential events with high erosive capacity given the generally dry regional climate. Such torrentiality could be associated with reinvigoration of the El Niño frequencies recorded in lake Pallcacocha (Moy et al., 2002). Although our PC1 peak predated the increase of El Niño frequency, the mismatch could be associated with the uncertainty of our age-depth model. From ~1700 to 1400 cal BP, samples were barren of cladocerans and Ti/Ca, Zr/Ca, and PC1 decreased (Figure 6), indicating low detrital input to the lake, and probably an almost complete desiccation of the water body. This period coincides with an abrupt onset of dry conditions that have been reported for the basin of Serdán Oriental at ~1600 cal BP (Bhattacharya et al., 2015). This same pattern has been reported for other sites at elevations above 1000 m asl (e.g. Metcalfe et al., 2010; Park et al., 2019), thus, indicating a regional drought.

5.4.3. ZONE III (209-164 CM BLF, ~1350-900 CAL BP)

A recovery of the water body was evidenced by the recolonization of cladocerans and the increase of indicators of detrital input, Ti/Ca, Zr/Ca, and PC1 (Figure 6). Whereas L. cf. louisi louisi suggests a relatively shallow lake, moderate concentrations of Liederobosmina sp. and Bosmina sp. suggest a large open water area. Three consecutive and relatively evenly spaced peaks of indicators of detrital input suggest a high-frequency mode of variability of the regional hydrological cycle. A similar mode of variability of successive centennial droughts for this time period has also been reported for the Yucatan Peninsula (Hodell et al., 2001), possibly implying solar activity as a common driver. Peaks of detrital materials (Ti/Ca and Zr/Ca) could also have been associated with torrential rains in a matrix of generally dry conditions. Indeed, absence of littoral taxa in most of the samples could have been caused by the massive input of sediments causing frequent disturbances to the littoral environment, impeding the establishment of littoral vegetation and, thus, preventing the establishment of a community associated with these areas.

The time period comprised by this zone roughly coincides with the Medieval Warm Period, a time when global temperatures were around one degree higher than those that characterized the 20th century (Mann et al., 2009). In Central Mexico, this time period was characterized by an intensification of the regional droughts (Metcalfe et al., 2010; Rodríguez-Ramírez et al., 2015; Bhattacharya et al., 2015; Park et al., 2019), however, the Metztitlán record shows conditions that seem more erosive. Such high erosion was probably associated with an intensification of the regional summer rains and therefore an accentuation of precipitation seasonality. Today, the main source of moisture for the Metztitlán region are the easterlies that bring humid air from the Gulf of Mexico (Douglas et al., 1993). During the MWP, sea surface temperature in the Gulf of Mexico was around 2 °C higher than modern (Richey et al., 2007), producing higher regional convective activity, and thus explaining the high erosive activity reflected in the sedimentary record of Lake Metztitlán.

5.4.4. ZONE IV (164-47 CM BLF, ~900-200 CAL BP)

The increase of all Cladocera species suggests higher water levels in the lake with a higher diversity of habitats. Substantial increases of Liederobosmina sp. and Bosmina sp. (Figure 5) suggest a relative deepening of the lake environment, as well as a widening of the lake surface. Regional records suggest an increase of regional moisture availability that coincides with the onset of this stage of the lake (Bhattacharya et al., 2015). The development of mature littoral areas is evidenced by permanent presence of both benthic and littoral Cladocera species. Indeed, historic accounts report that the lake occupied much of the southern areas that are today occupied by towns and extensive crops (Acuña, 2017). Relative flat and diminishing ratios of Ti/Ca, Zr/Ca and PC1 scores suggest that moisture availability was low but constant (Figures 3 and 6). Increasing Rb/Ti ratios suggest a progressive increase of finer material in the detrital input. Authigenic precipitation of carbonates is suggested by low values of the Ti/Ca and Zr/Ca ratios together with increasing ratios of Fe/Mn (Cohen, 2003). Overall, under generally dry conditions suggested by carbonate precipitation, moisture availability was enough for maintaining a large lake size that allowed the development of littoral areas. Most regional records show high variability through the time period comprised by this zone (Metcalfe et al., 2010; Cuna et al., 2014; Rodríguez-Ramírez et al., 2015; Bhattacharya et al., 2015; Park et al., 2019; Wogau et al., 2019). Regional environmental variability during this time interval was closely associated with anthropogenic activity, as well as with changes in moisture availability caused by the Little Ice Age (LIA) (Mann et al., 2009). In Metztitlán, the littoral community was established by the beginning of the LIA (~600 cal BP, Figure 6), suggesting that this anomaly in the region was associated with the enlargement of the water body.

Under regional colder conditions, evapotranspiration probably decreased, ameliorating the effects of the relatively drier climate. The most important event that took place through this stage was marked by a peak of indicators of terrigenic detrital input between ~450 and 400 cal BP (1500-1550 CE, Figure 6).

This peak is preceded by the Spörer solar minimum (Bard et al., 2000), which in the region was associated with dry conditions (Bhattacharya et al., 2011; Rodríguez-Ramírez et al., 2015) that were probably associated with a low vegetation vigor in the catchment basin of Lake Metztitlán. Thus, once dry conditions were overcome, there was a temporal invigoration of the hydrological cycle that strengthened the erosive capacity of rain, acting on a more erodible landscape.

5.4.5. ZONE V (47 CM BLF TO THE TOP, THE LAST ~200 YEARS).

The disappearance of almost all benthic and littoral taxa and rebounding of planktonic species (Figure 6) indicate a relatively deep and open lacustrine environment, with a weakly developed littoral zone. After the LIA, global temperatures rebounded and as a consequence, evapotranspiration rates probably increased, implying the retraction of the water body from the alluvial plain previously covered by water. General low moisture availability is reflected in low Ti/Ca and Zr/Ca ratios, with high Rb/Ti ratios suggesting a preferential deposition of fine materials in the coring site. As the valley was infilled by sediments, the influence of the mouth of the river on the coring site was probably progressively lower.

However, during the 19th and 20th centuries, the regional hydrology was highly modified by hydraulic infrastructure (Waitz, 1947). Although infiltration through the obstruction created by the landslide continued to be the main control of water level (Suter, 2004), man-made infrastructure probably played a major control on the deposition process that took place through this time period.

Indeed, besides natural process, the disappearance of littoral communities was probably also promoted by the occupation of the alluvial plain by productive activities.

6. Conclusion

Our multiproxy reconstruction based on the sedimentary record of Lake Metztitlán shows a history of high environmental variability through the last 3,000 years. The chronological model suggests that Lake Metztitlán has been in place for at least three millennia, shedding new light on the age of the ecosystem, formerly located at between ~1100 and 500 cal BP (Suter, 2004). Within this time frame, our results demonstrate that, although the regional climate has always been arid, it has been subjected to multiple forcings that have resulted in a highly dynamic system. Whereas the millennial scale trend in the regional climate does not seem to follow a specific pattern, the effects of decadal-to centennial drivers are evident.

Our results suggest that the local environmental evolution reflected by the cladoceran assemblages was highly conditioned by changes in the regional hydrological cycle. Changes in the regional precipitation regime conditioned the productivity and morphology of the water body, producing cycles of high erosion associated with lake infilling and development of littoral areas.

During the MWP, the lake was shallow, whereas the regional climate was characterized by successive peaks of drought and wet conditions. Contrastingly, the LIA was characterized by a relatively wet climate and an expansion of the lake surface, with a very wet period during the intermission between the solar minima of Maunder and Spörer. Overall, our record demonstrates the close association between local and regional environmental conditions. Through the last three millennia, the region of the Barranca de Metztitlán Biosphere Reserve has been submitted to a high environmental variability.

Acknowledgements

This research was funded by Consejo Nacional de Ciencia y Tecnología (grant No. 256406). We thank Luis Gerardo Martínez Jardines and Astrid Vázquez Salgado for their assistance with XRF analyses. We are grateful to Comisión Nacional de Áreas Naturales Protegidas (CONANP, Mexico) for granting us access to the study area.

References

Acuña, R., 2017, Relaciones Geográficas del Siglo XVI: Antequera tomo 2: México, Instituto de Investigaciones Antropológicas, 632 p.

Adamczuk, M., 2014, Niche separation by littoral-benthic Chydoridae (Cladocera, Crustacea) in a deep lake - potential drivers of their distribution and role in littoral-pelagic coupling: Journal of Limnology, 73, 490-501. https://doi.org/10.4081/jlimnol.2014.884

Ahnert, F., 1998, Introduction to Geomorphology: London, Arnold, 352 p.

Algeo, T.J., Liu, J., 2020, A re-assessment of elemental proxies for paleoredox analysis: Chemical Geology, 540, 119549. https://doi.org/10.1016/j.chemgeo.2020.119549

Bard, E., Raisbeck, G., Yiou, F., Jouzel, J., 2000, Solar irradiance during the last 1200 years based on cosmogenic nuclides: Tellus B, 52, 985-992. https://doi.org/10.3402/tellusb.v52i3.17080

Bays, J.S., Crisman, T.L., 1983, Zooplankton and trophic state relationships in Florida lakes: Canadian Journal of Fisheries and Aquatic Sciences, 40, 1813–1819. https://doi.org/10.1139/f83-210

Beaver, J.R., Renicker, T.R., Tausz, C.E., Vitanye B.T., 2018, Distribution of six taxa in the family Bosminidae Baird (Crustacea: Branchiopoda: Anomopoda) in the plankton of lakes and reservoirs within the continental United States, including expanded range of the invasive cladoceran Bosmina (Eubosmina) coregoni Baird: Zootaxa, 4407(4), 506-520. https://doi.org/10.11646/zootaxa.4407.4.3

Birabén, M., 1939, Los Cladóceros de la familia “Chydoridae”: Revista de la Sociedad Argentina de Ciencias Naturales, 17, 651-671.

Bhattacharya, T., Beach, T., Wahl, D., 2011, An analysis of modern pollen rain from the Maya lowlands of northern Belize: Review of Palaeobotany and Palynology, 164, 109-120. https://doi.org/10.1016/j.revpalbo.2010.11.010

Bhattacharya, T., Byrne, R., 2016, Late Holocene anthropogenic and climatic influences on the regional vegetation of Mexico’s Cuenca Oriental: Global and Planetary Change, 138, 56-69. http://dx.doi.org/10.1016/j.gloplacha.2015.12.005

Bhattacharya, T., Byrne, R., Böhnel, H., Wogau, K., Kienel, U., Ingram, B.L., Zimmerman, S. 2015, Cultural implications of late Holocene climate change in the Cuenca Oriental, Mexico: Proceedings of the National Academy of Sciences, 112, 1693-1698. https://doi.org/10.1073/pnas.1405653112

Blaauw, M., Christen, J.A., 2011, Flexible paleoclimate age-depth models using an autoregressive gamma process: Bayesian Analysis, 6(3), 457-474. https://doi.org/10.1214/11-BA618

Boës, X., Rydberg, J., Martinez-Cortizas, A., Bindler, R., Renberg, I., 2011, Evaluation of conservative lithogenic elements (Ti, Zr, Al, and Rb) to study anthropogenic element enrichments in lake sediments: Journal of Paleolimnology, 46, 75-87. https://doi.org/10.1007/s10933-011-9515-z

Bradley, R.S., 2015, Paleoclimatology: reconstructing climates of the Quaternary: Oxford, UK, Elsevier Inc., 675 p. https://doi.org/10.1016/s0074-6142(99)x8001-6

Bush, M.B., 2003, Ecology of a Changing Planet: New Jersey, Prentice Hall, 477 p.

Caballero‐Rodríguez, D., Lozano‐García, S., Correa‐Metrio, A., 2017, Vegetation assemblages of central Mexico through the late Quaternary: modern analogs and compositional turnover: Journal of vegetation Science, 28, 504-514. https://doi.org/10.1111/jvs.12515

Calvert, S.E., Pedersen, T.F., 2007, Chapter fourteen elemental proxies for palaeoclimatic and palaeoceanographic variability in marine sediments: interpretation and application: Developments in Marine Geology, 1, 567-644. https://doi.org/10.1016/S1572-5480(07)01019-6

Cane, M.A., 2005, The evolution of El Niño, past and future: Earth and Planetary Science Letters, 230, 227-240. https://doi.org/10.1016/j.epsl.2004.12.003

Carpenter, S.R., Caraco, N.F., Correll, D.L., Howarth, R.W., Sharpley, A.N., Smith, V.H., 1998, Nonpoint pollution of surface waters with phosphorus and nitrogen: Ecological applications, 8(3), 559-568. https://doi.org/10.1890/1051-0761(1998)008[0559:NPOSWW]2.0.CO;2

Carrasco-Velázquez, B.E., Martínez-Hernández, E., Ramírez-Arriaga, E., Solé Viñas, J., 2008, Estratigrafía de la formación Metztitlán del Plioceno (estado de Hidalgo, centro-este de México): Boletín de la Sociedad Geológica Mexicana, 60(1), 83-99. http://dx.doi.org/10.18268/BSGM2008v60n1a6

Cohen, A.S., 2003, Paleolimnology: The history and evolution of lake systems: New York, Oxford University Press,500 p.

Colinvaux, P., de Olivera, P.E., Moreno, P.J.E., 1999, Amazon Pollen Manual and Atlas: Amsterdam, Harwood Academic Publishers, 322p. https://doi.org/10.1201/9781482283600

Comisión Nacional de Áreas Naturales Protegidas (CONANP), 2003, Programa de Manejo Reserva de la Biósfera Barranca de Metztitlán: México, D.F., Dirección General de Manejo para la Conservación, 11 de marzo de 2003, 204 p.

Correa-Metrio, A., Vélez, M.I., Escobar, J., St-Jacques, J.-M., López-Pérez, M., Curtis, J., Cosford, J., 2016, Mid-elevation ecosystems of Panama: future uncertainties in light of past global climatic variability: Journal of Quaternary Science, 31(7), 731-740. https://doi.org/10.1002/jqs.2899

Cuna, E., Zawisza, E., Caballero, M., Ruiz-Fernández, A.C., Lozano-García, S., Alcocer, J., 2014, Environmental impacts of Little Ice Age cooling in central Mexico recorded in the sediments of a tropical alpine lake: Journal of Paleolimnology, 51, 1-14. https://doi.org/10.1007/s10933-013-9748-0

Delcourt, H.R., Delcourt, P.A., 1991, Quaternary Ecology: A paleoecological perspective: Cornwall, New York, Chapman & Hall, 242 p.

Douglas, M.W., Maddox, R.A., Howard, K., Reyes, S., 1993, The mexican monsoon: Journal of Climate, 6, 1665-1677. https://doi.org/10.1175/1520-0442(1993)006%3C1665:TMM%3E2.0.CO;2

Elías-Gutiérrez, M., Suárez-Morales, E., Gutiérrez-Aguirre, M., Silva-Briano, M., Granados-Ramírez, J.G., Garfias-Espejo, T., 2008, Cladocera y Copepoda de las aguas continentales de México. Guia Ilustrada: Ciudad de Mexico, México, Universidad Nacional Autónoma de México, 322 p.

Escobar, J., Serna, Y., Hoyos, N., Velez, M.I., Correa-Metrio, A., 2020, Why we need more paleolimnology studies in the tropics: Journal of Paleolimnology, 64, 47-53. https://doi.org/10.1007/s10933-020-00120-6

Esperon-Rodriguez, M., Beaumont, L.J., Lenoir, J., Baumgartner, J.B., McGowan, J., Correa-Metrio, A., Camac, J.S., 2019, Climate change threatens the most biodiverse regions of Mexico: Biological Conservation, 240, 108215. https://doi.org/10.1016/j.biocon.2019.108215

Franco‐Gaviria, F., Correa‐Metrio, A., Núñez‐Useche, F., Zawisza, E., Caballero, M., Prado, B., Wojewódka, M., Olivares, G., 2020, Millennial‐to‐centennial scale lake system development in the mountains of tropical Mexico: Boreas, 49(2), 363-374. https://doi.org/10.1111/bor.12430

Frey, D., 1986, Cladocera Analysis, in Berglund, B. (ed.), Handbook of Holocene Palaeoecology and Palaeohydrology: Chichester, UK, John Wiley & Sons, 677-692.

García Arizaga, M.T., Lugo Hubp, J., Palacios, D., 1996, La obturación de valles por procesos de ladera: el origen de la Vega de Metztitlán (México), en IV Reunión de Geomorfología: España, Sociedad Española de Geomorfología, 325-335.

Haug, G.H., Hughen, K.A., Sigman, D.M., Peterson, L.C., Röhl, U., 2001, Southward migration of the Intertropical Convergence Zone through the Holocene: Science, 293, 1304-1308. https://doi.org/10.1126/science.1059725

Hodell, D.A., Brenner, M., Curtis, J.H., Guilderson, T., 2001, Solar forcing of drought frequency in the Maya lowlands: Science, 292, 1367-1370. https://doi.org/10.1126/science.1057759

Hofmann, W., 1987, Cladocera in space and time: analysis of lake sediments: Hydrobiologia, 145, 315-321. https://doi.org/10.1007/BF02530293

Jenkin, P. M., 1934, Report on the Persy Sladen Expedition to some Rift Valley Lakes in Kenya in 1929. VI. Cladocera from the Rift Valley Lakes in Kenya: Annals and Magazine of Natural History, 13, 281-308.

Jolliffe, I.T., 1986, Principal Component Analysis: Springer-Verlag, 271 p. https://doi.org/10.1007/b98835

Korhola, A., 1999, Distribution patterns of Cladocera in subarctic Fennoscandian lakes and their potential in environmental reconstruction: Ecography, 22, 357-373. https://doi.org/10.1111/j.1600-0587.1999.tb00573.x

Korhola, A., Rautio, M., 2001, Cladocera and other branchiopod crustaceans, in Smol, J.P., Birks, J.B., Last, W.M. (eds.), Tracking environmental change using lake sediments. Vol 4: Zoological indicators: Dordrecht, Kluwer Academic Publishers, 5-41. https://doi.org/10.1007/0-306-47671-1_2

Kotov, A.A., Elías-Gutiérrez, M., 2004, Notes on Aloninae Dybowski & Grochowski, 1894 emend. Frey, 1967 (Cladocera: Anomopoda: Chydoridae): 2. Leydigia cf. striata Biraben, 1939 in South Mexico: Arthropoda Selecta, 13, 1-6.

Kotov, A.A., Elías-Gutiérrez, M., Guadalupe Nieto, M., 2003, Leydigia louisi louisi Jenkin, 1934 in the Neotropics, L. louisi mexicana n. subsp. in the Central Mexican highlands: Hydrobiologia, 510, 239-255. https://doi.org/10.1023/B:HYDR.0000008645.71534.81

Kylander, M.E., Ampel, L., Wohlfarth, B., Veres, D., 2011, High‐resolution X‐ray fluorescence core scanning analysis of Les Echets (France) sedimentary sequence: new insights from chemical proxies: Journal of Quaternary Science, 26(1), 109-117. https://doi.org/10.1002/jqs.1438

Löwemark, L., Chen, H.-F., Yang, T.-N., Kylander, M., Yu, E.-F., Hsu, Y.-W., Lee, T.-Q., Song, S.-R., Jarvis, S., 2011, Normalizing XRF-scanner data: a cautionary note on the interpretation of high-resolution records from organic-rich lakes: Journal of Asian Earth Sciences, 40, 1250-1256. https://doi.org/10.1016/j.jseaes.2010.06.002

Magaña, V.O., Vázquez, J.L., Pérez, J.L., Pérez, J.B., 2003, Impact of El Niño on precipitation in Mexico: Geofísica Internacional, 42(3), 313-330.

Mann, M.E., Zhang, Z., Rutherford, S., Bradley, R.S., Hughes, M.K., Shindell, D., Ammann, C., Faluvegi, G., Ni, F., 2009, Global signatures and dynamical origins of the Little Ice Age and Medieval Climate Anomaly: Science, 326, 1256-1260. https://doi.org/10.1126/science.1177303

Mayewski, P.A., Rohling, E.E., Stager, J.C., Karlén, W., Maasch, K.A., Meeker, L.D., Meyerson, E.A., Gasse, F., van Kreveld, S., Holmgren, K., Lee-Thorp, J., Rosqvist, G., Rack, F., Staubwasser, M., Schneider, R.R., Steig, E.J., 2004, Holocene climate variability: Quaternary Research, 62, 243-255. https://doi.org/10.1016/j.yqres.2004.07.001

Metcalfe, S.E., Jones, M.D., Davies, S.J., Noren, A., MacKenzie, A., 2010, Climate variability over the last two millennia in the North American Monsoon region, recorded in laminated lake sediments from Laguna de Juanacatlán, Mexico: The Holocene, 20(8), 1195-1206. https://doi.org/10.1177/0959683610371994

Moy, C.M., Seltzer, G.O., Rodbell, D.T., Anderson, D.M., 2002, Variability of El Niño/Southern Oscillation activity at millenial timescales during the Holocene epoch: Nature, 420: 162-165. https://doi.org/10.1038/nature01194

Oksanen, J., Blanchet, F.G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., Minchin, P.R., O’Hara, R.B., Simpson, G.L., Solymos, P., Stevens, M.H.H., Szoecs, E., Wagner, H., 2019, Vegan: Community Ecology Package, R package version 2.5-6: The R Project for Statistical Computing, available at <http://CRAN.R-project.org/package=vegan>, consulted on July 15, 2020.

Park, J., Byrne, R., Böhnel, H., 2019, Late holocene climate change in Central Mexico and the decline of Teotihuacan: Annals of the American Association of Geographers, 109(1), 104-120. https://doi.org/10.1080/24694452.2018.1488577

R Core Team, 2020, R: A Language and Environment for Statistical Computing: Vienna, Austria, R Foundation for Statistical Computing, available at <https://www.r-project.org/>, consulted on July 15, 2020.

Reimer, P.J., Austin, W.E.N., Bard, E., Bayliss, A., Blackwell, P.G., Ramsey, C.B., Butzin, M., Chen, H., Edwards, R.L., Friedrich, M., Grootes, P.M., Guilderson, T.P., Hajdas, I., Heaton, T.J., Hogg, A.G., Hughen, K.A., Kromer, B., Manning, S.W., Muscheler, R., Palmer, J.G., Pearson, C., van der Plicht, J., Reimer, R.W., Richards, D.A., Scott, E.M., Southon, J.R., Turney, C.S.M., Wacker, L., Adolphi, F., Büntgen, U., Capano, M., Fahrni, S.M., Fogtmann-Schulz, A., Friedrich, R., Köhler, P., Kudsk, S., Miyake, F., Olsen, J., Reinig, F., Sakamoto, M., Sookdeo, A., Talamo, S., 2020, The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP): Radiocarbon, 62(4), 725-757. https://doi.org/10.1017/RDC.2020.41

Richey, J.N., Poore, R.Z., Flower, B.P., Quinn, T.M., 2007, 1400 yr multiproxy record of climate variability from the northern Gulf of Mexico: Geology, 35(5), 423-426. https://doi.org/10.1130/G23507A.1

Rodríguez-Ramírez, A., Caballero, M., Roy, P., Ortega, B., Vázquez-Castro, G., Lozano-García, S., 2015, Climatic variability and human impact during the last 2000 years in western Mesoamerica: evidences of late Classic and Little Ice Age drought events: Climate of the Past, 11, 1239-1248. https://doi.org/10.5194/cp-11-1239-2015

Rose, C.W., 1985, Developments in soil erosion and deposition models, in Stewart, B.A. (ed.), Advances in soil science: New York, Springer, 1-63. https://doi.org/10.1007/978-1-4612-5088-3_1

Rothwell, R.G., Rack, F.R., 2006, New techniques in sediment core analysis: an introduction: Geological Society, London, Special Publications, 267, 1-29. https://doi.org/10.1144/GSL.SP.2006.267.01.01

Rzedowski, J., 2006, Vegetación de México: México, Comisión Nacional para el Conocimiento y Uso de la Biodiversidad, 504 p.

Sageman, B.B., Lyons, T., 2003, Geochemistry of fine-grained sediments and sedimentary rocks: Treatise on Geochemistry, 7, 115-158. https://doi.org/10.1016/B0-08-043751-6/07157-7

Sand-Jensen, K., Borum, J., 1991, Interactions among phytoplankton, periphyton, and macrophytes in temperate freshwaters and estuaries: Aquatic Botany, 41, 137-175. https://doi.org/10.1016/0304-3770(91)90042-4

Servicio Meteorológico Nacional, 2020, Normales climatológicas, Hidalgo, Estación Metztitlán (en linea): México, Comision Nacional del Agua, available at <https://smn.conagua.gob.mx/es/informacion-climatologica-por-estado?estado=hgo>, May 10, 2020.

Sigala, I., Caballero, M., Correa-Metrio, A., Lozano-García, S., Vázquez, G., Pérez, L., Zawisza, E., 2017, Basic limnology of 30 continental waterbodies of the Transmexican Volcanic Belt across climatic and environmental gradients: Boletín de la Sociedad Geológica Mexicana, 69(2), 313-370. http://dx.doi.org/10.18268/BSGM2017v69n2a3

Suter, M., 2004, A neotectonic-geomorphologic investigation of the prehistoric rock avalanche damming Laguna de Metztitlán (Hidalgo State, east-central Mexico): Revista Mexicana de Ciencias Geológicas, 21(3), 397-411.

Szeroczyńska, K., 2002, Human impact on lakes recorded in the remains of Cladocera (Crustacea): Quaternary International, 95-96, 165-174. https://doi.org/10.1016/S1040-6182(02)00037-X

Szeroczyńska, K., Sarmaja-Korjonen, K., 2007, Atlas of subfossil Cladocera from central and northern Europe: Poland, Friends of the Lower Vistula Society, 83 p.

Tribovillard, N., Algeo, T.J., Lyons, T., Riboulleau, A., 2006, Trace metals as paleoredox and paleoproductivity proxies: an update: Chemical Geology, 232, 12-32. https://doi.org/10.1016/j.chemgeo.2006.02.012

Tribovillard, N., Algeo, T. J., Lyons, T., & Riboulleau, A.,2006, Trace metals as paleoredox and paleoproductivity proxies: an update: Chemical Geology, 232(1-2), 12-32. https://doi.org/10.1016/j.chemgeo.2006.02.012

Vázquez-Molina, Y., Correa-Metrio, A., Zawisza, E., Franco-Gaviria, J.F., Pérez, L., Romero, F., Prado, B., Charqueño-Célis, F., Esperón-Rodríguez, M., 2016, Decoupled lake history and regional moisture availability in the middle elevations of tropical Mexico: Revista Mexicana de Ciencias Geológicas, 33(3), 355-364.

Waitz, P., 1947, Dos grandes derrumbes que causaron la formación de lagos, uno moderno en el Perú y otro antiguo en el Estado de Hidalgo: Ingeniería Hidráulica en México, 1, 145-160.

Wetzel, R.G., 2001, Limnology: lake and river ecosystems: San Diego, USA, Elsevier Academic Press, 1006 p.

Wetzel, R.G., Likens, G.E., 2000, Limnological analysis: Springer-Verlag New York, 429 p. https://doi.org/10.1007/978-1-4757-3250-4

Whiteside, M.C., Swindoll, M.R., 1988, Guidelines and limitations to cladoceran paleoecological interpretations: Palaeogeography, Palaeoclimatology, Palaeoecology, 62, 405-412. https://doi.org/10.1016/0031-0182(88)90065-X

Wogau, K.H., Arz, H.W., Böhnel, H.N., Nowaczyk, N.R.,Park, J., 2019, High resolution paleoclimate and paleoenvironmental reconstruction in the Northern Mesoamerican Frontier for Prehistory to Historical times: Quaternary Science Reviews, 226, 106001. https://doi.org/10.1016/j.quascirev.2019.106001

Wojewódka, M., Zawisza, E., Cohuo, S., Macario-Gonzáles, L., Schwalb, A., Zawiska, I., Pérez, L., 2016, Ecology of Cladocera species from Central America based on subfossil assemblages: Advances in Oceanography and Limnology, 7(2), 145-156. https://doi.org/10.4081/aiol.2016.6266

Wojewódka, M., Sinev, A.Y., Zawisza, E., 2020a, A guide to the identification of subfossil non-chydorid Cladocera (Crustacea: Branchiopoda) from lake sediments of Central America and the Yucatan Peninsula, Mexico: part I: Journal of Paleolimnology, 63, 269-282. https://doi.org/10.1007/s10933-020-00115-3

Wojewódka, M., Sinev, A.Y., Zawisza, E., Stańczak, J., 2020b, A guide to the identification of subfossil chydorid Cladocera (Crustacea: Branchiopoda) from lake sediments of Central America and the Yucatan Peninsula, Mexico: part II: Journal of Paleolimnology, 63, 37-64. https://doi.org/10.1007/s10933-019-00102-3

Manuscript received: March 16, 2020

Corrected manuscript received: August 25, 2020

Manuscript accepted: September 14, 2020

Peer Reviewing under the responsibility of Universidad Nacional Autónoma de México.

This is an open access article under the CC BY-NC-SA license(https://creativecommons.org/licenses/by-nc-sa/4.0/)