Long-Term Consequences of Water Pumping on the Ecosystem Functioning of Lake Sekšu, Latvia

: Cultural eutrophication, the process by which pollution due to human activity speeds up natural eutrophication, is a widespread and consequential issue. Here, we present the 85-year history of a small, initially Lobelia – Isoëtes dominated lake. The lake’s ecological deterioration was intensiﬁed by water pumping station activities when it received replenishment water for more than 10 years from a eutrophic lake through a pipe. In this study, we performed a paleolimnological assessment to determine how the lake’s ecosystem functioning changed over time. A multi-proxy (pollen, Cladocera, diatoms, and Chironomidae) approach was applied alongside a quantitative reconstruction of total phosphorus using diatom and hypolimnetic dissolved oxygen with chironomid-based transfer functions. The results of the biotic proxy were supplemented with a geochemical analysis. The results demonstrated signiﬁcant changes in the lake community’s structure, its sediment composition, and its redox conditions due to increased eutrophication, water level ﬂuctuations, and erosion. The additional nutrient load, particularly phosphorus, increased the abundance of planktonic eutrophic–hypereutrophic diatoms, the lake water’s transparency decreased, and hypolimnetic anoxia occurred. Cladocera, Chironomidae, and diatoms species indicated a community shift towards eutrophy, while the low trophy species were suppressed or disappeared.

temperature  in areas close to the Baltic Sea coast is +6.8 to +7.4 °C. July is the warmest month of the year, with an average air temperature +17.4 °C (average maximum +22.3 °C). The coldest month of the year is February, with an average air temperature of −3.7 °C. The annual rainfall is 692 mm. According to the data on climate change, the air temperature and precipitation in the area are increasing [13].
The lake is mostly surrounded by inland dune forest and dominated by pine growing on sandy soil. The southern part of the lake is dominated by birches growing on organic soils [14]. The lake has no runoff apart from a small, shallow ditch that transports humic substances to the lake. In this area, there is also a peat soil-based, meliorated forest [9,13].
According to the EU Water Framework Directive criteria (total phosphorus, Secchi depth, and chlorophyll a), the lake's current ecological status is good. The latest macrozoobenthos studies confirmed the lake's high biodiversity. Following Carlson's trophic state index, in terms of water transparency (SD), Lake Sekšu is a mesotrophic lake, while its chlorophyll a concentration (CA) indicates a eutrophic lake status. The average value of the indices corresponds to a eutrophic lake. Nevertheless, the lake still features low eutrophication. The main problem for such lakes is a loss of transparency [15].

Lake Sediment Coring
A 46 cm sediment core from the deepest part of Lake Sekšu was taken on 13th February 2019 using a Kayak/HTH gravity-type corer. The sediment core was divided in the field into 1 cm sections

Lake Sediment Coring
A 46 cm sediment core from the deepest part of Lake Sekšu was taken on 13th February 2019 using a Kayak/HTH gravity-type corer. The sediment core was divided in the field into 1 cm sections and stored in a cold room. In the laboratory, 1 cm 3 subsamples of fresh sediment from each section were taken for analysis of the biological (Cladocera, Chironomidae, diatoms, pollen) and spheroidal carbonaceous particles (SCPs). Dried material was used for 137 Cs/ 210 Pb dating, as well as chemical and physical analysis.

Core Chronology
A sediment core from Lake Sekšu was dated with 137 Cs and 210 Pb at the Geochronology Laboratory at the Gdańsk University according to the standard procedure [16]. The activity of 137 Cs and 226 Ra was determined directly by gamma-ray spectrometry. Gamma measurements were carried out using a HPGe well-type detector (GCW 2021, Canberra) coupled to a multi-channel analyzer and shielded by a 15 cm thick layer of lead. Counting efficiency was determined using reference materials (CBSS-2 for 137 Cs at 661.6 keV, and RGU-1 for 226 Ra via 214 Pb at 352 keV) with the same measurement geometry as the samples. The counting time for each sediment sample was 24 h.
The activity of total 210 Pb was determined indirectly by measuring 210 Po using alpha spectrometry. Dry and homogenized sediment samples of 0.2 g were transferred into Teflon containers, spiked with a 209 Po yield tracer, and digested with concentrated HNO 3 , HClO 4 , and HF at a temperature of 100 • C using a CEM Mars 6 microwave digestion system. After 24 h, the solution was transferred into a Teflon beaker, evaporated with 6 M HCl to dryness, and then dissolved in 0.5 M HCl. Polonium isotopes were spontaneously deposited within 4 h on silver discs. After deposition, the discs were analyzed for 210 Po and 209 Po using a 7200-04 APEX Alpha Analyst integrated alpha-spectroscopy system (Canberra) equipped with PIPS A450-18AM detectors. The samples were counted for 24 h. A certified mixed alpha source ( 234 U, 238 U, 239 Pu, and 241 Am; SRS 73833-121, Analytics, Atlanta, Georgia, USA) was used to check the detector counting efficiencies, which varied from 30.9% to 33.9% for the applied geometry.
In addition, an abundance of spheroidal carbonaceous particles (SCPs) was estimated throughout the sediment sequence following the methodology of Rose [17] and Alliksaar [18]. The analysis was performed in the Department of Geography, University of Latvia. According to the black carbon combustion continuum model of Hedges et al. [19] and Masiello [20], SCPs only form during industrial fuel combustion at high temperatures (greater than 1000 • C). Therefore, the peak followed the fuel combustion pattern: 1950-the rise of SCPs, 1982-the peak of SCPs, and 1991-the decrease of SCPs. The peak SCP emissions were previously established for Latvia at 1982 ± 10 years [21].
In the final step, we combined the results from the radionuclides and the SCP analyses to build an age-depth model using the Clam 2.2 deposition model [22] with a 95.4% confidence level in the R environment [23].

