quaternary Article Modern Analogue Approach Applied to High-Resolution Varved Sediments—A Synthesis for Lake Montcortès (Central Pyrenees) Teresa Vegas-Vilarrúbia 1, * , Valentí Rull 2 , Maria del Carmen Trapote 1,2 , Min Cao 3 , Antoni Rosell-Melé 3,4 , Teresa Buchaca 5 , Joan Gomà 1 , Pilar López 1 , Javier Sigró 6 , Elisabet Safont 1 , Núria Cañellas 2 , Sandra Garcés-Pastor 1 , Santiago Giralt 2 , Juan Pablo Corella 7 and Núria Pérez-Zanón 6 1 Faculty of Biology, Universitat de Barcelona, Av. Diagonal 643, 08028 Barcelona, Spain; mctrapote84@gmail.com (M.d.C.T.); jgoma@ub.edu (J.G.); marilopez@ub.edu (P.L.); esafont@ub.edu (E.S.); sgarcespastor@gmail.com (S.G.-P.) 2 Institute of Earth Sciences Jaume Almera (ICTJA-CSIC), C. Lluis Sole i Sabarís s/n, 08028 Barcelona, Spain; vrull@ictja.csic.es (V.R.); nuriacatcb@gmail.com (N.C.); sgiralt@ictja.csic.es (S.G.) 3 Institut de Ciència i Tecnologia Ambientals (ICTA-UAB), Universitat Autònoma de Barcelona, Bellaterra, 08193 Catalonia, Spain; smilecaomina@hotmail.com (M.C.); antoni.rosell@uab.cat (A.R.-M.) 4 Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, 08010 Catalonia, Spain 5 Centre for Advanced Studies of Blanes (CEAB-CSIC), Accés a la Cala St. Francesc 14, 17300 Blanes, Spain; buch@ceab.csic.es 6 Center for Climate Change (C3), Universitat Rovira i Virgili, C/ Joanot Martorell 15, Vila-seca, 43500 Tarragona, Spain; javier.sigro@urv.cat (J.S.); nuria.perez@bsc.es (N.P.-Z.) 7 Universite Grenoble Alpes, CNRS, IRD, Grenoble INP, IGE, 38000 Grenoble, France; pcorella@iqfe.csic.es * Correspondence: tvegas@ub.edu Academic Editors: Bernd Zolitschka and Matthew Peros   Received: 21 May 2019; Accepted: 19 December 2019; Published: 2 January 2020 Abstract: In Quaternary paleosciences, the rationale behind analogical inference presupposes that former processes can be explained by causes operating now, although their intensity and rates can vary through time. In this paper we synthesised the results of different modern analogue studies performed in a varved lake. We discuss their potential value to obtain best results from high resolution past records. Different biogeochemical contemporary processes revealed seasonality and year-to-year variability, e.g., calcite precipitation, lake oxygenation, production and deposition of pollen and phytoplankton growth. Fingerprints of the first two of these processes were clearly evidenced in the varve-sublayers and allow understanding related to past events. Pollen studies suggested the possibility of identifying and characterizing seasonal layers even in the absence of varves. Marker pigments in the water column were tightly associated with phytoplankton groups living today; most of them were identified in the sediment record as well. We observed that 50% of these marker pigments were destroyed between deposition and permanent burying. In another study, seasonality in the production/distribution of branched glycerol dialkyl glycerol tetraethers (brGDGTs) and derived temperature estimates were investigated in catchment soils and particles settling in the lake. The signatures of brGDGTs in depositional environments mainly were representative of stable conditions of soils in the catchment that last over decades; no brGDGTs seemed to be produced within the lake. The main contribution of this review is to show the advantages and limitations of a multiproxy modern-analogue approach in Lake Montcortès as a case study and proposing new working hypotheses for future research. Keywords: endogenic varves; calcite precipitation; pollen traps; meromixis; freshwater GDGTs; pigments; long-term ecology Quaternary 2020, 3, 1; doi:10.3390/quat3010001 www.mdpi.com/journal/quaternary Quaternary 2020, 3, 1 2 of 23 1. Introduction Analogy is a powerful form of reasoning that allows two domains to be considered as similar, based on their relationships [1–3]. The analogue approach permits the extension of knowledge of a base domain to a target past domain by virtue of their expected similarity and by means of analogical inference. This process involves reproducing structures and processes from the base to the target in which missing information is generated. Nonetheless, constraints must be placed on what information is to be carried over [4], in order to obtain inferences that are as realistic as possible. In Quaternary paleosciences, the rationale behind analogue inferences is based on postulating that past natural changes can be explained by causes operating today (principle of uniformitarianism), but that intensity and rates at which processes take place vary through time [5]. On the other hand, reconstruction and derived inferences for any proxy (natural recorders of variability) demand a thorough comprehension of the spatial and temporal frame in which processes and environmental conditions leading to the existence and persistence of each particular proxy make sense. Examples of modern analogue studies of proxies and their applications are widespread. In paleoecology, fossil assemblages (pollen, diatoms, pigments, ostracods etc.) are examined for coincidences with modern assemblages, and the direction and rate at which communities vary through space and time can be measured quantitatively by a wide array of numerical approaches and techniques [6–8]. For instance, in the Peruvian Andes, the pollen/altitude relationship in a forest transect was studied by means of correspondence analysis, and the good statistical correlation between altitude and the first component allowed building a transfer function with statistical reliability to infer altitude from pollen data [9]. In another case, a surface-sediment survey of pigments in 100 lakes (Northern Sweden) was combined with a reconstruction of sedimentary pigments from Lake Seukokjaure. Light regime and terrestrial supplies of 13 C-depleted carbon turned out to be important controls of primary production, whereas depth, elevation and conductivity discriminated between different phototrophic communities. The application of these results to the sediment record demonstrated the importance of terrestrial vegetation for the development of the lake, and of specific factors to control algal community composition, across the Holocene [10]. The modern analogue approach can be used to build models for climate change reconstruction. Central to these models is the development of transfer functions that, by definition, are “empirically derived equations for calculating quantitative estimates of past atmospheric or oceanic conditions from paleontological data” [11]. Transfer functions have been widely applied during the last decades. For example, in Patagonia, chironomids of surface sediment of 63 lakes were used to model the atmospheric annual average temperature. The good correlation between estimated and measured temperatures allowed applying the model to the chironomid record of Lake Potrok Aike. The reconstructed trends of temperature were consistent with other quantitative records of the Southern Hemisphere [12]. In the Pyrenees, present chrysophyte cysts distribution is influenced by altitude, and a transfer function was developed to assess the local “altitude anomaly” on a lake site throughout time. A warming trend was detected showing sub-millennial climatic variability from the early Holocene to 4 kiloyears (kyr) before present (BP) [13]. Sometimes, finding modern analogues is not straightforward. For example, a modern analogue study of the Otter lake helped understand why and how siderite (FeCO3 ) accumulated only in varved sediments older than 1200 calibrated years, even though limnological data did not evidence siderite precipitation. The authors combined modern data with siderite time series, abundances of major element in sediment and experimental results on siderite precipitation. They showed that the lake system alternated between ferruginous, meromictic phases with intense precipitation of siderite and holomictic phases with very low siderite precipitation [14]. Temporal resolution that can be achieved with non-laminated sediment is often constrained to more than decadal resolution. Reconstructions of environmental variability at high resolution scales using varved sediments have provided an unprecedented level of detail that needs to be properly extracted and interpreted [15]. This implies knowledge of fine-tuned mechanisms that need to be observed today at the same temporal resolution as they occur in past records [16]. Lakes with varved sediments are especially well suited for paleoecological and paleoclimatic research demanding a high Quaternary 2020, 3, 1 3 of 23 level of detail, as it is possible to downscale to annual and seasonal resolution [15,17]. The interpretative power of such high-resolution paleoenvironmental reconstructions depends on the availability of suitable modern analogues with the same temporal resolution [16]. With the advent of new techniques and proxies that allow acquiring an increased level of detail, modern analogue studies are gaining importance, a fact we want to highlight with this paper. We have taken advantage of the varved sediments of the Mediterranean Lake Montcortès and conducted several independent modern analogue studies. Our aim was to document the potential of the modern analogue approach not only by determining how similar modern and past communities are, but also by investigating different contemporary biological and physico-chemical processes that are expected to also have occurred in the past on annual and seasonal time scales. We also explored which factors and relationships between modern and subfossil proxies can limit or hinder the applicability of the modern analogue approach representing a potential source of error in analogue analysis. Such errors should be foreseen and minimised [6]. In this context, we examined how the different modern analogue approaches can help extract important cause-effect relationships at seasonal, annual and subdecadal time scale. We discussed, evaluated and synthesised the published results of four modern analogue studies carried out for Lake Montcortès: (1) Relationships between calcite precipitation, environmental variables and primary producers and their influences on the thickness and patterns of calcareous varve sublayers were examined as high-resolution paleoclimate indicators [18]. (2) Seasonal pollen sedimentation was compared with local meteorological variables to determine whether contemporary seasonal deposition is consistent with the two layers of the varve pattern archived in the sediments [19]. (3) Contemporary shifts in dissolved O2 , Fe, Mn and anaerobic phototrophic bacteria markers were assessed, in order to investigate whether mixing and oxygenation dynamics Lake Montcortès have been the same over the last five centuries [20]. (4) Distribution of branched glycerol dialkyl glycerol tetraethers (GDGTs) in settling particulate matter and catchment soils was analysed in order to determine whether branched glycerol dialkyl glycerol tetraethers (brGDGTs) are seasonally biased in environments where annual brGDGTs production may not be constant [21]. Additionally, we provided and discussed the first results from ongoing research that aimed at relating contemporary phytoplankton taxa in Lake Montcortès to their corresponding marker pigments, and then to subfossil pigment records at subdecadal time scales. 2. Environmental Settings Lake Montcortès (Figure 1) belongs to the southern flank of the Central Pyrenees (Catalonia, Spain) at 42◦ 190 N, 0◦ 590 E and 1027 m altitude (Figure 1). This karstic lake located on carbonated and evaporitic Mesozoic substrate is chiefly characterised by Triassic limestones, marls and evaporites, as well as Oligocene-carbonate conglomerates [22,23]. Climatic data collected at La Pobla de Segur meteorological station show that the annual average air temperature is 12.8 ◦ C (2.9 ◦ C in January to 23.2 ◦ C in July). Total annual precipitation is 669 mm (33.4 mm in February and 88.4 mm in May). Maximum and minimum annual air temperatures were 41 and −20 ◦ C, respectively. The vegetation of the catchment area is evergreen with deciduous oak forests (Quercus rotundifolia and Q. pubescens), conifer forests (Pinus nigra), pastures and crops. The vegetation belt at the shore is dominated by hygrophyte communities of Typha domingensis or Cladium mariscus, reed beds of Phragmites australis and communities of Carex riparia. Rush formations and grasslands grow on occasionally flooded soils [24]. The lake surface area is 0.14 km2 and its maximum water depth 32 m [25]. Only an ephemeral inlet exists, but water is drained by an outlet located at the northern shore. The lake’s sedimentary record shows thin, well preserved biogenic varves formed during the Late Holocene. The varve structure Quaternary 2020, 3, 1 4 of 23 alternates between sublayers of endogenic calcite and organic detritus with additional detrital layers and turbidites Quaternary 2020, 3, embedded x in the varve succession [22,26,27]. 4 of 23 Figure Figure1.1. (A) (A)Location Locationof ofLake LakeMontcortès Montcortèsininthe thePyrenean Pyreneanmountain mountainrange range(NE (NEIberian IberianPeninsula). Peninsula). (B) Bathymetric map with position of the cores and the traps. (B) Bathymetric map with position of the cores and the traps. 3.3.Materials Materialsand andMethods Methods Here,we Here, we provided provided aa summary summary ofof the the main main field field and and laboratory laboratory methods methods used. used. More More detailed detailed information can be found in corresponding publications and literature therein [18–20]. information can be found in corresponding publications and literature therein [18–20]. 3.1.Monitoring 3.1. Monitoringand andCoring Coring Thelimnological The limnological characteristics characteristics of of Lake Lake Montcortès Montcortèswere wereassessed assessedmonthly monthlyfromfromOctober October20132013to October 2015. Detailed vertical profiles of temperature (T), dissolved oxygen to October 2015. Detailed vertical profiles of temperature (T), dissolved oxygen (DO), electric (DO), electric conductivity (EC) and pH (EC) conductivity wereandobtained pH wereusing a [28,29]using obtained multiparametric datasonde Hydrolab a [28,29] multiparametric datasondeDS5.Hydrolab Water samples DS5. were obtained Water samples werewith obtained an adapted withsubmersible an adapted bilge pump at submersible three bilge pump depths in the at three epi-,inmeta- depths and the epi-, hypolimnion for alkalinity, cations, total phosphorus (TP), total nitrogen meta- and hypolimnion for alkalinity, cations, total phosphorus (TP), total nitrogen (TN), suspended(TN), suspended solids and molecular marker analyses (pigments and brGDGTs), and at two solids and molecular marker analyses (pigments and brGDGTs), and at two depths (epi- and depths (epi- and metalimnion) for phytoplanktonfor metalimnion) microscopic phytoplankton examination. microscopic Between 1.5 and 3Between examination. L of water 1.5was andfiltered through 3 L of water wasWhatman filtered GF/F glass-fibre filters for pigment analysis. These filters were protected through Whatman GF/F glass-fibre filters for pigment analysis. These filters were protected with aluminium foilwith and frozen to prevent aluminium foil anddegradation until laboratory frozen to prevent degradation analysis. until laboratory analysis. A set of sediment traps was deployed A set of sediment traps was deployed at 20 m at 20 m water waterdepth depthfrom from aa floating floating platform platformlocated locatedat at the deepest part of the lake, coinciding with the site where sediment the deepest part of the lake, coinciding with the site where sediment coring was performed. Totalcoring was performed. Total suspended solids, suspended solids,calcite calcitecrystals crystals [18], pollen [18], and spores pollen [19], marker and spores pigmentspigments [19], marker and settlingand particulate settling matter for brGDGT extractions [20] were collected in these traps, particulate matter for brGDGT extractions [20] were collected in these traps, with monthly with monthly or quarterly (three or months) frequency depending on the variable of interest. Additionally, replicated quarterly (three months) frequency depending on the variable of interest. Additionally, replicated topsoil samples from a depthsamples topsoil of 2–5 cm were from collected a depth monthly of 2–5 cm werefor brGDGT collected analysis monthlyatfor three sites within brGDGT analysis a 10 at m radius three in sites the catchment area. These samples were combined to reduce soil within a 10 m radius in the catchment area. These samples were combined to reduce soil heterogeneity. Two UWITEC gravity cores (MON12-3A-1G, 78 cm length, and MON12-2A-1G, 106 cm length) heterogeneity. wereTworetrieved UWITEC at the deepest gravity area(MON12-3A-1G, cores of the lake basin78 (30 cm mlength, depth) and at the same location of MON12-2A-1G, 106the cmsediment length) trap (Figure 1) [27]. The varve chronology was performed by double counting were retrieved at the deepest area of the lake basin (30 m depth) at the same location of the sediment in 14 overlapping thin sections. trap (FigureOnly 1% of 1) [27]. Thevarves varvewere interpolated chronology due to poorby was performed varve doublepreservation counting in [27]. Two independent 14 overlapping thin dating methods using AMS 14 C and 210 Pb radiometric dating techniques were previously achieved sections. Only 1% of varves were interpolated due to poor varve preservation [27]. Two independent in the Lake dating methodsMontcortès using AMS sedimentary 14C and 210 sequence [22,26,27]. Pb radiometric datingThe agreementwere techniques between these independent previously achieved in chronologies the supportssedimentary Lake Montcortès the robustness of the varve sequence chronology. [22,26,27]. The agreement between these independent chronologies supports the robustness of the varve chronology. Quaternary 2020, 3, 1 5 of 23 3.2. Laboratory Analyses The calcite saturation index (Ω) was estimated from carbonate and calcium concentrations. Ca2+ was measured with inductively coupled plasma atomic emission spectroscopy, and CO3 2− was calculated from alkalinity, pH and temperature with the CO2SYS_XLS v 2.1 program [28]. To confirm the presence of calcite crystals and to assess seasonal variations of their size and shape, trapped material from traps was analysed by scanning electron microscopy. Total suspended solids retained in the traps were analysed for total particulate carbon (TPC) and total particulate organic carbon (POC) with an elemental organic analyser, Thermo EA 1108. Diatom specimens of quarterly sediment traps were prepared, identified and counted using standard methods [29]. A minimum of 300 valves were counted per sample; identification was carried out following Krammer and Lange-Bertalot [30,31]. Epilimnetic and metalimnetic phytoplankton was settled using sedimentation chambers; quantification was performed applying Utermöhl’s technique for counting combined with biovolume conversion as mm3 L−1 . Analysis of samples for pollen analysis followed standard methods [32] after addition of exotic Lycopodium spores. Pollen sum included all pollen types except the superabundant Pinus and Quercus. See the original palynological work for more detailed information on chemical processing, counting and plotting/statistical procedures [11]. As in earlier studies [33–35], identification of pollen was based on our own general reference collection for the Pyrenees, enriched after detailed floristic and vegetation studies around Lake Montcortès area [24]. Geochemical data to trace high resolution oxic/anoxic events of the hypoliminon were obtained on core MON12-3A-1K by X-ray fluorescence (XRF) analysis (2000 A, 10–30 kV and 20–50 s measuring time) (0.2 mm resolution) using an XRF core-scanner (Avaatech Analytical X-Ray Technology, Netherlands). We obtained sub-annual resolution time series for Fe, Mn, Ca, Si, S, Ti and Br. The former five elements were normalised to Ti. The Fe/Mn ratio was used as the main redox tracer, together with marker pigments of sulphur bacteria, such as okenone (oken) and isorenieratene (isore), which gave information at the subdecadal level. A second core (MONT-0713-G05) was used for pigment analyses; sediment samples were obtained at sampling intervals of 5 mm, reaching subdecadal resolution, the frozen sediment samples were freeze dried before analysis. Pigment analyses were also performed on epi-, meta- and hypolimnetic water samples and sediment trap material collected in Whatman GF/F glass fibre filters. Marker pigments were extracted from frozen filters or freeze-dried sediment in 5 ml of 90% acetone. The extract was centrifuged, filtered and analysed with ultra-high-performance-liquid-chromatography following the method of [36] with some modifications. Pigments were identified by evaluating their retention times and absorption spectra against a library based on algae and photosynthetic bacterial cultures [37]. To isolate GDGTs, soil samples and filters from sediment traps were extracted monthly, evaporated and then fractionated over an aminopropyl silica column. The fractions containing the lipids were dried, re-dissolved and filtered. brGDGTs were analysed by high-performance liquid chromatography coupled with mass spectrometers and detected in selected ion monitoring modes; prior to extraction, a synthetic tetraether lipid with a structure typical of neutral archaeal membrane lipids was added to the sample. Further analyses in soils were performed for water content and total organic carbon (TOC). 4. Results 4.1. Biogeochemical Processes Involved in Varve Formation: Where, When, and How? Processes driving varve formation are remarkably variable and site-specific. Consequently, detailed knowledge of the local processes that promote and modulate particle flux dynamics to the sediment for each site becomes essential. This knowledge is a pre-requisite for proper paleoenvironmental and paleoclimatic interpretations of the sedimentary archive [17,38] and for building reliable varve-based chronologies. An independent absolute varve chronology was obtained from the composite sedimentary sequence of both sediment cores extending from 2012 back to Common Era (CE) 1347 [27]. The varves Quaternary 2020, 3, 1 6 of 23 Quaternary 2020, 3, x 6 of 23 of Lake Montcortès are formed by couplets of white calcite and brownish organic layers. The white back to Common Era (CE) 1347 [27]. The varves of Lake Montcortès are formed by couplets of white layer was assumed to have deposited during spring/summer and the dark layer in fall/winter [17,27]. calcite and brownish organic layers. The white layer was assumed to have deposited during Biogenic varveand spring/summer formation the darkis layer highlyincomplex as it[17,27]. fall/winter adds physical Biogenicand chemical varve processes formation to the is highly complex complex biological processes involved [15]. as it adds physical and chemical processes to the complex biological processes involved [15]. With With this this modern modern analogue analogue study study several several goals goals areare pursued pursued forfor Lake Lake Montcortès: Montcortès: (i)(i) to to identify identify thethemain main limnological limnological and andsedimentological sedimentologicalmechanisms mechanisms leading to varves leading formation; to varves (ii) to establish formation; (ii) to establish in which season(s) the deposition of individual calcite sub-laminae occurs and (iii)iftothe in which season(s) the deposition of individual calcite sub-laminae occurs and (iii) to determine obtained if determine results are coherent the obtained resultswith thecoherent are varve structure with theand previously varve structure proposed mechanisms and previously of their proposed formation [26]. mechanisms of their formation [26]. Elevated Elevated ΩΩ values values wererecorded were recordedininthetheepi- epi-andandmetalimnion metalimnionLake Lake Montcortès Montcortès fromfrom late late spring spring to toearly earlyfall, fall,suggesting suggestingthatthatcalcite calciteprecipitation precipitation mostly mostly occurred occurred during during thisthis period. period.TheThe deposition deposition rates were highest in summer when photosynthetic organisms were most rates were highest in summer when photosynthetic organisms were most active (Figure 2A,B).active (Figure 2A,B). Additionally, higher POC higher Additionally, values were POC observed values were during this period, observed during indicating thatindicating this period, the main origin that theof calcite crystals main origin ofwere newly calcite formed crystals wereautochthonous newly formed CaCO 3 precipitates autochthonous CaCO (Figure 2C). 3 precipitates (Figure 2C). Figure 2. Summary of monthly water column and sediment trap monitoring. (A) Water column Figure 2. Summary of monthly water column and sediment trap monitoring. (A) Water column temperature at 1 m resolution. (B) Total mass fluxes (TMF) of monthly sediment traps compared with temperature at 1 m resolution. (B) Total mass fluxes (TMF) of monthly sediment traps compared the phytoplankton biovolume of epi- and metalimnion. (C) Calcite fluxes versus organic matter (OM) with the phytoplankton biovolume of epi- and metalimnion. (C) Calcite fluxes versus organic matter fluxes of sediment trap material. (D) Well-formed blocky and rhombohedral calcite crystals, calcite (OM) fluxes of sediment trap material. (D) Well-formed blocky and rhombohedral calcite crystals, crystals and diatom frustules and one of the small centric diatoms found in sediment trap material. calcite crystals and diatom frustules and one of the small centric diatoms found in sediment trap Redrawn and modified from [18]. material. Redrawn and modified from [18]. In fact, the presence of blocky, polyhedral calcite crystals was confirmed in all sediment trap In fact,and samples theshowed presencea of blocky, clear polyhedral seasonal pattern,calcite being crystals smaller was whenconfirmed inspring formed in all sediment trap and summer samples and showed a clear seasonal pattern, being smaller when formed in spring and summer (Figure 2D). Deposition of smaller crystals was not only concurrent with Ω, but also with larger Quaternary 2020, 3, 1 7 of 23 (Figure 2D). Deposition of smaller crystals was not only concurrent with Ω, but also with larger fluxes of very tiny diatoms (~10 µm), suggesting that phytoplankton cells could be acting as calcite condensation nuclei under enhanced calcite Ω, pH and productivity conditions. Other sources of calcite input, such as detrital carbonates from the catchment, were unimportant for the studied period. Other seasonal events are precipitation of calcite and OM (Figure 2C). Both variables displayed maximum fluxes in summer to fall. However, calcite precipitation shows the lowest values from winter to spring and was exceeded by OM deposition suggesting a match between the timing of contemporary precipitation of calcite and OM deposition and the alternating couplets of white calcite and brownish organic sublayers of the sediment which would correspond to summer/fall and winter/spring, respectively. These results did not completely back previous paleolimnological inferences that assumed that the light calcite layer was formed in spring/summer and the organic layer was formed in fall/winter [26]. Because of the lack of local or regional modern analogue studies, the authors relied on available literature [17] to explain the observed mechanisms of varve formation. According to our results, the process of varve forming during the period of study would probably match one of the varve patterns described in the sediment record of Lake Montcortès, once deposited although this remains to be further confirmed. The seasonal calcite-crystal distribution observed in our quarterly sediment traps would produce calcite sub-layering with coarsening upward, i.e., fine-grained calcite crystals deposited in summer and coarse-grained calcite crystals in fall. This texture is often observed in the sedimentary sequence from CE 1350 to 1850 [26]. Modern analogue studies of crystal sizes in calcite sublayers may allow inferring environmental conditions at seasonal resolution, which would be of paleoclimatic significance. 4.2. Is Contemporary Seasonal Deposition of Pollen Consistent with a Two-Layer Varve Model? The formation of annually laminated lake sediments relies on seasonal deposition linked to climate and on the related flux of biotic and abiotic matter of multiple autochthonous and allochthonous sources to the sediment. Thus, pollen from the catchment and beyond arrives and is buried in lacustrine sediments. Pollen data have been widely used to reconstruct past climate change and have provided insights into the influence of intra-annual weather conditions to pollination patterns. We now argue that combining pollen data from sediment traps and varved lake sediments with modern meteorological data may produce long, high-resolution ecological time series. An example of what can be achieved with this approach are the continuous, long-term and high-resolution climate series obtained by associating paleoclimate data derived from tree rings and other similar proxies with instrumental climate measurements, at annual and seasonal resolutions [39]. To empirically support this hypothesis, this study is primarily intended to identify modern seasonal depositional patterns of pollen in sediment traps deployed in the varved Lake Montcortès that would be useful for interpreting subfossil pollen records from the same lake. Pollen maxima were recorded during spring, concurrently with the flowering season of the involved taxa, and preceded the temperature and precipitation maxima in summer by one month (Figure 3). Total pollen influx and taxonomic composition show a strong seasonal signal for the study period, peaking positively during spring/summer and negatively during fall/winter. The major components of the pollen assemblages, Pinus (pine) and Quercus (oak), dominate the pollen counts and match seasonal trends in temperature and precipitation, with Pinus almost tripling Quercus percentages in spring/summer but falling to similar levels in fall/winter. As far as the remaining species are concerned, most significant differences were observed in Plantago, Chenopodium, Typha/Sparganium, Cyperaceae, Fraxinus and Juniperus/Cupressus, which prevailed in the spring/summer assemblage, and Cannabis and Corylus, which were more abundant in the fall/winter assemblages (Figure 3). Pollen of Cannabis (hemp) prevailed during the fall. Cannabis is a cultivated plant with pollen being present and fairly abundant in the surroundings of Montcortès for the last 1200 years. Nonetheless, we still have not been able to locate the exact source of that pollen [33,34]. Quaternary 2020, 3, 1 8 of 23 Quaternary 2020, 3, x 8 of 23 Figure Figure 3. Summary of 3. Summary of palynological palynological results results showing showing thethe relationship relationship among among the the most most relevant relevant meteorological meteorological parameters (lower panel), the pollen influx of the superabundant types parameters (lower panel), the pollen influx of the superabundant pollen pollen(middle types panel) and the percentages of selected pollen types (upper panel). Average temperature (middle panel) and the percentages of selected pollen types (upper panel). Average temperature (Tm), relative humidity (Hm)humidity (Tm), relative and wind (Hm) velocity and(W) are represented wind velocity (W) byare lines, total precipitation represented by lines,(PTT) totalisprecipitation represented by bars. The predominant directions of the wind (Wd) are shown as circles (Sp: (PTT) is represented by bars. The predominant directions of the wind (Wd) are shown as circles spring, Sm: summer,(Sp: Fl: fall, Wt: winter). Pollen percentages were calculated excluding the abundant Pinus spring, Sm: summer, Fl: fall, Wt: winter). Pollen percentages were calculated excluding the abundant and Quercus. Taxa Pinusare andordered by Taxa Quercus. their respective are ordered flowering by their seasons (grey respective bands) from flowering bottom seasons (greytobands) top andfrom frombottom left to right. to top and from left to right. The flowering season of all species of the different genera present in[40] The flowering season of all species of the different genera present in the Montcortès region the has been considered. Montcortès Trees region [40] hasare shown been in green, Trees considered. shrubsare areshown in orange and herbs in green, are in shrubs areblue. Cultivated in orange and plants, such as Cerealia, Cannabis and families including many genera (Poaceae, Cyperaceae), are herbs are in blue. Cultivated plants, such as Cerealia, Cannabis and families including many genera located based on their pollen patterns due to the difficulty of establishing a definite flowering season. (Poaceae, Cyperaceae), are located based on their pollen patterns due to the difficulty of establishing Redrawn and modified form [19]. a definite flowering season. Redrawn and modified form [19]. Canonical correspondence analysis (CCA) shows the first two axes accounting for 70.7% of total Canonical correspondence analysis (CCA) shows the first two axes accounting for 70.7% of total variance. The strongest gradient coincides with axis 1 (56.8% of the total variance), which was highly variance. The strongest gradient coincides with axis 1 (56.8% of the total variance), which was highly and positively correlated with relative humidity and atmospheric pressure and negatively correlated and positively correlated with relative humidity and atmospheric pressure and negatively with wind velocity (Figure 4). Along this gradient, pollen samples became aligned from spring (left) to correlated with wind velocity (Figure 4). Along this gradient, pollen samples became aligned from winter (right) suggesting a seasonal succession with summer and fall occupying intermediate positions. spring (left) to winter (right) suggesting a seasonal succession with summer and fall occupying Pollen taxa were ordered according to the same gradient. The spring/summer group was highly intermediate positions. Pollen taxa were ordered according to the same gradient. The correlated with temperature, precipitation and wind velocity. The fall/winter group was also highly spring/summer group was highly correlated with temperature, precipitation and wind velocity. The correlated with wind direction from WSW. fall/winter group was also highly correlated with wind direction from WSW. Quaternary 2020, 3, 1 9 of 23 Quaternary 2020, 3, x 9 of 23 Figure Figure 4. Canonicalcorrespondence 4. Canonical correspondence analysis analysis (CCA) (CCA) biplot biplot using using the the scores scores of first of the the first two two axes axes accounting for 70.74% of the total variance. Samples are represented by red dots and pollen taxaand accounting for 70.74% of the total variance. Samples are represented by red dots (for pollen taxa (forsee abbreviations abbreviations seerepresented Figure 3) are Figure 3) arebyrepresented blue dots. by blue dots. Samples Samples are are identified identified seasonally as seasonally as follows: Sp—Spring, follows: Sp—Spring, Sm—Summer, Sm—Summer, Fl—Fall, Fl—Fall, Wt—Winter. Wt—Winter. Meteorological Meteorological parameters parameters are are abbreviated abbreviated as asin inFigure Figure3:3: Tm—average Tm—averagetemperature, temperature, Hm—average Hm—average relative relative humidity, humidity, Pm—average Pm—average pressure, W—average wind velocity, PPT—average precipitation, Wd—average wind pressure, W—average wind velocity, PPT—average precipitation, Wd—average wind direction. Pollen types: direction. Pollen Aln—Alnus, Art—Artemisia, types: Aln—Alnus, Art—Artemisia, Bux—Buxus, Bux—Buxus, Can—Cannabis, Can—Cannabis, Cer—Cerealia, Cer—Cerealia, Che—Chenopodium, Che—Chenopodium, Cor—Corylus, Cyp—Cyperaceae Cor—Corylus, (other), Fra—Fraxinus, Cyp—Cyperaceae Jun—Juniperus/Cupressus, (other), Fra—Fraxinus, Ole—Olea, Poa—Poaceae Jun—Juniperus/Cupressus, (other), Pin—Pinus, Ole—Olea, Pla—Plantago, Poa—Poaceae Que—Quercus, (other), Typ—Typha. Pin—Pinus, Redrawn and Pla—Plantago, modified form [19]. Que—Quercus, Typ—Typha. Redrawn and modified form [19]. Specific aspects of pollen sedimentation require further attention. For example, a lag in pollen Specific aspects of pollen sedimentation require further attention. For example, a lag in pollen sedimentation was observed between production and deposition throughout the year. This lag can sedimentation was observed between production and deposition throughout the year. This lag can have several causes; for example, the influence of internal water dynamics and resuspension, or the have several causes; for example, the influence of internal water dynamics and resuspension, or the fact that soils can be washed into the lake along several successive months after pollen deposition. fact that soils can be washed into the lake along several successive months after pollen deposition. With respect to lake internal processes, it is worth mentioning that strong thermal water column With respect to lake internal processes, it is worth mentioning that strong thermal water column stratification occurred from March to November 2014 coinciding with the main pollen production stratification occurred from March to November 2014 coinciding with the main pollen production period. This stable stratification could have slowed pollen sedimentation rates, whereas in winter period. This stable stratification could have slowed pollen sedimentation rates, whereas in winter vertical mixing of the water disrupted the thermal stability in January and February 2015, facilitating a vertical mixing of the water disrupted the thermal stability in January and February 2015, release of pollen to the lake bottom several months after their production. facilitating a release of pollen to the lake bottom several months after their production. The similarity of sedimentation patterns between Pinus and Quercus pollen was unexpected, The similarity of sedimentation patterns between Pinus and Quercus pollen was unexpected, because their respective pollen grains show important morphological differences that make their because their respective pollen grains show important morphological differences that make their capabilities for air suspension significantly distinct. However, once the pollen is submerged in the lake capabilities for air suspension significantly distinct. However, once the pollen is submerged in the water of Lake Montcortès, the settling of both types of pollen was quite similar even during summer, lake water of Lake Montcortès, the settling of both types of pollen was quite similar even during when the thermal stratification was stable. This finding suggests that internal lake dynamics neutralise summer, when the thermal stratification was stable. This finding suggests that internal lake the effect of such differences and that resuspension or catchment runoff may become more important dynamics neutralise the effect of such differences and that resuspension or catchment runoff may about pollen sedimentation. In general, the seasonal pattern in pollen deposition coincided well with become more important about pollen sedimentation. In general, the seasonal pattern in pollen the suggested two-layered varve model. In summary, pollen analysis was able to identify two well deposition coincided well with the suggested two-layered varve model. In summary, pollen differentiated modern assemblages, one corresponding to spring/summer and the other representing analysis was able to identify two well differentiated modern assemblages, one corresponding to fall/winter. This matches well to the seasonal pattern identified in sedimentological (varve) studies [41]. spring/summer and the other representing fall/winter. This matches well to the seasonal pattern Thus, utilising appropriate transfer functions supports the hypothesis that the combination of pollen identified in sedimentological (varve) studies [41]. Thus, utilising appropriate transfer functions data from sediment traps, varved lake sediments and modern meteorological data may produce long supports the hypothesis that the combination of pollen data from sediment traps, varved lake and high-resolution ecological time series. sediments and modern meteorological data may produce long and high-resolution ecological time series. Quaternary 2020, 3, 1 10 of 23 4.3. Revisiting Meromictic Lake Montcortès: Has the Mixing Regime Changed? Oxygenation of lake water relies on oxygen diffusion from the atmosphere, addition from autochthonous primary production and external inputs of oxygenated water. Lake mixing is the main mechanism, which distributes oxygen throughout the water column and is therefore central for lake health. Mixing regimes can shift with climate changes [42] and human activities [43]. Recently, dissolved oxygen depletion leading to hypoxia or even anoxia has become a matter of concern around the world, and there is a high potential for deterioration under the current global warming [44]. Furthermore, it is suggested that climate change will favour cyanobacteria, which form unpleasant blooms and produce toxins that are harmful to humans [45]. At the same time, obtaining time series of instrumental oxygen records that extend beyond a century is problematic, thus, reducing the possibility of recording long-term changes in oxygen variations that could be associated with climate or human influence. One seminal exception is Lake Zürich, where 50 year long-series of temperature and oxygen profiles corroborated that the increase in the thermal stratification period tends to favour reduction in the homeothermal brake and consequently favours the onset of hypoxia or anoxia [46,47]. Where concurrent series of oxygen and temperature are not available, the combination of proxies of oxic/anoxic conditions with high-resolution paleolimnological records can result in a useful surrogate to assess long-term changes in oxic/anoxic shifts through time. The varved sediment of LM appears to be highly suitable to provide such high-resolution data. Since redox processes depend on the biogeochemical conditions of lakes and the chemical characteristics of the involved compounds, we firstly conducted a modern analogue study consisting of a monthly survey of main redox indicators, i.e., DO, Fe, Mn and marker pigments [48–51], between October 2013 and October 2016. The same type of marker pigments was also examined in sediments: isorenieratene (Chlorobiaceae) as a tracer of euxinia, okenone (Chromatiaceae) as indicators of anoxia and oscillaxanthin (Oscillatoriales) as a general marker of Cyanobacteria that often points to increased eutrophic conditions. Furthermore, we used selected elemental ratios of the sedimentary record to survey shifts in the oxic/anoxic conditions of the lake over the last 500 years: bromine (Br) as a tracer of organic matter [52], the Ca/Ti and Si/Ti ratios as indicators of biologically mediated calcite and silica production, respectively, and the S/Ti and S/Fe ratios as indicators of the presence of sulphur compounds [53]. Lake Montcortès is reported to be meromictic for years [54]; nonetheless, we observed brief episodes of complete mixing during the winters of 2013 and 2014 that triggered precipitation of Fe and Mn oxides. With the onset of anoxic conditions, precipitated Mn and Fe oxides started to redissolve and sulfur bacteria grew in the hypolimnion under anoxic conditions, when Fe and Mn reached maximum concentrations. The oscillaxanthin was absent in the water column during the sampling period. Such conditions endured over the year 2016, which was a non-mixing year. These results confirmed that redox proxies were working well and provided the information required to properly interpret the sedimentary record. Applying this knowledge to the results obtained from sediment analyses, we were able to distinguish four prevailing situations on the basis of selected redox proxies: (A) years with abrupt and substantial sediment input (turbidites); (B) years with mixing and oxygenation of the water column; (C) years with strong stratification, anoxia, intense activity of sulfur bacteria and increased biomass production; (D) years showing stratification and anoxia, but relatively low biomass production (Figure 5). Interestingly, approximately 45.3% of the years were monomictic years, which mostly occurred in groups of several consecutive years. Meromictic years show a similar pattern. Most (A) years coincide with the climatic instability of CE 1850–1899. The (B) years were rather scattered but were best represented between CE 1820–1849. Most (D) years happened from CE 1500 to 1820, when human activities were locally most intense for the studied period. Almost all (C) years belong to the 20th century. More than 90% of the years with climatic instrumental records belonged to B and C. Quaternary 2020, 3, 1 11 of 23 Quaternary 2020, 3, x 11 of 23 Figure 5. Graphical Graphical summary summary of of the the modern modern analogue analogue study study about about oxygen oxygen shifts shifts and and inferred inferred results: results: current seasonal variation of mean values of Fe, Mn, O22 and and phototrophic phototrophic bacteria pigments okenone and isorenieratene (oken) and isorenieratene(isoren) (isoren)ininhypolimnetic hypolimnetic waters waters of Lake of Lake Montcortès; Montcortès; varved varved sediments sediments and and proxies proxies obtained; obtained; graphic graphic synthesis synthesis of of inferences inferences derivedfrom derived fromPrincipal PrincipalComponent Component Analyses; Analyses; decadal evolution evolution of ofconditions conditionsrepresented representedbybythe subsets the A,A,B, subsets B, C Cand D, D, and since CECE since 1500; barsbars 1500; show the show percentage of years that fall within a particular decade. NHTA: Northern Hemisphere the percentage of years that fall within a particular decade. NHTA: Northern Hemisphere temperature anomaly [39].anomaly temperature Redrawn andRedrawn [39]. modifiedandform [20]. form [20]. modified 4.4. Are GDGTs Promising Indicators of Seasonal Temperature Shifts? 4.4. Are GDGTs Promising Indicators of Seasonal Temperature Shifts? Branched glycerol dialkyl tetraethers (brGDGTs) are bacteria derived and widespread lipids of Branched glycerol dialkyl tetraethers (brGDGTs) are bacteria derived and widespread lipids of terrestrial and aquatic environments. The global distribution of brGDGTs in soils and peats has been terrestrial and aquatic environments. The global distribution of brGDGTs in soils and peats has been associated with past temperature, and they are therefore used as proxies of mean annual temperature associated with past temperature, and they are therefore used as proxies of mean annual (MAT proxies) in such environments [55]. However, in some regional studies, brGDGTs indices temperature (MAT proxies) in such environments [55]. However, in some regional studies, appeared to be influenced by other variables, e.g., precipitation, humidity and soil properties [56]. brGDGTs indices appeared to be influenced by other variables, e.g., precipitation, humidity and soil This proxy was originally calibrated against annual averages of environmental variables. However, properties [56]. This proxy was originally calibrated against annual averages of environmental a hypothesis arose that these proxy estimates are biased towards particular seasons, because the variables. However, a hypothesis arose that these proxy estimates are biased towards particular bioproduction of brGDGTs is enhanced under the most favourable conditions, which occur mostly seasons, because the bioproduction of brGDGTs is enhanced under the most favourable conditions, in summer. which occur mostly in summer. To answer whether soil brGDGT proxy estimates are seasonally biased, the brGDGT distributions To answer whether soil brGDGT proxy estimates are seasonally biased, the brGDGT and the brGDGT-derived MAT estimates were examined in settling particulate matter (monthly traps) distributions and the brGDGT-derived MAT estimates were examined in settling particulate◦ matter and surface soil samples from the catchment area of Lake Montcortès (MAT: −3.3 to 17.6 C). No (monthly traps) and surface soil samples from the catchment area of Lake Montcortès (MAT: −3.3 to clear-cut seasonal pattern of the brGDGT distribution has been found in soils (Figure 6a), probably 17.6 °C). No clear-cut seasonal pattern of the brGDGT distribution has been found in soils (Figure because of the slow turnover time in terrestrial environments, which is on timescales of decades 6a), probably because of the slow turnover time in terrestrial environments, which is on timescales of or longer. These results confirm previous findings from mid-latitude soils showing that brGDGTs’ decades or longer. These results confirm previous findings from mid-latitude soils showing that distribution and some of the brGDGT-derived proxy measurements are relatively stable through the brGDGTs’ distribution and some of the brGDGT-derived proxy measurements are relatively stable year [57]. However, we found that annual shifts in abundance of brGDGTs were controlled by the through the year [57]. However, we found that annual shifts in abundance of brGDGTs were variability of specific regional factors, i.e., soil humidity and their pH [13]. controlled by the variability of specific regional factors, i.e., soil humidity and their pH [13]. Quaternary 2020,3,3,1x Quaternary2020, 12 of 12 of 23 23 Figure 6. (A) Concentration of soil brGDGTs, (B) brGDGT flux of particles in the sediment trap at 20 m Figuredepth, water 6. (A)(C) Concentration brGDGT-based of soil brGDGTs,estimates temperature (B) brGDGT flux using theofMBT/CBT particles in the sediment trap at 20 (methylation/cyclisation) m water depth, (C) brGDGT-based temperature estimates using proxy in soils (black line, average from three soils in the catchment around the lake) and the MBT/CBT particles from (methylation/cyclisation) the sediment trap (blue line).proxy (D)in soils (black Fractional line, average abundance from three of brGDGT soils in the distribution catchment in soils around and sediment the lake) traps. and particles Redrawn from from and modified the sediment [21]. trap (blue line). (D) Fractional abundance of brGDGT distribution in soils and sediment traps. Redrawn and modified from [21]. The results record seasonal trends in sediment flux from the catchment area. In the case of The results particulate matter,record heavyseasonal rain is thetrends maininfactor sediment flux from influencing thethe catchment brGDGT area. Inresulting abundance, the case in of particulate seasonal matter, of variations heavy brGDGTrain flux is the main factor (Figures 3 and 6b).influencing the brGDGT abundance, resulting in seasonal variations In terms of brGDGT of modern flux (Figures analogues, this study 3 and 6b).evidence that the signatures of brGDGTs in offers depositional environments are representative of soils inevidence In terms of modern analogues, this study offers that theand the catchment, signatures of brGDGTs that the accuracy in of the depositional derived environments temperature estimatesarebased representative of soils in the catchment, on methylation/cyclisation (MBT/CBT) andproxies that theinaccuracy of the soils depends derived on temperature soil properties, that inestimates turn depend based on methylation/cyclisation on regional factors. Additionally,(MBT/CBT) those brGDGT proxies in soils signatures of depends soils in theon soil properties, catchment probably that in turn average represent depend environmental on regional factors. Additionally, conditions thoseorbrGDGT over decades longer, signatures and of soils means this observation in the that catchment any derivedprobably proxyrepresent average reconstructions canenvironmental be used only to conditions over infer variability decades in or longer, environmental and thisover variables observation the samemeans that any timescales. Thederived absenceproxy reconstructions of seasonality in brGDGT can proxies be used only the from to sediment infer variability traps Lake in Montcortès environmentalis duevariables to the lackover the same“intimescales. of significant The production situ” lacustrine absence of seasonality in brGDGT proxies from the sediment traps Lake Montcortès is due to the lack of Quaternary 2020, 3, 1 13 of 23 of brGDGTs (Figure 6c). The most straightforward explanation is that the main sources of brGDGTs are catchment soils, and that their non-seasonal signals are transferred to the settling particles, which have a distribution of brGDGTs that appear to be a weighted mean of soil signals (Figure 6c). In addition, similar patterns of brGDGT fractional abundance in soils and sediment traps confirm soil-related sources to the particles in sediment traps (Figure 6d). It remains unclear if downcore variability in brGDGTs at an annual scale will also capture changes derived from in situ production, and if derived brGDGT proxies will be suitable to build high-resolution MAT reconstructions. 4.5. From Contemporary Phytoplankton to Subfossil Pigments, What Can We Learn about Community Changes? In general, qualitative and quantitative relationships between modern climatic variables, aquatic primary producers and their pigments and subfossil marker pigments of the sediment have received little attention despite their potential to build reliable climatic proxies. To help fill this gap, we performed a modern analogue study to determine how phytoplankton and marker pigment information compare with the sedimentary record. For this purpose, we applied an analogue matching (AM) technique, which is concerned about identifying contemporary sites that most closely match the species assemblage identified in the past [7]. We assessed and compared the annual cycle of phytoplankton in the epi- and metalimnion and of marker pigments in the epi-, meta- and hypolimnion, and also in a sediment trap deployed at 20 m water depth. Redundancy Analysis (RDA) was applied to identify potential associations between the phytoplankton taxa identified by microscopy. Then, we performed a similarity analysis in order to examine expected bias between modern and subfossil samples (CE 1493–2013) in terms of presence or absence of common pigment markers. The total biovolume of phytoplankton depicted a regular annual cycle with growing peaks in early spring and summer and low biomass in winter (Figure 2B). Phytoplankton biovolume of the metalimnion was highest in summer (>6 × 106 µm3 /mL). Centric Bacillariophyceae and Chlorophytes follow one another as the dominant taxa in 2013–2014, whereas smaller centric Bacillariophyceae took over during 2014–2015 growing periods (Figure 7). The RDA analysis allowed to identify marker pigments more tightly associated with each taxonomic group (Figure 8). Chlorophyta, Dinophyta, Cyanoprokaryota (Cyanobacteria and Cyanophyta), centric Bacillariophyceae, Volvocales and Cryptophyta biovolumes explain a significant percentage of variance in marker pigment concentration, while the percentage variance explained by the other taxonomic groups was insignificant. Each taxonomic group was associated with specific marker pigments. The biovolume of Chlorophyta was associated with the marker pigments neoxanthin, lutein, chlorophyll-b and zeaxanthin, the biovolume of Dinophyta mainly with peridinin, diadinoxanthin, fucoxanthin, Chl-c1 and c2, although Chl-c2 was also partially associated with the biovolume of Cryptophyta. The biovolume of Cryptophyta was mainly associated with alloxanthin and a-carotene, both partially explained by pennate Bacillariophyceae. We know from the literature that some of the pigments associated with Dinophyta are also present in Bacillariophyceae (e.g., fucoxanthin, diadinoxanthin, Chl-c1 and Chl-c2), which is coherent with the observed RDA biplot association. Cyanobacteria biovolume was associated with myxoxanthophyll-like, while association with marker pigments of phototrophic bacteria is due to co-occurrence in the samples with these groups. These data are essential to identify pigments with a higher potential to be used as modern analogues in the sediment record. All this information is relevant to better identify pigments from sediment traps and use them as modern analogues for the sediment record following taxonomic affinity criteria and preservation. presence or absence of common pigment markers. The total biovolume of phytoplankton depicted a regular annual cycle with growing peaks in early spring and summer and low biomass in winter (Figure 2B). Phytoplankton biovolume of the metalimnion was highest in summer (>6 × 106 µm3/mL). Centric Bacillariophyceae and Chlorophytes follow one another as the dominant taxa in 2013–2014, whereas smaller centric Bacillariophyceae Quaternary 2020, 3, 1 14 of 23 took over during 2014–2015 growing periods (Figure 7). Quaternary 2020, 3, x 14 of 23 Figure 7. Taxonomic phytoplankton groups with their biovolume and relative dominance (%) in the Figure epi- 7. Taxonomic (upper graph) andphytoplankton groups the metalimnion with (lower theirfrom graph) biovolume Octoberand 2013relative dominance to September (%) in the 2015 (abscissa). epi- (upper graph) and the metalimnion (lower graph) from October 2013 to September 2015 The concentration of most representative marker pigments was calculated for the entire water (abscissa). sampling period (2013–2015), in order to compare their spatial variations between epi-, meta- and The RDAas analysis hypolimnion, well as withallowed to identifyin marker their deposition pigments the sediment more tightly trap (Figure associated Materials). S1, Supplementary with each taxonomic group (Figure 8). Chlorophyta, Dinophyta, Cyanoprokaryota Marker pigments of Bacillariophyceae (fucoxanthin: Figure S1c) and Dinophyta (diadinoxanthin: (Cyanobacteria and Cyanophyta), Figure S1d) werecentric produced Bacillariophyceae, Volvocales and deposited in higher andin Cryptophyta amounts spring and summer biovolumes explain a while Chlorophyta significant marker percentage pigments of variance (Chl-b, lutein and inzeaxanthin: marker pigment concentration, Figure S1e,f,h) werewhile the percentage produced variance mainly in fall, and explained by the other taxonomic groups was insignificant. Each taxonomic marker pigments of Cryptophyta (alloxanthin: Figure S1g) were produced throughout the year with group was associated withexception the specific of marker pigments. isThe fall. Zeaxanthin biovolume a marker pigmentof also Chlorophyta indicativewas associated with of Cyanobacteria andthe marker decoupled pigments neoxanthin, lutein, chlorophyll-b and zeaxanthin, the biovolume from parent marker pigments of Chlorophyta during summer, it indicates a higher contribution of Dinophyta mainly of with peridinin,during Cyanobacteria diadinoxanthin, summer months fucoxanthin, (FigureChl-c1 and c2,pigments S1). Marker although of Chl-c2 was also partially photosynthetic sulphur associated bacteria werewith only thefound biovolume of Cryptophyta. in the hypolimnion during The biovolume summer of (Figure and fall Cryptophyta was marker S1i,j).Most mainly associated with alloxanthin and a-carotene, both partially explained by pennate pigments from the water column were also identified in the sediment record. A list of the marker Bacillariophyceae. We knowco-occurring pigments from the literature that some in the water of the column andpigments sedimentassociated with in record is given Dinophyta are also present Table S1 (Supplementary in Bacillariophyceae Materials). (e.g., fucoxanthin, Their taxonomic affinities arediadinoxanthin, also shown and Chl-c1 informand Chl-c2), about whichbiological their likely is coherent with sources, the observed RDA biplot association. Cyanobacteria biovolume was which are consistent with the taxonomic phytoplankton groups of the epi- and metalimnion (Figure associated with 7). myxoxanthophyll-like, To while association gauge the similarity between with and water column marker pigments sediment of phototrophic samples we applied the bacteria is due Sorensen to index (S) on a similarity matrix built with binary data. This very simple similarity index provides a greatera co-occurrence in the samples with these groups. These data are essential to identify pigments with higher potential “weight” to be to species used astomodern common analogues the samples in those than to the sediment found inrecord. only one sample. According to S, All this information is relevant to better identify pigments from sediment traps and use them as modern analogues for the sediment record following taxonomic affinity criteria and preservation. Quaternary 2020, 3, 1 15 of 23 epi-, meta- and hypolimnetic samples resemble at ≥80%, whereas meta and hypolimnetic samples resemble each other relatively more (92%), likely because of lower oxidation rates during anoxic conditions of the hypolimnion. Interestingly, when compared with sediment samples of the top of the core (CE 2007–2013 E), S diminished abruptly to 0.5, that is, the bias was 50% (Figure 9). This decrease was followed by a slower decrease of S, until values of 30% (circa 1850) and its posterior stabilisation at values of 40–60% that lasted circa two centuries. A Venn diagram highlights that past communities only have a partial modern analogue (Figure 9) [5]. These results indicate that approximately 50% of the marker pigments in the water column were destroyed between deposition and permanent burying in the first centimetres of the sediment record, whereas the remaining 40–60% were accurately represented over the period of study. This pattern is roughly consistent with the three phases of sediment loss proposed by [48]: (1) rapid oxidation, enzymatic metabolism and digestion by herbivores through the water column (half-life T1/2 = days) would account for a breakdown of pigments while they sink, (2) slower post-depositional loss in surface sediments would result from structural rearrangements, release of labile compounds and further oxidation (T1/2 = years) and (3) once buried Quaternary pigment 2020, 3, x degradation would have continued at a very slow pace (T1/2 = centuries).15 of 23 Figure 8.8.Redundancy Figure Redundancy analysis (RDA) analysis between (RDA) phytoplankton between taxa (biovolume; phytoplankton mm3 ·L−1mm taxa (biovolume; ) and 3·Lmarker −1) and −1 pigments (concentration; µg·L ). Red arrows are explanatory variables (53% of the variance in pigment marker pigments (concentration; µg·L−1). Red arrows are explanatory variables (53% of the variance concentration was explained by phytoplankton taxa), blue arrows are response variables (pigments). in pigment concentration was explained by phytoplankton taxa), blue arrows are response variables Diadin: diadinoxanthin; Fucox: fucoxanthin; Chlc2: chlorophyll-c2; Allox: alloxanthin; a-carot: (pigments). Diadin: diadinoxanthin; Fucox: fucoxanthin; Chlc2: chlorophyll-c2; Allox: alloxanthin; a-carotene; Ph4_6: unknown phorbin; Phytin-a1,a2: phaeophythin-a1, a2; Phytin-b: phaeophythin-b; a-carot: a-carotene; Ph4_6: unknown phorbin; Phytin-a1,a2: phaeophythin-a1, a2; Phytin-b: Car12.5: unknown carotenoid; Myxo-like: myxoxanthophyll-like; BChl-a: bacteriochlorophyll-a; phaeophythin-b; Car12.5: unknown carotenoid; Myxo-like: myxoxanthophyll-like; BChl-a: BChl-e: bacteriochlorophyll-e; Oken: okenone; Isoren: isorenieratene; Neox: neoxanthin; Lut: lutein; bacteriochlorophyll-a; BChl-e: bacteriochlorophyll-e; Oken: okenone; Isoren: isorenieratene; Neox: Chlb: chlorophyll-b; Zeax: zeaxanthin; Bphytin-a: bacteriophaeophythin-a; Astax: astaxanthin; Diato: neoxanthin; Lut: lutein; Chlb: chlorophyll-b; Zeax: zeaxanthin; Bphytin-a: bacteriophaeophythin-a; diatoxanthin; Chla: chlorophyll-a; b-carot: b-carotene; EpimerA: Chl-a epimer; AllomerA: Chl-a Astax: astaxanthin; Diato: diatoxanthin; Chla: chlorophyll-a; b-carot: b-carotene; EpimerA: Chl-a allomer; Violax: violaxanthin; Chlide_a: chlorophyllide-a; Chlc1: chlorophyll-c1; Perid: peridinin. epimer; AllomerA: Chl-a allomer; Violax: violaxanthin; Chlide_a: chlorophyllide-a; Chlc1: chlorophyll-c1; Perid: peridinin. The concentration of most representative marker pigments was calculated for the entire water sampling period (2013–2015), in order to compare their spatial variations between epi-, meta- and hypolimnion, as well as with their deposition in the sediment trap (Figure S1, Supplementary Materials). Marker pigments of Bacillariophyceae (fucoxanthin: Figure S1c) and Dinophyta (diadinoxanthin: Figure S1d) were produced and deposited in higher amounts in spring and summer while Chlorophyta marker pigments (Chl-b, lutein and zeaxanthin: Figure S1e,f,h) were represented over the period of study. This pattern is roughly consistent with the three phases of sediment loss proposed by [58]: (1) rapid oxidation, enzymatic metabolism and digestion by herbivores through the water column (half-life T1/2 = days) would account for a breakdown of pigments while they sink, (2) slower post-depositional loss in surface sediments would result from structural rearrangements, release of labile compounds and further oxidation (T1/2 = years) and Quaternary 2020, 3, 1 3) 16 of 23 once buried pigment degradation would have continued at a very slow pace (T1/2 = centuries). Figure 9. Figure ComparisonofofSorensen’s 9. Comparison Sorensen’ssimilarity similarityindex index(S)(S)between between modern modern and and sediment sediment samples. samples. A A Venn diagram comparing modern and past pigment assemblages is shown in the upper Venn diagram comparing modern and past pigment assemblages is shown in the upper right corner. right corner. Epi-, meta-, Epi-, meta-, hypo- hypo-refer refertotocurrent currentin-depth in-depthlayers of of layers thethe water column. water column.Abcissa: dates Abcissa: inform dates about inform the about age age the of the of sediment samples, the sediment Ordinate: samples, S values. Ordinate: S values. Finally, Table S2 (Supplementary Materials) shows which marker pigments were found to be Finally, Table S2 (Supplementary Materials) shows which marker pigments were found to be exclusive from present-day samples (6) and which were found only in subfossil samples (7). Some of exclusive from present-day samples (6) and which were found only in subfossil samples (7). Some of the unspecific chlorophyll derivatives are produced under distinct environmental conditions and by the unspecific chlorophyll derivatives are produced under distinct environmental conditions and by different mechanisms with time, chiefly through losses of Mg+2 , herbivory, viral attack and enzymatic different mechanisms with time, chiefly through losses of Mg+2, herbivory, viral attack and catalysis [48]. Noteworthy, however, is the presence of four Cyanobacterial markers in the subfossil enzymatic catalysis [48]. Noteworthy, however, is the presence of four Cyanobacterial markers in record that were not present in modern samples, suggesting either a detection limit for the identification the subfossil record that were not present in modern samples, suggesting either a detection limit for of these pigments in the water column or that past lake conditions were favourable to a high diversity the identification of these pigments in the water column or that past lake conditions were favourable of Cyanobacterial groups including N -fixing. In fact, massive blooms of filamentous Planktothrix to a high diversity of Cyanobacterial2 groups including N2-fixing. In fact, massive blooms of rubescens De Candolle ex Gomont have been reported to have flourished in Lake Montcortès in the filamentous Planktothrix rubescens De Candolle ex Gomont have been reported to have flourished in 1970s [54]. However, based on microscopical counts and marker pigment detection, Cyanobacteria Lake Montcortès in the 1970s [54]. However, based on microscopical counts and marker pigment seem to constitute only a low percentage of the phytoplankton currently thriving at Lake Montcortès detection, Cyanobacteria seem to constitute only a low percentage of the phytoplankton currently (Figure 6), an issue that needs further confirmation. thriving at Lake Montcortès (Figure 6), an issue that needs further confirmation. 5. Discussion, Conclusions and Guidelines for Future Research 5. Discussion, Conclusions and Guidelines for Future Research To date, most paleoenvironmental reconstructions for Lake Montcortès have been performed To date, most paleoenvironmental reconstructions for Lake Montcortès have been performed using conventional methods and models, which were originally thought for lower resolution studies using conventional methods and models, which were originally thought for lower resolution studies and higher latitudes. They are therefore constrained to the supra-decadal level. However, it is and higher latitudes. They are therefore constrained to the supra-decadal level. However, it is clear clear that the sublayers of varves host invaluable, high resolution information that still has to been that the sublayers of varves host invaluable, high resolution information that still has to been retrieved. Our modern analogue studies provide detailed information about fine-tuned processes retrieved. Our modern analogue studies provide detailed information about fine-tuned processes that take place at annual and sub-annual scales and show how a suite of widely used proxies that take place at annual and sub-annual scales and show how a suite of widely used proxies respond at the same temporal resolution. We expect that each proxy is physically deposited in respond at the same temporal resolution. We expect that each proxy is physically deposited in a a varve sublayer in temporal correspondence with the processes it represents. Being able to put this exact match in evidence is expected to boost the interpretative power of future high-resolution paleoenvironmental reconstructions. The main contribution of this review is to show the advantages and limitations of a multiproxy modern-analogue approach in Lake Montcortès as a case study. This has not been published previously in any other papers. There exist other studies with a similar approach, e.g., lake Zabiński (Poland) [58] or Lake La Cruz (Spain) [59–62], but results in every case are different because sedimentary and biogeochemical processes differ among lakes. Obviously, results of one particular lake may not be directly applicable at other lakes. Our results demonstrate that current lacustrine processes of Lake Montcortès involving calcite precipitation, phytoplankton and bacterial growth, oxygenation of the water column and pollen production, have a clear seasonality and also a year-to-year variability, depending on subtle interannual climatic variations. Therefore, it is reasonable to expect finding the corresponding fingerprints of these processes at the varve (sublayer) levels. Interestingly, the brGDGT production in soils is not seasonal, while the particle flux into the lake is. The similarity of brGDGT distributions indicated the same sources of brGDGTs in soils and lakes. The brGDGT-estimated Quaternary 2020, 3, 1 17 of 23 temperature reflected the interannual average conditions from the soil records. In order to provide an integrated overview of the basic functioning of Lake Montcortès related to our modern analogue studies, a synthesis of the main processes and their spatial and temporal course is presented (Figure 10). This figure is intended to serve as a guide for high resolution paleoecological and paleoenvironmental reconstructions and derived inferences in lake Montcortès and other similar lakes. It shows how limnological variables and processes give raise to the proxies that will be used in these reconstructions, the transit of the latter from present to past conditions, and where in the varve structure would these proxies be found, according to the season where they are produced and settled down (time-lag = 0 is assumed). In this figure, the limnological variables under study are listed (column 1); moreover, the contemporary processes (column 4) that occur in the lake or lake catchment over the year (columns 2, 3). Columns 5 to 8 show the seasons when these processes occur (blue) and when they are more intense (dark blue). It is also indicated which varve sublayer is representative of each season. White and a dark sublayer correspond to summer and winter, respectively, dark-white transitions to spring and white-dark transitions to fall. Column 9 displays degradation and diagenetic mechanism that are expected to act upon the limnological variables inside the sediment and column 10 the derived proxies. Column 11 suggests which varve sublayers should be targeted to retrieve these proxies for high resolution reconstructions, accepting that one single varve contains all seasonal information needed, i.e., in correspondence with white and dark sublayers and their transitions. We also assume that highly precise sampling methods are available in every case. This is an especially sensitive issue in the case of Lake Montcortès where varves are relatively thin, e.g., total mean varve thickness averages 1.3 mm between the 14th and 20th century, with calcite sublayers 0.35 mm and organic sublayers 0.34 mm [26], and probably in other lakes as well. Note for example that by now, neither diatoms nor pollen or pigments could be extracted from a single varve and even less from a sublayer, due to methodological limitations. Quaternary 2020, 3, 1 18 of 23 Quaternary 2020, 3, x 18 of 23 Modern-analogue variables Catchment Lake Current process Season Transformation in the sediment Sedimentary proxy Expected main target (diagenesis) sublayers Sp (tr d-w) Sm (w) Fl (tr w-d) Wt (d) Endogenic Ca precipitation, primary Calcite saturation index x x production. Calcite deposition and accumulation. Ca/Ti, calcite laminations. w; tr w-d Authigenic formation of Large calcite crystals x calcite crystals, low TºC. Accumulation and growth. Larger calcite crystals. tr w-d; d Biogenic formation of Small calcite crystals calcite crystals, high TºC. Accumulation and growth. Smaller calcite crystals. tr d-w; w Lacustrine primary Particulate organic C (POC) x x production. Organic carbon accumulation/oxydation. Br, S/Ti w; tr w-d Phytoplankton growth Phytoplankton biovolume x and succession. Degradation. Subfossil pigments. w; tr w-d Diatom species x Diatom blooms. Accumulation of frustules. Subfossil diatoms. tr d-w O2 consumption and Oxidised/reduced mineral forms Oxygen x supply . O2 consumption. tr d-w; tr w-d Hypolimnetic oxidation , Fe/Mn, Fe/Ti, absence of O2, Mn,Fe x x mixing. Precipitates of Mn, Fe oxides/hydroxides. anaerobic pigments. d Hypolimnetic reduction, Fe/Mn, Mn/Ti, presence of O2, Mn,Fe x x anoxia. Redissolution of Mn, Fe oxides/hydroxides. anaerobic bacterial pigments. w; tr w-d POLLEN Pinus x Flowering season Partial degradation Pollen abundance tr d-w; w Plantago x Flowering season Partial degradation Pollen abundance w Chenopodium x Flowering season Partial degradation Pollen abundance w Typha/Sparganium x Flowering season Partial degradation Pollen abundance tr d-w; w Fraxinus x Flowering season Partial degradation Pollen abundance tr d-w Juniperus/Cupressus x Flowering season Partial degradation Pollen abundance tr w-d; d Cannabis x Flowering season Partial degradation Pollen abundance tr w-d; d Corylus x Flowering season Partial degradation Pollen abundance d Total pollen x Flowering season Partial degradation Pollen abundance tr d-w; w LIPIDS Bioproduction of brGDGT in soils x brGDGT brGDGT brGDGT in traps x brGDGT drainage Accumulation and degradation brGDGT tr w-d PRIMARY PRODUCERS Pigment production Production of neox, lut, chl-b, Chlorophyta biovolume x Degradation and oxidation zeax. w; tr w-d Pigment production Production of perid, diadinox, Dinophyta biovolume x Degradation and oxidation fucox, Chl-c1 and c2 . tr d-w; w; tr w-d; d Cryptophyta biovolume x Pigment production Degradation and oxidation Production of allox, a-carot, all sublayers Pigment production, Production of allox, a-carot, fucox, diadinox, Chlc1, Chlc2. Bacillariophyta x silicate consumption Degradation and oxidation Diatom frustules. tr d-w Cyanobacteria x Pigment production Degradation and oxidation Production of myxo-like. tr w+d Phototrophic bacteria x Pigment production Degradation and oxidation Production of oken, isoren. w; tr w-d FigureFigure 10. Synthesis of main 10. Synthesis modern of main modernanalogue analoguevariables variables and processesatatLake and processes Lake Montcortès, Montcortès, derived derived proxies proxies andexpected and their their expected target sublayers. target sublayers. White (w); White (w); dark dark (d): white-darktransitions (d): white-dark transitions (w-d (w-d tr) tr) and anddark-white dark-whitetransitions (d-w transitions tr) between (d-w sublayers. tr) between Sp: spring, sublayers. Sm: summer, Sp: spring, Fl: Fall, Wt: Sm: summer, winter. Fl: Fall, Wt: winter. Quaternary 2020, 3, 1 19 of 23 Finally, we provide some guidelines for future research at Lake Montcortès and similar lakes. i. Calcite: There is evidence that a temperature increase enhances the amount of calcite precipitation and flux to the sediment. This relationship may be an important factor, when considering the role of calcite as a sink of atmospheric CO2 in the context of global warming and needs further examination. On the other hand, if the size of endogenic calcite crystals varies seasonally, inferences at the sub-annual scale may be derived from calcite crystals contained in sublayers of varves. ii. These hypotheses encourage the execution of further modern analogue studies. Hence, three improvements are suggested: (1) increase the frequency and duration of sampling, in order to better document shifts in calcite precipitation as related to the current warming trend; (2) install a meteorological station in the lake’s catchment, in order to capture local temperature and other meteorological variables that condition water temperature; (3) perform a detailed survey of calcite crystal features retrieved from seasonal traps and from varve sublayers to precisely determine their origin and diagnostic characteristics. iii. Pollen: Two advantages of pollen studies are (1) the ability of identifying and characterizing seasonal layers even in the absence of varves and (2) the possibility of recording interannual variability and associated meteorological drivers. iv. Oxygen: Even though our approach cannot provide estimations of past oxygen concentrations, it is a useful tool to connect and compare the dynamics of past and present oxic and anoxic events with annual resolution. However, if past oxygen concentrations are to be inferred, a more specific modern analogue study is required. Perhaps monthly and multiannual monitoring of the evolution of hypolimnetic oxygen concentrations and a concurrently examination of the content of Fe and Mn oxides in surface sediment samples, catches mixing and meromictic years. It should then be possible to establish appropriate transfer functions to estimate dissolved oxygen from most related proxies. If successful, truly long-term oxygen time series linking estimated and instrumental data could be established. v. Cyanoprokariota assemblages lack a modern analogue, since they are nearly absent today, whereas they were abundant, diverse and pervasive until, at least, the 1970s [54]. This observation is surprising because, to our knowledge, there have been few changes in land use or lake exploitation that could have affected Lake Montcortès’ trophic conditions. In contrast, we expect Cyanoprokariota to play a prominent role during the coming warmer years, when Lake Montcortès becomes more intensely stratified. Cyanobacteria marker pigments (myxoxanthophyll and zeaxanthin) have been detected in the water column during monitoring, the lack of correspondence with microscopical counts may be due to the low cell size of species involved, but this will need further studies. Therefore, considering the known deleterious effects of Cyanobacteria on aquatic ecosystems and human health, it would be advisable to continue present-day monitoring. The results would aid in obtaining a long-term picture of cyanobacterial successions and importance and allow deriving future scenarios, in the context of global warming. vi. brGDGT production in the soils of the lake catchment show no seasonality. However, in sediment trap particles, the brGDGT flux presents a clear seasonality. The seasonal variations of flux rates confirm transport variations from the catchment rather than in situ production of brGDGT in the lacustrine environment. The MBT/CBT temperature estimates indicate that brGDGT signatures in the sediment traps show a mixed signature and originate from surrounding soils. These results suggest that any inferences made from the high-resolution record must be made with caution. vii. From a methodological point of view, cylindrical traps with seasonal or quarterly recovery yield coherent results in all cases, i.e., total suspended solids, calcite crystals, pollen and spores and particles for brGDGT extractions. However, it is recommended that additional samples of surface sediments at the end of sampling years should be taken, to assess post-depositional transformations of these proxies. In our case, surface sediment samples would have provided Quaternary 2020, 3, 1 20 of 23 complementary information, such as total fluxes to the sediment and subsequent accumulation, interactions at the sediment-water interface and early diagenesis. For example, in the case of calcite precipitation as related to varve formation, initial data of organic matter content may be needed to estimate the degree of sediment compaction with time [63] and for a full understanding of structure and thickness of varve sublayers at any depth. In the end, all this information allows a better connection and comparison between contemporary proxies and their modified version after being buried in the sediment, making the modern analogue technique more powerful. viii. Sometimes, modern analogue studies performed at a seasonal or annual frequency cannot be fully exploited for the interpretation of sedimentary records, because of a mismatch in resolution between present-day and past samples. This mismatch occurs when a relatively large weight or volume of sediment is needed to extract enough material for proxy analyses, because contiguous varves are too thin and/or must be joined into one sample. In doing so, the resolution of the samples decreases. With time and targeted effort, this disadvantage will hopefully be overcome with the advent of new technologies. In the field of marker pigments, hyperspectral image spectroscopy is progressing quickly and offers a non-destructive and inexpensive approach that permits high resolution. Hyperspectral image spectroscopy is being used in lake sediments, e.g., for concentrations of sedimentary bacteriopheophytin “a”, based upon diagnostic spectral properties with high spatial and temporal resolution [64]. A good match of the resolutions of modern and past studies is indispensable to successfully connect both and to better understand the evolution of ecosystems with time. Supplementary Materials: The Supplementary Materials is available online at http://www.mdpi.com/2571-550X/ 3/1/1/s1. Figure S1: Comparison of most representative modern marker pigments from the epi-, meta- and hypolimnion as well as their deposition in the sediment trap, on the basis of annual values of the period 2013–2015; Table S1: Marker pigments co-occurring in the water column and sediment record and their taxonomic affinities for lake Montcortès; Table S2: Marker pigments and their taxonomic affinities for Lake Montcortès. Author Contributions: T.V.-V. conceived and wrote the paper. T.B., V.R., M.d.C.T., M.C., A.R.-M., J.G., J.P.C. and P.L. contributed by refining the text and with data, ideas and figures. E.S., M.d.C.T., N.C., J.G., S.G.-P., P.L., S.G., J.P.C. and T.B. carried out fieldwork and laboratory analyses. J.S. and N.P.-Z. advised on climatic aspects and contributed with meteorological data. All co-authors reviewed the paper. All authors have read and agreed to the published version of the manuscript. Funding: Research funding was granted by the Spanish Ministry of Economy and Competitivity (MINECO/FEDER) with the projects MONT-500, ref. CGL2012-33665; GLOBALKARST, ref. CGL2009-08145), the Institute of Catalan Studies with the projects PIRIMOD and POLMONT, and the Catalan University and Research Management Agency (AGAUR) with the projects 2014 SGR 1207, 2017 SGR 1116). Acknowledgments: We acknowledge the Busseing Pallars Company, the Council of Baix Pallars and the cultural association Lo Vent de Port for their direct involvement in the project and their continuous support. Conflicts of Interest: The authors declare no conflict of interest. References 1. Gentner, D. Structure-mapping: A theoretical framework for analogy. Cogn. Sci. 1983, 7, 155–170. [CrossRef] 2. Gentner, D. The mechanisms of analogical learning. In Similarity and Analogical Reasoning; Vosniadou, S., Ortony, A., Eds.; Cambridge University Press: London, UK, 1989; pp. 199–224. 3. Hesse, M.B. Models and Analogies in Science; University of Notre Dame Press: Notre Dame, IN, USA, 1970; p. 184. 4. Markman, A.B. Constraints on analogical inference. Cogn. Sci. 1997, 21, 373–418. [CrossRef] 5. Delcourt, H.R.; Delcourt, P. The Paleoecological Perspective. In Quaternary Ecology: A Paleoecological Perspective; Chapman Hall: London, UK, 1991; p. 241. 6. Jackson, S.T.; Williams, J.P. Modern Analogs in Quaternary Paleoecology: Here Today, Gone Yesterday, Gone Tomorrow? Ann. Rev. Earth Planet. Sci. 2004, 32, 495–537. [CrossRef] 7. Simpson, G.L. Analogue Methods in Palaeolimnology. In Tracking Environmental Change Using Lake Sediments; Birks, H.J.B., Lotter, A.F., Juggins, S., Smol, J.P., Eds.; Springer Science Business Media: London, UK, 2012; Volume 5, pp. 495–522. Quaternary 2020, 3, 1 21 of 23 8. Legendre, P.; Birks, H.J.B. From classical to canonical ordination. In Tracking Environmental Change Using Lake Sediments; Smol, J.P., Birks, H.J.B., Last, W.M., Eds.; Kluwer: Dordrecht, The Netherlands, 2012; pp. 201–248. 9. Weng, C.; Bush, M.B.; Silman, M.R. An analysis of modern pollen rain on an elevational gradient in southern Peru. J. Trop. Ecol. 2004, 20, 113–124. [CrossRef] 10. Reuss, N.; Leavitt, P.R.; Hall, R.I.; Bigler, C.; Hammarlund, D. Development and application of sedimentary pigments for assessing effects of climatic and environmental changes on subarctic lakes in northern Sweden. J. Paleolimnol. 2010, 43, 149–169. [CrossRef] 11. Sachs, H.; Webb, T.; Clark, D.R. Paleoecological transfer funcions. Ann. Rev. Earth. Planet. Sci. 1977, 5, 159–178. [CrossRef] 12. Massaferro, J.; Larocque-Tobler, I. Using a newly developed chironomid transfer function for reconstructing mean annual air temperature at Lake Potrok Aike, Patagonia, Argentina. Ecol. Ind. 2013, 24, 201–210. [CrossRef] 13. Pla, S.; Catalan, J. Chrysophyte cysts from lake sediments reveal the submillennial winter/spring climate variability in the northwestern Mediterranean region throughout the Holocene. Clim. Dyn. 2005, 24, 263–278. [CrossRef] 14. Wittkop, C.; Teranes, J.; Lubenow, B.; Dean, W.E. Carbon-and oxygen-stable isotopic signatures of methanogenesis, temperature, and water column stratification in Holocene siderite varves. Chem. Geol. 2014, 11, 153–166. [CrossRef] 15. Zolitschka, B.; Francus, P.; Ojala, A.E.; Schimmelmann, A. Varves in lake sediments–A review. Quat. Sci. Rev. 2015, 117, 1–41. [CrossRef] 16. Rull, V. Time continuum and true long-term ecology: From theory to practice. Front. Ecol. Evol. 2014, 2, 1–7. [CrossRef] 17. Brauer, A. Annually laminated lake sediments and their paleoclimatic relevance. In Climate in Historical Time: Towards a Synthesis of Holocene Proxy Data and Climate Models; Fisher, H., Ed.; Springer: Berlin/Heidelberg, Germany, 2004; pp. 108–128. 18. Trapote, M.C.; Vegas-Vilarrúbia, T.; López, P.; Puche, E.; Gomà, J.; Buchaca, T.; Cañellas-Boltà, N.; Safont, E.; Corella, J.P.; Rull, V. Modern sedimentary analogues and integrated monitoring to understand varve formation in the Mediterranean Lake Montcortès (Central Pyrenees, Spain). Palaeogeogr. Palaeoclimatol. Palaeoecol. 2018, 496, 292–304. [CrossRef] 19. Rull, V.; Trapote, M.C.; Safont, E.; Cañellas-Bolta, N.; Pérez-Zanón, N.; Sigró, J.; Buchaca, T.; Vegas-Vilarrúbia, T. Seasonal patterns of pollen sedimentation in Lake Montcortès (Central Pyrenees) and potential applications to high-resolution paleoecology: A 2-year pilot study. J. Paleolimnol. 2017, 57, 95–108. [CrossRef] 20. Vegas-Vilarrúbia, T.; Corella, J.P.; Pérez-Zanón, N.; Buchaca, T.; Trapote, M.C.; López, P.; Sigró, J.; Rull, V. Historical shifts in oxygenation regime as recorded in the laminated sediments of lake Montcortès (Central Pyrenees) support hypoxia as a continental-scale phenomenon. Sci. Tot. Environ. 2018, 612, 1577–1592. [CrossRef] [PubMed] 21. Cao, M.; Rueda, G.; Rivas-Ruiz, P.; Trapote, M.C.; Henriksen, M.; Vegas-Vilarrúbia, T.; Rosell-Melé, A. Branched GDGT variability in sediments and soils from catchments with marked temperature seasonality. Org. Geochem. 2018, 122, 98–114. [CrossRef] 22. Corella, J.P.; Moreno, A.; Morellón, M.; Rull, V.; Giralt, S.; Rico, M.T.; Pérez-Sanz, A.; Valero-Garcés, B.L. Climate and human impact on a meromictic lake during the last 6000 years (Montcortès Lake, Central Pyrenees, Spain). J. Paleolimnol. 2011, 46, 351–367. [CrossRef] 23. Rosell, J. Geological Map of Spain and Report. Scale 1:50,000, Tremp Sheet (252); Instituto Tecnológico Geográfico de España (IGME): Madrid, Spain, 1994. 24. Mercadé, A.; Vigo, J.; Rull, V.; Vegas-Villarrúbia, T.; Garcés, S.; Lara, A.; Cañellas-Boltà, N. Vegetation and landscape around Lake Montcortès (Catalan pre-Pyrenees) as a tool for palaeoecological studies of lake sediments. Collect. Bot. 2013, 32, 87–101. [CrossRef] 25. Corella, J.P.; Benito, G.; Wilhelm, B.; Montoya, E.; Rull, V.; Vegas-Vilarrúbia, T.; Valero-Garcés, B.L. A millennium-long perspective of flood-related seasonal sediment yield in Mediterranean watersheds. Glob. Plan. Sci. 2019, 177, 127–140. [CrossRef] 26. Corella, J.P.; Brauer, A.; Mangili, C.; Rull, V.; Vegas-Vilarrúbia, T.; Morellón, M.; Valero-Garcés, B.L. The 1.5-ka varved record of Lake Montcortès (Southern Pyrenees, NE Spain). Quat. Res. 2012, 78, 323–332. [CrossRef] Quaternary 2020, 3, 1 22 of 23 27. Corella, J.P.; Benito, G.; Rodríguez-Lloveras, X.; Brauer, A.; Valero-Garcés, B.L. Annually resolved lake record of extreme hydro-meteorological events since AD 1347 in NE Iberian Peninsula. Quat. Sci. Rev. 2014, 93, 77–90. [CrossRef] 28. Pierrot, D.; Lewis, E.; Wallace, D.W.R. MS Excel Program Developed for CO2 System Calculations, ORNL/CDIAC-105; Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy: Oak Ridge, TN, USA, 2006. 29. Battarbee, R.W.; Jones, V.; Flower, R.J.; Cameron, N.G.; Bennion, H.; Carvalho, L.; Juggins, S. Diatoms. In Tracking Environmental Change Using Lake Sediments; Smol, J.P., Birks, H.J.B., Last, W.M., Eds.; Kluwer: Dordrecht, The Netherlands, 2001; Volume 3, pp. 155–202. 30. Krammer, K.; Lange-Bertalot, H. Bacillariophyceae. In Süßwasserflora von Mitteleuropa, Bänder 2,3,4/5; Gustav FischerVerlag: Stuttgart, Germany, 2000. 31. Lange-Bertalot, H. Diatoms of the Europe Inland Waters and Comparable Habitats; Gantner: Ruggell, Germany, 2001. 32. Bennett, K.D.; Willis, K.J. Pollen. In Tracking Environmental Change Using Lake Sediments; Smol, J.P., Birks, H.J.B., Last, W.M., Eds.; Kluwer: Dordrecht, The Netherlands, 2002; pp. 5–30. 33. Rull, V.; Vegas-Vilarrúbia, T. Preliminary report on a mid-19th century Cannabis pollen peak in NE Spain: Historical context and potential chronological significance. Holocene 2014, 24, 1378–1383. [CrossRef] 34. Rull, V.; Vegas-Vilarrúbia, T. Crops and weeds from the Estany de Montcortès catchment, Central Pyrenees, during the last millennium: A comparison of palynological and historical records. Veg. Hist. Archaeobot. 2015, 24, 699–710. [CrossRef] 35. Rull, V.; González-Sampériz, P.; Corella, J.P.; Morellón, M.; Giralt, S. Vegetation changes in the southern Pyrenean flank during the last millennium inrelation to climate and human activities: The Montcortès lacustrine record. J. Paleolimnol. 2011, 46, 387–404. [CrossRef] 36. Buchaca, T.; Catalan, J. Factors influencing the variability of pigments in the surface sediments of mountain lakes. Fresh. Biol. 2007, 52, 1365–1379. [CrossRef] 37. Buchaca, T. Pigments Indicadors: Estudi del Senyal en Estanys dels Pirineus i de la Seva Aplicació en Paleolimnologia; Institut d’Estudis Catalans: Catalans, Spain, 2009. 38. Leemann, A.; Niessen, F. Varve formation and climatic record in an Alpine proglacial lake: Calibrating annually-laminated sediments against hydrological and meteorological data. Holocene 1994, 4, 1–8. [CrossRef] 39. Mann, M.E.; Bradley, R.S.; Hughes, M.H. Northern Hemisphere Temperatures during the Past Millennium’ Inferences, Uncertainties, and Limitations. Geophs. Res. Lett. 1999, 26, 759–762. [CrossRef] 40. Bolòs, O.; Vigo, J.; Masalles, R.M.; Ninot, J.M. Flora Manual dels Països Catalans; Editorial Pòrtic: Barcelona, Spain, 2000. 41. Trapote, M.C.; Rull, V.; Giralt, S.; Corella, J.P.; Montoya, E.; Vegas-Vilarrúbia, T. High-resolution (sub-decadal) pollen analysis of varved sediments from Lake Montcortès (southern Pyrenean flank): A fine-tuned record of landscape dynamics and human impact during the last 500 years. Rev. Palaeobot. Palynol. 2018, 259, 207–222. [CrossRef] 42. Hakala, A. Meromixis as a part of lake evolution; observations and a revisedclassification of true meromictic lakes in Finland. Boreal Environ. Res. 2004, 9, 37–53. 43. Jellison, R.; Romero, J.; Melack, J.M. The onset of meromixis during restoration of Mono Lake, California: Unintended consequences of reducingwater diversions. Limnol. Oceanogr. 1998, 43, 706–711. [CrossRef] 44. Jenny, J.P.; Francus, P.; Normandeau, A.; Lapointe, F.; Perga, M.E.; Ojala, A.; Schimmelmann, A.; Zolitschka, B. Global spread of hypoxia in freshwater ecosystems during the last three centuries is caused by rising local human pressure. Glob. Chang. Biol. 2016, 22, 1481–1489. [CrossRef] 45. Elliot, J.A. Is the future blue-green? A review of the current model predictions of how climate change could affect pelagic freshwater cyanobacteria. Water Res. 2014, 46, 1364–1371. [CrossRef] [PubMed] 46. Jankowski, T.; Livingstone, D.M.; Bührer, H.; Forster, R.; Niederhauser, P. Consequences of the 2003 European heat wave for lake temperature profiles, thermal stability, and hypolimnetic oxygen depletion: Implications for a warmer world. Limnol. Oceanogr. 2006, 51, 815–819. [CrossRef] 47. Rempfer, J.; Livingstone, D.M.; Blodau, C.; Forster, R.; Niederhauser, P.; Kipfer, R. The effect of the exceptionally mild European winter of 2006–2007 on temperature and oxygen profiles in lakes in Switzerland: A foretaste of the future? Limnol. Oceanogr. 2010, 55, 2170–2180. [CrossRef] Quaternary 2020, 3, 1 23 of 23 48. Leavitt, P.R. A review of factors that regulate carotenoid and chlorophyll deposition and fossil pigment abundance. J. Paleolimnol. 1993, 9, 109–127. [CrossRef] 49. Villanueva, J.; Grimalt, J.O.; De Wit, R.; Keely, B.J.; Maxwell, J.R. Sources and transformations of chlorophylls and carotenoids in a monomictic sulphate-rich karstic lake environment. Org. Geochem. 1994, 22, 739–757. [CrossRef] 50. Hodgson, D.A.; Wright, S.W.; Tyler, P.A.; Davies, N. Analysis of fossil pigments from algae and bacteria in meromictic Lake Fidler, Tasmania, and its application to lake management. J. Paleolimnol. 1998, 19, 1–22. [CrossRef] 51. North, R.P.; North, R.L.; Livingstone, D.M.; Köster, O.; Kipfer, R. Long-term changes in hypoxia and soluble reactive phosphorus in the hypolimnion of a large temperate lake: Consequences of a climate regime shift. Glob. Chang. Biol. 2014, 20, 811–823. [CrossRef] 52. Moreno, A.; Giralt, S.; Valero-Garcés, B.; Sáez, A.; Bao, R.; Prego, R.; Pueyo, J.J.; González-Sampériz, P.; Taberner, C. A 14kyr record of the tropical Andes: The Lago Chungará sequence (18◦ S, northern Chilean Altiplano). Quat. Int. 2007, 161, 4–21. [CrossRef] 53. Davies, S.J.; Lamb, H.F.; Roberts, S.J. Micro-XRF core scanning in Palaeolimnology:recent developments. In Micro-XRF Studies of Sediment Cores: Applications of a Non-destructive Tool for the Environmental Sciences; Croudace, I.W., Rothwell, R.G., Eds.; Springer: Dordrecht, The Netherlands, 2015; pp. 189–226. 54. Camps, J.; Gonzalvo, I.; Güell, J.; López, P.; Tejero, A.; Toldrà, X.; Vallespinós, F.; Vicens, M. El lago de Montcortès, descripción de un ciclo anual. Oecologia. Aquat. 1976, 2, 99–110. 55. Weijers, J.W.H.; Schouten, S.; van den Donker, J.C.; Hopmans, E.C.; Sinninghe Damsté, J.S. Environmental controls on bacterial tetraether membrane lipid distribution in soils. Geoch. Cosmoch. Acta 2007, 71, 703–713. [CrossRef] 56. Dang, X.; Yang, H.; Naafs, B.D.A.; Pancost, R.D. Direct evidence of moisture control on the methylation of branched glycerol dialkyl glycerol tetraethers in semi-arid and arid soils. Geoch. Cosmoch. Acta. 2016, 189, 24–36. [CrossRef] 57. Weijers, J.W.H.; Bernhardt, B.; Peterse, F.; Werne, J.P.; Dungait, J.A.J.; Schouten, S.; Sinninghe Damsté, J.S. Absence of seasonal patterns in MBT-CBT indices in mid-latitude soils. Geoch. Cosmoch. Acta 2011, 75, 3179–3190. [CrossRef] 58. Bonk, A.; Tylmann, W.; Amann, B.; Enters, D.; Grosjean, M. Modern limnology and varve-formation processes in lake Zabińskie, Northeastern Poland: Comprehensive process studies as a key to understand the sediment record. J. Limnol. 2015, 74, 358–370. [CrossRef] 59. Miracle, M.R.; Camacho, A.; Juliâ, R.; Vicente, E. Sinking processes and their effect on the sedimentary record in the meromictic Lake La Cruz (Spain). Verh. Int. Verein. Limnol. 2000, 27, 1209–1213. [CrossRef] 60. Romero, L.; Camacho, A.; Vicente, E.; Miracle, M.R. Sedimentation patterns of photosynthetic bacteria based on pigment markers in meromictic Lake La Cruz (Spain): Paleolimnological implications. J. Paleolimnol. 2006, 35, 167–177. [CrossRef] 61. Rodrigo, M.A.; Vicente, E.; Miracle, M.R. Short-term calcite precipitation in the karstic meromictic Lake La Cruz (Cuenca, Spain). Verh. Int. Ver. Limnol. 1993, 25, 711–719. [CrossRef] 62. Rodrigo, M.A.; Vicente, E.; Miracle, M.R. The meromictic Lake La Cruz (Central Spain): Patterns of stratification. Aquat. Sci. 2001, 63, 406–416. [CrossRef] 63. Maier, D.; Rydberg, J.; Bigler, C.; Renberg, I. Compaction of recent varved lake sediments. GFF 2013, 135, 231–236. [CrossRef] 64. Butz, C.; Grosjean, M.; Fischer, D.; Wunderle, S.; Tylmann, W.; Rei, B. Hyperspectral imaging spectroscopy: A promising method for the biogeochemical analysis of lake sediments. J. Appl. Remote Sens. 2015, 9, 1–20. [CrossRef] © 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).