Trophic Responses of the Asian Reservoir to Long-Term Seasonal and Interannual Dynamic Monsoon

: The main objectives of the study were to determine the trophic response of the temperate reservoir to seasonal and interannual variabilities of monsoon inorganic solids and nutrients along the gradients of the morphologically complex Asian reservoir using long-term datasets between 2000–2018. Nutrient regime (total nitrogen—TN, total phosphorus—TP), total suspended solids (TSS), and chlorophyll-a (CHL-a) were primarily a ﬀ ected by an intensity of summer monsoon and the longitudinal structure of riverine (Rz), transitional (Tz), and lacustrine (Lz) zone. The reservoir is a nitrogen-rich system and the phosphorus content of the water was relatively low, and it had low mean N:P ratios ( < 40), implying a P-limiting system. The Lz was a highly P-limited zone in comparison to Rz and Tz zone during both drought (2015) and ﬂood year (2011). The TP content was higher in the mainstem (S3) than the embankment (S4 and S6) of the reservoir due to the monsoon river inputs of the nutrients. Nonparametric Mann–Kendall tests indicated that TP decreased over the long-term years in the Rz, while it did not show any trend in Tz, Lz, IT1, and IT2. TN showed an increasing trend in Rz, Tz, Lz, and IT2 except for IT1. The empirical regression model for chlorophyll nutrients showed that CHL-a had a strong positive relationship with TP (R 2 = 0.67, p < 0.01) than TN (R 2 = 0.06, p < 0.01), supporting the view that algal growth in lentic systems responds to TP enrichment and TP may provide a reliable basis for predicting algal biomass. The seasonality of CHL-a and TP showed a monomodal pattern and indicates that summer TP inﬂuences summer algal growth in Tz, Lz, and IT2. The water clarity (SD) of the reservoir was signiﬁcantly ( p < 0.01) inﬂuenced by TP (R 2 = 0.62), TSS (R 2 = 0.67), and CHL-a (R 2 = 0.68) rather than TN (R 2 = 0.10). The non-algal light attenuation coe ﬃ cient (K na ) was determined mainly by suspended solids and the monsoon hydrology. The trophic state was much higher when assessments were based on TSI (CHL-a) than on TSI (TP) and TSI (SD). TSI (CHL-a) indicated the eutrophic state of the reservoirs except for the zone of Lz during the premonsoon season. Analysis of trophic state index deviation (TSID) suggested that the blue-green algae dominated the algal community, and the e ﬀ ects of non-algal turbidity and zooplankton grazing were minor in the reservoir. p < 0.01), WT (r 0.55, p < 0.01), and CHL-a (r 0.51, p < 0.01) and was negatively correlated with SD (r = −0.57, p < 0.01). The TP was positively correlated with WT (r = 0.71, p < 0.01) and CHL-a (r = 0.71, p < 0.01) and was negatively correlated with SD (r = −0.67, p < 0.01). The inflow (r = 0.88, p < 0.01) and outflow (r = 0.84, p < 0.01) of the reservoir are highly regulated by the rainfall.


Introduction
Reservoirs are facing substantial water quality problems that lead to eutrophication and pose significant environmental, social, and economic impacts all over the world [1][2][3]. Rapid industrialization, urbanization, and intensive agricultural farming has accelerated eutrophication and has caused troublesome algal blooms and associated problems in the reservoirs [4,5]. Reservoir ecosystems typically have prominent longitudinal heterogeneities in water quality from the riverine zone (headwater) to Among them, Carlson provided a two-dimensional graphical approach for assessing the trophic state of a reservoir based on total phosphorus, chlorophyll, and Secchi depth [30]. The trophic state of the reservoir gradually decreased and showed dynamic behavior from the riverine zone to the lacustrine zone due to longitudinal patterns and seasonal hydrology [10]. This two-dimensional method is an effective tool for monitoring the long-term changes of pelagic structures and functions of the reservoir and designates valuable information for its management [31].
This study aimed to determine the effects of monsoons on the trophic response of the reservoir using the long-term limnological data set of 2000-2018. We also determined how algal chlorophyll is regulated by nutrient input and monsoon hydrology. The present research provides a means to compare variabilities of functional zonations between strong and weak monsoon years along the longitudinal gradients from the riverine to the lacustrine zone. In addition, we compared the nutrients, chlorophyll-a, and water clarity between the intake tower, located in the embayment, for drinking water supply and the mainstem site. To determine the trophic response of the reservoir, we analyzed seasonal and interannual variabilities of Secchi transparency, nutrient availability, algal CHL-a, and suspended solids in relation to non-algal light attenuation coefficients, and applied the approach of two-dimensional trophic state index deviation (TSID).

