1. Introduction
Southern South America is a key region for the understanding of Southern Hemisphere past environmental variability: it is the only land mass present at this latitude, containing continuous sedimentary records suitable for paleoenvironmental reconstructions. However, there is a limited number of suitable resolved paleorecords (Bonisegna et al., 2009), particularly those covering the last two millennia. This period is fundamental to place modern trends of observed biological proxy patterns, climate variables (mainly temperature and precipitation) and human impact records into a longer-term context to assess earth system sensitivity to natural and anthropogenic pressure. In the last few years, there has been an ongoing debate about placing the start of human activity in geological sequences. It has been proposed that a new epoch, the Anthropocene, has begun and efforts to formalize it are being carried out. It would be a new interval of geological time in which human influence on Earth and its geological record dominate over natural processes (Zalasiewicz et al., 2016). Comprehensive high-quality data sets and a sound understanding of past dynamics at different temporal scales are critical to reducing the uncertainty about the magnitude and impact of current and future climate variability and human activity. Tree rings and glacial records from the southern Patagonian Andes (Stine, 1994; Villalba, 1994; Masiokas et al., 2009a) indicate the occurrence of century-scale climatic fluctuations that can be associated with well-known climatic events in the Northern Hemisphere: the Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA). In the Northern Hemisphere, mean annual temperatures were warmer for the period from AD 800 to 1400 compared to the periods from AD 200 to 800 and from AD 1400 to 1900 (Jones and Mann, 2004). However, such climate events have not been well identified in southern South America (Moy et al., 2009). One reason is the sensitivity of proxies to detect the timing, magnitude, and nature of these events. Luckman and Villalba (2001) found a broad synchronicity in the initiation and timing of the LIA when studying glacier advances and treering records along the Andes. However, the study showed considerable variability among sites in both the relative magnitude and date of the events as well as the age of the LIA maximum position. Another possibility for the lack of MCA and LIA signals could be attributed to the paucity of multiproxy records analyzed for this time period in southern South America, which in turn impedes the proper characterization of these climate variations. Recent studies performed in extra-Andean lakes (Markgraf et al., 2003; Haberzettl et al., 2007; Fey et al., 2009; Kastner et al., 2010; Massaferro et al., 2013) revealed that multiproxy analyses of lacustrine sediments constitute a promising strategy for reconstructing past responses of these systems to climate and non-climate changes.
Bioproxy-based reconstructions from Patagonian lakes and peat bogs covering the last two millennia, focused on the dynamics of natural and anthropogenic forest disturbances (Szeicz et al., 1998, 2003; Huber and Markgraf, 2003). Only a few bioproxy-based reconstructions tackled climatic variability (Haberzettl et al., 2005; Mancini, 2007, 2009; Moy et al., 2008). Pollen records from Cabo Vírgenes collected from both a mire and a soil profile (Figure 1) suggested drier than present conditions before AD 800-1000, changing to more humid conditions during the last 700 years (Mancini, 2007). The multiproxy record of Laguna Potrok Aike (Figure 1) (Haberzettl et al., 2005) indicated arid conditions during AD 1250-1410 shifting to wetter and colder conditions until AD 1940. Since AD 1940, hydrological proxies point to persistent dry conditions, while for pollen spectra it is difficult to differentiate climatic variability from human impact (Haberzettl et al., 2005). These records are representative for continental lowlands at the southeastern tip of Patagonia, which are influenced by westerly winds as well as Atlantic and Antarctic circulation patterns. Toward the Andes Mountains and under a stronger influence of westerly winds, the isotopic record of Lago Guanaco (Figure 1) suggests high evaporation for the periods AD 1050-1350 and AD 1550-1900. In the same record, pollen trends between AD 1500 and 1800 were interpreted as indicative of higher than present precipitation (Moy et al., 2008; Moreno et al., 2009).
The aim of the present study is to reconstruct the environmental history of Laguna Las Vizcachas, located on the basaltic plateau Meseta de Las Vizcachas at 1100 masl in the Andean Patagonian steppe. High elevation lakes such as Laguna Las Vizcachas (Figure 1) are ideal for this type of study because signals of temperature fluctuations are often amplified by altitude (Fey et al., 2009). The environmental reconstruction is based on terrestrial (pollen) and aquatic (diatoms and chironomids) bioproxy data, complemented with sedimentological and geochemical records from the same sediment sequence (Fey et al., 2009). We aim to evaluate the extent of the sensitivity of these bioproxies to short-term climate variability and to human disturbances over the last millennia.
1.1. ENVIRONMENTAL SETTINGS
Laguna Las Vizcachas (50o42´ S, 71o59´ W) is located about 70 km east of the Andean mountain range, in the southern part of the Santa Cruz Province, Argentina (Figure 1). It is a glacial cirque lake situated about 1100 masl surrounded by the basaltic plateau Meseta de las Vizcachas, which reaches elevations of 1400 masl At present, the lake is 18.7 m deep and has a maximum length (north to south) of 1300 m and a maximum width of 500 m (west to east). It has an inflow in form of a waterfall and an outflow at the southern shore (Figure 1). A more detailed geological and geomorphological description of the catchment area was presented by Fey et al. (2009).
The climate is determined by the intensity of the southern westerly winds throughout the year. The Andean mountain range represents a barrier for humid air masses coming from the west Pacific, leading to a pronounced precipitation gradient from 3000 mm along the western slopes to 200 mm in the eastern lowlands. Meseta Las Vizcachas, close to the Andes and above 1000 masl, is located along this precipitation gradient. The annual rainfall at Laguna Las Vizcachas is ca. 300 mm, and the mean annual temperature is between 4 and 6.5°C. In July the average temperature is below 0°C, indicating precipitation in the form of snow is common during winter.
The dominant vegetation at this latitude is Festuca pallescens grass steppe which is present as isolated grass patches on highlands above 700 masl. This steppe is characterized by a high diversity of dwarf shrubs or cushions like Nassauvia darwinii, Mulinum microphyllum and Nardophyllum obtusifolium. Herbs are represented by Carex argentina, C. andina, Armeria elongata, Calceolaria lanceolada, Acaena pinnatifida, Relbunium richardianum, Polygala darwiniana, Colobanthus lycopodioides and Perezia recurvata (León et al., 1998, Oliva et al., 2001). At Laguna Las Vizcachas the presence of the cushion-plant Empetrum rubrum is meaningful. This species is common in cold highlands and slopes with long snow periods and at lower altitudes, on acidic soils with little nutrients and high availability of aluminum (Pisano, 1977; Roig et al., 1985).
2. Methods
2.1. SAMPLING, LITHOLOGY, CHRONOLOGY, AND PREVIOUS GEOCHEMICAL ANALYSIS
In March 2005, a sediment core (VIZ05/6) of 83 cm was recovered in the central part of the lake at 16 m depth with a modified ETH-gravity corer (Kelts et al., 1986) (Figure 1). The chronology, physical properties, mineralogy, stable isotopes, geochemistry, and general trends of diatoms of this core were studied by Fey et al. (2009). In the following paragraphs we briefly describe the methodology applied in this study.
The sediment sequence was divided into three lithological units: Unit A (83-70 cm), contains a plant macro-remain layer at the bottom. Unit B (70-25 cm) consists of many dark layers. Unit C (25-0 cm) presents several plant macro-remain layers at different levels (25-23 cm, 16.5-13 cm, 8.5-7.5 cm, 5.5-4.5 cm and 1.5-0.5 cm). Seven samples of fine fraction (<100 μm) of bulk sediment were dated at the Poznań Radiocarbon Laboratory, Poland. AMS 14C dates were calibrated with the southern hemisphere calibration curve (SHCal04, McCormac et al., 2004) using the software CALIB 5.0.2 (Stuiver and Reimer, 1993; Stuiver et al, 2005). Ages for sediments below the lowermost AMS 14C date were estimated by extrapolation of sedimentation rates. More details about the agedepth model were published in Fey et al. (2009). Volume specific magnetic susceptibility (κ) was measured with a Bartington F-sensor. The mineralogical composition was determined by standard powder X-ray diffraction (XRD) analyses (Philips X’Pert Pro MD equipped with an X’Celerator Detector Array). Concentrations of total organic carbon (TOC) and total nitrogen (TN), among others, were determined with a CNS elemental analyzer (EuroEA, Eurovector). Biogenic silica (BiSi) was analyzed by applying an alkaline digestion in autoclaves and subsequent detection by a continuous flow system with UV-VIS spectroscopy.
2.2. BIOPROXY ANALYSIS
A sediment sample of 4 cm3 was taken every centimeter along the core. Each sample was divided and approximately 1 cm3 of sediment was assigned for pollen, 1 cm3 for diatom analyses and 2 cm3 for chironomids. These proxies constitute the original data supplied for the current study.
Each pollen sample was submerged in KOH overnight to remove humic acids and then sieved through a 125 μm mesh. Only the fraction <125 μm was chemically treated. Exotic Lycopodium tablets of a known concentration were added in order to calculate pollen concentration (Stockmarr, 1971). Procedures included the treatment with 10% HCl for dissolution of carbonates, 40% HF for silica removal and acetolysis for cellulose removal (Faegri and Iversen, 1989). The residues were mounted on slides and counted under the microscope at 1000 × until a pollen sum of 300 per sample was reached. The analysis was performed every 4 cm of depth. Pollen percentage is based on the sum of all pollen types identified. Algal microfossils like Botryococcus and Pediastrum were also counted in the pollen slides. Reference pollen samples and the pollen floras of Heusser (1971) and Markgraf and D’Antoni (1978) were employed to facilitate pollen identification. Diatoms samples were treated following standard procedures, including heating with H2O2 (Battarbee, 1986) and mounting permanently onto microscope slides using Naphrax® mounting medium. A minimum of 400 valves were counted to calculate relative frequencies. Diatom taxonomy follows the standard literature (Krammer and Lange-Bertalot, 1986, 1988; Simonsen, 1987; Rumrich et al., 2000). A more detailed laboratory procedure is described by Fey et al. (2009). Chironomids samples of 2 g dry weight were prepared according to Massaferro and Brooks (2002). Sediment was deflocculated using 10% KOH and sieved through 100 and 200 μm meshes. Material retained by the meshes was hand-picked from a modified Bogorov counting tray under a dissecting microscope at magnifications of 25-40 ×. Individual heads were placed onto microscope slides and mounted using Hydro-Matrix® mounting medium. Identifications were made using a Nikon phase microscope at 400 × magnifications. Head capsules were identified with reference to available taxonomic literature (Cranston, 2000; Massaferro and Brooks, 2002) and the Patagonian subfossil chironomid taxonomic identification guide (Massaferro et al., 2013).
Bioproxies stratigrafic diagrams were elaborated with the C2 software version 1.6.6 (Juggins, 2010).
3. Results
Chronology was established from seven AMS 14C ages which gave a basal age of AD 410 (Fey et al., 2009). To facilitate temporal correlation and discussion of changes in the bioproxy records, stratigraphic zones were established following lithological units A (AD 400-750), B (AD 750- 1750) and C (AD 1750-2005), described in Fey et al. (2009). Owing to conspicuous changes in the bioproxies at ca. AD 1550, unit B was split into two subzones, B1 (AD 750-1550) and B2 (AD 1550-1750). In addition, geochemical data (TOC, TOC/TN, BiSi, Si/Ti) published (Fey et al.,2009) were combined with our data to contrast our interpretations.
3.1. POLLEN RECORD
A total of forty pollen types were identified in the sediment core. Major trends are shown in a simplified pollen diagram (Figure 2). The terrestrial vegetation around the lake is represented by the dominance of Poaceae with low values of Acaena, Asteraceae subf. Asteroideae and Empetrum rubrum. Before AD 1550 (units A and B1) the vegetation is dominated by a grass steppe represented by > 40% of Poaceae values in more than half of the total samples and ca. 3-5% of the other pollen types mentioned above. In zones B2 and C Poaceae declines fluctuating between 35 and 40%. Empetrum rubrum is relatively unchanged, although slight increases were recorded at the beginning of zones B2 and C, and in present times. After AD 1750 (zone C) and until the top of the record, Solanaceae was continuously recorded and Senecio, Rumex, Plantago and Ranunculus-type appeared. The extra-regional pollen of the Nothofagus dombeyi-type is present all along the sequence with percentages of ca. 30%.
3.2. DIATOM RECORD
A summary diagram (Figure 3) shows the diatom record dominated by small benthic forms oscillating between ca. 40% to 60% up to AD 1750 (zones A, B1 and B2) and decreasing to 20-40% in zone C. Subdominant planktonic taxa, Aulacoseira distans and Discostella stelligera seem to follow comparable patterns of change, showing an increasing trend in zone A followed by a decrease in zone B1 up to AD 1500, where they start to increase again reaching the highest values (ca. 20% D. stelligera and ca. 40% A. distans) at ca. AD 1700/1750. Both species show a gradual decrease from AD 1700/1750 up to the present. Fluctuating low values of epiphytic taxa Cocconeis placentula, Gomphonema gracile, G. sp. aff. subclavatum, and benthic Karayevia clevei are recorded from the bottom of the core up to AD 1750. At that time, Karayevia clevei increases. After AD 1750, there is an increase of Gomphonema sp. aff. subclavatum, G. acuminatum, and Epithemia adnata together with the appearance of benthic taxa such as Fragilaria capucina, F. crotonensis, Encyonema minutum, E. silesiacum, and Achnanthidium minutissimum that remain high until present days.
3.3. CHIRONOMID RECORD
Chironomids are only present in zones B1, B2 and C (Figure 4) with low abundances values (between 10 and 30 head capsules achieving the maximum values in the last 300 years). From the 6 morphotaxa recorded, Tanytarsini 1B (50% of the total number) is dominant, followed by Paratrissocladius accuminatus (25% of the total number) and Cricotopus spp. (8% of the total number). The remaining taxa (Smittia, Paratrichocladius, and Tanytarsini 1A) are present with abundances of < 5%.
4. Discussion
4.1. ENVIRONMENTAL RECONSTRUCTION
4.1.1. PRIOR TO AD 1550
The almost exclusive dominance of the small benthic diatoms from AD 400 up to 900 suggest long periods of ice cover and persistent cold air temperatures compared to the present (Figure 3). Indeed, small benthic diatoms succeed when the ice-free growing season is short (Hobbs et al., 2010). The present winter temperatures in the catchment area of Laguna Las Vizcachas are below 0°C. Thus it is highly probable that the lake surface was also covered by ice during past cold seasons. Likewise, the low percentages of planktonic diatoms Aulacoseira distans and Discostella stelligera may indicate low stability and mixing of the water-column probably related to strong winds (Michelutti et al., 2015). In Sterken et al. (2008) low abundance of these taxa was interpreted as a depletion of oxygen in the aquatic environment. Indeed, periods of prolonged ice cover on the lake and snow cover in the catchment area may lead to a reduction of the organic matter oxidation due to lake water stratification under ice cover and prolonged anoxic conditions at the sediment/water interface (Fey et al., 2009).
Between AD 800 and 1000 the increase of Aulacoseira distans coupled with the decrease of small benthic diatoms (Figure 3) suggests temperature amelioration. Poaceae values exceeding 40% from ca. AD 900 indicate improved grass development, hence wetter than present conditions until AD 1550. Ecological studies performed on Festuca pallescens (Poacea) steppes showed that the water availability in the soil, mainly in summer, is a critical factor controlling the establishment and development of these grasses (Bertiller et al., 1996). Therefore, vegetation trends seem to indicate that precipitation began to increase gradually. Between AD 1200 and 1400, high precipitation is suggested by the presence of epiphytic diatoms Cocconeis placentula, Gomphonema spp. and Karayevia clevei which are commonly related to fluvial water (Sterken et al., 2008). Titanium (Ti), calcium (Ca) and magnetic susceptibility values recorded by that time support our interpretations. They are linked to enhanced fluvial clastic input due to increasing precipitation according to Fey et al. (2009). Summarizing, between AD 400 and AD 1550 cold and mostly wet conditions prevailed in Las Vizcachas environment, although a warm spell may have intercepted this cold trend at ca. AD 800 to 1000. The most likely mechanism causing this surplus in precipitation can be associated with intensification of the westerly wind circulation (Fey et al., 2009).
Records from different natural archives, like relict tree and shrub stumps rooted in current lakes and marshes of Lago Cardiel (48°57´ S, 71°26´ W) and Lago Argentino (Figure 1; 50°28 ́S, 72°58 ́W; Stine, 1994) evidence a period of persistent warm and dry conditions around AD 1200 related to the MCA in lowlands located at similar latitudes (Stine, 1994; Haberzettl et al., 2005; Moy et al., 2009). Oxygen isotope-based reconstructions from Lago Guanaco in Chile (Figure 1; 51°52´ S, 72°52´ W) suggested high evaporation and then aridity, spanning the period AD 1050-1350 (Moy et al., 2009). Further to the east, the multiproxy record from Laguna Potrok Aike (Figure 1) displayed fluctuations indicating also arid conditions during the MCA, at ca. 1250-1410 (Haberzettl et al., 2005). Considering these evidences, Fey et al. (2009) proposed that intensified westerly winds might produce higher precipitation at Laguna Las Vizcachas but dry conditions at lower elevations, especially to the east. Our bioproxy records from Laguna Las Vizcachas support this hypothesis of wet conditions at western high altitudes during this time period.
4.1.2. FROM AD 1550 TO AD 1750
A decline of Poaceae pollen and slight increase of Empetrum rubrum values combined with the highest values of small benthic diatoms shortly before AD 1600 (figs. 2 and 3) point to colder conditions compared to the time period described above, probably related to the LIA. Colder temperatures with precipitation stored as snow for longer periods of time may have limited the water availability and restricted the development of grass, as shown by lower Poaceae percentages. In this scenario, Empetrum rubrum growth would have been favored. The synchronous appearance of cold chironomids from the Tribe Tanytarsini (Massaferro and Brooks, 2002) in the aquatic environment reinforced this assumption. Changes in pollen, chironomids, and diatoms of Laguna Las Vizcachas during the time interval AD 1500-1600, seem to predate the beginning of the LIA, which was established by the time of glacier advances reconstructed for the Southern Patagonian Ice Field ( AD 1600) (Luckman and Villalba, 2001; Masiokas et al., 2009a, b). However, the multiproxy reconstruction based on stable isotopes and pollen from Lago Guanaco in Chile (51° 52 ́S, 72°52 ́W) characterized the onset of the LIA period by an intensification of westerly winds and enhanced precipitation across the Andes at AD 1500 (Moy et al., 2008). This lake is located at around the same latitude as Laguna Las Vizcachas but on the other side of the Andes in the same climatic transition zone, where steep environmental gradients play an important role in defining the glacial patterns (Luckman and Villalba, 2001).
4.1.3. POST AD 1750, THE ANTHROPOCENE
The most noticeable change in all the proxies of Laguna Las Vizcachas took place in the transition between zone B2 and C, at AD 1750. This is indicated by the chironomid assemblage composed by semi-terrestrial and littoral morphotaxa (Smittia, Parapsectro cladius, Cricotopus) mixed with profundal morphotaxa such as Tanytarsini. Together with this, the appearance of Ranunculaceae pollen and the explosion of mainly benthic and epiphytic diatom taxa and the increase of macrophyte remains seem to indicate an expansion of littoral zones and, as a consequence, new habitats for the biota. This change is probably associated with the lowering of lake water levels reflecting a decline in water supply. Low magnetic susceptibility values (k) are consistent with this trend (Fig. 5). During this time, the maximum richness and abundance values of chironomids and diatoms are probably the results of a longer ice-free season suggesting better conditions for aquatic biota to colonize and grow. At AD 1800, shorter periods of ice cover and better conditions are also indicated by the increase of Poaceae and the decline of small benthic diatoms (Fig. 3). In addition, the highest values of Dicostella and Aulacoseira recorded at the beginning of zone C also suggest warmer conditions than before. Michelutti et al. (2015) attributed the rise of Dicostella and Aulacoseira to recent climate warming, a shorter duration of seasonal ice cover, and enhanced stability of the water-column thermal stratification during the growing season. The authors interpreted the increase in both planktonic taxa as ecological evidence for the onset of the Anthropocene, detectable in many lake sediments of the northern hemisphere. In our record, the presence of these taxa reflects a combined signal of global warming, reduced wind intensity, and lake level changes. Indeed, the appearance of epiphytic diatoms during this period (Figure 3) was ascribed to both, increased wind intensities and better conditions due to macrophyte growth (Fey et al., 2009).
4.2. LAKE PRODUCTIVITY
Before AD 1200, TOC/TN values indicate a mixture of aquatic and terrestrial plant material (Figure 5). This ratio is used as a proxy for the contribution of terrestrial plants into the sediments. Values greater than 20 are associated with allochthonous inputs. Our bioproxy records also show wet conditions by that time. Hence, it is plausible that a combined effect of allochthonous input with a consequent dilution of nutrients by the minerogenic input and/or changes in the organic matter preservation could have taken place, explaining to a certain extent the absence of chironomids by that time. From AD 1200 up to the present, increasing TOC and Si/Ti values give evidence of increasing aquatic productivity (Figure 5). TOC reflects the balance between organic matter production and decomposition, whereas Si/Ti ratio (Figure 5) is an indicator of in situ diatom and chrysophyte productivity as explained in Fey et al. (2009). This trend continues until the present; however, after ca. AD 1750, the establishment of the different biotic communities combined with the increase of TOC and decrease of TOC/TN points towards an increase of Laguna Las Vizcachas autochthonous productivity. These changes agree with the chronology of recent global warming and with the increase of human activities in the area that certainly have contributed to this intensification in the lake trophy.
4.3. EARLY EUROPEAN SIGNALS AND SUBSEQUENT PERMANENT SETTLEMENTS
Human activity near Laguna Las Vizcachas is evident from the appearance of Rumex in the pollen record occurring at ca. AD 1600, indicating an early presence of European settlers in the area. A first appearance of Rumex dated to ca. AD 1620 was also found in the pollen record of Rio Rubens, a bog located south of Laguna Las Vizcachas (Figure 1; 52°08´15”S-71°52´53”W; Huber and Markgraf, 2003). Given that Rumex is widely recognized as the stratigraphic marker of the first European settlements in southern Patagonia at the end of the 19th century (Heusser, 2003), the earlier appearance in two different pollen records is remarkable. A continuous presence of exotic Plantago is observed in the pollen record slightly before AD 1900. Some species of Plantago are considered indicators of the presence of Europeans (Heusser, 2003). Coupled with the Rumex record, Plantago pollen type along with the appearance of epiphytic taxa like Achnanthidium minutissimum and Epithemia adnata, in the diatom record, may indicate a permanent settlement of Europeans. Moreover, it can be related to the beginning of rural activities around the lake before AD 1900 and onwards, given that Achnanthidium minutissimum and Epithemia adnata are both associated with mild to moderate eutrophication levels of lakes (Hobbs et al., 2010).
5. Conclusion
The combined analysis of pollen, diatoms and chironomids, indicates climatic variations during the last ~1500 years at Laguna Las Vizcachas. The results also indicate that changes in the biotic communities after AD 1750 are response to human presence in the area. Assuming that water deficit is the most prominent climatic characteristic in Patagonia, Poaceae fluctuations are considered a proxy for air humidity, whereas fluctuations in the abundance of small benthic diatoms were interpreted as a proxy for the duration of lake ice-cover. Since AD 700, pollen and diatom taxa coupled with magnetic susceptibility values suggest a period of high precipitation lasting more than 500 years. During this time, temperatures were cold enough to maintain lake surface ice during cold seasons. Such a trend is opposite to reconstructions from Chilean lowland records at similar latitudes, with mimic the MCA dry period (Stine, 1994; Haberzettl et al., 2005; Moy et al., 2009). Pollen, diatom and chironomid records support the hypothesis proposed by Fey et al. (2009) of inverse precipitation regimes between western and eastern parts of extra-Andean southern Patagonia. The coldest conditions are recorded at ca. AD 1550 by the increase of Empetrum rubrum and small benthic diatoms. This cold signal registered in Laguna Las Vizcachas precedes the beginning of the LIA, recorded by glacier fluctuations in the Southern Patagonian Ice Field (Masiokas et al., 2009a, b). Before AD 1700, the climate was cold and wet. Therefore, it is likely that the lake experienced prolonged periods of ice cover enabling lake water stratification and subsequent anoxic conditions at the sediment/water interface. Such a situation could affect aquatic biota by reducing their abundances and favoring the development of species tolerant to low oxygen concentration. The decreasing trend of small benthic diatoms after AD 1700 marks a change towards warmer conditions. An expansion of littoral habitats due to a decrease in precipitation, and hence a lowering of the lake level explains in part the pronounced increase of epiphytic and benthic diatoms, the appearance of chironomids, and the increase of autochthonous productivity. The appearance in the sequence of Rumex and Plantago pollen coeval with Achnanthidium minutissimum and Epithemia adnata diatoms suggest that both, the early human presence and later human activities in the area also contributed to modify terrestrial and aquatic biotic composition at the lake watershed.
Our paleoecological data give evidence that climate and human impacts started in AD 1750 modifying the lake physical structure and biotic composition. Global warming of the 20th century has certainly enhanced the present lake trophic status causing ecological changes in the aquatic communities. Likewise, our study highlights the potential of the multiproxy approach to provide more accurate environmental reconstructions, especially regarding the timing and onset of human activities in conjunction with changing climates and how these two major drivers have interacted to modify aquatic ecosystems in the past. Disentangling both, climate and human driven impacts on aquatic ecosystems, is a key challenge for understanding the onset and nature of the Anthropocene (Mills et al 2017).