Seasonal and Long-Term Connections between Trophic Status, Sestonic Chlorophyll, Nutrients, Organic Matter, and Monsoon Rainfall in a Multipurpose Reservoir

: Due to rapid eutrophication, sustainable water quality management and supply are essential in drinking water sources and aquatic biota in large reservoirs. We evaluated the potentially crucial factors inﬂuencing the algal chlorophyll (CHL-a), nutrients, and the links between the rainfall and other vital elements in a large multipurpose reservoir (Yongdam Reservoir) during 2013–2019. We developed the empirical models on algal CHL-a, total phosphorus (TP), total nitrogen (TN), and TN:TP’s ambient ratios considering the maneuvering inﬂuence of Asian monsoons. The intensive rainfall during the monsoon months strongly impacted the nutrient regime and other vital factors. The seasonal patterns of algal CHL-a varied in response to the nutrient contents (TN, TP), suspended solids, and ambient N:P ratios along the longitudinal gradient. The conditional plot analysis, empirical modeling, and observations supported an overall P-limitation scenario, as was evident from the magnitude of N:P ratios (R 2 = 0.36, F = 24.9, p < 0.001). Furthermore, the reservoir’s trophic status alluded to the larger particles and blue-green algae during the monsoon and postmonsoon months. The correlation analysis, Mann–Kendall trend test, and principal component analysis illustrated compelling links between CHL-a, TP, and rainfall regime. The outcomes suggested the reservoir was primarily controlled by phosphorus limitation, with an increasing CHL-a tendency along with nitrogen dilution. However, a slight decline in phosphorus was also detected. The Yongdam Reservoir is under the threat of recurrent eutrophication events that could jeopardize this vital drinking water facility due to increasing agricultural and anthropic activities.


Introduction
The drinking water supply crises tremendously challenge sustainable human society expansion due to global anthropogenic pollution [1][2][3]. Reports abound on the human-caused disturbances at the catchment scale (e.g., by extensive industrialization, urbanization, intensive agriculture) in the lentic ecosystems, resulting in various water quality issues [4][5][6]. One of the most conspicuous and long-term issues is consistent nutrient enrichment (eutrophication), causing significant water quality degradations in the lentic ecosystems [7][8][9]. Such nutrient-related dynamics are primarily determined by the anthropogenic discharges combined with the unobserved internal processes [10,11]. Nutrient enrichment causes uncontrolled phytoplankton growth with subsequent harmful algal blooms (HABs) events. The resultant impacts cause problems such as oxygen depletion, pH fluctuations, causing community-level shifts in the aquatic organisms, toxins from certain phytoplankton species [12], drinking water quality, fishery catch, and recreation [13,14]. Therefore, monitoring water quality in large and multipurpose reservoirs is essential by in-The altered flow regime and water quality degradations in the reservoir impacted fish assemblages and the downstream river's trophic status [48]. Therefore, interlinked zonal water quality assessment of the reservoir is more critical for linked riverine ecosystems and safe drinking water use.
The Yongdam reservoir is one of the leading freshwater sources for the inhabitants of the west-coast region in South Korea. It plays a larger role as a multipurpose reservoir and is a unique water body in several ways, including it being the tributary reservoir, a hub of hydroelectric generation, and greatly influenced by the artificial upstream reservoir where rainfall strongly impacts the water flow. It is of grander importance to investigate the ecological and drinking water quality that is equally presenting a lentic and lotic type of ecosystems. Therefore, the scientific hypothesis of our study included whether the spatiotemporal water quality patterns in Yongdam Reservoir and their relationships are changed along the longitudinal zonation of the reservoir during study period. We also investigated the seasonal and annual water quality trends and inter-relationship between a large set of water quality parameters. Furthermore, we hypothesized the seasonal trophic state and trophic state index deviations (TSID) along with the linear regression models to determine the spatio-temporal relationship between nutrients and CHL-a level. Additionally, we applied multivariate analytical techniques such as principal component analysis (PCA) and the non-parametric Mann-Kendall trend (MKT) test to evaluate the long-term trend of the water quality parameters in the system.

Study Area
The Yongdam dam is located upstream (127 • 31 40" E, 35 • 56 30" N) in the Geum River, which flows through the Midwest Korean Peninsula (Figure 1). The reservoir was constructed as a part of the Yongdam dam for flood control (flood control capacity: 137 million m 3 /year), for water supply to Jeonju city (700,000 m 3 /d with a planned increase to 1,050,000 m 3 /d), and hydropower generation (24,400 kW). The principal morphometric features of the Yongdam reservoir show it as a narrow and deep reservoir with a surface area of about 37 km 2 , a maximum width of 1 km, and a maximum depth of about 70 m. The reservoir terrain slope ranges from 2% to 45%, while the average slope is 15%. The reservoir has a retention time of 318 days and is fed by a catchment area of approximately 930 km 2 with five tributaries. This constitutes 9.45% of the total basis area of the Geum River watershed. The reservoir catchment area's land-use patterns consist mainly of forest (68%), agricultural (30%), and urban regions of 1.2%.