Study Site Description
Daecheong Reservoir is a warm, monomictic man-made lake, located in the central part of the Korean peninsula and formed at the end of December in 1980 by the impounding of the Geum River. This multipurpose reservoir has great importance because it supplies 922,000 m 3 of drinking water per day to 2 million people dwelling in two megacities (Daejeon and Cheongju) and is also used for flood control and hydroelectric power generation [32]. The reservoir is morphologically complex and has a surface area of 6.8 × 10 7 m 2 and a volume of 14.3 × 10 8 m 3 , with a mean depth of 21.2 m and a maximum depth of 69 m as well as a mean water residence time of 180 days [10]. Regular sampling of water quality parameters from the surface water body had been collected from S1 and S2 (riverine zone; 36.3743, 127.639, and 36.362, 127.592), S3 (transitional zone; and 36.435, 127.553), S5 (lacustrine zone: dam site; 36.475, 127.484), S4, and S6 (intake tower for drinking water supply; 36.373, 127.553 and 36.511, 127.503) as shown in Figure 1. The longitudinal patterns of the reservoirs were divided into three zones: riverine zone (Rz); transitional zone (Tz); and lacustrine zone (Lz). The Rz reflects the headwater zone (upper reach that has a narrow and channelized basin with relatively higher flow rates), the Lz means the down lake near the dam region (lower part of the reservoir with lower flow rates) and the Tz reflects the intermediate characteristics between Rz and Tz (a relatively broader basin than Rz with intermediate flow rates). The major tributaries of the reservoir are the Seowha Stream, Bouchoung Stream, Soung Stream, Youngdong Stream, Whebuk Stream, and Janggae Stream [7]. Based on the land use pattern, S1 consists of 40% rural area + 50% crop farming + 10% livestock farming. In comparison, S2 consists of 90% urban household + 10% crop farming and directly influenced by the wastewater disposal plants located at the upper reaches of the Seowha stream [7]. In contrast, S3 has been composed of 45% crop farming + 55% forest area, and S4, S5, and S6 consist of 78% forest area and 22% rural area, respectively [7].

Analysis of Water Quality Parameters, Precipitation, and Hydrological Data
The monthly surface water quality parameters data (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017)(2018) were obtained from the Korean Ministry of Environment. Water samples were collected once per month at each site from the epilimnetic zone of the surface to 50 cm using a sampling bottle and stored in ice to prevent exposure from sunlight and were immediately transferred in the laboratory. Electrical conductivity (EC), dissolved oxygen (DO), water temperature (WT), and hydrogen ion concentration (pH) had been determined by a portable multi-parameter analyzed (YSI Sonde Model 6600) during the sampling time. Secchi depth (20 cm metal disk; SD) was also measured at the time of sample collection. Total phosphorus (TP) was determined by the ascorbic acid method after persulfate oxidation, which was standardized by the Ministry of Environment, Korea [33]. Total nitrogen (TN) was estimated using the UV spectrophotometric method after the potassium persulfate digestion [34,35]. Biological oxygen demand (BOD) and chemical oxygen demand (COD) were determined by the Eatson and Franson method, which was also standardized by the Ministry of Environment, Korea [35].
Chlorophyll-a (CHL-a) was measured using a spectrophotometer after extraction in acetone, which was also standardized by the Korean Ministry of Environment [33]. To ensure the accuracy of the data, TP, TN, and CHL-a analyses were carried out in duplicate. The monthly precipitation data, as well as inflow and outflow, were collected from the Korean Meteorological Administration and Korean Water Resource Corporation, respectively.

Analysis of Water Quality Parameters, Precipitation, and Hydrological Data
The monthly surface water quality parameters data (2000-2018) were obtained from the Korean Ministry of Environment. Water samples were collected once per month at each site from the epilimnetic zone of the surface to 50 cm using a sampling bottle and stored in ice to prevent exposure from sunlight and were immediately transferred in the laboratory. Electrical conductivity (EC), dissolved oxygen (DO), water temperature (WT), and hydrogen ion concentration (pH) had been determined by a portable multi-parameter analyzed (YSI Sonde Model 6600) during the sampling time. Secchi depth (20 cm metal disk; SD) was also measured at the time of sample collection. Total phosphorus (TP) was determined by the ascorbic acid method after persulfate oxidation, which was standardized by the Ministry of Environment, Korea [33]. Total nitrogen (TN) was estimated using the UV spectrophotometric method after the potassium persulfate digestion [34,35]. Biological oxygen demand (BOD) and chemical oxygen demand (COD) were determined by the Eatson and Franson method, which was also standardized by the Ministry of Environment, Korea [35]. Chlorophyll-a (CHL-a) was measured using a spectrophotometer after extraction in acetone, which was also standardized by the Korean Ministry of Environment [33]. To ensure the accuracy of the data, TP, TN, and CHL-a analyses were carried out in duplicate. The monthly precipitation data, as well as inflow and outflow, were collected from the Korean Meteorological Administration and Korean Water Resource Corporation, respectively.