Physical and Chemical Sediment Analyses
Sediment geochemical characteristics were determined using loss on ignition (LOI) combustion analysis and inductively coupled plasma-optical emission spectrometry (ICP-OES). Altogether, 45 subsamples at 1 cm intervals were analyzed, but only the topmost two sediment samples were merged into one sample representing depths from 0 to 2 cm. The sediment organic matter and carbonate content were investigated using the LOI method [24]. A measure of 0.1-0.2 g of fresh sub-sample sediment was weighed in a crucible and dried at 105 • C for 12 h; then, it was combusted at 550 • C for 4 h and, finally, at 950 • C for 2 h. Between each step, the samples were cooled with an exicator and weighed. The organic matter (OM) content was measured as the LOI from the combusted samples at 550 • C, and the carbonate matter (CM) was calculated as the difference between the LOI at 950 • C and the LOI at 550 • C multiplied by 1.36 [24]. Non-carbonate siliciclastic matter, here referred to as minerogenic matter (MM) content, was obtained by subtracting OM and CM from the total sample weight after final combustion. MM also includes biogenic silica (opal). However, siliciclastic matter here mainly represents terrigenous clastic matter.
A set of homogenized sub-samples was freeze dried and powdered using an agate mortar. The total carbon (TC) and total nitrogen (TN) contents were analyzed using thermal combustion elemental analyses (Element Analyzer Vario EL III) with an uncertainty of ±5%. Approximately 7-8 mg Water 2020, 12, 1459 5 of 21 of prepared material was weighed into tin cups for analyses that were performed at an accredited laboratory (ISO/IEC/17025).
For the ICP-OES analyses, 0.1-0.2 g of dry sediment powder was weighed in Teflon vials, and 5 mL of 65% HNO 3 was added. The vials were closed loosely, and digestion was carried out at 160 • C for 30 min. Gases and volatiles from the digestion were allowed to escape through the cap. Following digestion, 10 mL H 2 O was added to dilute the remaining acid, and the vials were weighed to determine the dilution factor. A measure of 1 mL of the resulting solution was then pipetted into 15 mL centrifuge tubes, and 9 mL 0.5 mol/l HNO 3 was added. These steps were performed to give the samples a final acid concentration of approximately 1 mol/l HNO 3 and to have the element concentrations in a suitable range for the ICP-OES analysis. The sulfur (S) concentration was determined by ICP-OES (Thermo Scientific iCAP 6000) using a concentric nebulizer. Certified control samples and blank samples were used to ensure the quality of the analytical process.

Pollen and Non-Pollen Palynomorphs
Samples of known volume for the pollen analysis were processed using standard procedures [25]. Known quantities of Lycopodium spores were added to each sample to allow the calculation of pollen concentrations [26]. At least 500 terrestrial pollen grains per sample were counted under a light microscope (400× magnification). Taxa were identified to the lowest possible taxonomic level using the reference collection at the Department of Geography at the University of Latvia along with published pollen keys [27]. The percentage of dry-land taxa was estimated using arboreal and non-arboreal pollen sums (excluding the sporomorphs of aquatic and wetland plants). Counts of spores were calculated as the percentages of the total sum of terrestrial pollen. Non-pollen palynomorphs were recorded throughout the pollen analysis and identified using the published literature listed in Miola [28], as well as from the descriptions of Sweeney [29] and Finsinger and Tinner [30]. Non-pollen palynomorphs were expressed as presence or concentrations. The pollen diagram was compiled using TILIA software [31].

Diatom Analysis
Samples were prepared according to the standard methods [32]. The material was treated with 10% HCl to remove calcareous matter, washed with distilled water, and then treated with 30% H 2 O 2 in a water bath to remove organic matter. The material was repeatedly washed with distilled water, and a known number of microspheres in a solution (concentration 8.02 × 10 6 microspheres/cm 3 ) was added to the diatom suspensions to estimate the diatom concentrations [33]. A few drops of diatom suspension were dried on a cover glass. At least 300 diatom valves per sample were analyzed using oil immersion at 1000× magnification under a light microscope. For identification, a selection of published keys was used [34][35][36][37][38][39][40][41].
The diatom ecological groups were determined using the OMNIDIA software (Version 4.2) [42]. Next, the resulting groups were distinguished according to Denys [43] and van Dam [44]. We considered the following indicator parameters: habitat category [43], dominant taxa (abundance over 2%) and preference for pH, saprobic level (OM contamination), and trophy [44]. The percentage diatom diagram was prepared with the Tilia software [31].
The total phosphorus (TP) concentration was reconstructed based on changes in the diatom species composition (DI-TP). The reconstruction was performed using the European Diatom Database (EDDI) in the ERNIE software [45]. A model based on inverse regression had a root mean square error of prediction (RMSEP) of 0.33 µg L −1 and a coefficient of determination (r 2 ) of 0.64. The reconstruction of the TP was based on the diatom taxa present at more than 2% abundance. The DI-TP was calculated using the combined TP dataset (derived from nine datasets with 347 samples in total), covering a TP range of 2-1189 µg L −1 , with a mean of 98.6 µg L −1 . The weighted averaging (WA) method with good empirical predictive ability was used [45]. The 46 fresh sediment samples were prepared in a laboratory according to the standard procedures [46], heated in 10% KOH, and sieved using a 38 µm mesh size. Microscope slides were prepared from 0.1 mL of each sample and examined with a light microscope under magnifications of ×100, ×200, and ×400. For each sample, 1-3 slides were scanned, and all skeletal elements (head shields, shells, and postabdomens) were counted until 70-100 individuals were found, which is regarded as an adequate number to characterize the assemblages [47]. Identification of the cladoceran remains was based on the key by Szeroczyńska and Sarmaja-Korjonen [46]. The stratigraphic diagrams presenting the results were prepared using C2 freeware [48].
The Cladocera composition is presented in the stratigraphic diagrams with percentage values. For trophic state reconstruction, we used changes in the percentage of species regarded as indicators for the eutrophic state, which, according to Flössner [49], are Alona rectangula and Chydorus cf. sphaericus. Lake water level changes were reconstructed by using the added percentage values of planktonic and littoral cladocerans.

