
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39294208
70408
10.1038/s41598-024-70408-4
Article
Identifying environmental impacts on planktonic algal proliferation and associated risks: a five-year observation study in Danjiangkou Reservoir, China
Nong Xizhi nongxizhi@gxu.edu.cn

12
Guan Xian 1
Chen Lihua 1
Wei Jiahua 2
Li Ronghui lironghui@gxu.edu.cn

1
1 https://ror.org/02c9qn167 grid.256609.e 0000 0001 2254 5798 School of Civil Engineering and Architecture, Guangxi University, Nanning, 530004 China
2 https://ror.org/03cve4549 grid.12527.33 0000 0001 0662 3178 State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing, 100084 China
16 9 2024
16 9 2024
2024
14 215685 7 2024
16 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.
Understanding the risks of planktonic algal proliferation and its environmental causes is crucial for protecting water quality and controlling ecological risks. Reservoirs, due to the characteristics of slow flow rates and long hydraulic retention times, are more prone to eutrophication and algal proliferation. Chlorophyll-a (Chl-a) serves as an indicator of planktonic algal biomass. Exploring the intricate interactions and driving mechanisms between Chl-a and the water environment, and the potential risks of algal blooms, is crucial for ensuring the ecological safety of reservoirs and the health of water users. This study focused on the Danjiangkou Reservoir (DJKR), the core water source of the Middle Route of the South-to-North Water Diversion Project of China (MRSNWDPC). The multivariate statistical methods and structural equation modeling were used to explore the relationships between chlorophyll-a (Chl-a) contents and water quality factors and understand the driving mechanisms affecting Chl-a variations. The Copula function and Bayesian theory were combined to analyze the risk of changes in Chl-a concentrations at Taocha (TC) station, which is the core water source intake point of the MRSNWDPC. The results showed that the factors driving planktonic algal proliferation were spatially heterogeneous. The main factors affecting Chl-a concentrations in Dan Reservoir (DR) were water physicochemical factors (water temperature, dissolved oxygen, pH value, and turbidity) with a total contribution rate of 60.18%, whereas those in Han Reservoir (HR) were nutrient factors (total nitrogen, total phosphorus, and ammonia nitrogen) with a total contribution rate of 73.58%. In TC, the main factors were water physicochemical factors (turbidity, pH, and water temperature) and nutrient factors (total phosphorus) with total contribution rates of 39.76% and 45.78%, respectively. When Chl-a concentrations in other areas of the DJKR ranged from the minimum to the uppermost quartile, the probabilities that Chl-a concentrations at the TC station exceeded 3.4 μg/L (the benchmark value of Chl-a for lakes in the central-eastern lake area of China) owing to the influence of these areas were all less than 10%. Thus, the risk of planktonic algal proliferation at the MRSNWDPC intake point is low. This study developed an integrated framework to investigate spatiotemporal changes in algal proliferation and their driving factors in reservoirs, which can be used to support water quality management in mega hydro projects.

Keywords

Multivariate statistical analysis
Structural equation modeling
Algal proliferation
Water quality
Reservoir
Subject terms

Environmental sciences
Environmental impact
http://dx.doi.org/10.13039/501100018571 Specific Research Project of Guangxi for Research Bases and Talents No.AD22035185 Nong Xizhi National Natural Science Foundation of ChinaNo.52309016 No.52079033 Nong Xizhi Li Ronghui Belt and Road Special Foundation of the National Key Laboratory of Water Disaster PreventionNo.2022nkms06 Nong Xizhi issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Due to the growing global population and socioeconomic development, water resource shortages and the balance between supply and demand are becoming increasingly acute, threatening sustainable development and ecosystem health/stability1,2. Water pollution is one of the main causes of water resource scarcity3. Globally, approximately 44% of domestic wastewater is discharged into the environment without safe treatment. The situation in developing countries is worse, wherein nearly 90% of sewage is discharged without treatment into natural water systems4. To address water shortage and pollution problems, many countries have constructed reservoirs and water diversion projects that provide reliable water sources for agriculture, industry, and urban development, improve living conditions, and play substantial roles in rational water resource use and protection3.

Reservoirs worldwide face numerous ecological threats5. Among these threats, eutrophication and the algal blooms it triggers are the most common ecological challenges facing reservoirs6. An algal bloom refers to the phenomenon in which phytoplankton proliferate and accumulate in quantities that exceed a threshold in an aquatic environment with sufficient nutrient salts7,8. Algal bloom events are more likely to occur in water bodies like lakes and reservoirs, where the long hydraulic retention time leads to the easy accumulation of nutrient8,9. As China contains over 20% of the global population and more than 97,000 reservoirs, water pollution and eutrophication risks affect the health and safety of a portion of the population. The safety of the water supply in Chinese reservoirs also substantially influences water quality and water resource management worldwide.

The Danjiangkou Reservoir (DJKR) is the largest artificial lake in Asia and representative of reservoirs in China, as it is the core water source of the Middle Route of the South-to-North Water Diversion Project of China (MRSNWDPC), which is the longest inter-basin water diversion project in the world. Since the commissioning of the MRSNWDPC in 2014, the DJKR has transferred more than 60.5 billion m3 of water to North China, making it the primary water source for 26 large and medium cities and directly benefiting more than 100 million people10. The water quality of the DJKR directly affects the efficiency and effectiveness of the MRSNWDPC and the drinking water safety in the receiving area; therefore, it is essential to ensure the water quality and safety of the DJKR.

Researchers have conducted many investigations on the water quality and phytoplankton distribution status in the DJKR. The water quality of the DJKR was not very good before the water diversion project was commissioned in 201411. Of which, the water quality of the DJKR deteriorated from Classes I or II to II or III from 2007 to 201012. After the water diversion project opened, the DJKR maintained good water quality conditions13. The phytoplankton population in the DJKR in 2009–2010 comprised diatoms and cyanobacteria, and the community structure was primarily influenced by conductivity14. To prepare for water diversion to the MRSNWDPC, the storage level in the DJKR was raised, after which the phytoplankton biomass decreased and cyanobacterial abundance increased15. After the water diversion project was commissioned, the proportion of diatoms in the DJKR decreased considerably, and the cyanobacteria and green algae became the dominant species from 2018 to 202016. The DJKR currently exhibits mesotrophic conditions17. Although the eutrophication level of the DJKR is not high. Still, as the core water source of the MRSNWDPC, it would jeopardize the water supply of North China and threaten over 100 million people in the adverse event of environmental risk. In the past years, some sections of the DJKR have also experienced severe nitrogen and phosphorus pollution18, nearing eutrophic conditions19. This indicates that the potential risk of algal proliferation does exist in the DJKR. Under this premise, assessing the potential risk of algal proliferation is essential. Therefore, it is crucial to understand the complex relationships between Chl-a concentrations and environmental factors in the DJKR to prevent the occurrence of algal blooms and safeguard water quality in both the DJKR and downstream areas.

Previous studies of the DJKR have shown that water quality between the different water bodies can influence each other due to their hydrologic connectivity13. Since the DJKR is the water source for the MRSNWDPC, the freshwater resources from the Dan and Han reservoirs pass through the Taocha (TC) station and are transported by the main canal to North China. Thus, the water quality at TC station can be affected by the sub-reservoirs, which in turn affects the water quality of the MRSNWDPC. Therefore, it is necessary to monitor the water quality at the TC intake of the MRSNWDPC to determine its algal proliferation potential risk under the influence of water exchange in the reservoir area.

Phytoplankton are affected by various environmental factors, including water temperature, pH value, nutrient contents, and their combined effects20–22. Under suitable environmental conditions (e.g., adequate nutrient contents, appropriate light and temperature conditions, and low water flow), algae proliferate rapidly23,24. This proliferation is also related to multiple environmental factors that interact with each other25–27. There have been numerous studies on the effects of environmental factors on phytoplankton in the DJKR28–30. These studies primarily used methods such as correlation analysis and Mantel tests, which can identify and describe the correlation between variables but are not effective in inferring causal relationships between multiple factors and planktonic algae. A more innovative approach is needed to explore the driving mechanisms of planktonic algae under the influence of multiple factors. Structural equation modeling (SEM) can be used to overcome this difficulty, as it can handle causal relationships between multiple factors31. SEM can also capture the direct effects between variables, as well as indirect effects generated by the mediating variables32. SEM visualizes the pathways through which environmental variables affect planktonic algae and has been widely applied to determine the mechanisms that influence the response and explanatory variables33,34.

The objectives of this study were to investigate the intricate interactions and driving mechanisms between Chl-a and environmental factors in the DJKR, and to assess the risk of algal blooms at the intake of the MRSNWDPC while considering the impacts of surrounding water areas. By employing structural equation modeling and multivariate statistical techniques, the research will provide a more comprehensive and detailed understanding of the interactions between Chl-a and environmental factors across the entire reservoir and at different stations, as well as the mutual influences of Chl-a across different stations, thus enhancing the understanding of algal proliferation mechanisms under multiple factors. However, due to limitations in data and numerical models, this study does not account for the effects of climate and hydrodynamic factors on algae, nor the impacts of pollutant inputs and biological migration under different hydrological conditions, which need to be addressed in future research.

Materials and methods

Study area and data collection

The DJKR (32° 36′–33° 48′ N, 110° 59′–111° 49′ E) straddles Hubei and Henan provinces and is the largest lake-type reservoir and protected drinking water source area in Asia, with a standard storage level of 170 m, storage capacity of 29.05 billion m3, and watershed covering 95,000 km2. The DJKR comprises the Dan (DR, Henan Province) and Han (HR, Hubei Province) reservoirs35 and is located in a low-mountainous and hilly region with a subtropical monsoon climate (i.e., with a warm spring, hot and rainy summer, cool fall, and cold and less rainy winter)36. The average annual temperature in the DJKR area ranges from 15 to 16 °C, and the annual precipitation ranges from 800 to 1000 mm, nearly 80% of which occurs from May to October. The main soil types in the study area are brown, yellow–brown, and yellow-cinnamon soils37, and the main land use types include woodlands, agricultural land (rice paddies and uplands), orchards, and wetlands38. Since water diversion via the MRSNWDPC officially began in December 2014, the DJKR has transferred more than 60.5 billion m3 of freshwater resources to North China. The DJKR is the primary water source for 26 large and medium cities, benefits more than 108 million people, and replenishes more than 9.5 billion m3 of ecological water in more than 50 rivers in northern China.