Trophic State Index Deviation and Non-Algal Light Attenuation Coefficient
The conventional criteria of trophic states based on total phosphorus, chlorophyll-a and Secchi depth had been used for the evaluation of eutrophication in the Daecheong Reservoir [36]. The range of average TSI for Oligotrophic is 30-40, Mesotrophic is 40-50, Eutrophic is 50-70, and Hypereutrophic is >70, respectively [36]. The trophic state index has been calculated using the following formulas [37] Deviations of the trophic state index, based on two-dimensional approaches, were inferred using the relationship: TSI (CHL-a) − TSI (SD) and TSI (CHL-a) − TSI (TP). This approach has also frequently been used to measure the degree of eutrophication and provides the limiting nutrient status in reservoirs [31]. Non-algal light attenuation coefficient (Kna) has been estimated by the following equation [38]:

Statistical Analysis
The Rosner outlier test and Mann-Kendall trend analysis of the water quality parameters had been done by ProUCL version 5.1. software [39]. Before empirical regression analysis, the concentration of TP, TN, CHL-a, TSS, and SD was converted into log 10 to boost normality. Sigma Plot 10.0 was used to perform the regression analysis of water quality parameters. Pearson correlation analysis of water quality and hydrological parameters were conducted by PAST software to know how the parameters are correlated in the reservoir [40]. This research considers the seasonal variability of water quality parameters. The terms of premonsoon (January-June), monsoon (July-August), postmonsoon (September-December), spring (March-May), summer (June-August), and fall (September-November) were used to define the temporal variations of water quality parameters and the reservoir trophic conditions.

Hydrology Pattern
Hydrology of the Korean reservoir is significantly influenced by the summer monsoon and showed distinct differences from 2000-2018 ( Figure 2). Annual precipitation in the reservoir during 2000-2018 was an average of 1115.10 mm, and half of the annual rainfall happened during monsoon season (July-August). In consequence, the rainfall was significantly greater in the monsoon season