Chironomidae Analysis
Standard methods were applied in the fossil Chironomidae analysis [50]. The wet sediment was sieved through mesh (100 µm), and the residue was examined under a stereomicroscope for larval head capsule extraction using a target counting sum of 50 per sample. The head capsules were mounted with Euparal on microscope slides for taxonomic identification following Brooks et al. [50] under a light microscope (400× magnification).

Core Chronology
The depth profile of 137 Cs activity is long and smooth, with only one wide maximum between a 28 and 38 cm sediment depth ( Figure 2). The sharp decrease below 41 cm indicates the presence of sediments older than 1950 below this depth. The lack of two independent peaks indicating global fallout during 1961-1965 and the Chernobyl peak in 1986 might suggest the deep mixing of surface sediments.
However, the excess 210 Pb activities decrease regularly with mass depth, which demonstrates no significant disturbance of the sediment column and relatively stable mass accumulation rates. Thus, we used the CFCS (Constant Flux Constant Sedimentation) model to calculate the mean value of the mass accumulation rate (MAR) for the entire core. The mean MAR value is 34.5 ± 2.2 mg/cm 2 /yr and allows us to estimate a maximum age of 84 ± 5 yrs (1935 ± 5) at a 45.5 cm sediment depth. The age of the depth interval with the greatest 137 Cs activities (37-32 cm) is 1956-1965 according to the CFCS model, which is consistent with global fallout history. This comparison suggests that the lack of a Chernobyl peak is related to 137 Cs migration within the sediment column rather than physical sediment mixing. Based on the SCP results, it was possible to define the years 1950 and 1991 with an error of +/− 10 years at depths 42 and 17 cm, respectively. The locations of these peaks are, with relative certainty, consistent with the 137 Cs and 210 Pb data.
The final age-depth model, including the radionuclide and SCP data, shows relatively stable sedimentation rates with only a moderate increase in the topmost part of the profile (Figure 3). The maximum age of sediment at a depth of 46 cm is 84 ± 9 yrs (95%). The mean sedimentation rate is 0.6 cm per year.  Based on the SCP results, it was possible to define the years 1950 and 1991 with an error of +/− 10 years at depths 42 and 17 cm, respectively. The locations of these peaks are, with relative certainty, consistent with the 137 Cs and 210 Pb data.

Sediment Composition
The final age-depth model, including the radionuclide and SCP data, shows relatively stable sedimentation rates with only a moderate increase in the topmost part of the profile ( Figure 3). The maximum age of sediment at a depth of 46 cm is 84 ± 9 yrs (95%). The mean sedimentation rate is 0.6 cm per year. Based on the SCP results, it was possible to define the years 1950 and 1991 with an error of +/− 10 years at depths 42 and 17 cm, respectively. The locations of these peaks are, with relative certainty, consistent with the 137 Cs and 210 Pb data.
The final age-depth model, including the radionuclide and SCP data, shows relatively stable sedimentation rates with only a moderate increase in the topmost part of the profile ( Figure 3). The maximum age of sediment at a depth of 46 cm is 84 ± 9 yrs (95%). The mean sedimentation rate is 0.6 cm per year.

Sediment Composition
The Lake Sekšu sediment composition is dominated by organic matter (50%-61%) with a significant component of minerogenic matter (35%-46%, Figure 4). The sediments contain 2%-5% carbonate matter. Based on sediment composition and the variation in element concentrations, the sedimentary data are divided into two zones, Zones I and II, which correspond to the approximate time period before 1950 and between 1950 and 2018. Within Zone II, short-lived events can be distinguished.
Water 2020, 12, x FOR PEER REVIEW 8 of 21 The Lake Sekšu sediment composition is dominated by organic matter (50%-61%) with a significant component of minerogenic matter (35%-46%, Figure 4). The sediments contain 2%-5% carbonate matter. Based on sediment composition and the variation in element concentrations, the sedimentary data are divided into two zones, Zones I and II, which correspond to the approximate time period before 1950 and between 1950 and 2018. Within Zone II, short-lived events can be distinguished. Zone I The relative organic matter (OM), total carbon (TC), and total nitrogen (TN) content show high values in this time period ( Figure 4). The OM is around 61%, while TC is around 32%, and TN is around 2.5%. The C/N value, expressed as the atomic ratio, is around 15.7. The sulfur (S) content is approximately 500 ppm. The relative minerogenic matter (MM) (36%) displays the lowest values, while approximately 3.5% is carbonate matter (CM).
Zone II At the boundaries of Zone I and Zone II, a large change occurs. The OM content shows a sudden drop from 61% to 51% accompanied by a fall in TC and TN content from 33% to 27% and from 2.45% to 1.95%, respectively. Simultaneously, the MM content increases from 36% to 45%. After this rapid change, a steady gradient emerges from 41 cm to the sediment surface, during which MM, OM, and TC gradually reset towards their Zone I value. However, TN increases over the same period to values far exceeding those of Zone I (approximately 2.8% at the sediment surface, Figure  4). Consequently, C/N displays a continuous gradient towards lower values throughout the entire core. Similarly, S content increases gradually throughout the entire core, with some pronounced variability within Zone II.

