Next Article in Journal
Impact of the October 2018 Storm Vaia on Coastal Boulders in the Northern Adriatic Sea
Next Article in Special Issue
North German Lowland Lakes Miss Ecological Water Quality Standards—A Lake Type Specific Analysis
Previous Article in Journal
A Rapid Assessment Method to Identify Potential Groundwater Flooding Hotspots as Sea Levels Rise in Coastal Cities
Previous Article in Special Issue
The Effect of Hydrological Connectivity on the Zooplankton Structure in Floodplain Lakes of a Regulated Large River (the Lower Vistula, Poland)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Quality Assessment of a Meromictic Lake Based on Physicochemical Parameters and Strontium Isotopes (87Sr/86Sr) Analysis: A Case Study of Lubińskie Lake (Western Poland)

by
Mateusz Zieliński
1,
Anna Szczucińska
2,
Mateusz Drożdżyński
3,
Marcin Frankowski
4 and
Andrzej Pukacz
5,*
1
Institute of Geoecology and Geoinformation, Adam Mickiewicz University, Bogumiła Krygowskiego 10, Pl-61-680 Poznań, Poland
2
Institute of Physical Geography and Environmental Planning, Adam Mickiewicz University, Bogumiła Krygowskiego 10, 61-680 Poznań, Poland
3
ul. Zbietka 8, PL-62-290 Mieścisko, Poland
4
Faculty of Chemistry, Adam Mickiewicz University Poznań, Poland, Uniwersytetu Poznańskiego 8, PL-61-614 Poznań, Poland
5
Collegium Polonicum, Adam Mickiewicz University, Kościuszki 1, PL-69-100 Słubice, Poland
*
Author to whom correspondence should be addressed.
Water 2019, 11(11), 2231; https://doi.org/10.3390/w11112231
Submission received: 19 September 2019 / Revised: 21 October 2019 / Accepted: 21 October 2019 / Published: 25 October 2019
(This article belongs to the Special Issue Water Quality of Freshwater Ecosystems in a Temperate Climate)

Abstract

:
In 2017, hydrochemical surveys of meromictic Lubińskie Lake (W Poland) and its water inflows were carried out. The lake experienced complete mixing in 2008 due to a series of orkan winds, and since 2015, intensifying worsening of water quality in the lake has been observed. Our aim was to determine the degree of transformation of Lubińskie Lake based on water chemistry and to identify the source of pollution of the lake using strontium isotopes (87Sr/86Sr) as a new chemical tracking tool. The physicochemical analysis confirmed the meromictic character of the lake. The comparison with previous studies (2003 and 2008) showed significant year-to-year differentiation, indicating intensifying eutrophication of the lake’s water, both in the epilimnion and the hypolimnion. Nine spring niches, directly supplying the lake, provide water with very high phosphorus and nitrogen concentrations (up to 10 kg of nitrogen and 0.9 kg of phosphorus daily). The strontium isotopes (87Sr/86Sr) analysis indicated that the lake’s water was supplied mostly by the springs, and recharge from deep aquifers is of secondary importance. Moreover, strontium isotope data and the relationship between Sr and Cl content support the finding that the high load of nutrients is of anthropogenic origin and reaches the lake through springs.

1. Introduction

Climate changes triggered by human activity have caused profound alterations in lake ecosystems, including imbalances in the water budget, eutrophication, and increased sediment resuspension, which lead to the worsening of water quality and the loss of habitats [1,2,3]. These negative changes have affected both large water bodies, such as the Great Lakes or Lake Chad [4,5], and small lakes in remote areas [6,7]. Anthropogenically driven eutrophication is considered as one of the most critical threats for water quality in lakes, resulting in plankton blooms, decreased water clarity, oxygen deficiencies, and finally, loss of habitats [8]. It is mostly induced by an increased supply of nutrients coming from manure and mineral fertilizers applied on crop fields. Nutrients reach the lake through ground water recharge and surface runoff [9]. They can also be delivered to the lake by sublacustrine springs supplying some lakes of early post-glacial areas [10,11]. This applies in particular to the springs whose alimentation area is located within the zone of strong anthropogenic influence [12].
In addition to external supply, an important source of nutrients is also internal loading from sediments [13]. This is particularly important in deep, stratified lakes, in which the decomposition of organic matter under anaerobic conditions causes the release of various substances, such as ammonia, sulfur, or heavy metals, which endanger the quality of water [14]. Thus, stratification and the associated cyclic mixing regime are important for the functioning of lake ecosystem and water quality as well.
The meromictic type is very special and rare with regards to natural lakes. Unlike dimictic lakes (a type of holomixis), in these lakes the deepest stratum of water (monimolimnion) remains isolated from the upper mixing stratum (mixolimnion) by a stratum called chemocline [15]. Hutchinson [16] distinguished three types of meromictic lakes: (i) Ectogenic meromictic lakes with meromixis caused by external events (e.g., the inflow of saline water), (ii) crenogenic meromictic lakes (inflow of saline spring water at a constant, low temperature), and (iii) biogenic meromixis, where differences in water density leading to limited mixing result from the decomposition of organic matter at the bottom. In all these types, stratification may be aided by the morphometric characteristics of the lake as well as shelter from the prevailing wind by surrounding hills or tall vegetation at the shore [9].
Due to permanent density stratification within the monimolimnion, there is no direct supply of oxygen from the atmosphere and the vertical transport of solutes is limited. As a result, sedimentary substances are trapped and accumulated for a long time in anoxic conditions. In holomictic lakes there is a systematic exchange between surface and bottom waters. In meromictic lakes hypolimnion water accumulates huge amounts of ammonia, carbon dioxide, and hydrogen sulfide, which may pose an internal threat to the stability of the ecosystem [17]. Consequently, if the meromixis is disturbed, substances accumulated in the monimolimnion for a long time, including substantial amounts of nutrients, can be rapidly released and cause negative changes within the whole ecosystem. This negative process may be accelerated by external supply from the catchment, e.g., via springs or surface runoffs [13].
Nowadays, water quality assessment of the lakes, mostly based on trophy-related parameters, has been replaced by ecological status assessment, whose main tools are bioindicators [18]. However, the parameters that directly determine water quality (i.e., visibility, nutrient concentration, and chlorine concentration) are still widely used and largely reflect the quality and functioning of the ecosystem as a whole. This is particularly important for assessing the degree and sources of pollution, given the simple and relatively quick analytical methods [19].
Numerous studies show that radiogenic strontium isotopes (87Sr/86Sr) have the potential to fingerprint anthropogenic additives (typically, mineral fertilizers, municipal sewage, and industrial emissions) in aquatic systems [20,21,22,23,24,25,26]. Moreover, strontium isotopes proved to be a reliable tool for the identification of water reservoirs and interactions between water bodies in various aquatic systems [27,28,29,30,31]. They are also a valuable tracer of water-rock interactions and are widely applied in weathering studies [32,33,34,35,36,37].
The strontium isotope method utilizes differences in relative concentration of two strontium isotopes, 87Sr, and 86Sr. As the amount of 87Sr increases with time due to the decay of rubidium isotope 87Rb, rocks of different age develop different 87Sr/86Sr ratios. Strontium is released from rocks by weathering into water and soil and becomes available for further uptake by plants and animals [38]. Since the strontium isotope composition of surface water results from the mixing of strontium derived from various natural (bedrock, soil, precipitation, sea spray) and anthropogenic (industrial emissions, sewage, mineral fertilizers) sources with different isotopic ratios, the 87Sr/86Sr ratio might be applied to trace strontium sources [39,40,41]. Although strontium isotopes are commonly applied in studies of lakes, works investigating meromictic lakes are scarce [42].
In 2004, Pełechaty et al. [43] described a new meromictic lake in mid-Western Poland: Lubińskie Lake. According to Pełechata et al. [44], this lake experienced complete mixing in 2008 due to a series of storm winds. Since then, intensifying worsening of water quality (i.e., increase in nutrient concentration, phytoplankton blooms, decrease of water visibility) in the lake has been observed. Additionally, in 2015, increased fertilization of fields, mainly with liquid manure, occurred within the catchment area.
Our aim was to determine the degree of transformation of Lubińskie Lake based on water chemistry and to identify the source of pollution of the lake using strontium isotopes (87Sr/86Sr) as a new chemical tracking tool. We hypothesized that complete mixing of the water column after storm winds triggered the eutrophication process, which was then reinforced by the inflow of nutrients from the catchment.