Sampling Sites
The monitoring sites explored for water quality evaluation are located at reasonable distances from each other and distinctly cover the longitudinal zones (riverine, transition and lacustrine) of Yongdam Reservoir (Figure 1). Site 1 (S1) was exactly located where the water intake tower is erected; therefore, it represents the lacustrine zone (Lz). Sites 2 and 3 (S2, S3) represented the transition zone (Tz), while the riverine zone (Rz) was mainly defined by site 4 (S4). The Yongdam reservoir is a dendritic-shaped reservoir that receives freshwater from various feeding stream channels.

Field Sampling and Rainfall Data
The reservoir's monthly water quality data collected during 2013-2019, were procured from the Korean Ministry of Environment. The water samples of the given 17 variables were collected once per month from the four disparate sites related to the reservoir's fluvial and lentic zone. We measured the water temperature (WT), conductivity (Cond), pH, dissolved oxygen (DO), and chlorophyll-a (CHL-a) in situ with multi-parameter water quality sensors (YSI Sonde Model 6600, YSI Incorporated Yellow Springs, OH, USA). Samples for nutrient analyses were stored in plastic water sample bottles. Rainfall and inflow data during 2013-2019 were obtained from the Korean Meteorological Administration.

Sampling Sites
The monitoring sites explored for water quality evaluation are located at reasonable distances from each other and distinctly cover the longitudinal zones (riverine, transition and lacustrine) of Yongdam Reservoir ( Figure 1). Site 1 (S1) was exactly located where the water intake tower is erected; therefore, it represents the lacustrine zone (Lz). Sites 2 and 3 (S2, S3) represented the transition zone (Tz), while the riverine zone (Rz) was mainly defined by site 4 (S4). The Yongdam reservoir is a dendritic-shaped reservoir that receives freshwater from various feeding stream channels.

Field Sampling and Rainfall Data
The reservoir's monthly water quality data collected during 2013-2019, were procured from the Korean Ministry of Environment. The water samples of the given 17 variables were collected once per month from the four disparate sites related to the reservoir's fluvial and lentic zone. We measured the water temperature (WT), conductivity (Cond), pH, dissolved oxygen (DO), and chlorophyll-a (CHL-a) in situ with multi-parameter water quality sensors (YSI Sonde Model 6600, YSI Incorporated Yellow Springs, OH, USA). Samples for nutrient analyses were stored in plastic water sample bottles. Rainfall and inflow data during 2013-2019 were obtained from the Korean Meteorological Administration.

Water Quality Analyses
Nutrients were analyzed by the standard methods of the Ministry of Environment, Korea (2000). Total phosphorus (TP) was determined using unfiltered water digested by the ascorbic acid method after persulfate oxidation. Orthophosphate (PO4-P) was investigated without digestion using filtered water based on the ascorbic acid method. Total dissolved phosphorus (TDP) was measured by the ascorbic acid method through filtered water digestion. Total nitrogen (TN) was determined by UV spectrophotometric method followed by a potassium sulfate digestion. Ammonia nitrogen (NH 4 -N) was investigated at 630 nm according to the phenate method after filtering the water through GF/C filters, and nitrate-nitrogen (NO 3 -N) was determined by the ion chromatography method after filtering the water. Organic matters indicators (BOD, COD, TOC), TSS, and chlorophyll-a (CHL-a) were measured by standard methods of APHA [49]. Total suspended solids (TSS) were determined after drying at 105 • C for one hour. Chlorophyll-a concentration was measured by using a spectrophotometer after extraction in ethanol [49]. The detailed water quality assessment methodology for all the tested parameters can be consulted from Atique and An [4].

Trophic Status Analysis
Natural and anthropogenic nutrient contributing factors primarily influence the trophic state of the lentic ecosystem. The trophic status was evaluated based on Carlson's trophic state index (TSI), which integrates the leading trophic state indicators in the lakes and reservoirs. The TSI was individually calculated by using the natural log