Pollen and Non-Pollen Palynomorphs
There is no significant change in surrounding vegetation before or after the pipe installation. The vegetation includes the stable dominance of pine (Pinus), birch (Betula), spruce (Picea), and alder (Alnus) over the studied period ( Figure 5). The presence of conifer stomata in the lake sediment supports local abundance of spruce and pine. Although, the surrounding landscape is forested, there is continuous evidence of human-activities in the vicinity. For instance, the presence of pollen of flax (Linum), rye (Secale cereale), barley (Hordeum), and wheat/oat (Triticum/Avena) is a direct indication of agricultural practices. Our results also underline that there have not been agricultural fields directly at the Sekšu lake shores, but only further away, as evidenced by the high (nearly 85-90 %) forest Zone I The relative organic matter (OM), total carbon (TC), and total nitrogen (TN) content show high values in this time period (Figure 4). The OM is around 61%, while TC is around 32%, and TN is around 2.5%. The C/N value, expressed as the atomic ratio, is around 15.7. The sulfur (S) content is approximately 500 ppm. The relative minerogenic matter (MM) (36%) displays the lowest values, while approximately 3.5% is carbonate matter (CM).
Zone II At the boundaries of Zone I and Zone II, a large change occurs. The OM content shows a sudden drop from 61% to 51% accompanied by a fall in TC and TN content from 33% to 27% and from 2.45% to 1.95%, respectively. Simultaneously, the MM content increases from 36% to 45%. After this rapid change, a steady gradient emerges from 41 cm to the sediment surface, during which MM, OM, and TC gradually reset towards their Zone I value. However, TN increases over the same period to values far exceeding those of Zone I (approximately 2.8% at the sediment surface, Figure 4). Consequently, C/N displays a continuous gradient towards lower values throughout the entire core. Similarly, S content increases gradually throughout the entire core, with some pronounced variability within Zone II.

Pollen and Non-Pollen Palynomorphs
There is no significant change in surrounding vegetation before or after the pipe installation. The vegetation includes the stable dominance of pine (Pinus), birch (Betula), spruce (Picea), and alder (Alnus) over the studied period ( Figure 5). The presence of conifer stomata in the lake sediment supports local abundance of spruce and pine. Although, the surrounding landscape is forested, there is continuous evidence of human-activities in the vicinity. For instance, the presence of pollen of flax (Linum), rye (Secale cereale), barley (Hordeum), and wheat/oat (Triticum/Avena) is a direct indication of agricultural practices. Our results also underline that there have not been agricultural fields directly Water 2020, 12, 1459 9 of 21 at the Sekšu lake shores, but only further away, as evidenced by the high (nearly 85-90%) forest pollen component in the landscape. It is interesting that the rye pollen accumulation rates above 1000 cm −2 year −1 ( Figure 6) point to cereal fields within a 2 km radius of Lake Sekšu [53,54]. Based on USSR topographic maps from the time period of 1941-1991, the cereal fields could be present in the northwest and east of the Sekšu Lake where pastoral activities have also been evident. It is, however, possible that the high concentration of rye pollen can be partially a result of soil erosion and water pumping (enlarging the pollen source area outside the watershed) into the Lake Sekšu. This assumption is further supported by a strong increase of fungi hyphae and corroded pollen grains.
Regarding in-lake vegetation, lake quillwort/Merlin's grass (Isoëtes lacustris) was recorded throughout the sediment sequence. The highest relative abundance of I. lacustris was recorded prior the pipe installation, after which values continuously declined and did not reach previous values. Sporadically, the pollen of water lilies (Nymphaeaceae), bulrush (Typha), spiked water-milfoil (Myriophyllum spicatum), pondweed (Potamogeton), and bur-reed (Sparganium) were found.
Water 2020, 12, x FOR PEER REVIEW 9 of 21 pollen component in the landscape. It is interesting that the rye pollen accumulation rates above 1000 cm −2 year −1 ( Figure 6) point to cereal fields within a 2 km radius of Lake Sekšu [53,54]. Based on USSR topographic maps from the time period of 1941-1991, the cereal fields could be present in the northwest and east of the Sekšu Lake where pastoral activities have also been evident. It is, however, possible that the high concentration of rye pollen can be partially a result of soil erosion and water pumping (enlarging the pollen source area outside the watershed) into the Lake Sekšu. This assumption is further supported by a strong increase of fungi hyphae and corroded pollen grains. Regarding in-lake vegetation, lake quillwort/Merlin's grass (Isoëtes lacustris) was recorded throughout the sediment sequence. The highest relative abundance of I. lacustris was recorded prior the pipe installation, after which values continuously declined and did not reach previous values. Sporadically, the pollen of water lilies (Nymphaeaceae), bulrush (Typha), spiked water-milfoil (Myriophyllum spicatum), pondweed (Potamogeton), and bur-reed (Sparganium) were found.

Diatom Analysis
The results show a medium to deplorable state for the diatom frustules, with numerous traces of destruction and dissolution. In total, 178 species of diatoms were identified. There were 47 pollen component in the landscape. It is interesting that the rye pollen accumulation rates above 1000 cm −2 year −1 ( Figure 6) point to cereal fields within a 2 km radius of Lake Sekšu [53,54]. Based on USSR topographic maps from the time period of 1941-1991, the cereal fields could be present in the northwest and east of the Sekšu Lake where pastoral activities have also been evident. It is, however, possible that the high concentration of rye pollen can be partially a result of soil erosion and water pumping (enlarging the pollen source area outside the watershed) into the Lake Sekšu. This assumption is further supported by a strong increase of fungi hyphae and corroded pollen grains. Regarding in-lake vegetation, lake quillwort/Merlin's grass (Isoëtes lacustris) was recorded throughout the sediment sequence. The highest relative abundance of I. lacustris was recorded prior the pipe installation, after which values continuously declined and did not reach previous values. Sporadically, the pollen of water lilies (Nymphaeaceae), bulrush (Typha), spiked water-milfoil (Myriophyllum spicatum), pondweed (Potamogeton), and bur-reed (Sparganium) were found.

Diatom Analysis
The results show a medium to deplorable state for the diatom frustules, with numerous traces of destruction and dissolution. In total, 178 species of diatoms were identified. There were 47