2. Study Object

Lubińskie Lake is a small lake with the surface area of 22.7 ha located in western Poland (52°18′32″ N, 14°54′25″ E), in Lubuskie Lakeland (Figure 1). Its catchment area only covers 5.5 km2. It is one of 21 lakes located in the post-glacial tunnel valley running by Gronów, Rzepin, and Torzym towns. The valley spreads NW-SE but in the southern part, its direction changes to W-E, so the basin of the water body has an arched shape. In the southern part of the lake, a bay has formed, separated from the main pool with a 60 m wide isthmus [45]. The maximum depth of the lake is 21.6 m [46,47], whereas the mean depth is 8.8 m [47]. This results in steep slopes of the lake bottom. Detailed parameters and morphometric indices of the lake are presented in Table 1. The west and north-east slopes are particularly steep. The slopes of the bay, however, are gentle, and its depth does not exceed 5 m. The north and west shores of the lake are directly adjacent to the edge of the valley and have the form of steep scarps, over a dozen meters high. The south and north sides are milder.
The bottom of the lake is mostly made of fine- or medium-grained sands, and at some points also of gravels. In the bay and in the north part of the lake, the bottom is mostly covered with gyttjas [45]. The lake basin cuts through alternating deposits of fluvial sands and gravels [48] whose thickness fluctuates between several and over a dozen meters, which facilitate the infiltration of water from the catchment. In the lower part, the basin probably also cuts through the deposits of post-glacial clay (Figure 2), whose thickness in the open pit located approximately 1 km north of the lake is up to 35 m. Above and below the layer of clay there are approximately 10 m thick deposits of fine- and medium-grained sands, whose waters directly feed the lake in the form of groundwater upwellings [48].
As evidenced on an orthophotomap made in 2010 [49], woodlands and farmlands prevail in the catchment, covering the area of 2.7 and 2.6 km2, respectively (Figure 3). The remaining area is covered by waters and the urban territory of Lubin town. The lake is surrounded by forests from all sides, but on the west and south side the forest buffer is very narrow, only between about a dozen and several dozens of meters from the lake shore. Woodlands predominantly occur in the eastern part of the catchment, while arable lands are located in the western part.

3. Materials and Methods

The lake was investigated in April and August 2017, both for comparison with the previous results in 2003 [43] and 2008 (unpublished data) and with regard to the potential meromictic nature of the lake. In both months, the electrolytic conductivity (EC) and pH were measured in the deepest part of the lake, using an Elmetron CX-701 meter (Elmetron Sp.j.). The Secchi disk was used for visibility. Oxygen concentration and temperature were measured in the vertical profile using an Elmetron CX-401 meter. In order to confirm the mixing type, additional measurements for oxygen and temperature were taken at the end of November 2017.
For further chemical analysis, the samples were collected in 1-l plastic bottles (one from the surface and one from the above-bottom layer—using an electric pump) and kept in a portable refrigerator. Then, the samples were kept in a refrigerator (at 4 °C) until the remaining chemical analyses were performed at the Department of Water and Soil Analysis, Faculty of Chemistry, Adam Mickiewicz University in Poznań.
Concentrations of total nitrogen (TN) were determined in the filtered samples by high-temperature combustion (HTC) using a Shimadzu TOC-L Total Organic Carbon analyzer (Shimadzu, Japan). For TN determination, the analyzer was equipped with a TNM-L TN unit (Shimadzu, Japan). The concentrations of NO2, NO3, and NH4+ ions were determined by the ion chromatography method with a conductivity detector (Shimadzu, Japan), using Dionex AS22 with AG22 columns. Total water hardness was determined by the versenate method. In order to determine Ca2+ and Mg2+ ion concentrations, a Shimadzu ion chromatograph with conductivity detection was applied, using Dionex IonPac SCS 1 Silica Cation Separator Column.
The inductively coupled plasma optical emission spectrometer (ICP-OES) ICPE-9820 (Shimadzu, Japan) was used for the qualitative and quantitative detection of total phosphorus (TP). Prior to the analysis, the Sigma-Aldrich (USA) containing 1000 mg L−1 P in 10% nitric acid (comprising hydrofluoric acid traces) was used for the calibration of ICP-OES. In order to preserve the standard/sample conditions, the matrix match method was used.
Soluble reactive phosphorus (SRP) in the water was determined by the molybdate method with ascorbic acid as a reducer. Prior to the analyses, the water samples were filtered through a rinsed 0.45 µm pore size filter (Macherey-Nagel, Germany).
Based on Secchi disc visibility and the concentration of TP in water, Carlson’s trophic state index (TSI) was calculated [50]. The TN and TP concentrations were used for the N:P ratio.
The same methodology as above was applied in the case of the analysis of physicochemical parameters of water in springs feeding the lake, also carried out in April and August 2017. Methods described by Pełechaty et al. [43] were applied with regard to the compared years 2003 and 2008.
The flow and discharge of springs flowing out of the lake were measured three times: On 15 March, 15 April, and 15 August, using an electromagnetic flow meter Valeport model 801 with measurement accuracy of ± 0.5% of reading plus zero stability < 0.005 m s−1.
Separate samples were taken for 87Sr/86Sr analysis of the dissolved strontium load. Additionally, rain water was collected at Słubice, about 30 km to the west, and groundwater was collected at Boczów (about 2.5 km to the north-east) in order to evaluate its contribution to the Sr budget of the lake. A total of 16 samples were taken. Grab samples were taken from the springs. Lake water was sampled with an electric pump from the deepest part of the lake and from the south-eastern bay (Figure 1). In the former, it was gathered at points from the surface up to the depth of 20 m in 4 m intervals, whereas in the latter, samples were collected from the surface and the depth of 4 m. The rain was caught on the roof of Collegium Polonicum of Adam Mickiewicz University to a precleaned rain collector. The obtained sample encompass rainwater from several rain events during two weeks. The groundwater was collected from water intakes. All the samples were taken unfiltered to precleaned polyethylene bottles (250 mL) after flushing the bottles with the sampled water. The samples were filtered through Rotilabo cellulose syringe filters (0.45 µm), acidified with Merck sub-boiled ultrapure HNO3 (1 mL of acid per 100 mL of water), and stored at 4 °C immediately after returning to the laboratory.
The analytical part of this work, including chemical separation of Sr and the measurements of Sr isotope ratios, were carried out in the Isotope Laboratory of the Adam Mickiewicz University in Poznań. About 15 mL of water was evaporated on a hot plate in PFA vials. Subsequently, a 340 µL portion of 2N HNO3 was added to the sample and left overnight to equilibrate. The Sr was separated from matrix elements on PFA columns filled with Sr.Spec resin using the miniaturized chromatographic technique described by Pin et al. [51]. Some modifications in the column size and concentration of reagents were introduced by Dopieralska [52]. Strontium was loaded with a TaCl5 activator on a single W filament and analysed in the dynamic collection mode on a Finnigan MAT 261 mass spectrometer. During the course of this study, the NBS 987 Sr standard yielded 87Sr/86Sr = 0.710210 ± 10 (2σ mean on thirteen analyses). The measured ratios were normalized to the nominal value of 0.710240 for the standard NBS 987. Total procedure blanks were below 83 pg.
Statistical analyses were performed using STATISTICA 12 (StatSoft, Kraków, Poland) software. The normality of distributions of the analyzed variables and the homoscedasticity of the samples were tested with the Shapiro-Wilk and Levene tests, respectively. The significance of differences was tested by ANOVA with the Kruskal-Wallis test. p < 0.05 was accepted as being statistically sound. The principal components analysis (PCA) was performed for the statistical analysis of the physicochemical parameters data set. Prior to this analysis, the data (except pH values) was subject to logarithmic transformation (log (1 + x)) to minimize the discrepancies between the empirical distribution of the variables and their theoretical normal variation. Out of the studied environmental variables, speciation forms of nutrients that were highly interrelated were reduced to TP and TN and included in the PCA. Additionally, to reduce the number of variables, a simple factor analysis was applied. Covariation was tested using non-parametric Spearman rank correlation. The significant levels were set to α = 0.05.