Physicochemical Properties of Water Quality Parameters
The water quality parameters of the Daecheong Reservoir showed zonal disparities and varied from Rz to Lz and IT1 to IT2 (

Physicochemical Properties of Water Quality Parameters
The water quality parameters of the Daecheong Reservoir showed zonal disparities and varied from Rz to Lz and IT1 to IT2 ( Table 1). The TP concentration was lower compared to TN in the reservoir. The highest mean TP (28.15 ± 22.62), TN (1.91 ± 0.54), CHL-a (10.6 ± 8.38), TSS (3.51 ± 3.37), BOD (1.27 ± 0.39), and EC (155.1 ± 39.17) concentrations were observed in Rz rather than the other zones (Tz, Lz, IT1, and IT2) due to external loads from the tributary streams and watersheds. The higher mean TN:TP ratio (144.3) was found in the dam site (Lz), which means it is highly P-limited compared to the other sites. The average Secchi depth was lower in the Rz (2.35 ± 1.06) zone than the other sites. The highest water clarity was observed in the Lz (3.75 ± 1.55). The higher mean COD level was found in IT1 (3.75 ± 0.76) rather than Rz (3.17 ± 0.53), Tz (3.02 ± 0.53), Lz (2.85 ± 0.49), and IT2 (3.09 ± 0.59). The maximum WT was detected in IT1 (31 • C) and was paralleled to other zones. The highest mean pH (8.02 ± 0.54) and DO (10.11 ± 2.24) levels were determined in IT2.
The annual variation of water quality parameters also showed discrepancies in the reservoir as shown in Supplementary Materials ( Figure S1). The annual concentration of nutrients (TP, TN), organic matters (BOD, COD), ions (EC), suspended solids (TSS), and algal chlorophyll-a (CHL-a) were primarily affected by the annual precipitation rate in the reservoir catchment area. Based on the yearly precipitation rate, the high flood and drought year has been observed. The highest annual precipitation occurred in 2011 (1603 mm), while the lowest rainfall happened in 2015 (690 mm). During the high flood year (2011), the mean and maximum concentrations of TP was higher, whereas it was lowest during the drought year (2015). The highest TN concentration was observed in 2016. It was followed in the same pattern for BOD and COD during the study period. The reservoir showed the highest P-limitation during the drought year (2015), meaning that the TN:TP ratio was maximum compared to other years. In 2006, the mean and maximum EC concentration was higher in comparison to other years. The average TSS concentration was higher during the flood year (2011), but the maximum TSS was observed in 2005. During the high flood year (2011), the average algal growth was at its peak. Table 1. Zonal variations of water quality parameters in the Daecheong Reservoir during 2000-2018 (TP: total phosphorus, TN: total nitrogen, CHL-a: chlorophyll-a, SD: Secchi depth, TSS: total suspended solids, BOD: biological oxygen demand, COD: chemical oxygen demand, EC: electrical conductivity, WT: water temperature, pH: hydrogen ion concertation, DO: dissolved oxygen, Rz: riverine zone, Tz: transitional zone, Lz: lacustrine zone, IT1: drinking water intake tower 1, IT2: drinking water intake tower 2, SD: standard deviation, Min: minimum, Max: maximum, and CV: coefficient of variation; monthly data from 2000-2018 were used in Rz-IT2 to make the summary statistics of water quality parameters).

Variations of Physicochemical Parameters during Drought and Flood Year and in Embankment and Mainstem
The average annual precipitation at the reservoir during 2000-2018 was 1115.10 mm. The significant differences in rainfall between 2011 and 2015 were observed during the monsoon season. Total precipitation during July-August was 262 mm in 2015 whereas it was 759 mm in 2011. The inflow and outflow of the reservoir are highly influenced by the rainfall pattern. The interannual differences of the inflow directly influenced the outflow, reservoir stage, and water residence time. The hydrological differences between the flood and drought year of the reservoir have a great impact on the water quality parameters (Figures 3 and 4). The concentration of total phosphorus was higher during the flood year in contrast to the drought year. The   I  II III IV V VI VII VIII XI X XI XII  I  II III IV V VI VII VIII XI X XI XII   I  II III IV V VI VII VIII XI X XI XII  I  II III IV V VI VII VIII XI X XI XII   I  II III IV V VI VII VIII XI X XI XII  I  II III IV V VI VII VIII XI X XI XII   I  II III IV V VI VII VIII XI X XI XII  I II III IV V VI VII VIII XI X XI XII

Pearson Correlation Analysis
The Pearson correlation analysis revealed a strong positive and negative relationship between the water quality parameters of the Daecheong Reservoir ( Table 2). The moderate positive correlation had been found among pH and water temperatures (r = 0.36, p < 0.01) while it was strongly correlated with rainfall (r = 0.74, p < 0.01). The dissolved oxygen levels decreased with an increasing BOD (r = −0.46, p < 0.01), COD (r = −0.53, p < 0.01), TSS (r = −0.54, p < 0.01), TP (r = −0.70, p < 0.01), WT (r = −0.82, p < 0.01), and CHL-a (r = −0.64, p < 0.01) in the reservoirs. The biological oxygen demand (BOD) of the reservoir had been increasing with TSS (r = 0.58, p < 0.01), TN (r = 0.59, p < 0.01), TP (r = 0.74, p < 0.01), WT (r = 0.54, p < 0.01), and CHL-a (r = 0.64, p < 0.01). The COD was positively correlated with TSS (r = 0.56, p < 0.01), WT (r = 0.69, p < 0.01), TP (r = 0.59, p < 0.01), and CHL-a (r = 0.64, p < 0.01). The TSS was positively correlated with TP (r = 0.85, p < 0.01), WT (r = 0.55, p < 0.01), and CHL-a (r = 0.51, p < 0.01) and was negatively correlated with SD (r = −0.57, p < 0.01). The TP was positively correlated with WT (r = 0.71, p < 0.01) and CHL-a (r = 0.71, p < 0.01) and was negatively correlated with SD (r = −0.67, p < 0.01). The inflow (r = 0.88, p < 0.01) and outflow (r = 0.84, p < 0.01) of the reservoir are highly regulated by the rainfall.  I II III IV V VI VII VIII IX X XI XII   I II III IV V VI VII VIII IX X XI XII  I II III IV V Figure S2). The mainstem contained higher concentrations of TP than the embankment. During the monsoon, the TP concentration was higher in comparison to other months. In the mainstem, the concentration of TN was lower during the winter period. The phosphorus limitation was higher in the mainstem during August in comparison to other months. In the embankment, the chlorophyll-a growth indicated an increasing pattern in August, September, and October. On the other hand, in the mainstem, the chlorophyll growth was at peak during August for S3 and September for S5. The water clarity was lower during August and September in the mainstem and embankment due to higher chlorophyll growth and nutrients.