Diatom Analysis
The results show a medium to deplorable state for the diatom frustules, with numerous traces of destruction and dissolution. In total, 178 species of diatoms were identified. There were 47 dominant taxa, whose share was more than 2% of the relative abundance. Based on the changes in species composition and the proportions between the ecological diatom groups, the data are divided into two zones, Zones I and II, which closely correspond to the geochemical zones (Figures 7 and 8). Species were classified according to their habitat category, preference for pH, saprobic state, and trophic state (Figure 7). The total phosphorous (TP) reconstruction presented values between 21.43 and 102.14 µg L −1 (Figure 11). trophic state (Figure 7). The total phosphorous (TP) reconstruction presented values between 21.43 and 102.14 μg L −1 (Figure 11).
Zone I This zone was characterized by the dominance of tychoplanktonic taxa, especially Staurosira construens and Achnanthidium minutissimum (Figure 8). In terms of pH, circumneutral taxa dominated. At the beginning of the phase, we observed an increase in alkaliphilous diatoms. The analysis of saprobic preferences revealed the domination of β-mesosaprobous diatoms. In terms of its trophic state, a high proportion of oligo to eutraphentic taxa was observed. The TP reconstruction remained low, between 21.43-57.66 μg L −1 (Figure 11).

Cladocera Analysis
Zone I This zone was characterized by the dominance of tychoplanktonic taxa, especially Staurosira construens and Achnanthidium minutissimum (Figure 8). In terms of pH, circumneutral taxa dominated. At the beginning of the phase, we observed an increase in alkaliphilous diatoms. The analysis of saprobic preferences revealed the domination of β-mesosaprobous diatoms. In terms of its trophic state, a high proportion of oligo to eutraphentic taxa was observed. The TP reconstruction remained low, between 21.43-57.66 μg L −1 (Figure 11).
Zone II Several short-term fluctuations in species composition were observed. Around 1950, the planktonic taxa (euplanktonic) increased by up to 65%, represented by Asterionella formosa, Aulacoseira ambigua, A. granulata, Diatoma tenuis, Fragilaria nanana, Cyclotella planktonica, Handmania bodanica, and Pantocsekiella comensis (Figure 8). The alkaliphilous species dominated the pH group, with an abundance between 38.7 and 76.5 R. Around 1950, the α-mesosaprobous species increased by 0.3%-17.0%. An increase in eutrophic and hypereutrophic taxa was also observed (Figure 7). The TP reconstruction values increased up to 102.14 μg L −1 (Figure 11).   Zone I This zone was characterized by the dominance of tychoplanktonic taxa, especially Staurosira construens and Achnanthidium minutissimum (Figure 8). In terms of pH, circumneutral taxa dominated. At the beginning of the phase, we observed an increase in alkaliphilous diatoms. The analysis of saprobic preferences revealed the domination of β-mesosaprobous diatoms. In terms of its trophic state, a high proportion of oligo to eutraphentic taxa was observed. The TP reconstruction remained low, between 21.43-57.66 µg L −1 (Figure 11).

Cladocera Analysis
The remains of 44 cladoceran species were found (Figure 9). Throughout the core, the dominant planktonic species were Bosmina (Eubosmina) spp. and Bosmina (Bosmina) longirostris, while among the littoral species, Chydorus cf. sphaericus, Alona quadrangularis and Alonella nana dominated. The overall species composition changes, as well as the indicator species appearance and continuation, allowed us to distinguish between Zones I and II, which correspond to the time before and after the pipe started to operate (Figure 9). The remains of 44 cladoceran species were found (Figure 9). Throughout the core, the dominant planktonic species were Bosmina (Eubosmina) spp. and Bosmina (Bosmina) longirostris, while among the littoral species, Chydorus cf. sphaericus, Alona quadrangularis and Alonella nana dominated. The overall species composition changes, as well as the indicator species appearance and continuation, allowed us to distinguish between Zones I and II, which correspond to the time before and after the pipe started to operate (Figure 9).
Zone I The littoral cladocerans comprised more than 50% of all species, among which Alona affinis, A. quadrangularis, and Alonella nana dominated. The Cladocera eutrophic indicator, which is based on the littoral species, showed the lowest values in the sequence. The rare species Ophryoxus gracilis was present almost continuously. Among the pelagic species, Bosmina (E.) spp. dominated, and Daphnia spp. remains were constantly present. Zone II During this time, several fluctuations in the species share and Cladocera indicator are observable. The pelagic species constitute more that 50% of all species until 2014 CE. However, an increase in the littoral taxa share was noted at the beginning of 1950, 1960, 1980, and ~1995. Next, Chydorus cf. sphaericus becomes the dominant littoral species. After 2014 CE, a clear decrease of pelagic species and an increase of littoral species occurred, and littoral A. rectangula started to dominate. The trophic indicator increased from around 1950 CE, but even more visibly from the ~1978 (Figure 9). The highest values of the indicator species were observed almost continuously from ~2010 until present.