4. Results

4.1. Lake Physio-Chemical Characteristics

The thermal and oxygen profiles (Figure 4) performed during the spring overturn and summer stagnation and repeated during the autumnal overturn indicate the meromictic character of Lubińskie Lake. In all seasons, the monimolimnion zone began at 13 m, which was indicated both by a constant temperature (approx. 4.5 °C) and by anoxia.
The analyzed water properties were typical for eutrophy (TSI = 61), with high nutrients concentration (TN = 3.81 mg N L−1 and TP = 0.13 mg P L−1, on average), high water mineralization (mean EC = 380 µS cm−1), and poor light conditions (mean SD visibility = 2 m; Table 2).
Comparison with the previous studies showed significant (p < 0.05) year-to-year differentiation for three parameters only: Hardness, TN, and NO3. For all of them, a significant increase in values was observed between 2003, 2008, and 2017. A similar pattern of increase was also observed for TP concentration and trophic indicators: TSI and N:P. The results indicate intensifying eutrophy, both in the epilimnion and hypolimnion waters. Interestingly, in the epilimnion, the N:P ratio increases with increasing TSI (phosphorus limitation), while in the hypolimnion, N:P decreases (nitrogen limitation). Although similar differentiation between the surface (epilimnion) and the above-bottom (hypolimnion) layer was observed in all the study years, the results from 2017 are the most divergent, indicating the highest mineralization and trophy.
The factor analysis showed that pH, EC, hardness, Ca2+, Mg2+, TN, and TP were the variables significantly differentiating the physicochemical composition of water within the studied periods. Based on these variables, the PCA analysis was performed (Figure 5). The first two main components explained almost 80% of the observed variance. The PCA output proved pH, hardness, Mg2+, TN, and TP (correlated with the first axis, r > 0.80) as well as EC (correlated with the second axis, r > 0.80) to be primarily responsible for the observed variance.

4.2. Site-to-Site Differentiation of Water Chemistry

The water analyses showed that springs provide water with a significantly lower temperature (from about 10 °C in winter to 12 °C in summer) and much higher mineralization (about 50 µS cm−1 higher in spring and about 100 µS cm−1 in summer) than in the lake’s water (Figure 6). This was also reflected in the values of total hardness and concentrations of Ca2+ and Mg2+. In the case of pH and oxygen, however, the lowest values were found in spring waters and the highest in the outflow waters.
The spring waters were characterized by much higher nutrient concentrations. This is particularly noticeable in the case of TP and SRP as well as TN and NO3-N concentrations (Figure 6). It is worth noting that the waters of the bay directly fed by springs (site P1) had higher values than the waters in site P2. Even higher values were found in outflow waters (site O). Total nitrogen concentrations in spring were almost twice as high (up to 8.2 mg L−1 on average) as in lake waters. The majority of this load was NO3-N, the highest values of which were found in spring S1, both in the spring and the summer (24.13 and 22.86 mg L−1, respectively). The exceptions were: NO2-N and Norg.-N, which was found to be almost non-existent in spring waters, and lower concentrations of NH4-N.

4.3. Water Inflow to the Lake

Groundwater upwellings occur both in the littoral zone and in the pelagic zone (Figure 1). They are fed from the shallow aquifer at the depth of 2 to 5 m below ground level and from the aquifer occurring at the depth below 10 m and isolated from the surface with post-glacial clay deposits. The gravity springs feeding the lake are permanent. In the littoral zone of the lake, nine spring niches were recorded, from which waters are directly supplied to the lake. The flow discharge rate of the four measured springs was in the range of 1 to 8 dm3 s−1 (Table 3). The other flow discharge rates were estimated at below 0.5 dm3 s−1. The total flow discharge rate of all the springs is on average 20 dm3 s−1. In the southern part of the lake, two zones of sublacustrine springs were also recorded. They are approximately 10–20 m from the shoreline. However, it is difficult to measure their flow discharge. It is estimated to be about 50% of the inflow from overground springs, but this value may be underestimated. The discharge rate of the stream flowing out of the lake in the spring, directly after the disappearance of the ice cover was 40 dm3 s−1, and in the summer, 21 dm3 s−1. The difference between the amount of water feeding from the springs and the outflow from the lake indicates additional feeding with groundwater from aquifer outcrops on the bottom of the lake.

4.4. Strontium Isotope Composition in Water

The profiles of 87Sr/86Sr and Sr content for Lubińskie Lake are shown in Figure 7, and Sr isotope and elemental characteristics for components of the hydrologic system of the lake are given in Table 4. Since the bay yields a uniform Sr isotope composition and Sr concentration throughout the water column, it will not be discussed in detail. The lake is characterized by a narrow range of 87Sr/86Sr values, between 0.710741 and 0.710798. Thus, it is isotopically homogeneous throughout the water column (Figure 7a). Similarly, Sr content of lake water is highly uniform in the profile and changes from 0.124 to 0.148 mg L−1 (Figure 7b). It is worth mentioning that the highest concentration of dissolved Sr occurs in the lower part of the water column. As can be expected, the stream draining the lake has the 87Sr/86Sr ratio close to that of surface lake water, but is enriched in Sr (Figure 7).
The springs show a wide range of 87Sr/86Sr values from 0.710412 to 0.711851 and relatively diverse Sr dissolved load (0.127–0.175 mg L−1). The strontium isotope composition of the groundwater falls within the range for the springs (87Sr/86Sr values of 0.710693–0.710906), yet it has a clearly higher Sr content (0.241–0.266 mg L−1). Finally, the rain is the least radiogenic component (with the lowest 87Sr/86Sr ratio) and yields 87Sr/86Sr value of 0.708879.

5. Discussion

5.1. Water Chemistry Analysis