Statistical Analysis
We used the Kolmogorov-Smirnov check to normality in the dataset. The water chemistry factors and nutrients (CHL-a, TN, TP, TN:TP) were log-transformed for a normal distribution. We determined the seasonal and spatial variations of the water quality variables using analysis of variance (ANOVA) with the Kruskal-Wallis test, a non-parameter method. We applied various multivariate analytical techniques such as principal component analysis (PCA) and Mann-Kendall test (MKT) to identify the multiple links between the water chemistry factors. PCA helped to identify the data variance and while MKT illustrated the prevalent trends during the study duration. We used the conditional plotting method to testify the limiting nutrient(s) with the help of the R program. Correlation analysis was also performed in R studio to identify the prevalent links between the inflow, rainfall, and selected water quality parameters. The trophic state index deviation (TSID) was performed to see the dominant trophic condition in the reservoir identified based on seasonal rainfall patterns. Sigma Plot v. 12.5 (Systat Software Inc., San Jose, CA, USA) was used to perform the ANOVA, linear regression models, while we used PAST (ver. 4.3) for the principal component analysis (PCA).

Role of Monsoon and Water Chemistry Variability
The seasonal rainfall patterns significantly influenced the critical water chemistry parameters, especially TP, TSS, TOC, and light regime ( Figure 2). The nutrients, organic matters, TSS, and CHL-a decreased along the longitudinal gradient from site 4 to site 1, showing significant variations, whereas TN:TP ratio increased significantly ( Table 1). The annual and maximum BOD and TOC values showed low organic matter pollution in the reservoir (Supplementary Figure S1). The seasonal variations in TSS, CHL-a, TOC, and TP noticeably increased from site 1 to site 4 ( Table 2). TN and its associated chemical forms increased from site 1 to site 4, while it was not significant (H = 4.5, p > 0.05). The averages and minimum TN values at different locations showed that the reservoir is a high-nitrogen content system like other temperate reservoirs [51]. Due to the lower TP value, the average of the ambient TN:TP ratios at each site were higher than different types of the reservoir, such as agricultural and estuarine reservoirs [4]. The minimum values of TN:TP was more than 20 at each site, indicating a P-limitation. The CHL-a mean values showed an increasing trend from site 1 (2.88 ± 1.10 µg/L) to site 4 (4.86 ± 1.73 µg/L), displaying mesotrophic to eutrophic nutrient enrichment conditions [50]. TP's maximum values showed that cyanobacterial risk reached approximately 80% at site 4 and 40% at locations 3 and 4.
Water 2021, 13, 1720 6 of 21 average of the ambient TN:TP ratios at each site were higher than different types of the reservoir, such as agricultural and estuarine reservoirs [4]. The minimum values of TN:TP was more than 20 at each site, indicating a P-limitation. The CHL-a mean values showed an increasing trend from site 1 (2.88 ± 1.10 µg/L) to site 4 (4.86 ± 1.73 µg/L), displaying mesotrophic to eutrophic nutrient enrichment conditions [50]. TP's maximum values showed that cyanobacterial risk reached approximately 80% at site 4 and 40% at locations 3 and 4. Such findings are also supported by several previous studies [14,52]. The summer monsoon strongly influenced the hydrological variables in the reservoir. A total of 41% of the total rainfall occurred during the monsoon months (July-August). Additionally, the precipitation during the monsoon months was greater than an average line of months (Supplementary Figure S2). As a result of intense rainfall, water quality parameters were altered significantly. Despite higher TOC levels recorded, the TSS, BOD, TN, and TP were higher during the intensive summer monsoon season than in the year. TOC is now increasingly used as an indirect indicator of organic pollution that is also linked with rainfall events, as exhibited through our findings. The predominant contributors of TOC include the internal lake processes as well as the surrounding environment. TSS variations showed a close similarity with the CHL-a. The initially lower TP concentration with subsequent increments during and post-rainfall periods indicated that P loading increased from the watershed under the impact of monsoon as well as higher inflow of solid contents [25]. TP's peak values (more than 30 µg/L) indicated that the algal growth with cyanobacteria dominant occurs in August and September, and the risk of cyanobacteria is more than 40% [52]. As presented in several previous studies, this study also highlights the dominating role of intensive rainfall events and overall rainfall intensity.
That is why the inclusion and impact of monsoon on nutrients, flow regime, organic matter, and algal CHL-a productivity remain crucial and always relevant.  The average TN:TP ratios significantly decreased from pre-monsoon (222 ± 93.2) to monsoon (214 ± 121.4) and post-monsoon (182 ± 72.3, H = 9.59, p < 0.05). The growth of CHL-a continuously increased from 2.16 ± 1.23 µg/L to 6.05 ± 4.38 µg/L between May and November. The increasing trend in CHL-a's higher values indicated that the risk of chlorophyll-a for drinking water repeatedly reached moderate levels (more than 10 µg/L) from August to November [53]. The findings showed the nutrient loadings from watersheds caused by extreme rainfall events that increased the algal CHL-a growth [21,[54][55][56]. The increasing TSS levels, BOD, TP, and CHL-a during the flood years declined, indicating intensive rainfall mediated dilution [20,46,57,58].