In this study, six representative water quality monitoring stations in the DJKR area were selected based on the Chinese national water quality monitoring programs, which covered two sub-reservoirs (i.e., the Dan and Han reservoirs). The Qingshan (QS), Longkoulanghe (LK), and Baqian (BQ) stations are located in the HR, and the Madeng (MD), Taochaqujie (QJ), and Taocha (TC) stations are located in the DR. The details of the water quality monitoring stations used in this study can be found in Fig. 1 and Table S1. The Changjiang River Scientific Research Institute of Changjiang Water Resources Commission collected water samples daily from March 2018 to February 2023. Nine water quality indicators were analyzed, i.e., water temperature (WT, °C), dissolved oxygen (DO, mg/L), pH, turbidity (Tur, NTU), total nitrogen (TN, mg/L), total phosphorus (TP, mg/L), ammonia nitrogen (NH3-N, mg/L), permanganate index (CODMn, mg/L), and Chl-a (μg/L). WT, pH, and DO were measured using a water quality multi-parameter monitor. NH3-N, TP, TN, and Chl-a were measured using a UV–visible spectrophotometer. Tur was measured using a portable turbidimeter, and CODMn was measured by an acid buret. Descriptions of the analytical methods and instrumentation can be found in the Supplementary Material (Table S2). The seasons were divided into spring (March to May), summer (June to August), autumn (September to November), and winter (December to February of the following year).Figure 1 Locations of water quality monitoring stations of the DJKR from 2018 to 2023 in this study.

Methodology

The methodological framework of this study was divided into three parts to identify factors driving Chl-a, determine their roles, and estimate the influences of water bodies in the study area on algal proliferation at the TC intake (Fig. 2). We first identified the main factors affecting Chl-a at each station using stepwise multiple regression analysis. We then used principal components analysis (PCA) to downscale the multiple nutrient and water physicochemical indicators into composite variables. Using PCA eliminated potential problems related to covariance and reduced the number of variables used in the SEM, simplifying the model. SEM was then used to model the nutrient, water physicochemical, and organic pollution factors influencing Chl-a to determine the mechanism(s) driving planktonic algal blooms in the study area. Finally, the Copula function was used to construct the joint distribution functions of Chl-a in the different water bodies and at the intake. Bayesian theory was applied to estimate algal proliferation risk at the TC water intake under the influence of the different water bodies. The data used for the above calculations and analysis were daily averages from March 2018 to February 2023.Figure 2 Methodological framework used in this study to identify environmental driving pathways and assess the risk of algal proliferation.

Stepwise multiple regression analysis

The goal of the stepwise multiple regression method is to select the most important variables from a large set of variables. We used this method to determine the main factors driving Chl-a concentrations at each station. The finalized model is shown in Eq. (1):1 Chl-a=a0+∑i=1naiXi

where a0 is a constant, ai represents the different coefficients, and Xi represents the different environmental factors. A log transformation (lg(x + 1)) was applied to pre-process the data before performing the stepwise multiple linear regressions, both of which were performed using the SPSS 24 software (USA, IBM).

PCA

PCA was applied to reduce the data dimensionality and improve its interpretability. PCA minimizes information loss by creating new uncorrelated principal components that sequentially maximize the variance39. The driving factors identified in the stepwise multiple regression analysis were categorized as water physicochemical, nutrient, and organic pollution factors according to a priori knowledge. The multiple water physicochemical and nutrient factors were downscaled using PCA. PCA was not applied for organic pollution downscaling, as only one indicator (CODMn) was used in this study. The principal components obtained from the PCA were referred to as water physicochemical (e.g., PChm 1 and PChm 2) and nutrient (e.g., Nutr 1 and Nutr 2) factors. If only one principal component remained after dimensionality reduction, it was referred to as PChm or Nutr. The PCA was performed using the SPSS 24 software (USA, IBM).

Each principal component is a linear combination of the original variables, with one principal component reflecting the characteristics of multiple variables. Before computing the principal components, data standardization is necessary. This is because different variables may have different units and scales, and variables with larger scales can disproportionately influence the covariance matrix, leading to significant bias in the PCA results. This step eliminates the effects of differing units and enhances the robustness of the results. When performing principal component analysis (PCA) using SPSS, the software automatically standardizes the data at the beginning of the calculation.

Structural equation modeling

SEM is a quantitative method based on statistical techniques that can be used to deal with multi-factor causal relationships31. SEM can reflect the direct influences of the variables, as well as the indirect influences generated by mediating variables32. The direct and indirect causal relationships between multiple variables are represented as paths (directional arrows) in an SEM diagram40. The topological structure comprises multiple variable nodes, and the paths between nodes reflect their driving mechanisms. SEM has been widely used to explore the driving mechanism between multiple variables41,42.

In this study, SEM was used to determine the environmental driving mechanism of Chl-a concentrations in the DJKR. This study utilized Confirmatory Factor Analysis (CFA). Based on prior knowledge, water physicochemical factors directly affect Chl-a and indirectly influence Chl-a through organic pollution factor and nutrients. An initial Structural Equation Modeling (SEM) model was constructed according to this theory, followed by path analysis to quantify the relationships between variables. Finally, model fit was evaluated using fit indices, and the model was refined to achieve higher explanatory power and goodness of fit.

Copula function and Bayesian theory

Relationships among hydrologic variables are complex, and hydrologic variables have various probability distributions. Thus, the Copula function is commonly used to determine the correlations between variables43. The Copula function divides a joint distribution into the marginal distribution and correlation structure, which are treated separately, preserve all information about the correlation, and describe the nonlinear and asymmetric correlations between variables44. According to Sklar's theorem, a binary joint distribution function with continuous marginal distributions F(x) and G(y) has a unique Copula function as follows:2 H(x,y)=CFx(x),Fy(y)

where H(x,y) is the bivariate joint distribution function of the random variables X and Y; F(x) is the marginal distribution function of the random variable X; G(y) is the marginal distribution function of the random variable Y; and C(F(x), G(y)) is the copula function that connects the marginal distribution functions F(x) and G(y).

In this study, the Copula function was used to construct bivariate joint distribution functions for the Chl-a concentrations between different monitoring stations. The conditional distribution functions for different conditions were obtained by applying Bayesian theory. Specifically, to calculate the probability that different Chl-a concentrations at the MD, QJ, QS, LK, and BQ (Chl-aMD, Chl-aQJ, Chl-aQS, Chl-aLK, and Chl-aBQ, respectively) stations will cause the Chl-a concentrations at the TC (Chl-aTC) station to change. The minimum–25%, 25–50%, 50–75%, and 75%–maximum Chl-a concentrations at the respective monitoring stations were used in this analysis. The conditional distribution function can be defined as follows:3 F(Chl-aTC>x|y2>Chl-aothers≥y1)=1-C(x,y1)-C(x,y2)FChl-aothers(y2)-FChl-aothers(y1)

where Chl-aTC represents the Chl-a monitored at the TC measurement station; Chl-aothers represents the Chl-a monitored at other stations excluding the TC station; x is a specific value of Chl-aTC; y1 and y2 are both specific values of Chl-aothers; C(x,y) is the copula function that captures the dependence between Chl-aTC and Chl-aothers; FChl-aothers is the marginal distribution function of Chl-aothers.

A Chl-a concentration of 3.4 μg/L is a benchmark nutrient concentration for lakes in central-eastern China, as defined by the Ministry of Ecology and Environment of the People’s Republic of China45. Algal proliferation does not jeopardize water body functions when Chl-a concentrations are below this benchmark value. The Copula function and Bayesian theory were calculated using the R software package. The R packages used are “copula” and “VineCopula”.

Ethical approval

All authors have read, understood, and have complied as applicable with the statement on "Ethical responsibilities of Authors" as found in the Instructions for Authors.

Results

Water quality characteristics

Table 1 presents a statistical summary of the water quality indicators and Chl-a concentrations at the six monitoring stations in the DJKR. Kruskal–Wallis tests indicate that all indicators varied significantly among the monitoring stations (P < 0.05). During the monitoring period, WT and DO exhibited significant seasonal characteristics (Fig. S1). DO concentrations were significantly higher in spring and winter than those in summer and autumn, with a minimum of 6.10 mg/L and a maximum of 14.40 mg/L. The average pH values at all stations exceeded 8.0, with a minimum average of 8.08 (at TC station) and a maximum average of 8.21 (at MD station), indicating that the water quality in the DJKR was predominantly neutral to alkaline. In addition, the seasonal variations in pH differed significantly (Fig. S1). The average Tur value at the QS station (4.21 NTU) was higher than that at the other stations. The minimum and maximum Tur values were also observed at the QS station (0.00 and 35.10 NTU, respectively). The concentrations of NH3-N and CODMn were low, all meeting the Chinese surface water quality standards for Class I or II. The minimum observed TP concentration was 0.003 mg/L, and ranged from 0.003 to 0.041 mg/L among all stations. Significant seasonal differences in TP concentrations were observed at all stations except the MD (Fig. S1). Table 1 Statistical water quality parameters at the six monitoring stations in the DJKR from 2018 to 2023 (Avg.: Average; S.D.: Standard deviation).