Trend Analysis of Water Quality Parameters
The trend analysis of water quality parameters was separately analyzed in the Rz, Tz, Lz, IT1, and IT2 using the nonparametric Mann-Kendall test (Table 3). Statistic S value for TP was −1717 (p < 0.04) in Rz, indicating that there was a decreasing trend, while TP did not show any trend of increasing and decreasing in Tz, Lz, IT1, and IT2. Statistic S was positive for TN in Rz (3875), Tz (6945), Lz (8862), and IT2 (7591) except for IT1 (−7491), representing an increasing trend. It was followed by the same pattern for TN:TP ratios and BOD and increased in Rz, Tz, Lz, and IT2. The COD revealed an increasing trend in Rz, Lz, IT1, and IT2 except for Tz. In Tz, the COD was showing no trend.  Table 3. Long-term trend analysis of water quality parameters in the Daecheong Reservoir during 2000-2018 (TP: total phosphorus, TN: total nitrogen, BOD: biological oxygen demand, COD: chemical oxygen demand, TSS: total suspended solids, EC: electrical conductivity, CHL-a: chlorophyll-a, SD: Secchi depth, Rz: riverine zone, Tz: transitional zone, Lz: lacustrine zone, IT1: drinking water intake tower 1, IT2: drinking water intake tower 2, D T : decreasing trend, I T : increasing trend, N T : no trend, monthly data from 2000-2018 were used to see the trend of water quality parameters separately from Rz to IT2).

Empirical Relationship of Chlorophyll-Nutrients
Empirical models based on log-transformed trophic parameters had been shown in Table 4. The reservoir's chlorophyll growth was more influenced by the nutrients in ambient water, and the growth of algal chlorophyll was better explained by the concentration of TP than by TN. In the empirical-based regression analysis model, when two predictors (TP and TN) are strongly correlated, collinearity problems arise, the meaning of slope is ambiguous, and coefficients tend to have a higher estimation of uncertainty. Our results showed that TP and TN in the Daecheong Reservoir are not strongly correlated. The regression analysis of TP and TN indicates that the R 2 is 0.21. The chlorophyll growth was highly determined by the concentrations of TP (R 2 = 0.69 RMSE = 0.14, p < 0.01) compared to TN (R 2 = 0.06, RMSE = 0.25, p < 0.01), indicating that the Daecheong Reservoir is a P-limited system. Moreover, conditional plot analysis was used to determine the limiting nutrient. The CHL-a-TP relationship is relatively stable in the Daecheong Reservoir (Supplementary Materials; Figure S3). The smooth line in each of the four panels shows a similar slope, which suggests that the effect of TP on CHL-a is consistent regardless of the level of TN, an indication of a phosphorus-limited reservoir. The CHL-a-TN relationship is less pronounced in the Daecheong Reservoir (Supplementary Materials; Figure S4). Additionally, the conditional plot shows the case of no interaction between TP and TN, which means the Daecheong Reservoir is phosphorus-limited. In Rz, spring TN elucidates 15% summer algal growth while summer TN explained a 13% fall in algal growth in Lz. In the Daecehong reservoir, the variation of summer algal growth is significantly controlled by summer TP in Tz (R 2 = 0.24), Lz (R 2 = 0.40), and IT2 (R 2 = 0.18).

Relations of Water Clarity to Other Water Quality Parameters
Regression analysis of log-transformed SD-TP, SD-TN, SD-TSS, and SD-CHL-a has been demonstrated in Figure 5. The water clarity of the reservoir was highly regulated by the ambient nutrients, suspended solids, and algal chlorophyll. The variations of water clarity were explained 62% by the concentration of total phosphorus (TP) while only 10% by the concentration of total nitrogen (TN). The concentrations of TSS (R 2 = 0.67, n = 72, p < 0.01) and CHL-a (R 2 = 0.68, n = 72, p < 0.01) directly influenced the water clarity of the Daecheong Reservoir.
Regression analysis of log-transformed SD-TP, SD-TN, SD-TSS, and SD-CHL-a has been demonstrated in Figure 5. The water clarity of the reservoir was highly regulated by the ambient nutrients, suspended solids, and algal chlorophyll. The variations of water clarity were explained 62% by the concentration of total phosphorus (TP) while only 10% by the concentration of total nitrogen (TN). The concentrations of TSS (R 2 = 0.67, n = 72, p < 0.01) and CHL-a (R 2 = 0.68, n = 72, p < 0.01) directly influenced the water clarity of the Daecheong Reservoir.