Long-Term Trends in the Water Quality Parameters
We used the non-parametric Mann-Kendall trend (MKT) test to check the long-term tendencies in the water quality factors in the Yongdam reservoir, and the results are showing in Table 3. The MKT analyses indicated an increasing trend in the organic matter (COD and TOC), while no trend was detected in the case of BOD. The pH, WT, and TSS alluded to an increasing trend during the study that could indicate that the changing weather conditions strongly impacted the reservoir. The aquatic biota is prone to changes because of such physical changes in the water chemistry. Table 3. Long-term trend analysis of limnological variables in the reservoir over the period 2013-2019 (S-Statistic, Z-Standardized variable, (+)-significantly upward trend, (0)-no trend, (−)-significantly downward trend, PP-particulate phosphate, TCB-total coliform bacteria). The EC, BOD, TN, and allied chemical species (TDN and NH 4 -N) displayed no trends, while NO 3 -N indicated to a negative trend. This declining trend in nitrate pointed towards the reduction in nitrate-containing fertilizers in the reservoir watershed. However, the increasing TP (TDP and PO 4 -P) allied chemical species alluded to a growing preference for phosphate fertilizers in the region. However, TP did not propose any trend based on MKT. The algal chlorophyll and TCB displayed an increasing overall trend in the reservoir that explained that the reservoir is subjected to various pressures, including impending water quality deterioration, especially human consumption. Our findings suggested that the key regulatory factors related to human water consumption are continuously increasing. The water quality managers need to advance to mitigate the looming water quality crises as the Yongdam reservoir is predominantly used as a human drinking water resource. It also alluded to the changing trends in fertilizer utilization among the agricultural community and shifting industrial production trends.

Spatio-Seasonal Links between Nutrients, Ambient Ratios, and Sestonic CHL-a
The spatial and seasonal empirical modeling illustrated algal CHL-a was regulated by TP in a moderately strong regression relationship rather than linked with the reservoir TN (Figure 3, Table 4). The TN:TP ambient ratios and TP demonstrated strong (R 2 < 0.70) negative links during the pre-monsoon, monsoon, and post-monsoon periods (Table 4). However, there were strong links between the ambient ratios and TN (R 2 = 0.76) during the monsoon. The algal CHL-a:TP showed the moderate-strong relationship to log-transformed CHL-a (n = 48, R 2 = 0.35, p < 0.001, F 1, 46 = 24.53). However, the monthly TN had weak positive correlation with CHL-a (n = 48, R 2 = 0.13, p < 0.05, F 1, 46 = 6.85). Water 2021, 13, x FOR PEER REVIEW 10 of 22    On the other hand, the regression trend in the log-transformed TN:TP ratios showed that the ratios declined linearly (n = 48, R 2 = 0.90, p < 0.001, F 1, 46 = 403.65) with an increase of TP, while there were no links with the TN (p > 0.05, F 1, 46 = 1.79). The TN:TP ratio dynamics explained 35% of the log-transformed CHL-a relationship, just like TP. This indicated the potential P-limitation scenario in the Yongdam reservoir. The reservoir shows a severe P-limitation scenario, which is also supported by the studies in the past [18,59,60].

Conditional Plotting
It is well-established that the lentic ecosystem's algal productivity is predominantly linked to the ambient nutrient (TP, TN) regime and is mainly contributed by agricultural and other anthropic actions [4,24]. We performed the conditional plotting using the R program, and the results showed that the TP concentrations better explained the algal growth (CHL-a) in the temperate reservoir (Yongdam) than that of TN ( Figure 4). This strongly suggested that the reservoir is firmly a P-limited ecosystem [20]. Moreover, conditional plotting has been used to identify limiting nutrients in aquatic systems [46,61]. When two predictors are strongly correlated (R 2 > 0.70), collinearity problems may emerge that could hinder the predominant nutrient's determination, causing the algal growth limitation. However, we performed the conditional plotting, and the findings have shown that TP and TN are poorly correlated (R 2 = 0.10) in the Yongdam Reservoir. In the present case, the conditional plots showed that the relationship between CHL-a and TP was relatively stable in Yongdam Reservoir, as indicated by the smooth lines on the four panels with similar slopes, which stated that the impact of TP on CHL-a was consistent regardless of the TN level, specifying a P-limited nutrient condition in the reservoir.
The conditional plot did not illustrate any overlapping interactions between TP and TN, which has further confirmed that Yongdam Reservoir is a P-limited system. The algal productivity is predominantly affected by the phosphorus regime in the reservoir. These findings are beneficial in designating the role of nutrients in algal productivity and are supported mainly by previously published research in the Asian region, predominantly impacted by the intensive summer monsoon rainfall. The conditional plot did not illustrate any overlapping interactions between TP and TN, which has further confirmed that Yongdam Reservoir is a P-limited system. The algal productivity is predominantly affected by the phosphorus regime in the reservoir. These findings are beneficial in designating the role of nutrients in algal productivity and are supported mainly by previously published research in the Asian region, predominantly impacted by the intensive summer monsoon rainfall.