Chironomidae Analysis
In the 45 sediment samples, 52 different chironomid taxa were observed. None of the taxa occurred in all the samples. The Polypedilum nubeculosum-type was the most frequent, as it was observed in 40 samples. The highest maximum abundance values belonged to Procladius (35.3%), P. nubeculosum-type (31.4%), and Heterotanytarsus apicalis (30.0%). The P. nubeculosum-type also had the highest mean abundance (9.3%), followed by the Psectrocladius sordidellus-type (7.9%) and the Tanytarsus pallidicornis-type (6.6%) (Figure 10). The hypolimnetic oxygen reconstruction showed values between 0.2 and 15.7 mg L −1 (Figure 11). Zone I The littoral cladocerans comprised more than 50% of all species, among which Alona affinis, A. quadrangularis, and Alonella nana dominated. The Cladocera eutrophic indicator, which is based on the littoral species, showed the lowest values in the sequence. The rare species Ophryoxus gracilis was present almost continuously. Among the pelagic species, Bosmina (E.) spp. dominated, and Daphnia spp. remains were constantly present.
Zone II During this time, several fluctuations in the species share and Cladocera indicator are observable. The pelagic species constitute more that 50% of all species until 2014 CE. However, an increase in the littoral taxa share was noted at the beginning of 1950, 1960, 1980, and~1995. Next, Chydorus cf. sphaericus becomes the dominant littoral species. After 2014 CE, a clear decrease of pelagic species and an increase of littoral species occurred, and littoral A. rectangula started to dominate. The trophic indicator increased from around 1950 CE, but even more visibly from the~1978 ( Figure 9). The highest values of the indicator species were observed almost continuously from~2010 until present.

Chironomidae Analysis
In the 45 sediment samples, 52 different chironomid taxa were observed. None of the taxa occurred in all the samples. The Polypedilum nubeculosum-type was the most frequent, as it was observed in 40 samples. The highest maximum abundance values belonged to Procladius (35.3%), P. nubeculosum-type (31.4%), and Heterotanytarsus apicalis (30.0%). The P. nubeculosum-type also had the highest mean abundance (9.3%), followed by the Psectrocladius sordidellus-type (7.9%) and the Tanytarsus pallidicornis-type (6.6%) (Figure 10). The hypolimnetic oxygen reconstruction showed values between 0.2 and 15.7 mg L −1 (Figure 11). Zone I The initial part of the sediment profile (before ~1950 CE) was dominated by oligotrophic taxa, such as H. apicalis and the Ablabesmyia monilis-type. The hypolimnetic oxygen reconstructed values remained high at ~10 mg L -1 .
Zone II In the middle portion of the sediment profile (1950-1990 CE), the P. sordidellus-type, which is a common species, was the most abundant. From 2000 CE onwards, Procladius, which prefers nutrient enriched waters, began to dominate. During the most recent years, the eutrophy indicating the Chironomus plumosus-type and Glyptotendipes pallens-type also significantly increased. From the 1950s until the 1980s, the oxygen values decreased to a level of ~4-8 mg L −1 . The hypolimnetic oxygen reconstructed values remained at ~4 mg L −1 from the late 1980s until 2010, after which the values further decreased to anoxic levels in recent years. Figure 11. Comparison of indicators used to reconstruct trophic state changes: chironomid-inferred hypolimnetic oxygen reconstruction (O2) mg L −1 , diatom-inferred total phosphorus reconstruction (TP) μg L −1 , total diatom flux (frustules × 10 6 /cm 2 /yr −1 ), Cladocera based eutrophy indicators, and water-level changes indicators, such as the euplanktonic diatoms (%) and pelagic (P) to littoral (L) Cladocera species ratio (%). Zone I The initial part of the sediment profile (before~1950 CE) was dominated by oligotrophic taxa, such as H. apicalis and the Ablabesmyia monilis-type. The hypolimnetic oxygen reconstructed values remained high at~10 mg L −1 .
Zone II In the middle portion of the sediment profile (1950-1990 CE), the P. sordidellus-type, which is a common species, was the most abundant. From 2000 CE onwards, Procladius, which prefers nutrient enriched waters, began to dominate. During the most recent years, the eutrophy indicating the Chironomus plumosus-type and Glyptotendipes pallens-type also significantly increased. From the 1950s until the 1980s, the oxygen values decreased to a level of~4-8 mg L −1 . The hypolimnetic oxygen reconstructed values remained at~4 mg L −1 from the late 1980s until 2010, after which the values further decreased to anoxic levels in recent years. Zone I The initial part of the sediment profile (before ~1950 CE) was dominated by oligotrophic taxa, such as H. apicalis and the Ablabesmyia monilis-type. The hypolimnetic oxygen reconstructed values remained high at ~10 mg L -1 .