Parameters	WT	DO	pH	Tur	NH3-N	TP	TN	CODMn	Chl-a	
°C	mg/L	–	NTU	mg/L	mg/L	mg/L	mg/L	μg/L	
Avg. ± S.D	MD	19.05 ± 6.61	9.69 ± 1.36	8.21 ± 0.31	3.41 ± 3.34	0.033 ± 0.010	0.014 ± 0.003	1.17 ± 0.15	1.97 ± 0.16	2.38 ± 1.17	
QJ	18.89 ± 6.75	9.43 ± 1.17	8.09 ± 0.35	3.00 ± 2.99	0.031 ± 0.009	0.015 ± 0.004	1.17 ± 0.13	1.90 ± 0.19	2.12 ± 0.80	
TC	18.76 ± 6.93	9.65 ± 1.31	8.08 ± 0.36	3.27 ± 3.28	0.031 ± 0.009	0.013 ± 0.004	1.17 ± 0.16	1.90 ± 0.18	2.15 ± 0.88	
QS	18.47 ± 6.69	9.69 ± 1.30	8.13 ± 0.31	4.21 ± 4.75	0.035 ± 0.020	0.016 ± 0.005	1.17 ± 0.16	1.92 ± 0.20	2.23 ± 1.29	
LK	19.06 ± 6.70	9.37 ± 1.09	8.09 ± 0.25	3.27 ± 2.89	0.036 ± 0.018	0.015 ± 0.003	1.18 ± 0.14	1.92 ± 0.13	2.07 ± 0.81	
BQ	19.13 ± 6.75	9.28 ± 1.23	8.12 ± 0.32	3.30 ± 2.90	0.035 ± 0.014	0.015 ± 0.004	1.15 ± 0.13	1.90 ± 0.18	2.17 ± 1.29	
Min.	MD	6.00	6.10	7.30	0.10	0.005	0.004	0.81	1.56	0.21	
QJ	6.80	6.80	7.00	0.18	0.009	0.004	0.86	1.36	0.20	
TC	4.90	6.10	6.96	0.30	0.012	0.003	0.73	1.43	0.78	
QS	6.30	6.15	6.90	0.00	0.003	0.004	0.64	1.08	0.27	
LK	4.40	7.10	7.00	0.10	0.012	0.007	0.88	1.45	0.20	
BQ	7.10	6.65	7.14	0.30	0.005	0.004	0.77	1.29	0.10	
Max.	MD	32.80	14.40	8.90	27.00	0.080	0.026	1.66	2.77	8.20	
QJ	32.50	12.50	8.81	26.80	0.082	0.034	1.61	2.71	5.06	
TC	33.40	12.70	9.00	30.20	0.073	0.023	1.61	2.37	6.96	
QS	33.30	13.71	9.31	35.10	0.173	0.041	1.82	2.74	11.60	
LK	32.50	12.00	8.70	23.80	0.147	0.027	1.61	2.46	5.46	
BQ	33.00	13.90	9.30	28.23	0.101	0.029	1.68	2.76	11.38	
Kruskal–Wallis test	P < 0.05	P < 0.05	P < 0.05	P < 0.05	P < 0.05	P < 0.05	P < 0.05	P < 0.05	P < 0.05	

Using the lowest or highest average values as examples, to illustrate the significant spatial and temporal differences for each indicator. The statistical analysis on the annual scale was provided in Table S3. Spatially, although WT passed the Kruskal–Wallis test (P = 0.41), the adjusted significance in pairwise comparisons was > 0.5. Annually, WT in 2020 was significantly lower than in other years. The remaining indicators exhibited significant differences both spatially and temporally (P < 0.5). DO was significantly lower at BQ compared to the other four stations except for LK, and significantly higher in 2020 than in other years. pH was significantly higher at MD compared to other stations, and significantly higher in 2022 than in other years. Tur was significantly higher at QS compared to other stations, and significantly higher in 2021 than in other years. NH3-N was significantly higher at LK compared to the other four stations except for BQ, and significantly higher in 2022 than in all other years except for 2021. TP was significantly higher at QS compared to the other five stations, and significantly higher in 2019 than in other years. TN was significantly lower at BQ compared to the other four stations except for MD, and significantly higher in 2022 than in all other years except for 2018. CODMn was significantly higher at MD compared to other stations, and significantly higher in 2021 than in other years.

Spatiotemporal distributions of Chl-a

Chl-a exhibited significant spatial and temporal variations (P < 0.5), with MD showing significantly higher values than other stations, and 2022 being significantly higher than other years. At the stations, the average Chl-a concentrations ranged from 2.07 to 2.38 μg/L in the entire reservoir, from 2.12 to 2.38 μg/L in the DR, and from 2.07 to 2.23 μg/L in the HR (Fig. 3a). The Chl-a concentrations in the DR were low (maximum of only 8.20 μg/L at MD station), whereas higher Chl-a concentrations were more often observed in the HR (QS and BQ stations reached 11.60 μg/L and 11.38 μg/L, respectively). Except for the QS station, significant seasonal differences in Chl-a concentrations were observed at all stations (Fig. 3b). Chl-a concentrations were significantly higher in spring than in other seasons, except at QS and MD stations. In addition, outlying high Chl-a concentrations were observed in the winter at the QS station and in the spring at the BQ station. Chl-a concentrations gradually increased from 2018 to 2023, and more outliers were observed in 2022 (Fig. 3c). The average spring Chl-a concentrations in 2018, 2019, 2020, and 2022 were higher than those during the rest of the year (Fig. 3d).Figure 3 Spatial (a, b) and seasonal (c, d) variations in chlorophyll-a (Chl-a) in the DJKR from 2018 to 2023.

Correlations between Chl-a and environmental factors

The Spearman correlation analysis (Fig. 4) indicates that WT was positively correlated with Chl-a only at the MD station in the DR (r = 0.18, P < 0.001) and did not show a significant correlation with Chl-a at QJ and TC stations. DO was less correlated with Chl-a, as it was negatively correlated at the MD station (r = − 0.08, P < 0.01), positively correlated at the QJ station (r = 0.14, P < 0.001), and uncorrelated at TC station. At all DR stations, both pH and Tur were highly and significantly positively correlated with Chl-a. Similarly, NH3-N, TN, and CODMn were all positively correlated with Chl-a but had smaller correlation coefficients (0.10 to 0.37) compared with those of pH value and Tur. In contrast, TP and Chl-a were significantly negatively correlated (r = − 0.29 to − 0.07, P < 0.01).Figure 4 Spearman correlation matrix of six water quality monitoring stations in the DJKR from 2018 to 2023.

In the HR, WT was negatively correlated with Chl-a only at BQ station (r = − 0.07, P < 0.01), and was not directly correlated with Chl-a at QS and LK stations. DO was weakly correlated with Chl-a, with positive correlations at LK and BQ stations (r = 0.10 to 0.24, P < 0.001) and no direct correlation at the QS station. Tur was highly positively correlated with Chl-a at all HR stations. Similarly, NH3-N and CODMn were both positively correlated with Chl-a but had smaller correlation coefficients (r = 0.06 to 0.28) than those between Tur and Chl-a. TP was significantly positively correlated with Chl-a at the QS station (r = 0.11, P < 0.001), negatively correlated at the LK station (r = − 0.12, P < 0.001), and not directly correlated at the BQ station. TN and Chl-a were significantly positively correlated at BQ and LK stations (r = 0.13, P < 0.001) and no correlated at the QS station.

Overall, the correlations between Chl-a and the water quality factors varied among the monitoring stations. The Chl-a concentrations exhibited low correlations with the individual factors. However, the correlations with most of the water quality factors suggest that changes in Chl-a concentrations were the result of the combined effects of multiple factors.

Key factor determination and downscaling

The stepwise multiple regression analysis results indicate that at least five key driving factors at each station explained the Chl-a concentration variations (Table S4). These results verify that the changes in Chl-a concentrations were synergistically driven by multiple water quality parameters. For instance, the key factors at MD and QJ stations were CODMn, WT, TP, DO, and pH; Tur, DO, TN, TP, WT, NH3-N, and CODMn, respectively.

To reduce the number of variables, the key drivers were identified using PCA prior to SEM calculation (Table S5 and Table 2). The results of the KMO and Bartlett’s sphericity tests were presented in Table S5. Some variables had KMO values below 0.5. Usually, the KMO value above 0.5 is typically required for conducting PCA. However, these variables with low KMO values all passed Bartlett’s Test of Sphericity, with significance levels less than 0.05, indicating that the correlation matrix significantly deviates from the identity matrix and is still suitable for PCA. Therefore, even though some variables had low KMO values, the principal components obtained through PCA are still meaningful. Table 2 Principal component analysis results for the monitoring stations and sub-reservoirs in the DJKR.

Stations/ Sub-area reservoirs	Factors	Factors of dimensionality reduction	Principal component	Eigenvalue	Percentage of explained variance (%)	Percentage of cumulative explained variance (%)	
MD	Water physicochemical indicators	WT, DO, and pH	1	1.826	60.868	60.868	
2	0.898	29.939	90.808	
3	0.276	9.192	100.000	
QJ	Water physicochemical indicators	Tur, DO, and WT	1	1.787	59.565	59.565	
2	0.966	32.193	91.758	
3	0.247	8.242	100.000	
Nutrients	TN, TP, and NH3-N	1	1.417	47.234	47.234	
2	1.016	33.861	81.094	
3	0.567	18.906	100.000	
TC	Water physicochemical indicators	Tur, pH, and WT	1	1.399	46.649	46.649	
2	0.944	31.481	78.130	
3	0.656	21.870	100.000	
DR	Water physicochemical indicators	WT, DO, Tur, and pH	1	1.955	49.887	49.887	
2	1.050	26.259	76.146	
3	0.816	20.407	96.553	
4	0.138	3.447	100.000	
Nutrients	TN, TP, and NH3-N	1	1.439	47.974	47.974	
2	0.981	32.692	80.665	
3	0.580	19.335	100.000	
QS	Water physicochemical indicators	pH, Tur, and DO	1	1.257	41.902	41.902	
2	0.968	32.272	74.174	
3	0.775	25.826	100.000	
Nutrients	NH3-N and TP	1	1.494	74.690	74.690	
2	0.506	25.310	100.000	
LK	Water physicochemical indicators	Tur, DO, and WT	1	1.605	53.510	53.510	
2	1.013	33.780	87.290	
3	0.381	12.710	100.000	
BQ	Water physicochemical indicators	DO, Tur, pH, and WT	1	1.650	41.251	41.251	
2	1.080	27.006	68.257	
3	0.823	20.573	88.830	
4	0.447	11.170	100.000	
Nutrients	TN and NH3-N	1	1.317	65.825	65.825	
2	0.683	34.175	100.000	
HR	Water physicochemical indicators	WT, DO, Tur, and pH	1	1.806	45.140	45.140	
2	1.074	26.862	72.003	
3	0.833	20.827	92.829	
4	0.287	7.171	100.000	
Nutrients	TN, TP, and NH3-N	1	1.430	47.672	47.672	
2	1.186	39.543	87.214	
3	0.384	12.786	100.000	

The results of the downscaling were presented in Table 2. At the MD station, the water physicochemical indicators were determined to be PC1. At the QJ station, the water physicochemical indicators were identified as PC1, and the nutrients were selected as PC1 and PC2. At the TC station, the water physicochemical indicators were identified as PC1. At the QS station, the water physicochemical indicators and nutrients were selected as PC1. At the LK station, the water physicochemical indicators were selected as PC1 and PC2. At the BQ station, the water physicochemical indicators were selected as PC1 and PC2, and the nutrients were identified as PC1.

Driving pathway identification of algal proliferation based on SEM