Non-Algal Light Attenuation
Regression analysis of Kna-TP, Kna-TN, Kna-TN:TP, Kna-TSS, and Kna-CHL-a has been presented in Figure 6. The ambient nutrient concentration (TP, TN) and suspended solids act as a putative light-limiting factor and affected the growth of algae in the waterbody. The variations of non-algal light attenuation has been explained by 38% by the concentration of TP while it was 28% by the concentration of TN. The concentrations of TSS were largely controlled by the monsoon rain and during the monsoon, the suspension from the sediment bed increased in the reservoir and impeded the light penetration. The non-algal light attenuation was largely explained by the TSS (R 2 = 0.52, n = 72, p < 0.01) in comparison to TN:TP (R 2 = 0.10, n = 72, p < 0.01) and CHL-a (R 2 = 0.04, n = 72, p < 0.01).

Non-Algal Light Attenuation
Regression analysis of Kna-TP, Kna-TN, Kna-TN:TP, Kna-TSS, and Kna-CHL-a has been presented in Figure 6. The ambient nutrient concentration (TP, TN) and suspended solids act as a putative light-limiting factor and affected the growth of algae in the waterbody. The variations of non-algal light attenuation has been explained by 38% by the concentration of TP while it was 28% by the concentration of TN. The concentrations of TSS were largely controlled by the monsoon rain and during the monsoon, the suspension from the sediment bed increased in the reservoir and impeded the light penetration. The non-algal light attenuation was largely explained by the TSS (R 2 = 0.52, n = 72, p < 0.01) in comparison to TN:TP (R 2 = 0.10, n = 72, p < 0.01) and CHL-a (R 2 = 0.04, n = 72, p < 0.01).

Trophic State Index and Its Deviation
The trophic state index of the Daecheong Reservoir based on TP, CHL-a, and SD had been estimated using a Carlson trophic model and it showed spatial and seasonal variations (Figure 7). During the premonsoon season, the lowest average TSI (TP) value was observed in Lz (36.94), indicating an oligotrophic condition, while Rz, Tz, IT1, and IT2 was in a mesotrophic condition.

Trophic State Index and Its Deviation
The trophic state index of the Daecheong Reservoir based on TP, CHL-a, and SD had been estimated using a Carlson trophic model and it showed spatial and seasonal variations (Figure 7). During the premonsoon season, the lowest average TSI (TP) value was observed in Lz (36.94), indicating an oligotrophic condition, while Rz, Tz, IT1, and IT2 was in a mesotrophic condition. During the monsoon season, the Rz and Tz were in the eutrophic condition based on TSI (TP) values because of nutrient runoff from the watershed. Based on mean TSI (CHL-a) throughout the year, the reservoir was in the eutrophic condition except in Lz during the premonsoon season. Average TSI (SD) values display a mesotrophic state in the reservoir during the premonsoon and postmonsoon seasons. During monsoon season based on average TSI (SD), Rz (57.07), and IT1 (50.42) they showed a eutrophic state.  The analysis of trophic state index deviation (TSID) indicated that blue-green algae were dominated in the reservoir during the premonsoon, monsoon, and postmonsoon seasons based on the relationship between TSI (CHL-a)-TSI (SD) and TSI (CHL-a)-TSI (TP) (Figure 8). During the monsoon and postmonsoon seasons, the concentrations of smaller particles, non-algal turbidity, and zooplankton grazing were minor in Rz, Tz, Lz IT1, and IT2. The yearly TSID analysis points out that, during the premonsoon season, the blue-green algae was dominated over the study period ( Figure  9). Whereas, during monsoon and postmonsoon seasons, the zooplankton grazing happened in the years of 2016 and 2018, respectively. The analysis of trophic state index deviation (TSID) indicated that blue-green algae were dominated in the reservoir during the premonsoon, monsoon, and postmonsoon seasons based on the relationship between TSI (CHL-a)-TSI (SD) and TSI (CHL-a)-TSI (TP) (Figure 8). During the monsoon and postmonsoon seasons, the concentrations of smaller particles, non-algal turbidity, and zooplankton grazing were minor in Rz, Tz, Lz IT1, and IT2. The yearly TSID analysis points out that, during the premonsoon season, the blue-green algae was dominated over the study period ( Figure 9). Whereas, during monsoon and postmonsoon seasons, the zooplankton grazing happened in the years of 2016 and 2018, respectively. Water 2020, 12, x FOR PEER REVIEW 20 of 26

Discussion
The need for the present study was aroused due to the presence of blue-green algae around the year as well as the fact that the reservoir is in a eutrophic state based on TSI (CHL-a) except for Lz during the premonsoon season. The Daecheong Reservoir is susceptible to nutrients and organic pollution due to urbanization, industrialization, and intensive agricultural farming [7]. The study of spatial and temporal variations of water chemistry, algal chlorophyll, and trophic state indexes has