Discussion
Lake Sekšu exemplifies how environmental changes have occurred, in this case, as a result of artificial water pumping station activities, and also shows how such changes at the trophic level can have a lasting influence on lake functioning. In 1953 CE, the pipe was built and for 12 years (1953-1965) pumped water to Lake Sekšu from the large eutrophic lake, Mazais Baltezers. Our results clearly show that this action persistently changed the lake's ecosystem by causing increased eutrophication.
A major cause of eutrophication is an increase in the nutrient supply, particularly phosphorus [55]. Besides excessive plant production, algal blooms, and anoxia, the lake community structure changes through cascading trophic effects and benthic-pelagic coupling [56,57]. Another stressor controlling ecosystem functioning is water level fluctuations, which are essential for aquatic-terrestrial shared boundary processes. The physical and biological effects are especially pronounced in the littoral zone and in shallow water ecosystems in general (e.g., erosion, sedimentation, habitat alterations, and biota changes) [58]. Furthermore, erosion is a source of enhanced nutrient supply, and a greater sediment load decreases water quality. An increase of suspended solids reduces water transparency and increases turbidity, which negatively affects the photosynthesis processes but may also reduce the density of fish and invertebrates [56,59].
The Time Before~1950 Before the pipe started to operate in 1953 CE, Lake Sekšu was a low productivity lake. However, the Cladocera community composition found in our research and in a study by Kuptsch [60] shows that already at the beginning of our core (~1935), the lake trophy had increased compared to 1924. This increase is indicated by a disappearance of the oligotrophy indicating Holopedium gibberum and the appearance of the increasing trophy indicator, Bosmina (Bosmina) longirostris [61][62][63][64][65].
The oldest available investigations from 1924 report a transparency of 5.8 m in the middle of June [60], i.e., almost throughout the lake. It was also reported that Lake Sekšu is one of the 14 lakes located in the vicinity of Latvia's capital city, Riga, that are Lobelia-Isoëtes population dominated lakes [66]. At the beginning of 20 century, there were at least three species typical for Lobelia-Isoëtes lakes found in Lake Sekšu-L. dortmanna, Isoëtes lacustris, and Juncus bulbosus [67]. Our research revealed the presence and high share of lake quillwort pollen Isoëtes lacustris, as well as Cladocera species Ophryoxus gracilis, around 1935 ( Figure 9). Ophyroxus gracilis is a species characteristic of deep, transparent lakes, indicating low/moderate trophy [64,68]. Flössner [49] considers O. gracilis a species typical for Lobelia-Isoëtes lakes.
The diatom-based total phosphorus (TP) reconstruction for Lake Sekšu indicates values between 25-50 µg L −1 before the pipe started to operate (Figure 11), which, according to the OECD (Organisation for Economic Cooperation and Development) classification, indicates a meso-to eutrophic type lake [69]. However, this reconstruction seems to overestimate the TP values when considering the other paleoindicators. This might be due to weak diatom assemblage analogues [70][71][72][73]. Nevertheless, the diatom-inference model closely reflected a significant trend in the measured TP, namely, declining TP values before pipe installation and increasing TP values after pipe installation.
The low trophy status is confirmed by the high share of oligotraphentic diatoms and dominance of Chironomidae oligotrophic taxa [74,75]. The well oxygenated hypolimnetic water (values 9-16 mg L −1 , Figure 11) was inferred by chironomid-based reconstruction and supports a low trophic status. The chironomid-inferred oxygen values remained at a high level (>8 mg L −1 ) from the beginning of the record until the mid-1940s, when a step change to a lower level (4-8 mg L −1 between the 1950s and early 1980s) occurred. The relative oxygen saturation exceeded 50% [56]. Lake Sekšu was slightly acidic to neutral at that time, indicated by the dominance of acidophilous and circumneutral diatoms [76][77][78][79].
The geochemical results support the conclusion of low trophy before the construction of the pipe. However, at that time, the sediment was dominated by organic matter (OM) (Figure 4) with a very high concentration of carbon compared to other lakes in the region [80][81][82][83]. The high carbon/nitrogen (C/N) values suggest a significant contribution of terrigenous OM. Atomic C/N values of phytoplankton vary between 4 and 10, while cellulose-rich vascular plants have C/N ratios of more than 20 [84,85].
Therefore, high C/N rates are expected for oligotrophic lakes with low autochthonous production, and the values observed here are consistent with those of the densely forested catchment of Lake Sekšu. The low water levels of this period due to water pumping, and the subsequent water level decreases as early as 1930, may also have led to the concentration of OM in the coring location due to erosion from exposed organic rich littoral sediments from the lake. The water level in 1930 was reported to be very shallow at 1 m a.s.l [86]. The dominance of littoral over planktonic Cladocera confirms this information [87]. In addition, we reconstructed the high hypolimnetic water's oxygen values (Figure 11), which suggest that the sediment was exposed for efficient bioturbation. Water mixing is also confirmed by the dominance of tychoplanktonic diatoms, i.e., random planktonic diatoms, which require turbulent waters to remain suspended within the photic zone [77,88,89].
The Time Between~1950-2018 The Lake Sekšu water level was artificially increased by pumping water from Lake Mazais Baltezers starting in 1953 CE in order to elevate the ground water level. Apart from the physical changes caused by an increased water level, this action also caused significant changes to the lake's ecology and functioning. The lake became eutrophic according to the diatom-inferred TP values around 100 µg L −1 [69]. During the period of artificial water pumping, until 1965, a high P concentration existed at the water column. After the pipe operation was terminated, the TP decreased to about 70 µg L −1 until the present, indicating that elevated trophy continued. According to our results, the additional nutrient load increased the presence of planktonic eutrophic-hypereutrophic diatoms (Figure 7), which, in turn, decreased water transparency. For instance, in August 2013, the Secchi depth was 2.1 m [15], while in 2019, in April, June, and October, it was 1.6 m, 1.6 m, and 1 m, respectively [90].
Further evidence of the rapid shift towards higher trophy during artificial water pumping is provided by the other biological and geochemical results. The Cladocera eutrophic indicator, based on the species living in the littoral zone (Alona rectangula and Chydorus cf. sphaericus), increased, thereby indicating eutrophic conditions until the most recent sample ( Figure 11). Simultaneously, the higher trophy in the pelagic zone is indicated by the increase of the Bosmina (B.) longirostris share. This species is associated with higher concentrations of TP [65,94].
The sudden change in the relative minerogenic matter and OM input coincides with the implementation of the pipeline and water pumping from Lake Mazais Baltezers to Lake Sekšu. The increased MM content of around 10% indicates substantial changes in the sediment sources. The pipe was installed tens of meters from today's lake shore, and flowing waters likely eroded the substrate-sandy soil that is prone to erosion. A comparable increase occurred in carbonate content, which also suggests intensified catchment erosion. The coeval decrease in OM accumulation, as well as carbon and nitrogen content, is likely a result of dilution by an excess MM supply (Figure 4). Despite the prominent relative changes in sediment composition, no evidence of significant increases in the total sedimentation rate are observable. This could be explained by a sudden change from the erosion of previous littoral sediments at a low-water level stage to channel erosion of the ditch, as well as erosion higher on the shore where no former lake sediments exist. Another indicator of the highly dynamic environment after the pipe started to operate is the increased share of corroded pollen grains [95]. The erosion of the shores is also indicated by the presence of herbivores (via coprophilous fungal spores), Fungi hyphae remains, and Glomus spores [54,96].
Pumping caused a significant elevation of the water level from 1 m a.s.l up to 4.5 m a.s.l [86], which is also reflected in the paleolimnological results. The higher water stand is expressed by the larger share of euplanktonic diatoms [97,98] and the dominance of pelagic over littoral Cladocera species [87] (Figures 9 and 11). However, the fluctuation of water level during the whole 12-year period can be deduced from the different share of the Cladocera pelagic/littoral taxa [87]. Both the increasing and decreasing water levels are related to increased catchment erosion [99][100][101]. The variable MM and carbonate matter (CM) content during the period of external water pumping indicate that fluctuating water levels enhanced catchment erosion in addition to channel erosion of the ditch from the pipe-end to the lake.
The C/N ratio (Figure 4) began to decrease soon after the major sedimentary change occurred, which suggests a steadily increasing proportion of phytoplankton in the total organic matter up to the present day. This is consistent with the gradual eutrophication of the lake, as indicated by the microfossil data. The sediment sulfur (S) profile also provides indirect evidence of eutrophication. Sulfur in lake sediments may be derived from natural processes, such as the weathering of sulfur-bearing rocks and the oxidation of organic sulfur in the catchment [102]. However, changes in S contents in recently deposited lake sediments generally indicate the balance between anthropogenic inputs and the rate of sulfate reduction leading to iron sulfide precipitation in sediments [103]. The steady increase in sedimentary S towards the present day suggests an additional input of sulfate from Lake Mazais Baltezers during the period of pumping, followed by enhanced iron sulfide precipitation in anoxic, OM-rich sediments during subsequent eutrophication. Eutrophication due to sewage water input from Lake Mazais Baltezers was suggested by Leinerte [10] (p.2) and is further supported by the high abundance of α-mesosaprobous diatoms [104,105]. This diatom group flourishes under the lower oxygen saturation in the water column [44].
In the middle part of the sediment record, representing the early 1980s, a threshold change can be observed in the chironomid-inferred hypolimnetic oxygen values (Figure 11). At this point, the values decreased from the previous level of 4-8 mg L −1 (between~1945 and 1980) to a constant level of 4 mg L −1 , which prevailed until the 2010s. This is supported by the rapid decrease in hypolimnetic oxygen levels reconstructed from the chironomid data and recent measurements in April and June 2019, when the hypolimnetic oxygen concentration was 4.0 mg L −1 and 0.1 mg L −1 , respectively [90]. Jansons [9] (p. 27) and Zarina [15] (p.41) reported an amount of oxygen close to zero or zero in the hypolimnion of the deepest part of the basin, which agrees well with our reconstruction. Oxygen depletion was likely also favored by the water depth increase from <3 meters to 8 meters. Deeper waters reduced wind-driven mixing and enabled efficient thermal stratification.
From around 2000 CE, the water level lowered, reaching 2.5 m a.s.l. at present. This might have increased sedimentary nutrient release, further favoring eutrophication [58]. This higher productivity is indicated by the increase of the Chironomidae species Procladius, which prefers nutrient-rich waters [106] (Figure 10), but is also evidenced by the higher values of diatom-inferred TP ( Figure 11). The lowering of the water level inferred from the Cladocera pelagic/littoral species ratio was also confirmed by the local authority [107]. The water level decrease and water trophy increase resulted in oxygen depletion in the hypolimnion, which was revealed by the chironomid reconstruction and recent measurements. Another threshold change based on the chironomid reconstruction occurred during the most recent decade, when the values decreased from a level of 4 mg L −1 (between thẽ 1940s and 2010s) to <2 mg L −1 (Figure 11). Such a decrease in oxygen concentration is ecologically highly significant for chironomid larvae [108].
The lake continues to be significantly overgrown with water plants, which is evidenced by the slight increase in the abundance of the Cladocera phytophilous species, Alonella exigua, and A. excisa [49,109,110]. More sand and mud associated species (Rhynchotalona falcata, Leydigia spp., and Ilyocryptus sp.) have joined the existing community since 1924 [58]. The lake surface pH in 2019 varied within a season from 7.9 to 8.4 [88], which is consistent with previous studies [9,15].
Due to inappropriate water resource management, many lakes in Riga and its vicinity have lost their high ecological quality. The pumping of water from the eutrophic Lake Mazais Baltezers into the low trophic Lake Sekšu increased nutrient inflow and pushed the lake over the threshold, causing eutrophication. Presently, Lake Sekšu continues to be overgrown with macrophytes such as Phragmites australis, Nuphar sp., and Nymphaea sp. [15]. Phytoplankton blooms can be observed in the lake. In 1995/1996, the majority of the cyanobacterial mass contained toxic species [111], while in 2002, Balode et al. [112] found no toxic species.
The "Baltezers" water pumping station is still in use. However, its contribution to the drinking water supply is not as significant as it was in previous years because other pumping stations were established, and additional resources for water supply were found [8].

Conclusions
The pumping of water from the eutrophic lake changed the water level and contributed to faster eutrophication in Lake Sekšu. Before the pipe began operations, the lake was shallow with low productivity and could be defined as a Lobelia-Isoëtes population-dominated biotope. The core studied represents the time period~1935-2018. In 1953, the water level was artificially increased by pumping water from another lake. The sudden nutrient load increased the abundance of planktonic eutrophic-hypereutrophic diatoms, causing the lake water transparency to decrease and inducing hypolimnetic anoxia. The deterioration of the lake's ecological state, as indicated by the results of the biological analysis, was also confirmed by the geochemical results. Cladocera species indicating low trophy became suppressed or disappeared.
Reducing external nutrient loading is known as an important approach to improve the environmental state in small shallow lakes suffering from eutrophication. Nevertheless, in Lake Sekšu, after the water replenishment activities were terminated, the disturbance effect caused by the pumping continued to occur over decades, slowly changing assemblages of keystone species and the foodweb. This highlights the long legacy effects of eutrophication in lacustrine systems. Artificial restoration measures may, therefore, be considered as a strategy to accelerate the recovery of Lake Sekšu.