Different SEMs were constructed and their fits were evaluated using Chl-a and driving factors after downscaling. The evaluation criteria are presented in Table S6. The goodness-of-fit results indicate that all of the SEMs passed the corresponding criteria evaluations (χ2/df < 3, CFI > 0.9, GFI > 0.9, AGFI > 0.9, RMSEA < 0.05, and SRMR < 0.05). Therefore, the models were well-fitted and met the requirements of the study. The SEMs for each station are shown in Fig. 5. All insignificant pathways were removed from the SEMs, and only the significant relationships (P < 0.05) were retained. The factor loading of a pathway between two variables is the direct effect, whereas the indirect effect is the product of the factor loadings on all associated pathways (Table 3 and Fig. 6). The total effect between two variables is the sum of the direct and indirect effects (Fig. 7).Figure 5 Environmental factors affecting Chl-a variations at different water quality monitoring stations and sub-reservoirs based on structural equation models (SEMs) (Note: The variable at the beginning of each arrow is the independent variable, while the variable at the arrow endpoint is the response variable for that independent variable. Data on the arrows are the factor loadings for the variable pairs. (a) MD, (b) QJ, and (c) TC stations are located in the Dan Reservoir (DR). (d) Overall conditions in the DR. (e) QS, (f) LK, and (g) BQ stations are located in the Han Reservoir (HR). (h) Overall conditions in the HR).

Table 3 Driving effects of Chl-a variations at different monitoring stations and sub-reservoirs in the DJKR.

Stations/ Sub-area reservoirs	Variables	Direct effect	Indirect effect	Total effect	
MD	PChm	0.15	0.01	0.16	
TP	− 0.12	–	− 0.12	
CODMn	0.29	0.02	0.31	
QJ	PChm	–	− 0.01	− 0.01	
Nutr 1	0.11	–	0.11	
Nutr 2	0.26	− 0.01	0.25	
CODMn	0.05	0.11	0.16	
TC	PChm	0.27	0.06	0.33	
TP	− 0.38	–	− 0.38	
CODMn	0.23	−0.11	0.12	
DR	PChm 1	0.19	–	0.19	
PChm 2	− 0.40	− 0.09	− 0.49	
Nutr	− 0.27	–	− 0.27	
CODMn	0.16	0.02	0.18	
QS	PChm	0.18	− 0.07	0.11	
Nutr	0.26	–	0.26	
CODMn	− 0.14	0.03	− 0.11	
LK	PChm 1	–	0.02	0.02	
PChm 2	− 0.38	− 0.03	− 0.41	
TP	− 0.11	–	− 0.11	
CODMn	0.15	0.02	0.17	
BQ	PChm 1	− 0.20	− 0.03	− 0.23	
PChm 2	–	0.02	0.02	
Nutr	0.20	–	0.20	
CODMn	0.20	0.02	0.22	
HR	PChm 1	− 0.11	0.02	− 0.09	
PChm 2	–	− 0.01	− 0.01	
Nutr 1	0.12	–	0.12	
Nutr 2	0.27	–	0.27	
CODMn	–	0.04	0.04	

Figure 6 Direct, indirect, and total effects of driving factors on Chl-a variations at the water quality monitoring stations and sub-reservoirs based on the SEMs (Note: direct, indirect, and total effects were calculated using standardized path coefficients and are represented by red, blue, and green bars, respectively).

Figure 7 Relative contributions of different environmental driving factors to Chl-a variations at the water quality monitoring stations and sub-reservoirs (Note: a, b, c, e, f, and g are monitoring stations. d and h are sub-reservoirs).

In the DR, the SEMs for MD and TC stations had the same driving pathways and key factors, indicating that PChm and CODMn directly affected Chl-a concentrations and simultaneously had indirect effects on Chl-a via the mediator variable TP. However, the environmental driving pathway at QJ station was relatively complex, wherein the CODMn had direct and indirect effect on Chl-a by influencing two nutrient variables. In the HR, the relationships between the driving factors and Chl-a were similar to those in the DR (i.e., PChm and CODMn had both direct and indirect effects on Chl-a by influencing mediator nutrient variables). Figure 5 shows the overall scenarios in the DR and HR, respectively. The main difference between the sub-reservoirs was that CODMn in the DR had a direct effect on Chl-a, whereas that in the HR did not directly affect Chl-a.

In the DR (Fig. 6), the total effects of PChm, TP, and CODMn on Chl-a at the MD station were 0.16, − 0.12, and 0.31, respectively. The total effects of PChm, Nutr 1, Nutr 1, and CODMn on Chl-a at the QJ station were − 0.01, 0.11, 0.25, and 0.16, respectively. The total effects of PChm, TP, and CODMn on Chl-a at the TC station were 0.33, − 0.38, and 0.12, respectively. In the HR, the total effects of PChm, Nutr, and CODMn on Chl-a at the QS station were 0.11, 0.26, and − 0.11, respectively. The total effects of PChm 1, PChm 2, TP, and CODMn on Chl-a at the LK station were 0.02, − 0.41, − 0.11, and 0.17, respectively. The total effects of PChm 1, PChm 2, Nutr, and CODMn on Chl-a at the BQ station were − 0.23, 0.02, 0.20, and 0.22, respectively.

In the DR, the factors driving Chl-a were PChm 1, PChm 2, Nutr, and CODMn, which contributed 16.81%, 43.36%, 23.89%, and 15.93%% to Chl-a variations, respectively. Among them, water physicochemical indicators (PChm 1 and PChm 2) were the dominant factors (total contribution of 60.18%) driving Chl-a variations in the DR. In the HR, the driving factors of Chla were PChm 1, PChm 2, Nutr 1, Nutr 2, and CODMn, which contributed 16.98%, 1.89%, 22.64%, 50.94%, and 7.55% to Chl-a variations, respectively. Of which, nutrient factors (Nutr 1 and Nutr 2) were the dominant factors (total contribution of 73.58%) driving Chl-a variations in the HR.

Chl-a variation risks at TC station

The effects of the different sub-reservoirs on Chl-aTC concentrations were identical, i.e., the probability of exceeding a specific Chl-aTC concentration increased with increasing Chl-a concentrations in other parts of the DJKR (Fig. 8). Chl-a at MD station ranged from 0.21 to 8.20 μg/L, and the probabilities of Chl-aTC exceeding 3.4 μg/L in the four quartiles were 0.72%, 2.53%, 8.40%, and 22.92%, respectively. The probabilities of Chl-aTC exceeding 3.4 μg/L when the QJ station Chl-a was under the same conditions were 0.53%, 2.25%, 8.30%, and 23.77%, respectively. The probabilities of Chl-aTC exceeding 3.4 μg/L when the QS station Chl-a was under the same conditions were 0.72%, 2.70%, 8.79%, and 23.65%, respectively. The probabilities of Chl-aTC exceeding 3.4 μg/L when the LK station Chl-a was under the same conditions were 0.58%, 2.29%, 8.50%, and 24.32%, respectively. The probabilities of Chl-aTC exceeding 3.4 μg/L when the BQ station Chl-a was under the same conditions were 2.09%, 4.91%, 9.43%, and 17.56%, respectively.Figure 8 Conditional risk probability of Chl-a concentration variations in the TC station under the influence of other stations (Note: a, b, c, d, and e represent the effects of Chl-a concentration variations in MD, QJ, QS, LK, and BQ stations on TC station respectively).

The probabilities that TC station would be affected by an algal bloom were largest for the uppermost quartile to maximum Chl-a in the sub-reservoirs, which ranged from 17 to 25%. For the rest of the concentration ranges, the probabilities of algal bloom risk at TC station did not exceed 10%. The probability of actual Chl-a concentrations exceeding 3.4 μg/L at TC station during 2018–2023 was 7.17%, which is similar to predictions for the sub-reservoirs when the Chl-a was in the 50%–75% quartile range (8.30–9.43%), which were less than 10%.

Discussion

Chl-a concentrations were relatively lower in the DR, particularly at TC station, whereas higher concentrations were more frequently observed in the HR, which is consistent with the findings of previous studies46. The algal bloom events signify that the planktonic algal biomass has peaked, and the Chl-a concentration is an important indicator for the degree of algal bloom and eutrophication. Therefore, comparing Chl-a maxima can reflect the intensity and scale of potential algal blooms in the two sub-reservoirs. Chl-a concentrations at QS and BQ stations exceeded 10 μg/L (11.60 and 11.38 μg/L, respectively). Thus, the algae were more aggregated in the HR and the risk of algal blooms was higher than that in the DR. This spatial difference in Chl-a concentration may be related to the differing lengths of the reservoir backwaters and the reservoir shapes. The DR is a lake-type reservoir with an expansive water surface and short shoreline, whereas the HR area is a canyon-type reservoir with a long and narrow reservoir surface, many tributaries, and zigzag-shaped banks that receive more exogenous nutrients from surface runoff14. Therefore, water in the HR is more prone to algal proliferation and blooms. The high Chl-a concentrations observed at QS station may be related to higher TP and NH3-N inputs from the Han River46. Chl-a concentrations at BQ station were also higher, which can be attributed to nutrient retention at the dam and lower flow rates, which resulted in Chl-a enrichment47. Therefore, the tributary inflow and dam in the HR should be considered when conducting algal risk management in the DJKR.

Algal proliferation requires a suitable environment, including optimal temperature, pH value, and light/radiation conditions, as well as appropriate dissolved oxygen concentration status, organic pollutant concentrations, and primary nutrients (carbon, nitrogen, and phosphorus) loadings for proliferation and survival48–51.

The composite variable Nutr reflected information from multiple nutrient indicators52. Nitrogen and phosphorus are key nutrients that affect algal proliferation53. Different TN and TP levels affect spatiotemporal phytoplankton dynamics and species succession54, and the phytoplankton biomass exhibits spatiotemporal heterogeneity with changing TN and TP concentrations55. A previous case study has shown that changes in nitrogen and phosphorus affected microcystis blooms in Isahaya Bay Reservoir, Japan, and additional nutrients provided by precipitation events promoted algal proliferation and led to algal blooms in the reservoir56. The water quality in the DJKR is generally good13; however, the TN exceeds the standard and is between Class III and IV57. TN concentrations recorded during the monitoring period were consistent with those obtained by previous investigations, suggesting that the nitrogen salts in the DJKR have always been adequate. Higher TN or TP concentrations lead to higher algal densities when they exceed the optimal range for algal development58. Thus, the adequate TN concentrations were responsible for the significant positive correlation observed between TN and Chl-a in this study. Additionally, algal biomass may be limited by TN concentrations less than 1.0 mg/L or TP concentrations less than 0.025 mg/L59. TP concentrations were below 0.025 mg/L during the monitoring period, which may have limited algal proliferation. Previous studies have determined that phosphorus is one of the main factors that influences planktonic algae in the DJKR10. Among the nutrient factors at MD, TC, and LK stations, the only selected factor was TP. Based on the SEMs for these stations, TP had a negative effect on Chl-a, whereas the SEMs for the other stations indicate that multiple trophic factors had a positive effect on Chl-a. This indicates that the synergistic effects of TP, TN, and NH3-N in this study were beneficial for algal proliferation. NH3-N is the preferred inorganic nitrogen utilized by planktonic algae60. A combination of NH3-N and phosphorus also promotes phytoplankton growth61,62. The synergistic effects of TN, TP, and NH3-N also showed positive impacts on Chl-a in Xiannv Lake, China63, which is consistent with what we observed in the HR.