The conducted research showed that Lubińskie Lake is currently in the state of higher trophy, which is indicated by the concentrations of nutrients and low Secchi disk visibility. The comparison with previous years of research shows clearly worsening of the water quality. In 2003, despite high nutrient concentrations, the lake maintained the clean water state (visibility over 4 m) and was classified as a charophyte lake [43,53].
The results of the physicochemical and algological tests showed that the meromictic character of the lake maintained since 2003 was disturbed between January and March 2008 [44]. At the end of February 2008, the waters in the lake mixed completely. At that time, a series of storm winds went across the western part of Poland (“Paula”, “Zizi”, and “Emma”). The complete restoration of the meromictic conditions lasted until June 2008 [54] (Pukacz et al. 2010). Although the complete mixing of lake waters lasted only a few months (which is confirmed by measurements carried out in the subsequent years), it probably contributed to the release of nutrients and other substances (e.g., methane and hydrogen sulfide) deposited for years in the monimolimnion [55]. This resulted in gradual worsening of the condition of the whole ecosystem. One of the visible effects was the worsening of light conditions and gradual disappearance of submerged vegetation (Pukacz 2008–2018, unpublished data).
Although meromictic stratification was restored in the lake, further increase in nutrients concentration was observed. The study showed that an important source of phosphorus and nitrogen (especially in the form NO3-N and NH4-N) supplied to the lake is the spring waters. Based on the average concentrations and water supply from the springs, it can be estimated that each day, approximately 10 kg of nitrogen and 0.9 kg of phosphorus can be introduced to the lake. However, these values may be underestimated by approximately 25%, because they do not take into account the load supplied with discharges for which the volume was estimated, and the load supplied through the springs in the bottom of the lake.
The increased supply of those elements must be associated with the change in the use of land in the catchment area of Lubińskie Lake. Since 2015, in the fields located in the SE part of the catchment, corn, and rapeseed have been cultivated, which are not only fertilized with mineral fertilizers but also, in the spring, with organic fertilizers (own observations). Taking into account the directions of inflow of groundwater and the very high concentrations of nutrients in spring waters, it is a very probable cause of the increased inflow of nutrients to the lake.
Interestingly, the values of the studied parameters in the SE bay, where most springs are located, were not significantly higher than in the central basin. This may result from the small surface area of the lake and from feeding it with inflows located in the slope of the basin. Feeding the lake through springs located in the bottom of the lake basin may have an impact, not only on the lower quality of the lake waters, but also on disturbed stratification and crenogenic meromixis [56].

5.2. Strontium Isotopes Perspective

Strontium isotope composition coupled with Sr content enables the identification of reservoirs controlling the 87Sr/86Sr ratio of water and mixing processes involved [27,36,57]. There is a strong correlation between 87Sr/86Sr and 1/Sr for the springs dataset, which points to mixing between the two Sr reservoirs (Figure 8). The first reservoir is characterized by a relatively radiogenic 87Sr/86Sr ratio and moderately low 1/Sr (high Sr concentration), whereas the second reservoir has a less radiogenic 87Sr/86Sr ratio and higher 1/Sr (lower Sr content). Since the springs emerge from shallow Quaternary aquifers, the first reservoir is represented by silicates, mostly glacial deposits. It controls the isotopic composition of the S1 spring. Since the investigated springs are charged by atmospheric waters [58], we identify the second reservoir with rainwater, which strongly influences the Sr isotope budget of the S2-S4 springs. As shown by Raiber et al. [59], rain can be an important source of Sr to shallow groundwater. Strontium isotope composition of rainwater at Słubice (87Sr/86Sr value of 0.708879) is less radiogenic than that of modern seawater (87Sr/86Sr about 0.70917; [60]) being the source of rainwater, likely due to the addition of atmospheric dust. Inland deviation of precipitation 87Sr/86Sr ratio is a well-documented phenomenon [59,61,62]. Typically, rainwater is not rich in Sr [59,61,63]. Zieliński et al. [25] provide Sr concentration of about 0.001 mg L−1 for rain and snow at Poznań, about 140 km east of Lubińskie Lake. However, the admixture of dust elevates the Sr content [59]. Therefore, Sr dissolved load of tens of µg L−1 for the rain at Słubice can reasonably be assumed.
The role of rain as one of the major sources of Sr for the springs is supported by a clear relationship between Cl and Sr in the springs dataset (Figure 9), indicating their common origin. Chloride ions are a useful atmospheric input reference to the hydrosystem, because in the dissolved load of the hydrosystem they originate from seasalt and saltrock dissolution [64,65].
However, they can also result from various human activities (domestic sewage, fertilizers, and deicing salts) [66]. Since the Lubińskie Lake drainage basin is free of evaporates, they do not participate in the weathering reactions. Consequently, the general trend of relationships between Cl and Sr can be explained by the contribution of Cl from the atmospheric and anthropogenic sources. Since crop fields occur in hinterland of the spring area, the latter is likely represented by mineral fertilizers. This assumption is supported by the high load of nitrates (NO3) carried by the springs (Figure 6). Additionally, sewage leaking from surrounding villages may also provide some load of Cl.
Recently, Zieliński et al. [25] provided evidence for the common lowering of the present-day 87Sr/86Sr values of Noteć River water in western Poland due to the impact of fertilizers used in agriculture. The estimated scale of the lowering was not higher than 0.0006 to 0.0007 of the observed 87Sr/86Sr ratios in the downstream portion of the Noteć River. Given that the input of fertilizers-related Sr to the springs waters was recognized, the lowering of the lake water isotopic ratio should be assumed.
Figure 8 shows that data points for Lubińskie Lake in general plot in the springs domain on the 87Sr/86Sr versus 1/Sr diagram; thus, crucial contribution of spring-derived Sr is identified. Surprisingly, spring S1, the second largest yet richest in Sr, has a negligible impact on the Sr isotopic composition of the lake. Instead, springs S2–S4 with lower flow rates and Sr contents seem to contribute to it the most. This suggests additional input of shallow groundwater with 87Sr/86Sr ratios similar to these springs below the water table due to the lack of other springs along the lake’s shore. Indeed, we found several seepages at the bay’s bottom.
If the isotopic composition of Lubińskie Lake was solely related to the springs-derived Sr, data points for the lake water would lie on paths between the springs. Most of the data points, however, are moved aside the mixing trends, mainly towards lower Sr concentrations but without noticeable changes in 87Sr/86Sr, and evidence the contribution of less concentrated waters. Apart from rainwater supplied to Lubińskie Lake through the springs, as inferred above, it also reaches the lake directly. Therefore, it is likely responsible for the dilution of the lake’s Sr dissolved load.
The input of dust-derived Sr can account for differences in Sr content between Lubińskie Lake and the stream originating from it. Graustein and Armstrong (1983) [39] show that coniferous foliage is more effective in dust trapping than deciduous foliage. As the stream flows in a coniferous forest surrounding the northern fringe of the lake, dust deposited on the foliage and then incorporated into the soil by throughfall seems to charge soil waters and hence the stream by additional Sr.
Although the lake basin cuts deeply into the aquifer (Figure 2), and the recharge of the lake by groundwater arguably occurs, its impact on Sr isotope composition of the lake water seems negligible. Nonetheless, this impression is only apparent, because 87Sr/86Sr ratios of the groundwater fall within the 87Sr/86Sr range for the springs, so the contribution of the former will be masked. The input of groundwater-derived Sr to Lubińskie Lake is suggested by the rise of Sr concentration several meters above the bottom, which is consistent with the lateral input of groundwater enriched in Sr and its downward descend. Due to the meromictic character of the lake, it is, however, limited to its deepest part.