Discussion
The need for the present study was aroused due to the presence of blue-green algae around the year as well as the fact that the reservoir is in a eutrophic state based on TSI (CHL-a) except for Lz during the premonsoon season. The Daecheong Reservoir is susceptible to nutrients and organic pollution due to urbanization, industrialization, and intensive agricultural farming [7]. The study of spatial and temporal variations of water chemistry, algal chlorophyll, and trophic state indexes has great significance and helps to characterize the water quality problems in the Daecheong Reservoir for sustainable water resource management.
The pattern of the spatial and temporal distribution of water quality parameters showed heterogeneity and was regulated by the hydrological parameters, zonations of the reservoir, land-use patterns, and tributary streams [8,9,24,41]. Previous studies on lentic systems have found that reservoirs receive most external phosphorus and nitrogen loads from their watersheds and tributaries [4,20,24,42]. Therefore, the water quality of the tributaries determines the reservoir water quality [41]. The chemical conditions of the riverine and transitional zones of the Daecheon Reservoir are controlled by the Geum River, Soung Stream, Seowha Stream, Janggae Stream, Whebuk Stream, and wastewater disposal plants. The concentrations of TP and TN were higher in Rz and Tz than in Lz, IT1, and IT2 due to the external loads from the tributary streams and watersheds. Reservoir studies have shown that increased external loads of TP and TN are responsible for the higher algal production and decreased water clarity [24,43,44]. Our present study indicated that the higher algal growth and lower SD was observed in Rz. Reservoir water quality varied from the riverine zone to the lacustrine zone, which is similar to our present findings [5,6].
Summer monsoon was the key determinant for the nutrient runoff from the catchment and algal growth in the reservoir [10]. During the flood year (2011), the concentrations of TP, and TSS were higher, which means that the intense monsoon brings the nutrients and suspended solids to the reservoir. It directly affects the algal growth of the reservoir. In the lacustrine zone, the P-limitation was severe during both drought and flood year. The hydraulic residence time has a significant impact on ionic concentration and chlorophyll growth in the riverine zone. The concentration of EC and CHL-a was lower during July in the riverine zone due to the dilution effects and shorter residence time.
The sudden increase of flow rate, TP, and TSS in summer has unique features for the seasonal changes of aquatic environments in Korean streams and reservoirs [24]. The nutrient concentrations and chlorophyll growth were higher during the monsoon period in the mainstem rather than in the embankment of the reservoir. The phosphorus loadings in the reservoir are larger during summer rather than other seasons and are expected to bloom in the reservoir [20].
Long-term trend analysis of water quality parameters indicated that TP was showing a decreasing trend in Rz and did not show any trend in other parts of the reservoir. The external load of TP may be affected by the upper Yongdam dam in the Geum river watershed. It can alter the trophic functions in the Daecheong reservoir. Ignole and An found that the Yongdam dam slightly alters pH, water temperature, and EC in the Daecheong Reservoir systems [5]. This finding is provisional and should be validated through further research. One remarkable characteristic of the Daecheong reservoir is that TN concentrations were higher than 1 mgL −1 , which means it is a nitrogen-rich system. All over the world, most temperate lakes are nitrogen-rich systems [20,45]. The trend analysis of TN was showing an increasing trend in the reservoir except for IT1. The ionic (EC) and organic content (BOD and COD) of the Daecheong reservoir were showing an increasing trend, which is similar to the previous study of Ignole and An (2016) [5]. We found a downward trend for algal growth in the Daecheong Reservoir except for IT1. The present findings concurred with some previous studies in Korean reservoirs [5,46].
An empirical model based on CHL-a-TP and CHL-a:TN presented that the reservoir mean TP determined CHL-a rather than TN. Regression analysis of CHL-a vs. TP showed a strong linear fit, and TP explained 67% of CHL-a variation in the Daecheong Reservoir. The slope (0.98) of the present model was greater than the slope (0.96) of the previous CHL-a-TP model of 59 Korean reservoirs [4]. The value of R 2 (0.67) in this reservoir was lower than the R 2 values from the earlier nationwide model (0.77) and a single lake model (0.88) [10]. Despite slight differences of slope and R 2 value, the CHL-a-TP relationship reported here supports the view that algal growth in lentic systems respond to P enrichment, and that TP may provide a reliable basis for predicting the average algal abundance [4,25,46]. Previous studies of temperate and subtropical reservoirs have suggested that summer CHL-a is a function of spring or summer TP [3,25]. The present research observed that summer TP regulates the summer algal growth in Tz (R 2 = 0.24), Lz (R 2 = 0.40), and IT1 (R 2 = 0.18). This scenario had been observed frequently in Asian reservoirs [3,28,29]. In contrast, spring TP influences the summer algal chlorophyll in North American reservoirs [26,27]. Moreover, conditional plotting is one of the most useful tools to determine the limiting nutrients in the aquatic system [47]. Conditional plotting analysis suggested that the Daecheong Reservoir is a P-limited system. Secchi depth (SD) and non-algal light attenuation coefficient (Kna) are used widely to measure underwater light availability in lentic systems [2,23]. Light availability is essential for the algal growth in aquatic systems [48]. Light availability in lentic systems was controlled by the water column depth, nutrients (TP, TN), suspended solids (TSS), and algal biomass [2,23,42]. The water clarity of the reservoir is largely regulated by the TP (R 2 = 0.62), TSS (R 2 = 0.67), and CHL-a (R 2 = 0.68) than TN (R 2 = 0.10). The Secchi depth in the reservoir less than 1 m linked with impaired water quality problems due to organic turbidity and suspended solids [23]. The present study did not show such a low transparency. The light attenuation coefficient (Kna) was positively correlated with TP, TN, TSS, and CHL-a and was largely regulated by TSS (R 2 = 0.52) rather than TP (R 2 = 0.38), TN (R 2 = 0.28), and CHL-a (R 2 = 0.04). The positive relationship between Kna and TP in the Daecheong Reservoir supports the previous research conducted in Mark Twain Lake in Missouri, USA [42]. The empirical relationships among SD, Kna, TP, TSS, and CHL-a should be further examined to determine light availability in the reservoirs.
The blooming of blue-green algae is frequently associated with eutrophic conditions. Analysis of TSID indicated that the blue-green algae were dominated in the reservoir. The non-algal turbidity had happened during the monsoon and postmonsoon period due to the surface runoff from the watershed. A little bit of zooplankton grazing was observed during the study period. From 2000-2018, during the three seasons, the blue-green algae were subjugated except in 2016 and 2018 in the monsoon and postmonsoon seasons, respectively. During the premonsoon, monsoon, and postmonsoon seasons, the mean TSI (CHL-a) of the reservoir was in eutrophic conditions, indicating that algal blooms happened except for Lz. According to the USEPA, the trophic state of the water should be in oligotrophic-to-mesotrophic conditions for drinking purposes [49]. The intake towers of the reservoir are facing substantial bloom problems, giving the water supply unfavorable conditions for the citizens. Moreover, the negative correlation of DO with TSS, WT, TP, BOD, and COD suggested that the organic pollution of the reservoir is persistent, and that chemical, ecological, and environmental degradation had been happening [6].

