Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 11979 (2022)
651
1
1
Metrics details
The ecosystem of the Naroch Lakes (Belarus) includes three water bodies, Lake Batorino, Lake Myastro and Lake Naroch. These lakes have a common catchment area. At the end of the 80 s, the ecosystem of the Naroch Lakes underwent a transformation, during which the nutrient load on the catchment area decreased, and the concentration of phosphorus as a limiting factor in these water bodies decreased significantly. At the same time, the Naroch Lakes were exposed to zebra mussel (Dreissena polymorpha Pallas) invasion. In the mid-90 s, the biological and hydrochemical characteristics of the ecosystem of the Naroch Lakes stabilized. We show here that complex processes associated with the transformation of the lake ecosystem and affecting both trophic interactions in the Naroch Lakes and the influence of environmental factors on them can be represented using a single parameter, the phase-locking index, PLI. In this case, a statistically significant numerical value of PLI characterizes the phase synchronization of two processes, oscillations of the concentration of total phosphorus, TP, and oscillations of the concentration of chlorophyll, Chl. We show that the phase synchronization of these processes occurs only after the stabilization of the ecosystem of the Naroch Lakes. In the period preceding the transformation of the lake ecosystem, there was no synchronization. Numerical evaluation of PLI as a holistic parameter allows us to characterize the transformation of the lake ecosystem as a whole without resorting to study of complex interactions of various factors involved in this transformation.
The study of the phosphorus load, as well as the factors determining the dynamics of phosphorus in reservoirs, became of paramount importance in aquatic ecology after it was established in the middle of the twentieth century that phosphorus significantly limits the development of phytoplankton in freshwater reservoirs and, accordingly, determines their potential productivity1,2,3,4,5,6,7. Despite the long-lasting history of studying the effect of phosphorus on water bodies, a number of issues still remain not completely obvious, and the problems of eutrophication, which worsened in the middle of the twentieth century, remain relevant to the present7,8,9,10. All this determines the continuing interest in the study of phosphorus dynamics in various types of reservoirs.
Phosphorus exists in different forms in water. It can be dissolved, bound to particles of soil and other materials or contained within living or decaying plants and animals. In addition, phosphorus can enter lakes due to runoff from farmland or fertilized lawns10. Total phosphorus (TP) represents all the forms of phosphorus listed above combined. Relationship of phosphorus with chlorophyll (Chl) allows TP to be used to determine lake trophic state11.
Chlorophyll a, essential for plant photosynthesis, is regularly measured12,13 as an indicator of total biomass. The concentration of Chl is one of the key indicators to assessing the productivity and ecological state of aquatic ecosystems11,14.
In many lakes, TP is mainly associated with phytoplankton9,15,16,17,18,19. As a result, variations in phytoplankton biomass would be associated with changes in both Chl and TP, which may lead to a correlation between the two. However, in addition to phosphorus, nitrogen is sometimes a limiting nutrient for algal growth7,20,21,22. The Chl-TP relationships also depend on a number of other factors, such as trophic interactions in the phytoplankton-zooplankton chain and variations in the concentration of dissolved organic carbon10. The combination of these factors can result in a complex, poorly predictable character of variations in the Chl-TP relationship.
Here, we compare the Chl-TP relationships before and after the transformation of the ecosystem of the Naroch Lakes, which occurred as a result of both a biological invasion and a decrease in the nutrient load. As a quantitative characteristic of such a transformation, we used the magnitude of the phase synchronization of Chl and TP oscillations. Numerical estimates of synchronization of oscillatory processes are widely used in the analysis of the dynamics of natural, technical and social systems23. In this paper, the phase-locking index (PLI)24 acts as a synchronization measure. PLI has gained widespread application in medical research25. At the same time, as to our knowledge, PLI has rarely been used to study ecological processes (see25 as an example). We show here that the phase-locking index (PLI), which characterizes the phase synchronization of irregular Chl and TP oscillations, can be considered as the indicator of the transformation of the lake ecosystem.
We investigated the Chl-TP relationships in the Naroch Lake system, which is situated in the Northwestern Belarus (Fig. 1). It includes three lakes that are interconnected by channels, Lake Batorino, Lake Myastro and Lake Naroch, and are all polymictic, i.e., intensively mixed (Table 1).
Geographical location of the Naroch Lakes. The map was created using ArcMap 10.6.1 (https://www.arcgis.com/index.html). The coordinates of the sampling sites are marked (see also the “Methods” section).
The Naroch Lakes were monitored intensively during the last fifty years. In the late 1970s, the increased nutrient load on the catchment area was followed by eutrophication of the Naroch Lakes. Hereafter, two external factors had a significant impact on the ecosystem of the Naroch lakes: (1) implementation since 1981 of the State Program of integrated use and protection of water and land resources of the Lake Naroch basin, which, in fact, was the beginning of the period of re-oligotrophication; (2) the introduction and mass distribution of the water filtering mollusk Dreissena polymorpha Pallas in all three water bodies since the early 1990s.
In the early 80 s, a General Scheme for the integrated use and protection of water and land resources of the lake basin was adopted. A number of environmental measures were carried out in the catchment area of the Naroch Lakes. Farmland was repurposed, the water protection zone was planted, the remaining livestock farms were removed from the catchment area and improved, and a bypass collector was built for collecting, subsequent purification and discharge of wastewater outside the Naroch Lakes catchment area. As a result of the measures carried out, the phosphorus load from the catchment area on Lake Naroch decreased by 40–48%, on Lake Myastro by 22–27%, and on Lake Batorino by 34–36%27.
In addition, the Naroch Lakes were affected by the invasion of zebra mussel Dreissena polymorpha Pallas28,29,30. In Lake Myastro, Dreissena was first noted in 1984, and the first Dreissena finds in Lake Batorino date back to about the same time34. In the following years, there was an abrupt growth of the Dreissena population. In 1990, the average biomass of Dreissena (± standard error) in Lake Naroch was 1.5 ± 0.6 g/m2, and by 1993 was 99 ± 30 g/m2. In Lake Myastro, the average biomass of Dreissena in 1993 was 402 ± 187 g/m2, and in Lake Batorino was 79 ± 13 g/m2. Considering that within two years after the first detection in lakes, the abundance of Dreissena usually reaches values close to the maximum28, it can be assumed that since 1991 in Lake Naroch and a little earlier in Lake Myastro and Lake Batorino, Dreissena began to significantly affect ecosystem processes. The stabilization of biological and hydrochemical characteristics of these water bodies came in the mid-1990s29,31,32.
More information is available in the “Methods” section.
Chl and TP oscillations illustrate the transformation in the Naroch lake system from a period characterized by the high nutrient load (period I) to period characterized with low nutrient load and dreissenids (period II) (Fig. 2). Although the oscillations remain irregular throughout the entire observation period, their amplitude decreased significantly following the transition from period I to period II.
Time series representing the Chl (Cchloro) and TP (Cphosph) concentration oscillations: (a) Lake Batorino, (b) Lake Myastro, (c) Lake Naroch. Vertical dashed dotted lines show the boundaries of the division of the studied data into the period (I) of high anthropogenic load (1979–1988) and the period (II) of low anthropogenic load (1993–2018). Short-lived (1989–1992) transients are not analyzed here.
The time series of monthly chlorophyll and total phosphorus were correlated for Lake Myastro and Lake Batorino in both period I and II; however, the Chl oscillations in Lake Naroch had no significant correlation with the TP oscillations for either time period (Table 2). A comparison of these results will be discussed in context of the estimates of the phase synchronization (i.e., the phase-locking index, PLI) of the Chl and TP oscillations.
To estimate the phase synchronization of two oscillatory processes, it is necessary to remember that the numerical value of PLI can depend on the moment of time (tn) to which the initial phase value equal to 2πn (n = 0, 1, 2, …) is assigned (see Methods for more details). The values of PLI (for 1979–1988) for each of the reservoirs under study are shown in Fig. 3. As can be seen from Fig. 3, the PLI values are inside the PLI distributions that were obtained during testing using surrogate data, i.e. inside the distributions obtained as a result of random shuffling of the initial Chl and TP time series; the shuffled time series are independent and out of sync34. During period I (1979–1988), there was no synchronization detected between the oscillations of Chl and TP concentrations in any of the Naroch system lakes for either moments tn corresponding to the maximum or minimum values of Chl and TP time series (Fig. 3).
Values of PLI (shown by *) for 1979–1989, and the corresponding distributions of PLI values for the surrogate data resulted from 1000 random shuffles of the intial Chl and TP time series for each of the Naroch Lakes: (a) Lake Batorino (PLI = 0.26), (b) Lake Myastro (PLI = 0.18), (c) Lake Naroch (PLI = 0.04).
The results of the analysis of the phase relations between the Chl and TP time series for the values of tn (period II 1993–2018), which are associated with the maximum values of these time series (Fig. 4) and with the minimum values (Fig. 5), suggest that the oscillations of Chl and TP concentrations in every lake of the Naroch system are actually synchronized in phase. Note that although the PLI values (Fig. 4) lie within the surrogate data distributions, nevertheless, the values of the corresponding significance values (k = 96% for Lake Batorino and k = 98% for Lake Naroch) indicate phase synchronization of the dynamics of Chl and TP concentrations in these water bodies. Note that k = 100% corresponds to such numerical values of PLI that are outside the distribution of surrogate data. However, in Lake Myastro (provided that the moments tn correspond to the maxima of the Chl and TP time series) there is no such synchronization (since k = 86% < 95%). Nonetheless, the moments tn correspond to the minima of the Chl and TP time series in Lake Myastro are synchronized in phase (PLI = 0.30; k = 99% see Fig. 5).
Values of PLI (shown by *) for 1993–2018, and the corresponding distributions of PLI values for the surrogate data resulted from 1000 random shuffles of the initial Chl and TP time series for each of the Naroch Lakes: (a) Lake Batorino (PLI = 0.22; k = 96%), (b) Lake Myastro (PLI = 0.20; k = 86%), (c) Lake Naroch (PLI = 0.26; k = 98%). Here the tn moments correspond to the maximum values of these time series.
Values of PLI (shown by *) for 1993–2018, and the corresponding distributions of PLI values for the surrogate data resulted from 1000 random shuffles of the initial Chl and TP time series for each of the Naroch Lakes: (a) Lake Batorino (PLI = 0.17; k = 77%), (b) Lake Myastro (PLI = 0.30; k = 99%), (c) Lake Naroch (PLI = 0.18; k = 86%). Here the tn moments correspond to the minimum values of these time series.
The results of the analysis of the phase relations between the Chl and TP time series for the values of tn, which are associated with the minimum values of these time series (see Fig. 5), suggest that the oscillations of Chl and TP concentrations in Lake Myastro, which turned out to be out of sync in Fig. 4, are actually synchronized in phase.
In general, the results of the analysis of the phase synchronization of the time series of Chl concentrations and the time series of TP concentrations, taking into account the values of tn, which correspond to both the minimum and maximum values of the concentrations of Chl and TP, allow us to conclude that these time series are phase synchronized in 1993–2018.
Many studies have been devoted to the problem of the occurrence of oscillations in ecological systems (see, for example36,37. In particular, Chl and TP oscillations can be driven by a plethora of factors, such as trophic interactions, input from the catchment area, temperature fluctuations. The subject of this paper, however, is not the process of variations in the concentrations of Chl and TP over time, but the change in the synchrony of their oscillations following the biological invasion and reduction of the phosphorus load.
The time series that we present here (Fig. 2) clearly show that the environmental measures that were carried out in the catchment area of the Naroch Lakes, as well as the invasion of zebra mussel Dreissena polymorpha Pallas, were followed by significant changes in the nature of Chl and TP oscillations. These significant changes, however, were not evident when examining the correlation between oscillations in Chl and TP concentrations. Indeed, such correlations took place in Lake Batorino and Lake Myastro both before the invasion of mussels and environmental measures (period I) and after (period II). At the same time, in Lake Naroch, the correlation between oscillations in Chl and TP concentrations was practically absent both during period I and during period II (Table 2).
In search of such a numerical parameter that would be able to reflect the changes characterizing the transition from period I to period II (Fig. 2), we turned to the estimation of the phase synchronization24 (but see also38) of Chl and TP oscillations. Note that, unlike the calculation of the correlation coefficient, the calculation of phase synchronization does not use the amplitude of oscillations and is indifferent to the trend. We used numerical estimates of the phase-locking index, PLI (see Methods), in order to find out the effect of changes in nutrient load on the phase synchronization of irregular oscillations (Fig. 2) of Chl and TP, concentrations in the Naroch Lakes. This analysis approach identified that during period I of high nutrient load (1979—1988), the oscillations of Chl and TP concentrations were not synchronized in phase (Fig. 3). Following a significant decrease in the nutrient load in 1993–2018 and the invasion of zebra mussel Dreissena polymorpha, there was a noticeable phase synchronization of Chl and TP oscillations (Figs. 4 and 5). It is important to note that when determining the level of phase synchronization, the PLI value must necessarily be compared with an array of surrogate data.
While TP as a limiting nutrient is known to impact Chl concentrations, there are also other factors that could be affecting the oscillations of Chl, such as temperature, exposure to sunlight, and trophic interactions. The occurrence of the Chl-TP synchronization during period II may be due to redistribution of the relative contribution of these factors, aggravated by influence of the invasion of zebra mussel Dreissena polymorpha. When zebra mussels invade, the plankton biomass decreases while the abundance of benthophage fish increases, and primary production is more affected by higher trophic levels than in ecosystems without zebra mussels31,39. Zebra mussel can become a determining factor in the phosphorus availability40 and thereby contribute to a greater dependence of phytoplankton development on phosphorus and reduce the influence of other factors. It was also shown that the transformation of the ecosystem of the Naroch Lakes has also affected the structure of phytoplankton communities. In particular, in Lake Naroch and Lake Myastro, the biomass of diatoms decreased during this transformation41. In this context, it is worth noting that the phase-locking index, PLI, that characterizes the phase synchronization of the processes under study, reflects the impact of the plethora of biological and hydrological factors on the relationship between the Chl dynamics and TP variations, and therefore can be considered as a holistic parameter.
Along with phosphorus, water temperature has a significant effect on the vital activity of phytoplankton9,42,43,45. In this regard, we did a comparative assessment of the relationship using PLIs between temperature and Chl oscillations in the Naroch Lake ecosystem before and after the transformation.
Oscillations of Chl and variations in the temperature averaged over the water column remained correlated with each other in both period I and period II despite the drastic changes that the ecosystem of the Naroch Lakes has undergone (Table 3). Additionally, the Naroch Lakes are highly susceptible to wind mixing, thus are considered polimictic (Table 1) water bodies. In these water bodies, there were significant correlations between the average temperature and the temperature oscillations at different depths26.
The results of estimates of the PLI values in period I and period II are shown in Table 4. It is noteworthy that, unlike Chl and TP oscillations, synchronization between which was absent during period I (Fig. 3) and was manifested as statistically significant PLI values only during period II (Figs. 4 and 5), the Chl and water temperature oscillations demonstrated phase synchronization in both period I and period II (Table 4), and thus their synchronization could not be used as an indicator of the transformation in the ecosystem of the Naroch Lakes. It remains in this study that the PLI characterizing the synchronization of Chl and TP oscillations is a key indicator of the lake transformation. However, the mechanisms behind their synchrony require further study.
Up to now, the phase synchronization of ecological processes remains an insufficiently investigated phenomenon. In this context, it is time to open the gates for a more comprehensive approach to the study of phase synchronization of ecological processes, which involves the use of both field monitoring data analysis and mathematical modeling methods to cover the relationships (including synchronization) between processes that are overlooked in environmental studies.
Water samples were collected monthly at specific monitoring points of the pelagic zones of the Naroch Lakes during the vegetative season (from May to October) using two-liter Ruttner sampler. The samples were collected from The samples were collected from six different depths (0.5, 3, 6, 8, 12 and 16 m) in Lake Naroch (54° 53′ 10.56″ N, 26° 43′ 12.12″ E), from four depths (0.5, 4, 7 and 9 m) in Lake Myastro (54° 52′ 0.91″ N, 26° 52′ 49.86″ E), and from three depths (0.5, 3 and 5 m) in Lake Batorino (54° 50′ 47.94″ N, 26° 58′ 3.36″ E). The water samples from all depths were mixed. Volumes of samples collected from different depths were proportional to the total estimated volume of the water in this horizon in the lake (according to bathymetry tests). Temperature was measured at all tested depths using a mercury deep-water thermometer with a scale resolution of 0.1 °C.
The suspended matter for determination of chlorophyll content (without correction for the presence of pheopigments) was collected on the nuclear membrane filters with a pore diameter of 1.5 μm. The analysis of chlorophyll was carried out by standard spectrophotometry. Pigments were extracted using 90% acetone. The chlorophyll concentration was calculated as described in14,46.
The concentration of the total phosphorus was determined after mineralization of a raw water sample with potassium persulfate in acid medium on water bath47,48.
Spearman’s Rank Correlation Coefficient was used as nonparametric measure of correlation strength. The Spearman correlation between two variables is equal to the Pearson correlation between the rank values of those two variables:
(rho =frac{cov(M,S)}{sigma (M)sigma (S)}), where (cov(M,S)) is the covariance of the rank variables, (sigma (M)) and (sigma (S)) are the standard deviations of the rank variables.
The significance was tested using t-test: (t=rho sqrt{frac{n-2}{1-{rho }^{2}}}) with n − 2 degrees of freedom under the null hypothesis.
For irregular oscillations typical of the Chl and TP dynamics (Fig. 2) the phase of oscillations can be defined as the function23.
where t is time, tk (k = n, n + 1) is the k-th point in time at which the oscillation under study reaches its maximum/minimum.
To assess the degree of synchronization of two oscillatory processes, a measure of synchronization, the phase-locking index (PLI) was suggested24. It is defined as
where N is the number of measurements, and Δφ is the phase difference between oscillatory processes. PLI is restricted to the interval [0, 1] and reaches 1 iff the time series are strictly synchronized whereas for unsynchronized time series (i.e. for a uniform distribution of Δφ) PLI = 0.
In real data neither of these extreme values can be observed, but values between 0 and 1 are typical. Statistical significance testing must be done to establish whether a PLI value resulted from the analysis of phase relation between time series indicates a real dynamical coupling between the processes under study. Testing with surrogate data35 allows estimating how much synchronized the processes are.
The studied time series contained an insignificant amount of missing data (see Supporting Materials S1). The missing values were imputed by seasonally splitted missing value imputation (with interpolation) using the impute TS package49,50. Imputation algorithm splits the times series into seasons and afterwards performs imputation separately for each of the resulting time series datasets (each containing the data for one specific season). The time series obtained as a result of the imputation are shown in Fig. 2.
All data generated or analysed during this study are included in this published article [and its supplementary information files].
Sakamoto, M. Primary production by phytoplankton community in some Japanese lakes and its dependence on lake depth. Archiv für Hydrobilogie. 62, 1–28 (1966).
Google Scholar
Vollenweider, R. A. Scientific fundamentals of the eutrophication of lakes and flowing waters, with particular reference to nitrogen and phosphorus as factors in eutrophication (Organisation for Economic Co-operation and Development, 1968).
Google Scholar
Edmondson, W. T. Phosphorus, nitrogen, and algae in Lake Washington after diversion of sewage. Science 169, 690–691 (1970).
Article ADS CAS Google Scholar
Dillon, P. J. & Rigler, F. H. The phosphorus-chlorophyll relationship in lakes. Limnol. Oceanogr. 19, 767–773 (1974).
Article ADS CAS Google Scholar
Jones, J. R. & Bachmann, R. W. Prediction of phosphorus and chlorophyll levels in lakes. J. Water Pollut. Control Feder. 48, 2176–2182 (1976).
CAS Google Scholar
Schindler, D. W. Evolution of phosphorus limitation in lakes. Science 195, 260–262 (1977).
Article ADS CAS Google Scholar
Filstrup, C. T. & Downing, J. A. Relationship of chlorophyll to phosphorus and nitrogen in nutrient-rich lakes. Inland Waters. 7, 385–400 (2017).
Article CAS Google Scholar
Schindler, D. W. Recent advances in the understanding and management of eutrophication. Limnol. Oceanogr. 51, 356–363 (2006).
Article ADS Google Scholar
Quinlan, R. et al. Relationships of total phosphorus and chlorophyll in lakes worldwide. Limnol. Oceanogr. 66, 392–404 (2020).
Article ADS Google Scholar
Yuan, L. L. & Jones, J. R. Rethinking phosphorus–chlorophyll relationships in lakes. Limnol. Oceanogr. 65, 1847–1857 (2020).
Article ADS CAS Google Scholar
Carlson, R. E. A trophic state index for lakes. Limnol. Oceanogr. 11, 361–369 (1977).
Article ADS Google Scholar
Neveux, J. et al. Comparison of chlorophyll and phaeopigment determinations by spectrophotometric, fluorometric, spectrofluorometric and HPLC methods. Mar. Microb. Food Webs 4, 217–238 (1990).
Google Scholar
Lampert, W. & Sommer, U. Limnoecology (Oxford University, 2007).
Google Scholar
Kovalevskaya, R. Z., Zhukava, H. A. & Adamovich, B. V. Modification of the method of spectrophotometric determination of chlorophyll a in the suspended matter of water bodies. J. Appl. Spectrosc. 87, 72–78 (2020).
Article ADS CAS Google Scholar
Søndergaard, M., Lauridsen, T. L., Johansson, L. S. & Jeppesen, E. Nitrogen or phosphorus limitation in lakes and its impact on phytoplankton biomass and submerged macrophyte cover. Hydrobiologia 795, 35–48 (2017).
Article Google Scholar
Søndergaard, M., Jensen, J. P., Jeppesen, E. & Møller. P. H. Seasonal dynamics in the concentrations and retention of phosphorus in shallow Danish lakes after reduced loading. Aquat. Ecosyst. Health Manag. 5(1), 19–29 (2002).
Magumba, D., Atsushi, M., Michiko, T., Akira, K. & Masao, K. Relationships between Chlorophyll-a, phosphorus and nitrogen as fundamentals for controlling phytoplankton biomass in lakes. Environ. Control. Biol. 51(4), 179–185 (2013).
Article CAS Google Scholar
Smith, V. H. & Shapiro, J. Chlorophyll-phosphorus relations in individual lakes. Their importance to lake restoration strategies. Environ. Sci. Technol. 15(4), 444–451 (1981).
Pothoven, S. A. & Vanderploeg, H. A. Seasonal patterns for Secchi depth, chlorophyll a, total phosphorus, and nutrient limitation differ between nearshore and offshore in Lake Michigan. J. Great Lakes Res. 46, 519–527 (2020).
Article CAS Google Scholar
Søndergaard, M. & Jeppesen, E. Lake Søbygaard, Denmark: phosphorus dynamics during the first 35 years after an external loading reduction. In: Internal Phosphorus Loading: Causes, Case Studies, and Management (ed. Steinman, A.D. & Spears, B. M.) 285–299 (J. Ross, Plantation, 2020).
Guildford, S. J. & Hecky, R. E. Total nitrogen, total phosphorus, and nutrient limitation in lakes and oceans: Is there a common relationship?. Limnol. Oceanogr. 45, 1213–1223 (2000).
Article ADS CAS Google Scholar
Jones, J.R. et al. Nutrients, seston, and transparency of Missouri reservoirs and oxbow lakes. An analysis of regional limnology. Lake Reser. Manag. 24, 155–180 (2008).
Pikovsky, A., Rosenblum, M. & Kurths, J. Synchronization. A universal concept in nonlinear sciences (Cambridge University, 2001).
Book Google Scholar
Kuramoto, Y. Chemical Oscillations, Waves, and Turbulence (Springer, 1984).
Book Google Scholar
Sazonov, A. V. et al. An investigation of the phase locking index for measuring of interdependency of cortical source signals recorded in the EEG. Biol. Cybern. 100, 129–146 (2009).
Article Google Scholar
Medvinsky, A. B. et al. Temperature as a factor affecting fluctuations and predictability of the abundance of lake bacterioplankton. Ecol. Complex. 32, 90–98 (2017).
Article Google Scholar
Zhukova, T. V. & Ostapenya, A. P. Estimation of efficiency of nature protection measures in water catchment area of the Naroch lakes. Natural Resources. 3, 68–73 (2000) ((in Russian)).
Google Scholar
Burlakova, L. E., Karatayev, A. Y. & Padilla, D. K. Changes in the distribution and abundance of Dreissena polymorpha within lakes through time. Hydrobiologia 571, 133–146 (2006).
Article Google Scholar
Ostapenia, A. P. et al. Bentification of lake ecosystem: causes, mechanisms, possible consequences, prospect for future research. Trudy BGU. 7, 135–148 (2012) ((in Russian)).
Google Scholar
Karatayev, A.Y., Burlakova, L.E. & Padilla, D.K. Impacts of Zebra Mussels on aquatic communities and their role as ecosystem engineers. In: Leppäkoski, E., Gollasch, S., Olenin, S. (eds) Invasive Aquatic Species of Europe. Distribution, Impacts and Management (Springer, Dordrecht, 2002).
Adamovich, B. V. et al. The divergence of chlorophyll dynamics in the Naroch Lakes. Biophysics 60, 632–638 (2015).
Article CAS Google Scholar
Zhukova, T. V. et al. Long-term dynamics of suspended matter in Naroch Lakes: Trend or intervation. Inland Water Biol. 10, 250–257 (2017).
Article Google Scholar
Adamovich, B. V. et al. Eutrophication, oligotrophication, and benthiphication in Naroch Lakes: 40 years of monitoring. J. Siber. Federal Univ. Biol. 10, 379–394 (2017).
Article Google Scholar
Ostapenya A.P. et al. Ecological passport of Lake Myastro (EcoMir, Minsk, 1994) (in Russian).
Kantz, H. & Schreiber, T. Nonlinear time series analysis (Cambridge University, 1997).
MATH Google Scholar
Kot, M. Elements of mathematical ecology (Cambridge University, 2001).
Book Google Scholar
Turchin, P. Complex population dynamics. A Theoretical/Empirical Synthesis (Princeton University, Princeton, 2003).
MATH Google Scholar
Cazelles, B. & Stone, L. Detection of imperfect population synchrony in an uncertain world. J. Anim. Ecol. 72, 953–968 (2003).
Article Google Scholar
Karatayev, A. Y., Burlakova, L. & Padilla, D. K. The effects of Dreissena polymorpha (Pallas) invasion on aquatic communities in Eastern Europe. J. Shellfish Res. 16, 187–203 (1997).
Google Scholar
Lia, J. et al. Benthic invaders control the phosphorus cycle in the world’s largest freshwater ecosystem. PNAS 118(6), e2008223118. https://doi.org/10.1073/pnas.2008223118 (2021).
Article CAS Google Scholar
Mikheyeva, T. M. et al. The dynamics of freshwater phytoplankton stability in the Naroch Lakes (Belarus). Ecol. Ind. 81, 481–490 (2017).
Article Google Scholar
Harris, P. H. Phytoplankton ecology. Structure, functioning and flucttuation (Chapman & Hall, London, New York, 1986).
Jeppesen, E., Jensen, J. P., Søndergaard, M. & Lauridsen, T. L. Response of fish and plankton to nutrient loading reduction in eight shallow Danish lakes with special emphasis on seasonal dynamics. Freshw. Biol. 50, 1616–1627 (2005).
Article CAS Google Scholar
Nezlin, N.P. & Li, B-L. Time-series analysis of remote-sensed chlorophyll and environmental factors in the Santa Monica–San Pedro Basin off Southern California. J. Mar. Syst. 39, 185–202 (2003).
French, T. D. & Petticrew, E. L. Chlorophyll a seasonality in four shallow eutrophic lakes (northern British Columbia, Canada) and the critical roles of internal phosphorus loading and temperature. Hydrobiologia 575, 285–299 (2007).
Article CAS Google Scholar
SCOR-UNESCO Working Group no. 17. Determination of photosynthetic pigments in sea-water. Monographs on Oceanologic Methodology 9–18 (UNESSCO, Paris, 1966).
Semenov, A. D. Guide on the chemical analysis of continental surface waters (Gidrometeoizdat, 1977) ((in Russian)).
Google Scholar
Wetzel, R. G. & Likens, G. E. Limnological analysis (Springer, 2000).
Book Google Scholar
Steffen, M. & Bartz-Beielstein, T. imputeTS: time series missing value imputation in R. R J. 9(1), 207–218 (2017).
Article Google Scholar
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, 2020).
Download references
We would like to thank the anonymous reviewers for reading and commenting on an earlier version of this manuscript. The authors thank the researchers and technicians of Belarusian State University Naroch Biological Station who helped to collect and produce the data presented here. The authors thank A. Novikov from National Park “Narochansky” who helped to create the map presented on Fig. 1. This work was carried out within the framework of State Task no. 075-00381-21-00 and State Task of Belarus (GPNI GR no. 20212317), with the financial support of The Russian Foundation for Basic Research, project no. 20-54-00017, and The Belarusian Republican Foundation for Fundamental Research.
This article was funded by Russian Foundation for Basic Research (20-54-00017), Belarusian Republican Foundation for Fundamental Research, Ministry of Education of the Republic of Belarus.
Institute of Theoretical and Experimental Biophysics, Russian Academy of Sciences, Pushchino, 142290, Russia
Alexey V. Rusakov, Nailya I. Nurieva & Alexander B. Medvinsky
Belarusian State University, 220030, Minsk, Belarus
Boris V. Adamovich, Raisa Z. Kovalevskaya, Tamara M. Mikheyeva, Hanna A. Zhukava, Yulia K. Veres & Tatyana V. Zhukova
Moscow State University of Psychology and Education, Moscow, 127051, Russia
Nataly P. Radchikova
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
B.A. is the co-first author. B.A., R.K., T.M., H.Z., Y.V., T.Z. carried out the monitoring of the Naroch Lakes, performed field sampling and the water chemistry analysis. A.R., B.A., N.N., N.R., A.M. conducted the analysis of the time series, A.M., B.A., A.R., N.N., N.R. wrote the paper.
Correspondence to Alexander B. Medvinsky.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Rusakov, A.V., Adamovich, B.V., Nurieva, N.I. et al. Phase synchronization of chlorophyll and total phosphorus oscillations as an indicator of the transformation of a lake ecosystem. Sci Rep 12, 11979 (2022). https://doi.org/10.1038/s41598-022-16111-8
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-16111-8
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Scientific Reports (2023)
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.