6. Conclusions

The results of the study indicate that Lubińskie Lake has regained its meromictic character after complete mixing in 2008. That sudden event probably led to the release of nutrients from the monimolimnion layer, which resulted in an increase in trophy and deterioration of water quality in subsequent years. An additional, significant factor affecting the lake’s water quality is the nutrient inflow from the catchment. Although the lake is insulated with a belt of trees, water polluted with fertilizers infiltrates into the first aquifer and enters the lake through springs. The difference between the lake supply from springs and the outflow from the lake indicates lake charging with groundwater from underwater outcrops of aquifers. That was also confirmed by strontium isotopes (87Sr/86Sr) analysis. It indicates that water is supplied to Lubińskie Lake mostly by the springs, and recharge from deep aquifers is of secondary importance. Moreover, strontium isotope data and the relationship between Sr and Cl content support finding that high load of nutrients is of anthropogenic origin and reach the lake through springs. We show that the 87Sr/86Sr composition of lake water supplemented by major ion chemistry is a valuable tool in the investigation of meromictic lakes.

Author Contributions

Conceptualization, A.P. and A.S.; Methodology and investigation, M.Z., A.S., M.D., M.F., A.P.; Data analysis, M.Z., A.S., M.D., M.F., A.P.; Writing—original draft preparation, M.Z., A.S., A.P.; Writing—review and editing, M.Z., A.S., A.P.

Funding

The APC was founded by Adam Mickiewicz University and the research received no external funding.

Acknowledgments