Co-Relating Water Chemistry Parameters
We determined the strength of Pearson's correlation coefficient (r) as moderately strong (r ≥ 0.50-≤0.69) to strong (r ≥ 0.7) linear associations between the water quality parameters, which mainly indicated a moderate to strong correlation among the most critical parameters ( Figure 5). The correlation coefficient values between ≥0.3-≤0.49 were designated as weak linear relationship while less the 0.3 was considered as little to no linear relationship. For instance, BOD and COD showed a strong positive correlation (r ≥ 0.7) with TSS, TOC, TP, and allied chemical species, and TN, while a strong negative relationship to the ambient ratios of TN:TP (p < 0.001, r > 0.70, respectively). The COD showed strong positive correlations with BOD, TSS, TOC, TP, total dissolved phosphorus (TDP).

Co-Relating Water Chemistry Parameters
We determined the strength of Pearson's correlation coefficient (r) as moderately strong (r ≥ 0.50-≤0.69) to strong (r ≥ 0.7) linear associations between the water quality parameters, which mainly indicated a moderate to strong correlation among the most critical parameters ( Figure 5). The correlation coefficient values between ≥0.3-≤0.49 were designated as weak linear relationship while less the 0.3 was considered as little to no linear relationship. For instance, BOD and COD showed a strong positive correlation (r ≥ 0.7) with TSS, TOC, TP, and allied chemical species, and TN, while a strong negative relationship to the ambient ratios of TN:TP (p < 0.001, r > 0.70, respectively). The COD showed strong positive correlations with BOD, TSS, TOC, TP, total dissolved phosphorus (TDP). The other organic matter (TOC) measure manifested a robust positive relationship with BOD, COD, TSS, TP, TDP, and PO 4 -P. On the other hand, it illustrated a robust negative link with TN:TP. The strongly negative linear relationship between DO and WT suggested the dynamic shifts in reservoir DO that could be linked to the substantial temperature fluctuations or diurnal photosynthesis [62]. On the other hand, pH displayed weak to moderately strong positive linear relationship with most of the parameters with the COD (r = 0.55) showing moderately strong links. EC displayed the highest links with the TN (r = 0.60) and TDN (r = 0.62) while slightly lower links with TSS (r = 0.58). This alluded to the internal loading of ionic content from the alluvial sediments, potentially enriched with high nutrient contents [16,19].
The other organic matter (TOC) measure manifested a robust positive relationship with BOD, COD, TSS, TP, TDP, and PO4-P. On the other hand, it illustrated a robust negative link with TN:TP. The strongly negative linear relationship between DO and WT suggested the dynamic shifts in reservoir DO that could be linked to the substantial temperature fluctuations or diurnal photosynthesis [62]. On the other hand, pH displayed weak to moderately strong positive linear relationship with most of the parameters with the COD (r = 0.55) showing moderately strong links. EC displayed the highest links with the TN (r = 0.60) and TDN (r = 0.62) while slightly lower links with TSS (r = 0.58). This alluded to the internal loading of ionic content from the alluvial sediments, potentially enriched with high nutrient contents [16,19]. The organic pollution indicators (BOD and COD) are regularly used as signs of organic pollution on a global scale [63]. However, TOC has also been suggested as an alternative indicator of the organic matter prevalence in the freshwater ecosystems, especially those impacted by the WWTPs. The strong correlation between BOD, COD, TOC, and TSS indicated the presence of non-biodegradable organic pollutants in the reservoir that could The organic pollution indicators (BOD and COD) are regularly used as signs of organic pollution on a global scale [63]. However, TOC has also been suggested as an alternative indicator of the organic matter prevalence in the freshwater ecosystems, especially those impacted by the WWTPs. The strong correlation between BOD, COD, TOC, and TSS indicated the presence of non-biodegradable organic pollutants in the reservoir that could seriously damage the sustainable water quality of the Yongdam Reservoir. This also indicated that the runoff water mediated by the torrential monsoon currents and increasingly inflowing organic pollutants from the WWTPs could be a considerable problem in the future [18,23]. The strong connections between solids regime and nutrients indicated rich inflows of the adsorbed phosphorus and other dissolved nutrients (TN, TP, and allied chemical species) that could be linked with high monsoon rainfall events in the reservoir watershed [46,48].
The WT showed a strong influence on the TCB (r = 0.87), while algal CHL-a demonstrated a moderately strong (r = 0.65) relationship, thereby compelling links between water temperature, primary productivity, and bacterial communities. Furthermore, the sestonic CHL-a illustrated a moderately strong (r ≥ 0.50-≤0.69) association with BOD, COD, TSS, TOC, and TCB. Similarly, the precipitation regime and inflow displayed a relatively strong positive influence on the presence of TCB. However, the DO level demonstrated a strong negative (r = −0.82) with TCB, which indicated the bacterial communities are likely to populate in the places where the water quality is rapidly worsening due to chemical pollution. This will further endanger the aquatic life and the human health due to the continuous water consumption. On the other hand, the moderately strong relationship between TCB and water inflow (r = 0.59) and rainfall intensity (r = 0.57) indicated that the pollutants are contributed by variety of factors including human fecal matter, livestock, and other sources.
The strong positive (r = 0.93) correlation between the rainfall intensity and inflow alluded to the runoff water as the reservoir's predominant water resource. This also indicated that the reservoir was more vulnerable to natural events of intensive rainfall as it can potentially play as a connecting factor between variety of pollutants and nutrients present in the reservoir watershed [18]. TN and TDN indicated a weak or no correlation with most of the parameters. However, TP and allied chemical species displayed stronger (r ≥ 0.7) associations with each other and the TN:TP ratios, TOC, TSS, and weak links with the WT and DO. The ambient TN:TP ratios pointed to the more robust but negative (r ≥ 0.7) associations with TOC, BOD, COD, TP, and its other chemical forms, including particulate phosphorus. Moreover, the other rations, i.e., CHL-a:TP and CHL-a:TN, did not show any significant associations with the preponderance of the chemical water quality parameters. The primary productivity indicator (CHL-a), however, showed no explicitly strong associations with any of the critical factors; however, it showed moderately strong (r ≥ 0.50-≤0.69) with TP rather than with TN. Therefore, it is vital to mention that Pearson's correlation revealed the TP as the more convincing limiting nutrient in the Yongdam reservoir.
The moderately strong connections between algal CHL-a and TP regime strongly indicated that TP makes the greatest contributions to explain the algal CHL-a variations in the reservoir like most of the other multipurpose reservoirs in South Korea. The increasing human activities including resorting to the intensive agriculture practices, livestock farming, industrial actions are all connected to the reservoir through the rich inflows of the water currents after rainfall events during monsoon [47,55,58]. The other sources of water quality deterioration and strong connections between leading water quality parameters comprised of the WWTPs, industrial effluents, domestic releases. All these factors are strongly revealed by our current correlation assessment in the Yongdam Reservoir.