Nutrient concentrations influence algal proliferation, but interactions also occur with water physicochemical indicators64. The correlation analysis indicates that WT was not well correlated with Chl-a, which may have an indirect effect on Chl-a by modulating other factors65. However, the stepwise multiple regression analysis results indicate that WT was a selected driving factor at all stations except for QS station. Therefore, we hypothesized that other factors may have more significant direct effects on Chl-a than WT66. Previous studies have shown that DO affects phytoplankton biomass in the DJKR67. The correlation analysis indicates that WT and DO are highly negatively correlated, but together they reflect the primary environmental conditions required for algal proliferation. Therefore, DO was considered to be a key driver of algal proliferation in this study. Unlike common phenomena, this study found a low correlation between DO and Chla, while TP was generally negatively correlated with Chla. It may be related to the phytoplankton community in DJKR. Previous studies have shown that the dominants of the DJKR are Cyanophyta, Chlorophyta, and Bacillariophyta10. The DO and TP are significantly negatively correlated with Cyanophyta and Chlorophyta, significantly positively correlated with Bacillariophyta, and exhibit a negative and weaker correlation with the phytoplankton diversity index10. Chl-a is an indicator of the biomass of all planktonic algae in the water. When considering the entire phytoplankton community, the effects on individual communities may be neutralized, resulting in the correlations between DO, TP, and Chla as observed in this study.

Tur also affects algal proliferation, as low turbidity can provide suitable light conditions68. When an algal bloom event occurs, it reduces the water clarity and increases turbidity, which occurred at the QS station in the autumn (both Chl-a and Tur concentrations were relatively high). Tur recorded in this study were highly positively correlated with Chl-a, and was selected in the stepwise multiple regression models at all monitoring stations except MD. generally as the first or second factor. Therefore, among the physicochemical factors affecting Chl-a concentration variations in the DJKR, Tur was one of the key driving factors69.

As one of the driving factors, pH value is one of the components of the water physicochemical factor (PChm) in the SEMS established for the MD, TC, and QS stations. The PChm arrow in the model points to Chl-a, which reflects the significant and positive active role of pH in influencing Chl-a. During the monitoring period, the mean pH value range was 8.08 to 8.21, which is within the optimal range. The effects of pH and algae are reciprocal in that suitable pH facilitates algal blooms, while photosynthesis by the algae absorbs CO2 in the water, thereby changing the pH value70,71. Therefore, pH value was included as a key driving factor.

The SEM results in this study showed that CODMn was one of the key factors driving algal proliferation in the DJKR. Previous studies showed that the water quality of the MRSNWDPC was good, with low CODMn concentration and a significant correlation between CODMn and planktonic algal cell density10,72. Some algal proliferation decreases with increasing organic pollutant concentrations, whereas some algae survive and degrade at lower concentrations, producing nutrient sources in the degradation products that can promote algal growth73. Specifically, at low organic matter concentrations, organic matter can be used as a direct nutrient source and degradation product to promote algal growth by increasing enzymatic activity and lipid peroxidation in algal cells74. At higher organic matter concentrations, algal growth is inhibited by damage to cell membranes and reduced enzymatic activity75. In addition, planktonic algal growth also releases organic matter, which affects organic pollutant concentrations76.

Hydrodynamic conditions are likely to be a significant factor influencing algal growth. Research indicates that the DJKR experiences notable thermal stratification, with the thermocline situated between 2 and 15 m below the surface, averaging around 10 m in thickness, and the stratification influences the distribution of nitrogen and phosphorus, thereby impacting the community structure of phytoplankton77. Heavy rainfall can severely disrupt this stratification, resulting in increased dissolved oxygen levels that inhibit nutrient release from sediments and consequently reduce phytoplankton biomass78. Additionally, the reservoir’s phytoplankton diversity index is higher in spring and winter and lowest in autumn28, likely due to variations in flow velocity and volume. Historical data reveal that the monthly average of DJKR inflow is highest from July to October, peaking at 2410 m3/s. The intake flow velocity and volume are greatest in autumn and lowest in winter, with annual averages of 0.76 m/s and 203.48 m3/h, respectively28. High flow velocity conditions negatively affect algal growth79. If hydrodynamic conditions are incorporated as a factor in constructing a structural equation model, the initial model should illustrate that hydrodynamic factors directly influence physicochemical factors, nutrient factors, organic pollution factors, and phytoplankton, and indirectly influence phytoplankton through physicochemical, nutrient, and organic pollution factors. The causal pathways among the remaining three factors should remain consistent with those observed in this study. A thorough elucidation of the driving mechanisms of planktonic algae necessitates considering additional factors, which will be addressed in future research.

Conclusions

This study developed a methodological framework for analyzing spatiotemporal distributions of algal proliferation and quantifying the driving mechanism using multivariate statistical techniques and structural equation models. The Danjiangkou Reservoir, i.e., the core water source area of the South-to-North Water Diversion Project of China, was taken as a study case. The risks of Chl-a concentration changes at the TC station were analyzed using the Copula function and Bayesian theory. The main conclusions are as follows:Chl-a concentration distributions exhibited significant spatiotemporal variations in the DJKR, and Chl-a was relatively higher in spring than in other seasons. The factors driving algal proliferation differed spatially. Water physicochemical factors and CODMn had both direct and indirect effects on Chl-a by influencing mediator nutrient variables, whereas nutrient factors had direct effects on Chl-a variations.

In the DR, the dominant factor driving Chl-a was the water physicochemical factors, whereas nutrients were the dominant factor in the Han reservoir. When Chl-a concentrations in other areas of the DJKR ranged from the minimum to the upper quartile, the probabilities of the TC station Chl-a exceeded 3.4 μg/L (owing to the influences of other parts of the DJKR) were all less than 10%. Thus, the risk of planktonic algal proliferation at the MRSNWDPC intake is currently low.

Due to limitations in data and numerical models, this study does not account for the effects of climate and hydrodynamic factors on algae, nor the impacts of pollutant input and biological migration under different hydrological conditions. Future research should comprehensively consider the influence of climate/hydrodynamic factors on algae proliferation, establish a mathematical model of mechanism, and quantitatively analyze the environmental impact of pollutant input and algae migration under different hydrological conditions so as to make a more reliable and informative comparative analysis of the relevant statistical models and achievements.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70408-4.

Acknowledgements

The authors would like to thank the editors and the reviewers for their valuable suggestions and contributions which significantly improved this article, and would also need to acknowledge the Construction and Administration Bureau of the Middle Route of the South-to-North Water Diversion Project of China that supported the data collection.

Author contributions

Dr. Xizhi Nong: Data collection, Investigation, Conceptualization, Formal analysis, Writing-Original Draft, Draft review & editing, Funding Acquisition; Mr. Xian Guan: Methodology, Software, Writing-Original Draft; Prof. Lihua Chen: Formal analysis, Draft review & editing; Prof. Jiahua Wei: Supervision, Draft review & editing; Prof. Ronghui Li: Funding Acquisition, Resources, Draft review & editing.

Funding

This research was funded by the Specific Research Project of Guangxi for Research Bases and Talents (No. AD22035185), the National Natural Science Foundation of China (Nos. 52309016 & 52079033), the Belt and Road Special Foundation of the National Key Laboratory of Water Disaster Prevention (No. 2022nkms06).

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Consent to publication