We wish to thank Z. Belka (Poznań) for the constructive criticism of the MS as well as comments of two anonymous reviewers and journal editor help to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Evans, D.O.; Nicholls, K.H.; Allen, Y.C.; McMutry, M.J. Historical land use, phosphorus loading, and loss of fish habitat in Lake Simcoe, Canada. Can. J. Fish. Aquat Sci. 1996, 53 (Suppl. 1), 194–218. [Google Scholar] [CrossRef]
  2. Callisto, M.; Molozzi, J.; Barbosa, J.L.E. Eutrophication of lakes. In Eutrophication: Causes, Consequences and Control, 1st ed.; Ansari, A.A., Gill, S.S., Eds.; Springer: Dordecht, The Netherlands, 2014; pp. 55–71. [Google Scholar]
  3. Hijosa-Valsero, M.; Bécares, E.; Fernández-Aláez, C.; Fernández-Aláez, M.; Mayo, R.; Jiménez, J.J. Chemical pollution in inland shallow lakes in the Mediterranean region (NW Spain): PAHs, insecticides and herbicides in water and sediments. Sci. Total Environ. 2016, 544, 797–810. [Google Scholar] [CrossRef] [PubMed]
  4. Gao, H.; Bohn, T.J.; Podest, E.; McDonald, K.C.; Lettenmaier, D.P. On the causes of the shrinking of Lake Chad. Environ. Res. Lett. 2011, 6, 034021. [Google Scholar] [CrossRef] [Green Version]
  5. Gronewold, A.D.; Stow, C.A. Water loss from the Great Lakes. Science 2014, 343, 1084–1085. [Google Scholar] [CrossRef]
  6. Kang, W.; Chen, G.; Wang, J.; Huang, L.; Wang, L.; Li, R.; Hu, K.; Liu, Y.; Tao, J.; Blais, J.M.; et al. Assessing the impact of long-term changes in climate and atmospheric deposition on a shallow alpine lake from southeast Tibet. Sci. Total Environ. 2019, 650, 713–724. [Google Scholar] [CrossRef]
  7. Wan, D.; Mao, X.; Jin, Z.; Song, L.; Yang, J.; Yang, H. Sedimentary biogeochemical record in Lake Gonghai: Implications for recent lake changes in relatively remote areas of China. Sci. Total Environ. 2019, 649, 929–937. [Google Scholar] [CrossRef]
  8. Welch, E.B.; Jacoby, J.M. Pollutant Effects in Freshwater: Applied Limnology, 3rd ed.; CRC Press: London, UK, 2004; pp. 1–520. [Google Scholar]
  9. Wetzel, R.G. Limnology, Lake and River Ecosystems, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2001; pp. 1–1006. [Google Scholar]
  10. Choiński, A.; Ilyin, L.; Marszelewski, W.; Ptak, M. Lakes supplied by springs: Selected examples. Limnol. Rev. 2008, 8, 145–150. [Google Scholar]
  11. Cieśliński, R.; Piekarz, J.; Zieliński, M. Groundwater supply of lakes: The case of Lake Raduńskie Górne (northern Poland, Kashubian Lake District). Hydrol. Sci. J. 2016, 61, 2427–2434. [Google Scholar] [CrossRef]
  12. Scott, J.T.; Marcarelli, A.M. Cyanobacteria in freshwater benthic environments. In The Ecology of Cyanobacteria; Whitton, B.A., Ed.; Springer: Durham, UK, 2012; pp. 271–289. [Google Scholar]
  13. Dodds, W.K. Freshwater ecology. In Concepts and Environmental Applications, 1st ed.; Academic Press: San Diego, CA, USA, 2002; pp. 1–591. [Google Scholar]
  14. Last, W.M.; Smol, J.P. An Introduction to Physical and Geochemical methods Used in Paleolimnology. In Tracking Environmental Change Using Lake Sediments; Last, W.M., Smol, J.P., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2001; pp. 1–504. [Google Scholar]
  15. Lampert, W.; Sommer, U. Ecology of Inland Waters; PWN Scientific Press: Warsaw, Poland, 2001; pp. 1–415. [Google Scholar]
  16. Hutchinson, G.E. A contribution to the limnology of arid regions: Primarily founded on observations made in the Lahontan basin. Trans. Conn. Acad. Arts Sci. 1937, 33, 47–132. [Google Scholar]
  17. Göltenboth, F.; Timotius, K.H.; Milan, P.P.; Margraf, J. Ecology of Insular Southeast Asia; Elsevier: Amsterdam, The Netherlands, 2006; pp. 1–568. [Google Scholar]
  18. Bis, B. Assessing the Ecological Status Assessment of Freshwaters. In Freshwater Ecosystems in Europe–An Educational Approach; Voreadou, C., Ed.; Natural History Museum of Crete, Selena Press: Heraklion, Greece, 2008; pp. 56–69. [Google Scholar]
  19. Imteaz, M.A.; Asaeda, T.; Lockington, D.A. Modelling the effects of inflow parameters on lake water quality. Environ. Model. Assess. 2003, 8, 63–70. [Google Scholar] [CrossRef]
  20. Böhlke, J.K.; Horan, M. Strontium isotope geochemistry of groundwaters and streams affected by agriculture, Locust Grove, MD. Appl. Geochem. 2000, 15, 599–609. [Google Scholar] [CrossRef]
  21. Hosono, T.; Nakano, T.; Igeta, A.; Tayasu, I.; Tanaka, T.; Yachi, S. Impact of fertilizer on a small watershed of Lake Biwa: Use of sulfur and strontium in environmental diagnosis. Sci. Total Environ. 2007, 384, 342–354. [Google Scholar] [CrossRef] [PubMed]
  22. Chapman, E.C.; Capo, R.C.; Stewart, B.W.; Hedin, R.S.; Weaver, T.J.; Edenborn, H.M. Strontium isotope quantification of siderite, brine and acid mine drainage contributions to abandoned gas well discharges in the Appalachian Plateau. Appl. Geochem. 2013, 31, 109–118. [Google Scholar] [CrossRef]
  23. Vengosh, A.; Jackson, R.B.; Warner, N.; Darrah, T.H.; Kondash, A. A critical review of the risk to water resources from unconventional shale gas development and hydraulic fracturing in the United States. Environ. Sci. Technol. 2014, 48, 8334–8348. [Google Scholar] [CrossRef] [PubMed]
  24. Sherman, L.S.; Blum, J.D.; Dvonch, J.T.; Gratz, L.E.; Landis, M.S. The use of Pb, Sr, and Hg isotopes in Great Lakes precipitation as a tool for pollution source attribution. Sci. Total Environ. 2015, 502, 362–374. [Google Scholar] [CrossRef]
  25. Zieliński, M.; Dopieralska, J.; Belka, Z.; Walczak, A.; Siepak, M.; Jakubowicz, M. Sr isotope tracing of multiple water sources in a complex river system, Noteć River, central Poland. Sci. Total Environ. 2016, 548, 307–316. [Google Scholar] [CrossRef]
  26. Zieliński, M.; Dopieralska, J.; Belka, Z.; Walczak, A.; Siepak, M.; Jakubowicz, M. Strontium isotope identification of water mixing and recharge sources in a river system (Oder River, central Europe): A quantitative approach. Hydrol. Process. 2018, 32, 2597–2611. [Google Scholar] [CrossRef]
  27. Négrel, P.; Fouillac, C.; Brach, M. A strontium isotopic study of mineral and surface waters from the Cézallier (Massif Central, France): Implications for mixing processes in areas of disseminated emergences of mineral waters. Chem. Geol. 1997, 135, 89–101. [Google Scholar] [CrossRef]
  28. Négrel, P.; Casanova, J.; Aranyossy, J.P. Strontium isotope systematics used to decipher the origin of groundwaters sampled from granitoids: The Vienne Case (France). Chem. Geol. 2001, 177, 287–308. [Google Scholar] [CrossRef]
  29. Négrel, P.; Petelet-Giraud, E.; Barbier, J.; Gautier, E. Surface water-groundwater interactions in an alluvial plain: Chemical and isotopic systematics. J. Hydrol. 2003, 277, 248–267. [Google Scholar] [CrossRef]
  30. Ojiambo, S.B.; Lyons, W.B.; Welch, K.A.; Poreda, R.J.; Johannesson, K.H. Strontium isotopes and rare earth elements as tracers of groundwater-lake water interactions, Lake Naivasha, Kenya. Appl. Geochem. 2003, 18, 1789–1805. [Google Scholar] [CrossRef]
  31. Flusche, M.A.; Seltzer, G.; Rodbell, D.; Siegel, D.; Samson, S. Constraining water sources and hydrologic processes from the isotopic analysis of water and dissolved strontium, Lake Junin, Peru. J. Hydrol. 2005, 312, 1–13. [Google Scholar] [CrossRef]
  32. Blum, J.D.; Gazis, C.A.; Jacobson, A.D.; Chamberlain, C.P. Carbonate versus silicate weathering in the Raikhot watershed within the High Himalayan Crystalline Series. Geology 1998, 26, 411–414. [Google Scholar] [CrossRef]
  33. Jacobson, A.D.; Blum, J.D.; Walter, L.M. Reconciling the elemental and Sr isotope composition of Himalayan weathering fluxes: Insight from the carbonate geochemistry of stream water. Geochim. Cosmochim. Acta 2002, 66, 3417–3429. [Google Scholar] [CrossRef]
  34. Jin, Z.; Wang, S.; Zhang, F.; Shi, Y. Weathering, Sr fluxes, and controls on water chemistry in the Lake Qinghai catchment, NE Tibetan Plateau. Earth Surf. Process. 2010, 35, 1057–1070. [Google Scholar] [CrossRef]
  35. Lyons, W.B.; Welch, K.A.; Priscu, J.C.; Tranter, M.; Royston-Bishop, G. Source of Lake Vostok cations constrained with strontium isotopes. Front. Earth Sci. 2016, 4, 78. [Google Scholar] [CrossRef]
  36. Zieliński, M.; Dopieralska, J.; Belka, Z.; Walczak, A.; Siepak, M.; Jakubowicz, M. The strontium isotope budget of the Warta River (Poland): Between silicate and carbonate weathering, and anthropogenic pressure. Appl. Geochem. 2017, 81, 1–11. [Google Scholar] [CrossRef]
  37. Négrel, P.; Pauwels, H.; Chabaux, F. Characterizing multiple water-rock interactions in the critical zone through Sr-isotope tracing of surface and groundwater. Appl. Geochem. 2018, 93, 102–112. [Google Scholar] [CrossRef]
  38. Capo, R.C.; Stewart, B.W.; Chadwick, O.A. Strontium isotopes as tracers of ecosystem processes: Theory and methods. Geoderma 1998, 82, 197–225. [Google Scholar] [CrossRef]
  39. Graustein, W.C.; Armstrong, R.L. The use of strontium-87/strontium-86 ratios to measure atmospheric transport into forested watersheds. Science 1983, 219, 289–292. [Google Scholar] [CrossRef]
  40. Åberg, G. The use of natural strontium isotopes as tracers in environmental studies. Water Air Soil Poll. 1995, 79, 309–322. [Google Scholar] [CrossRef]
  41. Bain, D.C.; Midwood, A.J.; Miller, J.D. Strontium isotope ratios in streams and the effect of flow rate in relation to weathering in catchments. Catena 1998, 32, 143–151. [Google Scholar] [CrossRef]
  42. Faure, G.; Mensing, T.M. The Transantarctic Mountains. Rocks, Ice, Meteorites and Water; Springer: Dordrecht, The Netherlands, 2010; pp. 1–832. [Google Scholar]
  43. Pełechaty, M.; Pukacz, A.; Pełechata, A. Diversity of micro- and macrophyte communities in the context of the habitat conditions of a meromictic lake on Lubuskie Lakeland. Limnol. Rev. 2004, 4, 209–214. [Google Scholar]
  44. Pełechata, A.; Pełechaty, M.; Pukacz, A. Complete mixing of a meromictic lake as a result of extreme weather events and its impact on phytoplankton community. In Proceedings of the 23rd Congres of Polish Hydrobiologists, Koszalin, Poland, 8–12 September 2015. [Google Scholar]
  45. Zieleniewski, W. Fisherman report for waters of circumference I.5.8. In The Fishing Region of the Lubińskie Lake on an Unnamed Watercourse Flowing into the Rzepia River-No. 1 in the Ilanka River Basin in the Lower Oder and Western Catchment Water Region; Polish Fishermans Association: Zielona Góra, Poland, 2007. [Google Scholar]
  46. Choiński, A. Outline of Polish Physical Limnology; Adam Mickiewicz University Press: Poznań, Poland, 1995; pp. 1–298. (In Polish) [Google Scholar]
  47. Jańczak, J. Atlas of Polish Lakes; Bogucki Press: Poznań, Poland, 1996; pp. 1–256. [Google Scholar]
  48. Stupnicka, E. Regional Geology of Poland; University of Warsaw Press: Warsaw, Poland, 2007; pp. 1–346. (In Polish) [Google Scholar]
  49. Main Office of Geodesy and Cartography. Available online: http://www.gugik.gov.pl/ (accessed on 15 September 2019).
  50. Carlson, R.E. A trophic state index for lakes. Limnol. Oceanogr. 1977, 22, 361–369. [Google Scholar] [CrossRef]
  51. Pin, C.; Briot, D.; Bassin, C.; Poitrasson, F. Concomitant separation of strontium and samarium-neodymium for isotopic analysis in silicate samples, based on specific extraction chromatography. Anal. Chim. Acta 1994, 298, 209–222. [Google Scholar] [CrossRef]
  52. Dopieralska, J. Neodymium Isotopic Composition of Conodonts as a Palaeoceanographic Proxy in the Variscan Oceanic System. Ph.D. Thesis, Justus-Liebig-University, Giessen, Germany, 2003. [Google Scholar]
  53. Pełechata, A.; Pełechaty, M.; Pukacz, A. An attempt to the trophic status assessment of the lakes of Lubuskie Lakeland. Limnol. Rev. 2006, 6, 239–246. [Google Scholar]
  54. Pukacz, A.; Pełechaty, M.; Pełechata, A.; Sękowska, N. The influence of winter windstorms on the disorder of meromixia of Lubińskie Lake. In Proceedings of the 14th Polish Limnologic Conference Natural and Anthropogenic Transformations of the Lakes, Szczecin-Stare Drawsko, Poland, 20–23 September 2010. (In Polish). [Google Scholar]
  55. Boehrer, B.; von Rohden, C.; Schultze, M. Physical features of Meromictic Lakes: Stratification and circulation. In Ecology of Meromictic Lakes; Gulati, R.D., Zadereev, E., Degermendzhy, A.G., Eds.; Springer: Dordrecht, The Netherlands, 2017; pp. 1–405. [Google Scholar]
  56. Hakala, A. Meromixis as a part of lake evolution—Observations and a revised classification of true meromictic lakes in Finland. Boreal Environ. Res. 2004, 9, 37–53. [Google Scholar]
  57. Shand, P.; Darbyshire, D.P.F.; Love, A.J.; Edmunds, W.M. Sr isotopes in natural waters: Applications to source characterisation and water-rock interaction in contrasting landscapes. Appl. Geochem. 2009, 24, 574–586. [Google Scholar] [CrossRef]
  58. Szczucińska, A. Spatial distribution and hydrochemistry of springs and seepage springs in the Lubuska Upland of western Poland. Hydrol. Res. 2014, 45, 379–390. [Google Scholar] [CrossRef]
  59. Raiber, M.; Webb, J.A.; Bennetts, D.A. Strontium isotopes as tracers to delineate aquifer interactions and the influence of rainfall in the basalt plains of southeastern Australia. J. Hydrol. 2009, 367, 188–199. [Google Scholar] [CrossRef]
  60. McArthur, J.M.; Howarth, R.J.; Bailey, T.R. Strontium isotope stratigraphy: LOWESS version 3: Best fit to the marine Sr-isotope curve for 0-509 Ma and accompanying look-up table for deriving numerical age. J. Geol. 2001, 109, 155–170. [Google Scholar] [CrossRef]
  61. Négrel, P.; Guerrot, C.; Millot, R. Chemical and strontium isotope characterization of rainwater in France: Influence of sources and hydrogeochemical implications. Isot. Environ. Health S 2007, 43, 179–196. [Google Scholar] [CrossRef] [PubMed]
  62. Xu, Z.; Han, G. Chemical and strontium isotope characterization of rainwater in Beijing, China. Atmos. Environ. 2009, 43, 1954–1961. [Google Scholar] [CrossRef]
  63. Nakano, T.; Morohashi, S.; Yasuda, H.; Sakai, M.; Aizawa, S.; Shichi, K.; Morisawa, T.; Takahashi, M.; Sanada, M.; et al. Determination of seasonal and regional variation in the provenance of dissolved cations in rain in Japan based on Sr and Pb isotopes. Atmos. Environ. 2006, 40, 7409–7420. [Google Scholar] [CrossRef]
  64. Stallard, R.F.; Edmond, J.M. Geochemistry of the Amazon, 1. Precipitation chemistry and the marine contribution to the dissolved load at the time of peak discharge. J. Geophys. Res. 1981, 86, 9844–9858. [Google Scholar] [CrossRef]
  65. Meybeck, M. Atmospheric inputs and river transport of dissolved substances. In Proceedings of the International Association of Hydrological Sciences Conference Dissolved load of rivers and surface water quantity/quality relationships, Hamburg, Germany, 15–27 August 1983; pp. 173–192. [Google Scholar]
  66. Sherwood, W.C. Chloride loading in the South Fork of the Shenandoah River, Virginia, U.S.A. Environ. Geol. Water S 1989, 14, 99–106. [Google Scholar] [CrossRef]