Predominant Links between Water Chemistry Parameters
The principal component analysis is widely practiced in environmental quality data to reduce large datasets into smaller components (or principal components) with the most negligible information loss. The statistics method helps interpret the relationship among the measured parameters in large datasets to compellingly display the assortments to enhance our understanding of the parameters. The eigenvalues helped to select the principal components (PC) displaying most of the water quality variations. The loading values of the water quality variables were considered strong (>0.75), moderate (0.75-0.50), and weak (0.50-0.30), which also illustrated the strength of variation in the vital parameters. The first four PCs displayed 88.3% variance, with the first PC explaining 48.31% and the second as 16.68% (Table 5, Figure 6). Most of the water chemistry factors, including the nutrient and organic matter regime, displayed strong to moderately strong positive (>0.75, 0.75-0.50) loadings in the PC 1, alluding to the most significant PC. For instance, the WT, BOD, COD, TSS, TOC, TP, TDP, PO 4 -P, and CHL-a, alluded to strongly positive, while the DO and TN:TP revealed strong negative loadings. We also noticed moderate positive loading values of EC, pH, particulate phosphate (PP), TCB, weak positive loadings of TN, TDN, CHL-a:TN, inflow rate, and precipitation. The PC 1 indicated that the phosphorus regime and WT strongly influenced primary productivity growth (CHL-a). Consequently, the DO level decline with increasing organic matter loads, PP, and bacterial populations. The second PC interpreted 16.7% of the variance and displayed a strongly favorable loading of NO 3 -N only. However, moderate positive loadings of EC, TN, TDN, precipitation, and TN's weak positive loadings with TP and water inflow. Besides, the TOC showed moderate negative loading. It maintained intense rainfall, reducing organic matter and continuously transporting the reservoir's nitrogen contents, mostly NO 3 -N loading. Moreover, the negative loadings of TDP, PO 4 -P, NH 4 -N in the PC 2 were probably related to the ionic dilution by the intensive monsoon rainfall as the reservoir inflow is strongly impacted the high rainfall events. The fourth PC elucidated 8.1% of the total variance, with the EC, CHL-a, CHL-a:TP, and CHL-a:TP showing moderate positive loadings, whereas rainfall, inflow, PO 4 -P illustrated the weak negative loadings, indicating the dilution impact of intensive monsoon rainfall on the primary productivity, and an overall positive effect on the phosphorus loadings in the reservoir watershed [64]. Water 2021, 13, x FOR PEER REVIEW 16 of 22 Figure 6. Biplot on the first two principal components (color meaning: blue-site 1, brown-site 2, green-site 3, purple-site 4).
The PC 1 indicated that the phosphorus regime and WT strongly influenced primary productivity growth (CHL-a). Consequently, the DO level decline with increasing organic matter loads, PP, and bacterial populations. The second PC interpreted 16.7% of the variance and displayed a strongly favorable loading of NO3-N only. However, moderate positive loadings of EC, TN, TDN, precipitation, and TN's weak positive loadings with TP and water inflow. Besides, the TOC showed moderate negative loading. It maintained intense rainfall, reducing organic matter and continuously transporting the reservoir's nitrogen contents, mostly NO3-N loading. Moreover, the negative loadings of TDP, PO4-P, NH4-N in the PC 2 were probably related to the ionic dilution by the intensive monsoon rainfall as the reservoir inflow is strongly impacted the high rainfall events. The fourth PC elucidated 8.1% of the total variance, with the EC, CHL-a, CHL-a:TP, and CHL-a:TP showing moderate positive loadings, whereas rainfall, inflow, PO4-P illustrated the weak negative loadings, indicating the dilution impact of intensive monsoon rainfall on the primary productivity, and an overall positive effect on the phosphorus loadings in the reservoir watershed [64].