The authors authorize the publication of the article in this journal.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Greve P Global assessment of water challenges under uncertainty in water scarcity projections Nat. Sustain. 2018 1 486 494 10.1038/s41893-018-0134-9
Greve, P. et al. Global assessment of water challenges under uncertainty in water scarcity projections. Nat. Sustain. 1, 486–494. 10.1038/s41893-018-0134-9 (2018).
2. Li B Li J An G Zhao C Wang C Long-term and strong suppression against Microcystis growth and microcystin-release by luteolin continuous-release microsphere: Optimal construction, characterization, effects and proteomic mechanisms Water Res. 2021 10.1016/j.watres.2021.117448 34986458
Li, B., Li, J., An, G., Zhao, C. & Wang, C. Long-term and strong suppression against Microcystis growth and microcystin-release by luteolin continuous-release microsphere: Optimal construction, characterization, effects and proteomic mechanisms. Water Res.10.1016/j.watres.2021.117448 (2021).34986458
3. Boretti A Rosa L Reassessing the projections of the world water development report NPJ Clean Water 2019 10.1038/s41545-019-0039-9
Boretti, A. & Rosa, L. Reassessing the projections of the world water development report. NPJ Clean Water10.1038/s41545-019-0039-9 (2019).
4. Koncagül, E. & Connor, R. The United Nations World Water Development Report 2023: partnerships and cooperation for water; facts, figures and action examples. (2023).
5. Shen H Cai Q Zhang M Spatial gradient and seasonal variation of trophic status in a large water supply reservoir for the South-to-North Water Diversion Project China J. Freshwater Ecol. 2015 30 249 261 10.1080/02705060.2014.935748
Shen, H., Cai, Q. & Zhang, M. Spatial gradient and seasonal variation of trophic status in a large water supply reservoir for the South-to-North Water Diversion Project China. J. Freshwater Ecol. 30, 249–261. 10.1080/02705060.2014.935748 (2015).
6. Zhang X Xiong J Cheng J Yao Z Chen Y Application of fuzzy mathematics for evaluation of eutrophication in danjiangkou reservoir Environ. Monit. China 2017 33 99 105
Zhang, X., Xiong, J., Cheng, J., Yao, Z. & Chen, Y. Application of fuzzy mathematics for evaluation of eutrophication in danjiangkou reservoir. Environ. Monit. China 33, 99–105 (2017).
7. Voeroesmarty CJ Global threats to human water security and river biodiversity Nature 2010 467 555 561 10.1038/nature09440 20882010
Voeroesmarty, C. J. et al. Global threats to human water security and river biodiversity. Nature 467, 555–561. 10.1038/nature09440 (2010).20882010
8. Xu H Environmental controls of harmful cyanobacterial blooms in Chinese inland waters Harmful Algae 2021 10.1016/j.hal.2021.102127 34887007
Xu, H. et al. Environmental controls of harmful cyanobacterial blooms in Chinese inland waters. Harmful Algae10.1016/j.hal.2021.102127 (2021).34887007
9. Montefiore LR Downstream nutrient concentrations depend on watershed inputs more than reservoir releases in a highly engineered watershed Water Resour. Res. 2024 10.1029/2023wr035590
Montefiore, L. R. et al. Downstream nutrient concentrations depend on watershed inputs more than reservoir releases in a highly engineered watershed. Water Resour. Res.10.1029/2023wr035590 (2024).
10. Xiao Y Spatio-temporal variation of phytoplankton community and its relationship with environmental factors in Danjiangkou Reservoir J. Lake Sci. 2023 35 821 832 10.18307/2023.0306
Xiao, Y. et al. Spatio-temporal variation of phytoplankton community and its relationship with environmental factors in Danjiangkou Reservoir. J. Lake Sci. 35, 821–832 (2023).
11. Ma F A Bayesian method for comprehensive water quality evaluation of the Danjiangkou Reservoir water source area, for the middle route of the south-to-north water diversion project in China Front. Earth Sci. 2014 8 242 250 10.1007/s11707-013-0395-6
Ma, F. et al. A Bayesian method for comprehensive water quality evaluation of the Danjiangkou Reservoir water source area, for the middle route of the south-to-north water diversion project in China. Front. Earth Sci. 8, 242–250. 10.1007/s11707-013-0395-6 (2014).
12. Tan X Xia XL Li SY Zhang QF Water quality characteristics and integrated assessment based on multistep correlation analysis in the danjiangkou reservoir, China. J. Environ. Inf. 2015 25 60 70 10.3808/jei.201500296
Tan, X., Xia, X. L., Li, S. Y. & Zhang, Q. F. Water quality characteristics and integrated assessment based on multistep correlation analysis in the danjiangkou reservoir, China.. J. Environ. Inf. 25, 60–70. 10.3808/jei.201500296 (2015).
13. Zhang C Nong X Shao D Chen L An integrated risk assessment framework using information theory-based coupling methods for basin-scale water quality management: A case study in the Danjiangkou Reservoir Basin, China Sci Total Environ 2023 10.1016/j.scitotenv.2023.163731 38163615
Zhang, C., Nong, X., Shao, D. & Chen, L. An integrated risk assessment framework using information theory-based coupling methods for basin-scale water quality management: A case study in the Danjiangkou Reservoir Basin, China. Sci Total Environ10.1016/j.scitotenv.2023.163731 (2023).38163615
14. Tan X Xia X-L Cheng X-L Zhang Q-F Temporal and spatial pattern of phytoplankton community and its biodiversity indices in the Danjiangkou Reservoir Huanjing Kexue 2011 32 2875 2882 22279895
Tan, X., Xia, X.-L., Cheng, X.-L. & Zhang, Q.-F. Temporal and spatial pattern of phytoplankton community and its biodiversity indices in the Danjiangkou Reservoir. Huanjing Kexue 32, 2875–2882 (2011).22279895
15. Pan Y Effects of water level increase on Phytoplankton assemblages in a drinking water reservoir Water 2018 10.3390/w10030256
Pan, Y. et al. Effects of water level increase on Phytoplankton assemblages in a drinking water reservoir. Water10.3390/w10030256 (2018).
16. Xiong M Spatio-temporal characteristics of the phytoplankton community and assessment of fish productivity in the Danjiangkou Reservoir, the water source for the South-to-North Water Diversion Project, China J. Fish. Sci. China 2021 28 715 727
Xiong, M. et al. Spatio-temporal characteristics of the phytoplankton community and assessment of fish productivity in the Danjiangkou Reservoir, the water source for the South-to-North Water Diversion Project, China. J. Fish. Sci. China 28, 715–727 (2021).
17. Liu J Tang Y Tian K Study on evaluation of nutrition degree and eutrophication prevention and control in Danjiangkou reservoir China Resour. Comprehens. Util. 2022 40 189 191
Liu, J., Tang, Y. & Tian, K. Study on evaluation of nutrition degree and eutrophication prevention and control in Danjiangkou reservoir. China Resour. Comprehens. Util. 40, 189–191 (2022).
18. Pan X Distribution characteristics of carbon, nitrogen, and phosphorus and pollution load estimation of sediments in Danjiangkou Reservoir, China Water 2023 10.3390/w15162885
Pan, X. et al. Distribution characteristics of carbon, nitrogen, and phosphorus and pollution load estimation of sediments in Danjiangkou Reservoir, China. Water10.3390/w15162885 (2023).
19. Wan Y Zhang L Fu X Jin H Assessment of eutrophication degree and its prevention in Danjiangkou Reservoir J. Beijing Normal Univ. Nat. Sci. 2020 56 275 281
Wan, Y., Zhang, L., Fu, X. & Jin, H. Assessment of eutrophication degree and its prevention in Danjiangkou Reservoir. J. Beijing Normal Univ. Nat. Sci. 56, 275–281 (2020).
20. Sun R Wei J Zhang S Pei H The dynamic changes in phytoplankton and environmental factors within Dongping Lake (China) before and after the South-to-North Water Diversion Project Environ. Res. 2024 10.1016/j.envres.2024.118138 39288827
Sun, R., Wei, J., Zhang, S. & Pei, H. The dynamic changes in phytoplankton and environmental factors within Dongping Lake (China) before and after the South-to-North Water Diversion Project. Environ. Res.10.1016/j.envres.2024.118138 (2024).39288827
21. Liang L Deng Y Wang W Zhou S Zhang L Influences of lower pH on phytoplankton growth in alkaline lakes after water transfer: Insights from a coupled hydrodynamic–algal ecological model and experimental analysis Environ. Res. 2024 257 119201 10.1016/j.envres.2024.119201 38782337
Liang, L., Deng, Y., Wang, W., Zhou, S. & Zhang, L. Influences of lower pH on phytoplankton growth in alkaline lakes after water transfer: Insights from a coupled hydrodynamic–algal ecological model and experimental analysis. Environ. Res. 257, 119201 (2024).38782337
22. Dory F Interaction between temperature and nutrients: How does the phytoplankton community cope with climate change? Sci. Total Environ. 2024 10.1016/j.scitotenv.2023.167566 37802360
Dory, F. et al. Interaction between temperature and nutrients: How does the phytoplankton community cope with climate change?. Sci. Total Environ.10.1016/j.scitotenv.2023.167566 (2024).37802360
23. Çamdevyren H Demyr N Kanik A Keskyn S Use of principal component scores in multiple linear regression models for prediction of Chlorophyll-a in reservoirs Ecol. Modell. 2005 181 581 589 10.1016/j.ecolmodel.2004.06.043
Çamdevyren, H., Demyr, N., Kanik, A. & Keskyn, S. Use of principal component scores in multiple linear regression models for prediction of Chlorophyll-a in reservoirs. Ecol. Modell. 181, 581–589. 10.1016/j.ecolmodel.2004.06.043 (2005).
24. Paerl HW Controlling harmful cyanobacterial blooms in a hyper-eutrophic lake (Lake Taihu, China): The need for a dual nutrient (N & P) management strategy Water Res. 2011 45 1973 1983 10.1016/j.watres.2010.09.018 20934736
Paerl, H. W. et al. Controlling harmful cyanobacterial blooms in a hyper-eutrophic lake (Lake Taihu, China): The need for a dual nutrient (N & P) management strategy. Water Res. 45, 1973–1983. 10.1016/j.watres.2010.09.018 (2011).20934736
25. Nelson NG Munoz-Carpena R Phlips E Parameter uncertainty drives important incongruities between simulated chlorophyll-a and phytoplankton functional group dynamics in a mechanistic management model Environ. Modell. Software 2020 10.1016/j.envsoft.2020.104708
Nelson, N. G., Munoz-Carpena, R. & Phlips, E. Parameter uncertainty drives important incongruities between simulated chlorophyll-a and phytoplankton functional group dynamics in a mechanistic management model. Environ. Modell. Software10.1016/j.envsoft.2020.104708 (2020).
26. Alvarez S Lupi F Solis D Thomas M Valuing provision scenarios of coastal ecosystem services: The case of boat ramp closures due to harmful algae blooms in Florida Water 2019 10.3390/w11061250
Alvarez, S., Lupi, F., Solis, D. & Thomas, M. Valuing provision scenarios of coastal ecosystem services: The case of boat ramp closures due to harmful algae blooms in Florida. Water10.3390/w11061250 (2019).
27. Michalsen, M. M., Feist, S. M., VanZomeren, C. M., Seiter-Moser, J. M. & Boss, A. N. USACE freshwater harmful algal bloom research and development initiative. (2024).
28. Yan X Hydrologic and physicochemical factors co-drive seasonal changes of phytoplankton during dynamic water diversion processes in the Danjiangkou Reservoir J. Lake Sci. 2021 33 1350 1363 10.18307/2021.0505
Yan, X. et al. Hydrologic and physicochemical factors co-drive seasonal changes of phytoplankton during dynamic water diversion processes in the Danjiangkou Reservoir. J. Lake Sci. 33, 1350–1363 (2021).
29. Zhang C High-throughput sequencing of diatom community, its spatial and temporal variation and interrelationships with physicochemical factors in Danjiangkou Reservoir China. Water 2022 14 1609 10.3390/w14101609
Zhang, C. et al. High-throughput sequencing of diatom community, its spatial and temporal variation and interrelationships with physicochemical factors in Danjiangkou Reservoir. China. Water 14, 1609 (2022).
30. Zhu Y Land use impact on water quality and phytoplankton community structure in Danjiangkou Reservoir Diversity 2024 16 275 10.3390/d16050275
Zhu, Y. et al. Land use impact on water quality and phytoplankton community structure in Danjiangkou Reservoir. Diversity 16, 275 (2024).
31. Hayes AF Montoya AK Rockwood NJ The analysis of mechanisms and their contingencies: Process versus structural equation modeling Australas. Mark. J. 2017 25 76 81 10.1016/j.ausmj.2017.02.001
Hayes, A. F., Montoya, A. K. & Rockwood, N. J. The analysis of mechanisms and their contingencies: Process versus structural equation modeling. Australas. Mark. J. 25, 76–81 (2017).
32. Vinzi, V. E., Trinchera, L. & Amato, S. PLS path modeling: from foundations to recent developments and open issues for model assessment and improvement. Handbook of partial least squares: Concepts, methods and applications, 47–82 (2010).
33. Chen F Cropland carbon stocks driven by soil characteristics, rainfall and elevation Sci. Total Environ. 2023 10.1016/j.scitotenv.2022.160602 38163608
Chen, F. et al. Cropland carbon stocks driven by soil characteristics, rainfall and elevation. Sci. Total Environ.10.1016/j.scitotenv.2022.160602 (2023).38163608
34. Liu Q Qiao J Li M Huang M Spatiotemporal heterogeneity of ecosystem service interactions and their drivers at different spatial scales in the Yellow River Basin Sci. Total Environ. 2024 10.1016/j.scitotenv.2023.168486 39288879
Liu, Q., Qiao, J., Li, M. & Huang, M. Spatiotemporal heterogeneity of ecosystem service interactions and their drivers at different spatial scales in the Yellow River Basin. Sci. Total Environ.10.1016/j.scitotenv.2023.168486 (2024).39288879
35. Sheng Y Wang ZQ Feng X Potential ecological risk and zoning control strategies for heavy metals in soils surrounding core water sources: A case study from Danjiangkou Reservoir, China Ecotoxicol. Environ. Saf. 2023 10.1016/j.ecoenv.2023.114610 38000301
Sheng, Y., Wang, Z. Q. & Feng, X. Potential ecological risk and zoning control strategies for heavy metals in soils surrounding core water sources: A case study from Danjiangkou Reservoir, China. Ecotoxicol. Environ. Saf.10.1016/j.ecoenv.2023.114610 (2023).38000301
36. Guo XM Zhang QM Zhao TQ Jin C Fluxes, characteristics and influence on the aquatic environment of inorganic nitrogen deposition in the Danjiangkou reservoir Ecotoxicol. Environ. Saf. 2022 10.1016/j.ecoenv.2022.113814 38321688
Guo, X. M., Zhang, Q. M., Zhao, T. Q. & Jin, C. Fluxes, characteristics and influence on the aquatic environment of inorganic nitrogen deposition in the Danjiangkou reservoir. Ecotoxicol. Environ. Saf.10.1016/j.ecoenv.2022.113814 (2022).38321688
37. Zhang QZ The research on riverine hydrochemistry and controlling factors in the Danjiangkou Reservoir J. Radioanal. Nucl. Chem. 2020 324 507 519 10.1007/s10967-020-07096-7
Zhang, Q. Z. et al. The research on riverine hydrochemistry and controlling factors in the Danjiangkou Reservoir. J. Radioanal. Nucl. Chem. 324, 507–519. 10.1007/s10967-020-07096-7 (2020).
38. Li ZG Synergistic impacts of landscape, soil, and environmental factors on the spatial distribution of soil aggregates stability in the Danjiangkou reservoir area Catena 2024 10.1016/j.catena.2024.107840
Li, Z. G. et al. Synergistic impacts of landscape, soil, and environmental factors on the spatial distribution of soil aggregates stability in the Danjiangkou reservoir area. Catena10.1016/j.catena.2024.107840 (2024).
39. Jolliffe IT Cadima J Principal component analysis: A review and recent developments Philos. Trans. R. Soc. A 2016 10.1098/rsta.2015.0202
Jolliffe, I. T. & Cadima, J. Principal component analysis: A review and recent developments. Philos. Trans. R. Soc. A10.1098/rsta.2015.0202 (2016).
40. Grace JB Structural equation modeling and natural systems 2006 Cambridge Cambridge University Press
Grace, J. B. Structural equation modeling and natural systems (Cambridge University Press, Cambridge, 2006).
41. Yang J Dynamic distribution and driving mechanisms of antibiotic resistance genes in a human-intensive watershed Water Res 2022 10.1016/j.watres.2022.118841 36603307
Yang, J. et al. Dynamic distribution and driving mechanisms of antibiotic resistance genes in a human-intensive watershed. Water Res10.1016/j.watres.2022.118841 (2022).36603307
42. Lu K Driving mechanism of water replenishment on DOM composition and eutrophic status changes of lake in arid and semi-arid regions of loess area Sci. Total Environ. 2023 10.1016/j.scitotenv.2023.165609 38163598
Lu, K. et al. Driving mechanism of water replenishment on DOM composition and eutrophic status changes of lake in arid and semi-arid regions of loess area. Sci. Total Environ.10.1016/j.scitotenv.2023.165609 (2023).38163598
43. Yan B Pan Z Xue Y Duan M Modeling dependence and correlation in hydrological calculation J. Hydraul. Eng. 2017 48 1039 1046
Yan, B., Pan, Z., Xue, Y. & Duan, M. Modeling dependence and correlation in hydrological calculation. J. Hydraul. Eng. 48, 1039–1046 (2017).
44. Guo S Yan B Xiao Y Fang B Zhang N Multivariate hydrological analysis and estimation J. China Hydrol. 2008 28 3 1 7
Guo, S., Yan, B., Xiao, Y., Fang, B. & Zhang, N. Multivariate hydrological analysis and estimation. J. China Hydrol. 28(3), 1–7 (2008).
45. MEEPRC. Announcement on the Publication of the National Ecological Benchmark “Lake Nutrient Benchmark - East Central Lakes (Total Phosphorus, Total Nitrogen, Chlorophyll a)” (2020 Edition) and its Technical Report, https://www.mee.gov.cn/xxgk2018/xxgk/xxgk01/202012/t20201230_815561.html (2020).
46. Guo S Spatiotemporal characteristics of layered chlorophyll-a concentration and influencing factors in Danjiangkou Reservoir J. Lake Sci. 2021 33 366 376 10.18307/2021.0206
Guo, S. et al. Spatiotemporal characteristics of layered chlorophyll-a concentration and influencing factors in Danjiangkou Reservoir. J. Lake Sci. 33, 366–376 (2021).
47. Luo Y-F Spatial and temporal distribution of chlorophyll a and its relationship to algae and environmental factors in Aha reservoir Environ. Sci. 2017 38 4151 4159 10.13227/j.hjkx.201703177
Luo, Y.-F. et al. Spatial and temporal distribution of chlorophyll a and its relationship to algae and environmental factors in Aha reservoir. Environ. Sci. 38, 4151–4159. 10.13227/j.hjkx.201703177 (2017).
48. Liu C Effects of pH on growth and species changes of algae in freshwater J. Agro-environ. Sci. 2005 24 294 298
Liu, C. et al. Effects of pH on growth and species changes of algae in freshwater. J. Agro-environ. Sci. 24, 294–298 (2005).
49. Xu L Feng P Sun D Li F Numerical simulation for the effect of temperature on the algae growth J. Safety Environ. 2013 13 76 81
Xu, L., Feng, P., Sun, D. & Li, F. Numerical simulation for the effect of temperature on the algae growth. J. Safety Environ. 13, 76–81 (2013).
50. Tarabih OM Dang TD Paudel R Arias ME Lake operation optimization of nutrient exports: Application of phosphorus control in the largest subtropical lake in the United States Environ. Modell. Softw. 2023 10.1016/j.envsoft.2022.105603
Tarabih, O. M., Dang, T. D., Paudel, R. & Arias, M. E. Lake operation optimization of nutrient exports: Application of phosphorus control in the largest subtropical lake in the United States. Environ. Modell. Softw.10.1016/j.envsoft.2022.105603 (2023).
51. Alvarez S Brown CE Diaz MG Leary H Solis D Non-linear impacts of harmful algae blooms on the coastal tourism economy J. Environ. Manage. 2024 10.1016/j.jenvman.2023.119811 38581895
Alvarez, S., Brown, C. E., Diaz, M. G., Leary, H. & Solis, D. Non-linear impacts of harmful algae blooms on the coastal tourism economy. J. Environ. Manage.10.1016/j.jenvman.2023.119811 (2024).38581895
52. O'Leary H Alvarez S Bahja F What's in a name? Political and economic concepts differ in social media references to harmful algae blooms J Environ Manage 2024 10.1016/j.jenvman.2024.120799 39213843
O’Leary, H., Alvarez, S. & Bahja, F. What’s in a name? Political and economic concepts differ in social media references to harmful algae blooms. J Environ Manage10.1016/j.jenvman.2024.120799 (2024).39213843
53. Jiang M Nakano S-I The crucial influence of trophic status on the relative requirement of nitrogen to phosphorus for phytoplankton growth Water Res 2022 10.1016/j.watres.2022.118868 36603309
Jiang, M. & Nakano, S.-I. The crucial influence of trophic status on the relative requirement of nitrogen to phosphorus for phytoplankton growth. Water Res10.1016/j.watres.2022.118868 (2022).36603309
54. Nelson NG revealing biotic and abiotic controls of harmful algal blooms in a shallow subtropical lake through statistical machine learning Environ. Sci. Technol. 2018 52 3527 3535 10.1021/acs.est.7b05884 29478313
Nelson, N. G. et al. revealing biotic and abiotic controls of harmful algal blooms in a shallow subtropical lake through statistical machine learning. Environ. Sci. Technol. 52, 3527–3535. 10.1021/acs.est.7b05884 (2018).29478313
55. Liu X Spatiotemporal dynamics of succession and growth limitation of phytoplankton for nutrients and light in a large shallow lake Water Res 2021 10.1016/j.watres.2021.116910 35032895
Liu, X. et al. Spatiotemporal dynamics of succession and growth limitation of phytoplankton for nutrients and light in a large shallow lake. Water Res10.1016/j.watres.2021.116910 (2021).35032895
56. Umehara A Komorita T Takahashi T Tsutsumi H Estimation of production and sedimentation of cyanobacterial toxins (microcystin) based on nutrient budgets in the reservoir of Isahaya Bay, Japan Ecotoxicol. Environ. Saf. 2019 10.1016/j.ecoenv.2019.109477 31369939
Umehara, A., Komorita, T., Takahashi, T. & Tsutsumi, H. Estimation of production and sedimentation of cyanobacterial toxins (microcystin) based on nutrient budgets in the reservoir of Isahaya Bay, Japan. Ecotoxicol. Environ. Saf.10.1016/j.ecoenv.2019.109477 (2019).31369939
57. Dong G Hu Z Liu X Fu Y Zhang W Spatio-temporal variation of total nitrogen and ammonia nitrogen in the water source of the middle route of the South-to-North Water Diversion Project Water 2020 12 2615 10.3390/w12092615
Dong, G., Hu, Z., Liu, X., Fu, Y. & Zhang, W. Spatio-temporal variation of total nitrogen and ammonia nitrogen in the water source of the middle route of the South-to-North Water Diversion Project. Water 12, 2615 (2020).
58. Hu Y Zhang J Huang J Hu S Environmental drivers and aquatic ecosystem assessment of periphytic algae at inflow rivers in six lakes over the Yangtze River Basin Water 2022 14 2184 10.3390/w14142184
Hu, Y., Zhang, J., Huang, J. & Hu, S. Environmental drivers and aquatic ecosystem assessment of periphytic algae at inflow rivers in six lakes over the Yangtze River Basin. Water 14, 2184 (2022).
59. Zhu G Characteristics of diatom blooms in a reservoir-water supply area and the countermeasures in Taihu Basin, China J. Lake Sci. 2016 28 9 21 10.18307/2016.0102
Zhu, G. et al. Characteristics of diatom blooms in a reservoir-water supply area and the countermeasures in Taihu Basin, China. J. Lake Sci. 28, 9–21 (2016).
60. Herndon J Cochlan WP Nitrogen utilization by the raphidophyte Heterosigma akashiwo: Growth and uptake kinetics in laboratory cultures Harmful Algae 2007 6 260 270 10.1016/j.hal.2006.08.006
Herndon, J. & Cochlan, W. P. Nitrogen utilization by the raphidophyte Heterosigma akashiwo: Growth and uptake kinetics in laboratory cultures. Harmful Algae 6, 260–270. 10.1016/j.hal.2006.08.006 (2007).
61. Zeng Q Qin L Bao L Li Y Li X Critical nutrient thresholds needed to control eutrophication and synergistic interactions between phosphorus and different nitrogen sources Environ. Sci. Pollut. Res. 2016 23 21008 21019 10.1007/s11356-016-7321-x
Zeng, Q., Qin, L., Bao, L., Li, Y. & Li, X. Critical nutrient thresholds needed to control eutrophication and synergistic interactions between phosphorus and different nitrogen sources. Environ. Sci. Pollut. Res. 23, 21008–21019. 10.1007/s11356-016-7321-x (2016).
62. Zacharias Q Kaplan D The past, present, and potential future of phosphorus management in the Florida Everglades Restor. Ecol. 2023 10.1111/rec.13799
Zacharias, Q. & Kaplan, D. The past, present, and potential future of phosphorus management in the Florida Everglades. Restor. Ecol.10.1111/rec.13799 (2023).
63. Xia W Spatio-temporal dynamics of dissolved oxygen and its influencing factors in Lake Xiannv Jiangxi, China J. Lake Sci. 2023 35 874 885 10.18307/2023.0330
Xia, W. et al. Spatio-temporal dynamics of dissolved oxygen and its influencing factors in Lake Xiannv Jiangxi, China. J. Lake Sci. 35, 874–885 (2023).
64. Li NX Linking water environmental factors and the local watershed landscape to the chlorophyll a concentration in reservoir bays Sci. Total Environ. 2021 10.1016/j.scitotenv.2020.143617 34979230
Li, N. X. et al. Linking water environmental factors and the local watershed landscape to the chlorophyll a concentration in reservoir bays. Sci. Total Environ.10.1016/j.scitotenv.2020.143617 (2021).34979230
65. Lu Y Vertical distribution rules and factors influencing phytoplankton in front of a drinking water reservoir outlet Sci. Total Environ. 2023 10.1016/j.scitotenv.2023.166512 38163598
Lu, Y. et al. Vertical distribution rules and factors influencing phytoplankton in front of a drinking water reservoir outlet. Sci. Total Environ.10.1016/j.scitotenv.2023.166512 (2023).38163598
66. Shangpei XUN Yuanjian Y Binfang HE Wei FAN Hongqun Z Analysis of relationship between spring water temperature and chlorophyll-a in Lake Cha-ohu Sci. Limnol. Sin. 2011 23 767 772
Shangpei, X. U. N., Yuanjian, Y., Binfang, H. E., Wei, F. A. N. & Hongqun, Z. Analysis of relationship between spring water temperature and chlorophyll-a in Lake Cha-ohu. Sci. Limnol. Sin. 23, 767–772 (2011).
67. Wang Y Chen L Niu Y Yu H Luo M Spatio-temporal variation in phytoplankton community and its influencing factors in Danjiangkou Reservoir J. Lake Sci. 2016 28 1057 1065 10.18307/2016.0516
Wang, Y., Chen, L., Niu, Y., Yu, H. & Luo, M. Spatio-temporal variation in phytoplankton community and its influencing factors in Danjiangkou Reservoir. J. Lake Sci. 28, 1057–1065 (2016).
68. Zhou C Temporal and spatial distribution of environmental factors and Phytoplankton during algal bloom season in Pengxi river, Three Gorges Reservoir Environ. Sci. 2016 37 873 883
Zhou, C. et al. Temporal and spatial distribution of environmental factors and Phytoplankton during algal bloom season in Pengxi river, Three Gorges Reservoir. Environ. Sci. 37, 873–883 (2016).
69. Tang H Liu P Wu J Liao H He A Phytoplankton functional groups and their response to water physiochemical factors in Hongchaojiang reservoir in Guangxi Province J. Hydroecol. 2022 43 85 91
Tang, H., Liu, P., Wu, J., Liao, H. & He, A. Phytoplankton functional groups and their response to water physiochemical factors in Hongchaojiang reservoir in Guangxi Province. J. Hydroecol. 43, 85–91 (2022).
70. Hai XU Zhaopu LIU Lan Y Linzhang Y Effect of pH on growth of several freshwater algae Environ. Sci. Technol. 2009 32 27 30
Hai, X. U., Zhaopu, L. I. U., Lan, Y. & Linzhang, Y. Effect of pH on growth of several freshwater algae. Environ. Sci. Technol. 32, 27–30 (2009).
71. Tarabih OM Effects of the spatial distribution of best management practices for watershed wide nutrient load reduction Ecol. Eng. 2024 10.1016/j.ecoleng.2024.107211
Tarabih, O. M. et al. Effects of the spatial distribution of best management practices for watershed wide nutrient load reduction. Ecol. Eng.10.1016/j.ecoleng.2024.107211 (2024).
72. Nong X Shao D Shang Y Liang J Analysis of spatio-temporal variation in phytoplankton and its relationship with water quality parameters in the South-to-North Water Diversion Project of China Environ. Monit. Assess. 2021 10.1007/s10661-021-09391-6 34424412
Nong, X., Shao, D., Shang, Y. & Liang, J. Analysis of spatio-temporal variation in phytoplankton and its relationship with water quality parameters in the South-to-North Water Diversion Project of China. Environ. Monit. Assess.10.1007/s10661-021-09391-6 (2021).34424412
73. Echeveste P Agusti S Dachs J Cell size dependent toxicity thresholds of polycyclic aromatic hydrocarbons to natural and cultured phytoplankton populations Environ. Pollut. 2010 158 299 307 10.1016/j.envpol.2009.07.006 19631432
Echeveste, P., Agusti, S. & Dachs, J. Cell size dependent toxicity thresholds of polycyclic aromatic hydrocarbons to natural and cultured phytoplankton populations. Environ. Pollut. 158, 299–307. 10.1016/j.envpol.2009.07.006 (2010).19631432
74. Dang TD Modeling temporal and spatial variations of biogeochemical processes in a large subtropical lake: Assessing alternative solutions to algal blooms in Lake Okeechobee Florida J. Hydrol. Reg. Stud. 2023 10.1016/j.ejrh.2023.101441
Dang, T. D. et al. Modeling temporal and spatial variations of biogeochemical processes in a large subtropical lake: Assessing alternative solutions to algal blooms in Lake Okeechobee Florida. J. Hydrol. Reg. Stud.10.1016/j.ejrh.2023.101441 (2023).
75. Shen H Zhou P Advance in the studies on effect of environmental organic pollutants on the algae growth Acta Hydrobiol. Sin. 2002 26 529 535 10.3724/issn1000-3207-2002-5-529-v
Shen, H. & Zhou, P. Advance in the studies on effect of environmental organic pollutants on the algae growth. Acta Hydrobiol. Sin. 26, 529–535 (2002).
76. Yu H A method for quantifying the contribution of algal sources to CODMn in water bodies based on ecological chemometrics and its potential applications J. Environ. Chem. Eng. 2024 10.1016/j.jece.2024.111943
Yu, H. et al. A method for quantifying the contribution of algal sources to CODMn in water bodies based on ecological chemometrics and its potential applications. J. Environ. Chem. Eng.10.1016/j.jece.2024.111943 (2024).
77. Wang H Phytoplankton communities’ response to thermal stratification and changing environmental conditions in a deep-water reservoir: stochastic and deterministic processes Sustainability 2024 16 3058 10.3390/su16073058
Wang, H. et al. Phytoplankton communities’ response to thermal stratification and changing environmental conditions in a deep-water reservoir: stochastic and deterministic processes. Sustainability 16, 3058 (2024).
78. Li X Huang T Ma W Sun X Zhang H Effects of rainfall patterns on water quality in a stratified reservoir subject to eutrophication: Implications for management Sci. Total Environ. 2015 521 27 36 10.1016/j.scitotenv.2015.03.062 25828409
Li, X., Huang, T., Ma, W., Sun, X. & Zhang, H. Effects of rainfall patterns on water quality in a stratified reservoir subject to eutrophication: Implications for management. Sci. Total Environ. 521, 27–36. 10.1016/j.scitotenv.2015.03.062 (2015).25828409
79. Rigosi A Rueda FJ Hydraulic control of short-term successional changes in the phytoplankton assemblage in stratified reservoirs Ecol. Eng. 2012 44 216 226 10.1016/j.ecoleng.2012.04.012
Rigosi, A. & Rueda, F. J. Hydraulic control of short-term successional changes in the phytoplankton assemblage in stratified reservoirs. Ecol. Eng. 44, 216–226. 10.1016/j.ecoleng.2012.04.012 (2012).