Figure 1. Location of Lubińskie Lake and the points of water samples intake from: The lake (W1 and W2), outflow (O), springs (S1–S4), and wells P2 and P4 (A,C). The line A–A’ shows the line of hydrogeological cross-section (B).
Figure 1. Location of Lubińskie Lake and the points of water samples intake from: The lake (W1 and W2), outflow (O), springs (S1–S4), and wells P2 and P4 (A,C). The line A–A’ shows the line of hydrogeological cross-section (B).
Water 11 02231 g001
Figure 2. Hydrogeological cross-section of Lake Lubińskie drainage basin along the line Boczów—Lubińskie Lake—Mierczany.
Figure 2. Hydrogeological cross-section of Lake Lubińskie drainage basin along the line Boczów—Lubińskie Lake—Mierczany.
Water 11 02231 g002
Figure 3. Structure of land use in the catchment of Lubińskie Lake based on the orthophotomap.
Figure 3. Structure of land use in the catchment of Lubińskie Lake based on the orthophotomap.
Water 11 02231 g003
Figure 4. Vertical variability of temperature and oxygen concentration in Lubińskie Lake studied in 2017.
Figure 4. Vertical variability of temperature and oxygen concentration in Lubińskie Lake studied in 2017.
Water 11 02231 g004
Figure 5. The principal components analysis (PCA) for water analyses in Lubińskie Lake performed in 2003 [43], 2008 (unpublished data), and 2017—differentiation across the selected hydrochemical characteristics (preceded by factor analysis). Gray color indicates samples taken in August, and black—those taken in April.
Figure 5. The principal components analysis (PCA) for water analyses in Lubińskie Lake performed in 2003 [43], 2008 (unpublished data), and 2017—differentiation across the selected hydrochemical characteristics (preceded by factor analysis). Gray color indicates samples taken in August, and black—those taken in April.
Water 11 02231 g005
Figure 6. Differentiation of water properties between the studied sites. P1—pelagial 1, P2—pelagial 2, O—outflow, S—spring, * — mean values for all springs. Green bars—April, yellow bars—August.
Figure 6. Differentiation of water properties between the studied sites. P1—pelagial 1, P2—pelagial 2, O—outflow, S—spring, * — mean values for all springs. Green bars—April, yellow bars—August.
Water 11 02231 g006
Figure 7. Depth profiles for (a) 87Sr/86Sr and (b) Sr content of the lake water. Data for the springs, the stream, and the groundwater is also shown for comparison. Explanations for panel (a) also refer to panel (b).
Figure 7. Depth profiles for (a) 87Sr/86Sr and (b) Sr content of the lake water. Data for the springs, the stream, and the groundwater is also shown for comparison. Explanations for panel (a) also refer to panel (b).
Water 11 02231 g007
Figure 8. 87Sr/86Sr versus 1/Sr mixing diagram of the Lubińskie Lake hydrologic system. Note a binary mixing trend for the springs dataset and the deviation of data points for lake water from mixing paths between the springs.
Figure 8. 87Sr/86Sr versus 1/Sr mixing diagram of the Lubińskie Lake hydrologic system. Note a binary mixing trend for the springs dataset and the deviation of data points for lake water from mixing paths between the springs.
Water 11 02231 g008
Figure 9. Relationship between Cl and Sr in samples collected from the springs and Lubińskie Lake. Chlorine content in the lake was measured only for surface and near bottom waters (see Table 4).
Figure 9. Relationship between Cl and Sr in samples collected from the springs and Lubińskie Lake. Chlorine content in the lake was measured only for surface and near bottom waters (see Table 4).
Water 11 02231 g009
Table 1. Parameters and morphometric indices of Lubińskie Lake (according to Jańczak [47]).
Table 1. Parameters and morphometric indices of Lubińskie Lake (according to Jańczak [47]).
Parameters and IndicesUnitValue
Surface areaha22.7
Maximum lengthm1130
Maximum widthm270
Mean widthm200
Length ratio-4.2
Shoreline lengthm2950
Shoreline development-1.75
Volume103 m31997.6
Maximum depthm21.6
Mean depthm8.8
Relative depth *-0.0453
Depth ratio-0.41
Exposure ratio-2.6
* Maximum depth (m)/surface area (km2).
Table 2. Physicochemical characteristic of Lubińskie Lake epilimnion (E) and hypolimnion (H) water over three years of study: 2003 [43], 2008 (unpublished data) and 2017. The mean value ± standard deviation is given for each year.
Table 2. Physicochemical characteristic of Lubińskie Lake epilimnion (E) and hypolimnion (H) water over three years of study: 2003 [43], 2008 (unpublished data) and 2017. The mean value ± standard deviation is given for each year.
Parameter(unit)200320082017
EHEHEH
SD(m)3.2 ± 0.5 2.5 ± 0.9 2 ± 0.1
Temp.(°C)14.4 ± 10.74.4 ± 0.113.4 ± 9.84.3 ± 0.815.3 ± 12.74 ± 0.5
EC(µS cm−1)416 ± 11.3452 ± 19.8442 ± 15.6485 ± 28.3380 ± 7.1460.5 ± 0.7
pH(mg L−1 )8.4 ± 0.67.5 ± 08.6 ± 0.37.7 ± 0.68.4 ± 0.27.2 ± 0
O2(mg L−1)14 ± 2.20 ± 012.2 ± 1.52.8 ± 3.98.2 ± 0.60.2 ± 0.1
Ca2+(mg L−1)70 ± 8.0676.4 ± 7.0777.3 ± 13.7267.45 ± 0.4974.16 ± 2.7690.26 ± 1.75
Mg2+(mg L−1)4.3 ± 05.05 ± 2.191.91 ± 1.333.29 ± 0.793.54 ± 0.156.81 ± 0.11
Hard.(dH)10.8 ± 1.1312.2 ± 09.05 ± 1.4910.2 ± 0.1412 ± 0.1412.8 ± 0.14
TP(mg L−1)0.08 ± 0.0250.199 ± 0.1060.088 ± 0.0710.834 ± 0.0690.13 ± 0.020.998 ± 0.04
TN(mg L−1)1.524 ± 1.5614.357 ± 1.0081.964 ± 0.314.008 ± 3.0293.812 ± 0.1735.844 ± 0.047
Norg.-N(mg L−1)1.252 ± 1.1773.194 ± 0.5041.335 ± 0.0642.605 ± 1.5633.029 ± 0.0333.804 ± 0.157
NO3-N(mg L−1)0.1 ± 0.1410 ± 00.199 ± 0.1770.148 ± 0.2090.467 ± 0.1950.656 ± 0.109
NO2-N(mg L−1)0.002 ± 0.0030.004 ± 0.0050.005 ± 0.0040.003 ± 0.0040.001 ± 0.0010 ± 0
NH4-N(mg L−1)0.17 ± 0.241.16 ± 0.5090.426 ± 0.5551.253 ± 1.6790.316 ± 0.0131.384 ± 0.218
SRP(mg L−1)0.011 ± 0.0160.127 ± 0.0510.036 ± 0.0090.126 ± 0.0960.089 ± 0.0120.858 ± 0.084
N:P * 6.818.515.87.833.96.0
TSI 52.859.561
* For epilimnion water only.
Table 3. Efficiency of the springs and discharge of the outflow from the lake in 2017.
Table 3. Efficiency of the springs and discharge of the outflow from the lake in 2017.
DischargeQ (dm3 s−1)
Date15.0312.0402.08
Spring (S1)6.04.04.0
Spring (S2)8.08.07.0
Spring (S3)1.05.04.0
Spring (S4)1.02.01.0
Outflow (O)29.040.021.0
Table 4. Strontium isotope composition (87Sr/86Sr) and Sr content of components of the Lubińskie Lake hydrologic system.
Table 4. Strontium isotope composition (87Sr/86Sr) and Sr content of components of the Lubińskie Lake hydrologic system.
Component87Sr/86SrSr Content (mg L−1)
SpringSummerSpringSummer
Lake water-open lake
0 m0.710798 ± 100.710790 ± 100.1370.124
4 m0.710764 ± 100.710784 ± 90.1300.129
8 m0.710770 ± 100.710748 ± 110.1310.129
12 m0.710762 ± 100.710744 ± 80.1310.128
16 m0.710791 ± 120.710775 ± 90.1370.148
20 m0.710741 ± 90.710776 ± 100.1360.135
Lake water–bay
0 m0.710770 ± 10ND0.133ND
4 m0.710771 ± 10ND0.132ND
Stream
0.710765 ± 100.710741 ± 90.1430.135
Springs
S10.711851 ± 100.711710 ± 100.1750.144
S20.710746 ± 100.710525 ± 100.1350.127
S30.710463 ± 110.710412 ± 90.1290.133
S40.710813 ± 110.710844 ± 160.1420.132
Groundwater
P2ND0.710693 ± 10ND0.266
P40.710922 ± 100.710906 ± 100.2270.241
Rain
0.708879 ± 11NDNDND
ND—no data.