Conclusions
The Daecheong Reservoir has undergone eutrophication conditions, and the nutrient availability of the reservoir had enough for algal growth. The nutrient and organic content of the reservoir showed disparities from headwater to the dam site and was largely influenced by the external load from tributaries and watersheds. The summer monsoon was also a responsible factor for nutrient dispersion from Rz to Lz and regulated the algal growth in the reservoir. The chlorophyll growth was elevated during both drought and flood year. The algal growth was highly regulated by the TP rather than the TN. The reservoir showed a P-limited system. Analysis of the trophic state index showed that cyanobacterial blooms were dominated in the reservoir throughout the study period. To prevent the outbreak of cyanobacterial blooms and to improve water quality as a source of drinking water, the internal and external loading of TP and TN must be declined from tributaries and watersheds. The present study recommends that ecotechnology or biomanipulation could be used in the Daecheong Reservoir for improving the water quality and to manage the eutrophication.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2073-4441/12/7/2066/s1, Figure S1: Annual variation of chemical water quality parameters in Daecheong Reservoir from 2000-2018. TP: total phosphorus; TN: total nitrogen; EC: electrical conductivity; BOD: biological oxygen demand; COD: chemical oxygen demand; TSS: total suspended solids; CHL-a: chlorophyll-a; and TPR: total precipitation; red lines in the box indicate the mean value. Figure S2: Monthly variation (2000-2018) of physicochemical parameters in embankment (S4 and S6) and mainstem (S3 and S5) of Daecheong Reservoir (TP: total phosphorus, TN: total nitrogen, CHL-a: chlorophyll-a, SD: Secchi depth, and I-XII: January-December). Figure S3: The relationship between CHL-a and TP is plotted conditional on the range of TN. Figure S4: The relationship between CHL-a and TN is plotted conditional on the range of TP.