Trophic State Index Deviation (TSID)
We performed the site-based and seasonal trend analysis on the trophic state index (TSI) of TP, CHL-a, and SD. The graphical presentation highlights the predominant

Trophic State Index Deviation (TSID)
We performed the site-based and seasonal trend analysis on the trophic state index (TSI) of TP, CHL-a, and SD. The graphical presentation highlights the predominant trophic condition in the Yongdam reservoir (Figure 7). The TP trophic state indicated increasing eutrophication tendency during the monsoon period, which alluded to higher TP inflows with the runoff inflowing. The TP showed an expanding trend from oligotrophic to hypereutrophic from sites 1 to 4. The algal productivity (CHL-a) displayed clear eutrophication signs at site 4, indicating an increase in algal productivity along the reservoir gradient. However, the TSI (SD) demonstrated a mean mesotrophic state at all the study sites.
trophic condition in the Yongdam reservoir (Figure 7). The TP trophic state indicated increasing eutrophication tendency during the monsoon period, which alluded to higher TP inflows with the runoff inflowing. The TP showed an expanding trend from oligotrophic to hypereutrophic from sites 1 to 4. The algal productivity (CHL-a) displayed clear eutrophication signs at site 4, indicating an increase in algal productivity along the reservoir gradient. However, the TSI (SD) demonstrated a mean mesotrophic state at all the study sites. We plotted the two-dimensional seasonal response of the TSI factors (TP, CHL-a, and TP) given by Carlson [50] as TSI (CHL-a)-TSI (SD) along the x-axis while TSI (CHL-a)-TSI (TP) along the y-axis. The post-monsoon trophic status mainly deviated from the algal productivity (blue-green algae) and larger particles ( Figure 8). There was an overall tendency of larger particles presence and the blue-green algae, indicating a more substantial summer monsoon influence on the reservoir productivity. The seasonal deviations suggested that the reservoir's mainstream displayed a minimal to no zooplankton grazing activity during the year. However, the Yongdam reservoir showed a propensity of nonalgal light limitation during the pre-monsoon season, with a slight P-limitation. Slighter chances of zooplankton grazing occurred during the monsoon season as well. The burgeoning blue-green algal tendencies during the pre-monsoon, monsoon, and post-monsoon seasons mostly alluded to the higher chances of eutrophication due to the increasing We plotted the two-dimensional seasonal response of the TSI factors (TP, CHL-a, and TP) given by Carlson [50] as TSI (CHL-a)-TSI (SD) along the x-axis while TSI (CHL-a)-TSI (TP) along the y-axis. The post-monsoon trophic status mainly deviated from the algal productivity (blue-green algae) and larger particles ( Figure 8). There was an overall tendency of larger particles presence and the blue-green algae, indicating a more substantial summer monsoon influence on the reservoir productivity. The seasonal deviations suggested that the reservoir's mainstream displayed a minimal to no zooplankton grazing activity during the year. However, the Yongdam reservoir showed a propensity of non-algal light limitation during the pre-monsoon season, with a slight P-limitation. Slighter chances of zooplankton grazing occurred during the monsoon season as well. The burgeoning blue-green algal tendencies during the pre-monsoon, monsoon, and post-monsoon seasons mostly alluded to the higher chances of eutrophication due to the increasing nutrient enrichment. Moreover, such seasonal trends also revealed the potential impact of WRT and the fluctuating water volume.  The TSI is a multidimensional approach mainly applied to the lentic ecosystems, their physiognomic, and the essential energy transfer dynamics. The leading conclusions specify that no particular nutrient could be held responsible for the reservoirs' trophic state [4,50]. The water clarity and algal biomass are inversely linked based on the grazing principle, while phosphorus emerges as the crucial nutrient contributing factor in this condition [65]. Consistent nutrient enrichment from diffused sources is causing multiple and widespread problems, including the risks of recurrent eutrophication events in various parts of the world across the lentic, lotic and coastal ecosystems. Severe eutrophication actively drives the production of HABs that are strongly capable of producing harmful and toxic compounds that jeopardize the drinking and ecological water quality. As agriculture is the leading contributor of nutrients, nutrient reduction measures should help mitigate the recurrent nutrient enrichment if the reduction measures are adequately implemented. This study adequately supports this notion as sustained in the case of the seasonal trophic state. However, to specify the most critical and critical limiting nutrient, a bioassay is conducted to reveal the suitable answers to such queries and understand the underlying basis of a specific nutrient.

Conclusions and Management Strategies
The effective management and sustainable conservation of drinking water facilities are of imperative importance during the modern age. South Korea is rapidly turning towards intensive agricultural activities after urbanization and industrial success. As a result, the temperate reservoirs used for human drinking are under stress of nutrient enrichment and concomitant eutrophication. We performed this study in one of the largest The TSI is a multidimensional approach mainly applied to the lentic ecosystems, their physiognomic, and the essential energy transfer dynamics. The leading conclusions specify that no particular nutrient could be held responsible for the reservoirs' trophic state [4,50]. The water clarity and algal biomass are inversely linked based on the grazing principle, while phosphorus emerges as the crucial nutrient contributing factor in this condition [65]. Consistent nutrient enrichment from diffused sources is causing multiple and widespread problems, including the risks of recurrent eutrophication events in various parts of the world across the lentic, lotic and coastal ecosystems. Severe eutrophication actively drives the production of HABs that are strongly capable of producing harmful and toxic compounds that jeopardize the drinking and ecological water quality. As agriculture is the leading contributor of nutrients, nutrient reduction measures should help mitigate the recurrent nutrient enrichment if the reduction measures are adequately implemented. This study adequately supports this notion as sustained in the case of the seasonal trophic state. However, to specify the most critical and critical limiting nutrient, a bioassay is conducted to reveal the suitable answers to such queries and understand the underlying basis of a specific nutrient.

Conclusions and Management Strategies
The effective management and sustainable conservation of drinking water facilities are of imperative importance during the modern age. South Korea is rapidly turning towards intensive agricultural activities after urbanization and industrial success. As a result, the temperate reservoirs used for human drinking are under stress of nutrient enrichment and concomitant eutrophication. We performed this study in one of the largest drinking water temperate reservoirs in South Korea that spanned seven years (2013-2019). The outcomes have denoted that the monsoon rainfall regime strongly induced the water quality changes in the Yongdam reservoir. The intensive rainfall events primarily regulated the inflow and water residence time. We found that the potential factors influencing the algal chlorophyll (CHL-a) could be designated as the intensive rainfall patterns, influx, and phosphorus yielding nutrients. Various types of analyses corroborated each other, confirming the TP as the limiting factor. The PCA displayed that the reservoir water quality is primarily modified by the monsoon and organic matter closely linked to each other and the inflowing solids. The empirical models also supported the P-limitation scenarios along with the trophic state index deviation (TSID). The longitudinal gradient of the riverinetransition-lacustrine zones similarly displayed the variation patterns to the longitudinal zonation. Thus, the spatial and seasonal fluctuations between the algal CHL-a, nutrients, and TSS significantly accounted for the ambient water quality and nutrient enrichment status deviations.
Reservoir water quality's most common management issues include algal blooms, turbidity, nutrient enrichment, organic matter inflow, persistent large aquatic plants, water level fluctuations, and species shifts. Therefore, the most critical aspect of reservoir water quality management is the nutrient regulation and control of anthropic activities in the reservoir watershed, contributing to high nutrients and organic matter. The gradual reduction in TP loads and external nutrients control is the most commonly used measure in this case. Another successful option could be the diversion of nutrients discharge as well as strict monitoring of the WWTPs effluents along with other industrial releases. On the other hand, in-lake treatment methods (aeration and dredging) could also help alleviate the persistent water quality issues. However, such methods are mainly eutrophication targeting and offer only short-term treatment options. As lakes and reservoirs act as the sink for pollutants and inflowing sediments, active management of anthropic activities in the surrounding watersheds could offer long-term and comprehensive management options in the large lakes and reservoirs.

Data Availability Statement:
The data may be available upon request to the corresponding author, subject to the approval.