Share and Cite

MDPI and ACS Style

Zieliński, M.; Szczucińska, A.; Drożdżyński, M.; Frankowski, M.; Pukacz, A. Water Quality Assessment of a Meromictic Lake Based on Physicochemical Parameters and Strontium Isotopes (87Sr/86Sr) Analysis: A Case Study of Lubińskie Lake (Western Poland). Water 2019, 11, 2231. https://doi.org/10.3390/w11112231

AMA Style

Zieliński M, Szczucińska A, Drożdżyński M, Frankowski M, Pukacz A. Water Quality Assessment of a Meromictic Lake Based on Physicochemical Parameters and Strontium Isotopes (87Sr/86Sr) Analysis: A Case Study of Lubińskie Lake (Western Poland). Water. 2019; 11(11):2231. https://doi.org/10.3390/w11112231

Chicago/Turabian Style

Zieliński, Mateusz, Anna Szczucińska, Mateusz Drożdżyński, Marcin Frankowski, and Andrzej Pukacz. 2019. "Water Quality Assessment of a Meromictic Lake Based on Physicochemical Parameters and Strontium Isotopes (87Sr/86Sr) Analysis: A Case Study of Lubińskie Lake (Western Poland)" Water 11, no. 11: 2231. https://doi.org/10.3390/w11112231

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop