Monitoring the Spatial and Temporal Variations in The Water Surface and Floating Algal Bloom Areas in Dongting Lake Using a Long-Term MODIS Image Time Series

: Signiﬁcant water quality changes have been observed in the Dongting Lake region due to environmental changes and the strong inﬂuence of human activities. To protect and manage Dongting Lake, the long-term dynamics of the water surface and algal bloom areas were systematically analyzed and quantiﬁed for the ﬁrst time based on 17 years of Moderate Resolution Imaging Spectroradiometer (MODIS) observations. The traditional methods (index-based threshold algorithms) were optimized by a dynamic learning neural network (DL-NN) to extract and identify the water surface area and algal bloom area while reducing the extraction complexity and improving the extraction accuracy. The extraction accuracy exceeded 94.5% for the water and algal bloom areas, and the analysis showed decreases in the algal bloom and water surface areas from 2001–2017. Additionally, the variations in the water surface and algal bloom areas are greatly a ﬀ ected by human activities and climatic factors. The results of these analyses can help us better monitor human contamination in Dongting Lake and take measures to control the water quality during certain periods, which is crucial for future management. Moreover, the traditional methods optimized by the DL-NN used in this study can be extended to other inland lakes to assess and monitor long-term temporal and spatial variations in algal bloom areas and can also be used to acquire baseline information for future assessments of the water quality of lakes.


Introduction
Lakes, which not only enable shipping and supply water resources for human consumption but also improve the local microclimate and ecological environment, represent one of the most important components of the hydrosphere [1,2]. However, lakes are very sensitive to climate change and human factors [3][4][5]. As a result, due to intensive human activities and climate change in recent years, many lakes have decreased in size [6][7][8][9], and as a consequence, nearly half of all

Study Area
Dongting Lake (111°40′ E ~ 113°10′ E, 28°30′ N ~ 30°20′ N), the second largest freshwater lake in China, has a drainage area of approximately 262,800 km 2 and an average water depth of 6.39 m. The lake is located in the middle and lower reaches of the Yangtze River Basin in northeastern Hunan Province (Figure 1). The topography of the land surrounding the lake is high in the west and low in the east, and the lake can be divided into eastern, western and southern sections. Dongting Lake is adjacent to the Yangtze River to the north and four rivers (the Xiangjiang River, Zi River, Yuan River, and Li River) to the south, all of which run through Chenglingji and part of Yueyang and finally drain into the Yangtze River. The Dongting Lake Basin is one of the most industrialized regions in China and is surrounded by Yueyang, Yiyang, Yuanjiang, Changde and other important cities in Hunan Province, as well as many prefecture-level counties and cities with strong economic development, such as Xiangyin and The Dongting Lake Basin is one of the most industrialized regions in China and is surrounded by Yueyang, Yiyang, Yuanjiang, Changde and other important cities in Hunan Province, as well as many prefecture-level counties and cities with strong economic development, such as Xiangyin and Nan Counties. The population density is high in this area (about 22.24 million people in 2017), and economic development has been rapid, which has affected the water surface area and quality of Dongting Lake (Bi et al., 2018). Consequently, the ecological environment of Dongting Lake has deteriorated, and the eutrophication trend of the lake has worsened. In particular, algal bloom outbreaks occurred in the eastern part of Dongting Lake in [2005][2006]. According to statistics from the State Environmental Protection Administration, the overall water quality of Dongting Lake reached grade V in 2005 (National Surface Water Quality Standard GB3838-2002), where a higher grade corresponds to a worse water quality. At the end of 2007, Hunan Province strengthened the water quality protection standards for Dongting Lake and investigated the enterprises around the lake; as a result, more than 200 businesses were penalized to prevent them from introducing additional pollution into the lake. Subsequently, the water quality of Dongting Lake has gradually improved. Nevertheless, although the water quality of Dongting Lake reached grade IV following these actions (GB3838-2002), the water pollution in Dongting Lake increased from 2007 to 2008. Then, at the end of 2014, the Dongting Lake water quality further deteriorated coincident with increasing eutrophication, reaching a water quality worse than grade V (GB3838-2002). In particular, the phosphorus content in the water greatly exceeded the national standards ( GB3838-2002). Therefore, to protect the local water quality and ensure the safe supply of drinking water for the surrounding residents, it is necessary to monitor the changes in water quality in real time.

MODIS and Landsat Data
MODIS is an earth observation system (EOS) instrument that has 36 bands, and the Landsat Enhanced Thematic Mapper Plus (ETM+) sensor has 8 bands with a high spatial resolution. MODIS and ETM+ data are particularly useful because of their global coverage, radiometric resolution, dynamic ranges, and free availability.
In this study, MODIS L1B data were used primarily to extract the water surface and algal bloom areas for Dongting Lake, while 30-m-resolution Landsat 7 ETM+ data were used to validate the MODIS observations. MODIS data acquired by both the Terra satellite and the Aqua satellite over the period from January 2001 to February 2018 were obtained from https://ladsweb.modaps.eosdis.nasa.gov/, and Landsat 7 ETM+ data were obtained from https://earthexplorer.usgs.gov/. Due to the 16-day revisit frequency of Landsat 7, only a few Landsat scenes with 0-20% cloud cover over Dongting Lake were obtained.
We selected 1070 images with low cloud cover over the Dongting Lake area in spring (March, April, and May), summer (June, July, and August), autumn (September, October, and November), and winter (December, January, and February) from 2001-2017; the dataset included 1059 MODIS images from both Terra and Aqua and 11 Landsat 7 ETM+ images.

In Situ Data Acquisition
Water level data were used to verify the accuracy of the extracted water surface areas because the surface area and water level are highly correlated. These data comprising daily measured water level data from 2001 to 2017 were derived from the Chenglingji hydrological station. In this study, water level data corresponding to the dates of the remote sensing images are the primary source for validation.
Precipitation and temperature data (2001-2017) from eight stations (Yueyang, Huarong, Anxiang, Yuanjiang, Xiangyin, Guluo, Jianli, and Nanxian) around Dongting Lake were acquired from the China Meteorological Administration (http://www.cma.gov.cn/2011qxfw/2011qsjgx/). The location of each Remote Sens. 2020, 12, 3622 5 of 31 meteorological station is shown in Figure 1. These data were used to analyze the effects of precipitation and temperature on the water surface and algal bloom areas in Dongting Lake. The temperature and precipitation data used in this study are the average values from these stations.
Field investigations were performed monthly at predefined sites across Lake Dongting from 2001 to 2017 to obtain (1) manual outlines of the area where known aquatic macrophytes exist based on field observations and (2) field observations of nutrients (including total nitrogen (TN) and total phosphorus (TP)) and of chlorophyll-a (Chl-a concentration ≥0.01 mg/L, at which algal blooms accumulate on the surface of the water). These data were used to modify the extraction results of the Landsat 7 ETM+ and MODIS images because the algal blooms depended significantly on variations in the TN and TP concentrations [79].

Methods
The overall methodology is illustrated in Figure 2. We performed Rayleigh correction for the MODIS and Landsat data to generate Rayleigh-corrected reflectance (Rrc); then, traditional methods were used to extract and identify the extents of the surface water and algal blooms. More specifically, time series of the FAI, cyanobacteria and macrophyte index (CMI), and MNDWI products were generated, and the maximum gradient method was used to obtain the threshold employed to extract and identify them. Second, we used the DL-NN to optimize the traditional methods to extract and identify the water surface and algal bloom areas more concisely and accurately. After verifying the extraction accuracy, we analyzed the dynamic changes in the water area and algal blooms throughout Dongting Lake. More details are given in the following sections.
Remote Sens. 2020, 12, x FOR PEER REVIEW 5 of 32 Guanjiapu) around Dongting Lake were acquired from the Changjiang (Yangtze) Sediment Bulletin at http://www.cjw.gov.cn/zwzc/bmgb/. These data were used to analyze the causes of the changes in the water surface area. Field investigations were performed monthly at predefined sites across Lake Dongting from 2001 to 2017 to obtain (1) manual outlines of the area where known aquatic macrophytes exist based on field observations and (2) field observations of nutrients (including total nitrogen (TN) and total phosphorus (TP)) and of chlorophyll-a (Chl-a concentration ≥0.01 mg/L, at which algal blooms accumulate on the surface of the water). These data were used to modify the extraction results of the Landsat 7 ETM+ and MODIS images because the algal blooms depended significantly on variations in the TN and TP concentrations [79].

Methods
The overall methodology is illustrated in Figure 2. We performed Rayleigh correction for the MODIS and Landsat data to generate Rayleigh-corrected reflectance (Rrc); then, traditional methods were used to extract and identify the extents of the surface water and algal blooms. More specifically, time series of the FAI, cyanobacteria and macrophyte index (CMI), and MNDWI products were generated, and the maximum gradient method was used to obtain the threshold employed to extract and identify them. Second, we used the DL-NN to optimize the traditional methods to extract and identify the water surface and algal bloom areas more concisely and accurately. After verifying the extraction accuracy, we analyzed the dynamic changes in the water area and algal blooms throughout Dongting Lake. More details are given in the following sections.

Remote Sensing Data Processing
The MODIS L1B data were destriped using the National Oceanic and Atmospheric Administration (NOAA) destriping code [80,81]. Gaseous absorption and Rayleigh scattering effects were corrected using computer software provided by the MODIS Rapid Response Team based mainly on the 6S Radiative Transfer Model [82]; then, reprojection and geometric correction were performed. The Landsat ETM+ data were georeferenced with an error of less than 0.5 pixels using a  The MODIS L1B data were destriped using the National Oceanic and Atmospheric Administration (NOAA) destriping code [80,81]. Gaseous absorption and Rayleigh scattering effects were corrected using computer software provided by the MODIS Rapid Response Team based mainly on the 6S Radiative Transfer Model [82]; then, reprojection and geometric correction were performed. The Landsat ETM+ data were georeferenced with an error of less than 0.5 pixels using a 1:50,000 topographic map, and the radiative transfer calculations based on the second simulation of the satellite signal in the solar spectrum [82] were used for the atmospheric correction. Regions with clouds were manually outlined in a small number of MODIS and Landsat images, which were then excluded from further analysis.

MNDWI-Based Threshold Algorithm for Detecting Water Pixels
In numerous studies, the MNDWI has been demonstrated to be an effective and high-precision index for extracting the area of surface water [23]. The MNDWI was applied to the MODIS data in this study because it extracts water bodies more accurately than the NDWI [23,83,84]. Noise from built-up land, soil and vegetation can be reduced using the shortwave infrared (SWIR) and green bands, thereby enabling water boundaries to be accurately identified. Generally, the MNDWI values of water areas are positive, whereas the MNDWI values of soil, vegetation, and built-up land are negative. Therefore, the appropriate extent of Dongting Lake can be calculated by selecting an appropriate threshold for the MNDWI. Some binary MNDWI images can improperly delineate the boundaries of areas where clouds and shadows exist; in such cases, manual editing and merging are used to obtain the correct lake boundary. The MNDWI is expressed as follows: where R (1650) is the reflectance in the SWIR band, which corresponds to the 6th MODIS band and the 5th Landsat 7 ETM+ band, and R (560) is the reflectance in the green band, which corresponds to the 4th MODIS band and the 2nd Landsat 7 ETM+ band.

FAI and CMI-Based Threshold Algorithm for Detecting Algal Blooms
Water strongly absorbs solar light at red, NIR, and SWIR wavelengths; hence, in these bands, the reflectance of water is very low (close to zero). However, algal blooms and other organisms, such as aquatic macrophytes, on the surfaces of water bodies have higher reflectance in the NIR band than at other wavelengths and thus can be easily distinguished from the surrounding waters [85,86]. The FAI (Equation (2)) is defined as the difference between the Rayleigh-corrected reflectance (Rrc) in the NIR band (859 nm) and the calculated NIR value based on the linear interpolation of the Rayleigh-corrected reflectances in the red (645 nm) and SWIR (1240 nm) bands [87]. The FAI was used to identify algal blooms according to the differences in spectral characteristics between the algal blooms and surrounding surface water. Please note that Dongting Lake contains abundant aquatic macrophytes during April to September; in areas that are known to feature aquatic macrophytes, the FAI is actually an indicator of these macrophytes rather than of algal blooms. Considering that there are significant differences in spectral characteristics between these two organisms at blue, green, and SWIR wavelengths [85,86], the CMI (Equation (3)) was used to distinguish them. To ensure the accuracy of the CMI-based identification, we manually outlined the area that is known to feature aquatic macrophytes based on field observations and an aquatic macrophyte map of Dongting Lake to modify the aquatic macrophyte extraction results. This approach was used to exclude pixels representing aquatic macrophytes in the FAI images.

Threshold Determination
The key step in distinguishing intense algal blooms from water and other objects is to determine an appropriate threshold value in the imagery. An experienced analyst can select the most appropriate threshold by using trial and error together with visual interpretation, but this process is influenced by subjective factors, which can have a certain impact on the extraction of the target. Therefore, gradients and statistics of the FAI (CMI or MNDWI) were used to obtain the thresholds for the FAI (CMI or MNDWI). A pixel gradient was defined as the sum of the absolute value differences between each pixel in the image and the adjacent pixels in a 3 × 3 pixel window. The equation can be formulated as follows: where g is a pixel gradient in the MNDWI, CMI or FAI image; X i is the ith pixel value around the center pixel of a 3 × 3 pixel window; and X c is the center pixel of a 3 × 3 pixel window. This method can generate a gradient image for each FAI (CMI or MNDWI) image. However, the FAI values of particles and thick algae scum in a water body are also different from the FAI values in the surrounding region. Therefore, before generating the FAI gradient image, FAI values less than −0.01 were used to remove pixels containing water with high concentrations of particles, and FAI values greater than 0.02 were used to remove pixels with thick algae scum [63,64]. For each gradient image, the maximum gradient value was determined, and the FAI (CMI or MNDWI) value corresponding to the maximum gradient was obtained. Theoretically, this value is the threshold for extracting the target from the corresponding image. However, to rule out the effects of outliers, the maximum gradient value was not obtained from a single pixel but from a group of pixels around the interface between the target and surrounding pixels. This group of pixels includes the pixels associated with the maximum gradient mode determined by the gradient image histogram. The mean value of this group of pixels in the FAI (CMI or MNDWI) image was chosen as the threshold value to separate the target from the nontarget. Therefore, a different threshold was used for each MNDWI image to extract water bodies. However, due to the small accumulation areas of algal blooms or aquatic macrophytes in some images, the number of pixels pooled in the histogram was too low to generate meaningful statistical values. Consequently, the threshold obtained with this method was not useful for identifying algal blooms or aquatic macrophytes. In this study, two methods were used to obtain the FAI and CMI thresholds for the entire image series. (1) The threshold was obtained based on an individual image when the algal bloom or aquatic macrophyte patches in the image were large. The mean FAI or CMI value of the maximum gradient of pixels was used as the threshold for extracting the target. (2) The threshold was obtained based on multi-image statistics when the algal bloom or aquatic macrophyte patches in the image were small. The mean and standard deviation of the FAI (CMI) values in all other images were calculated, and the difference between the mean value and the standard deviation was used as the threshold for extracting algal blooms (aquatic macrophytes).

DL-NN
The DL-NN is one of the most popular pattern recognition methods using forward and backward propagation, thereby processing multidimensional image data with an input layer, an output layer, hidden layers, and activation functions [88]. The hidden layers extract the internal features of inputs (these features can be either the actual input to the system or the output from other neurons in the Remote Sens. 2020, 12, 3622 8 of 31 preceding layers) by a weighted summation of the input vector and learnable weights plus the bias term [89]. Then, each hidden layer is equipped with activation functions to provide nonlinear variations in the elements. Furthermore, activation functions can take many forms. The most common activation function is the nonlinear sigmoid function given by where w represents the weights (parameters), ς is the bias term, and x represents the input characteristic data for a previous layer. The DL-NN output layer uses the mean square error (MSE) as a loss function, and the functional equation for a layer is expressed as where y i is the target of the ith data in a training batch, and y i is the retrieved value given by the DL-NN. Forward and backward propagation constitutes the main process of dynamic learning optimization. During the forward propagation phase, training patterns are input and transmitted to the output layer through a hidden layer transformation. The output value in the forward propagation stage is calculated and compared with the target value to obtain the corresponding error. During the backward propagation phase, after computing the gradient vector of the loss function, the gradient descent algorithm is used to properly adjust the weight vector according to the gradient and learning rate so that the DL-NN output value of the training sample is as close as possible to the target value [85]. In this study, we used the adaptive moment estimation (Adam) algorithm as the gradient descent algorithm in the DL-NN backward propagation stage [90].

Training and Testing the DL-NN
The extents of the MNDWI of water bodies, FAI of algal blooms, and CMI of aquatic macrophytes were first acquired based on the results of the traditional methods. The MNDWI values of water bodies, FAI values of algal blooms, and CMI values of aquatic macrophytes at different times and places were then collected as sample-dependent variables as inputs for the DL-NN. Simultaneously, the MODIS L1B data were used as an independent variable to train and test the DL-NN. Then, the DL-NN was utilized to extract features from the MODIS images to obtain the MNDWI, FAI and CMI values of targets. To ensure the proper implementation of the DL-NN algorithm, the DL-NN input samples comprised the Rayleigh-corrected reflectance bands and their linear variation. The input samples of each index (MNDWI, FAI and CMI) covered a minimum of 61,231 sets (MNDWI-121490, FAI-61231, CMI-82657) for reflecting the site-and season-specific water bodies, algal blooms, and aquatic macrophytes. The input samples were randomly divided into a training dataset accounting for 70% of the total data and a test dataset consisting of 30% of the total data. Considering that different targets have different responses to remote sensing reflectance and vegetation index, explanatory variables with strong responses to water bodies, algal blooms, and aquatic macrophytes were chosen for the prediction of this indexes [63,85,86]. The DL-NN employed to obtain the MNDWI was trained using Rrc(1640), Rrc(555), Rrc(859), NDVI, and (Rrc(859)−Rrc(555))/(Rrc(859) + Rrc(555)); while that utilized to calculate the FAI was trained by Rrc  (2130). After continuous trials, the DL-NN that obtained the MNDWI exhibited the best accuracy when the number of hidden layers was four and the number of hidden nodes was 300; the mean error in the MNDWI was 4.32 × 10 −3 , and the standard deviation of the fit was approximately 4.18 × 10 −3 . The DL-NN used to obtain the FAI was the most accurate when the number of hidden layers was three and the number of hidden nodes was 100; the mean error in the FAI was 2.72 × 10 −4 , and the standard deviation of the fit was approximately 1.77 × 10 −4 . Finally, the DL-NN utilized to obtain the CMI was the most accurate when the number of hidden layers was three and the number of hidden nodes was 300; the mean error in the CMI was 4.27 × 10 −4 , and the standard deviation of the fit was approximately 2.49 × 10 −4 . The details regarding the training and testing datasets can be found in Tables 1-3.

Accuracy Assessment
To verify the accuracy of the extracted water surface and algal bloom areas for Dongting Lake by the two methods, we used concurrent higher-resolution Landsat 7 ETM + images combined with ground measurements. First, the conventional index-based threshold algorithms were used to extract the water surface and algal bloom areas from 11 concurrent Landsat 7 ETM + images. Second, visual interpretations and field investigations were combined to modify the extraction results from the Landsat 7 ETM + images. Finally, the modified results were used to validate the extraction results of the MODIS data. Simultaneously, we compared the MODIS extraction results with historical survey data to further verify our results.
The accuracy of the extracted water surface and algal bloom areas was assessed using the root mean square error (RMSE), which is used to measure the deviation between the observed and true values, the relative error (RE) between the actual and observed values, and the Pearson correlation coefficient (r), which is widely used to measure the closeness between two variables, namely, X and Y.
where S M is the extracted area, and S T is the observed area.
where n is the sample size,ŷ i is the ratio of the extracted area to the total area of Dongting Lake, and y i is the ratio of the observed area to the total area of Dongting Lake.
where n is the sample size, x i is the extracted area, y i is the observed area, and x and y are the means of x and y, respectively.

Comparison with Concurrent Landsat 7 ETM + Observations
We calculated the RE, RMSE and r values between the extracted area and observed area. The extraction accuracies (1-RE) of the conventional methods for the water surface and algal bloom areas were 93.12% and 94.11%, respectively, and the corresponding RMSE values were 4.32% and 0.76%, respectively. The extraction accuracies of the traditional methods optimized by the DL-NN for extracting the areas of surface water and algal blooms were greater than 94.53% and 94.66%, respectively, and the corresponding RMSE values were 3.60% and 0.71%, respectively. In addition, the r values between the water surface and algal bloom areas with the optimized methods were superior to those with the traditional methods, as determined by the fit between the extracted area and observed area ( Figure 3). The optimized methods can improve the extraction accuracy of surface water and algal blooms. Most importantly, the DL-NN overcomes the shortcomings of traditional methods that need to take a global threshold value for each image. Therefore, in this paper, the extraction results from the optimized methods were used to analyze the changes in the water surface and algal bloom areas in Dongting Lake from 2001 to 2017.

Comparison with Water Level Data
Water level time series from in situ stations were used to further verify the water extraction accuracy of the traditional methods optimized by the DL-NN in Dongting Lake because the surface area and water level are highly correlated. The dates of the water level data correspond to the dates of the MODIS images. The comparisons indicated very good agreement between the water levels and surface areas ( Figure 4). We also computed the r values between the water surface area and water level in Figure 4, which displays the resulting scatter plot of the correlations (r = 0.934). These correlations further indicate that the results of the optimized methods for the extraction of surface water are accurate.

Comparison with Water Level Data
Water level time series from in situ stations were used to further verify the water extraction accuracy of the traditional methods optimized by the DL-NN in Dongting Lake because the surface area and water level are highly correlated. The dates of the water level data correspond to the dates of the MODIS images. The comparisons indicated very good agreement between the water levels and surface areas (Figure 4). We also computed the r values between the water surface area and water level in Figure 4, which displays the resulting scatter plot of the correlations (r = 0.934). These correlations further indicate that the results of the optimized methods for the extraction of surface water are accurate. area and water level are highly correlated. The dates of the water level data correspond to the dates of the MODIS images. The comparisons indicated very good agreement between the water levels and surface areas (Figure 4). We also computed the r values between the water surface area and water level in Figure 4, which displays the resulting scatter plot of the correlations (r = 0.934). These correlations further indicate that the results of the optimized methods for the extraction of surface water are accurate.   The water surface area of Dongting Lake changed seasonally; the lake area was small in spring and winter and larger in spring (approximately 1000 km 2 ) than in winter (550-650 km 2 ). The lake area was large in both summer and autumn but peaked in summer; the largest area of 2114.32 km 2 occurred in the summer of 2002. In addition, the surface area of Dongting Lake in spring, summer, autumn, and winter showed a decreasing trend from 2001-2017; the largest rate of decrease was in summer, and the summer variation rate reached −16.49 km 2 /a, followed by a rate of −14.93 km 2 /a in autumn. The rates of decrease were small in spring and winter at −4.05 km 2 /a and −0.88 km 2 /a, respectively. These characteristics are related to precipitation, temperature conditions, the operation of the Three Gorges Reservoir (TGR), and the main flood season of the Yangtze River in the Dongting Lake region [73,91].

Analysis of the Interannual and Seasonal Changes in the Water Surface Area
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 32 The seasonal and interannual changes in the water area of Dongting Lake from 2001 to 2017 are shown in Figure 5. The water area of the lake displayed a declining trend between 2001 and 2017, with a minimum in 2011 at an area of 888.17 km 2 ; the maximum was reached in 2002 with an area of 1249.43 km 2 . The water surface area of Dongting Lake changed seasonally; the lake area was small in spring and winter and larger in spring (approximately 1000 km 2 ) than in winter (550-650 km 2 ). The lake area was large in both summer and autumn but peaked in summer; the largest area of 2114.32 km 2 occurred in the summer of 2002. In addition, the surface area of Dongting Lake in spring, summer, autumn, and winter showed a decreasing trend from 2001-2017; the largest rate of decrease was in summer, and the summer variation rate reached −16.49 km 2 /a, followed by a rate of −14.93 km 2 /a in autumn. The rates of decrease were small in spring and winter at −4.05 km 2 /a and −0.88 km 2 /a, respectively. These characteristics are related to precipitation, temperature conditions, the operation of the Three Gorges Reservoir (TGR), and the main flood season of the Yangtze River in the Dongting Lake region [73,91].

Analysis of the Monthly Spatial Changes
The spatial characteristics of the water surface distribution in Dongting Lake significantly differed over a 1-year period, and the variation characteristics of the water surface exhibited obvious regional differences (Figure 6). The water surface area of eastern Dongting Lake was large in summer and autumn, while only a few ribbons of water remained in winter and spring, when the variation in the water surface area was the largest. The variations in the water surface area in southern Dongting Lake and western Dongting Lake were smaller than those in eastern Dongting

Analysis of the Monthly Spatial Changes
The spatial characteristics of the water surface distribution in Dongting Lake significantly differed over a 1-year period, and the variation characteristics of the water surface exhibited obvious regional differences ( Figure 6). The water surface area of eastern Dongting Lake was large in summer and autumn, while only a few ribbons of water remained in winter and spring, when the variation in the water surface area was the largest. The variations in the water surface area in southern Dongting Lake and western Dongting Lake were smaller than those in eastern Dongting Lake. Similarly, the water surface areas were large in summer and autumn and small in winter and spring. Evidently, the water surface area was dynamically changing. In the dry season (December-February), Dongting Lake lost nearly 53.94% of its surface, the water area was small, and only a few ribbons of water remained (with an area of 598.48 km 2 ). In the rainy season (May-October), the water area was large with an average of 1299.23 km 2 . In addition, the variations in the water area of Dongting Lake in the dry season were opposite to those in the rainy season, and the spatial distributions of the changes in the water area were similar in each lake section. The largest dynamic spatiotemporal changes in the water surface area occurred in eastern Dongting Lake, and the seasonal hydrofluctuation belts were distributed mainly in the southern and western regions of eastern Dongting Lake with an average area of 358.53 km 2 . The second largest seasonal hydrofluctuation belt was in southern Dongting Lake with an average area of 239.86 km 2 , and a seasonal hydrofluctuation belt was distributed primarily in the southern portion of southern Dongting Lake. The area of the seasonal hydrofluctuation belt was the smallest in western Dongting Lake, with an average area of 130.83 km 2 , where the belt was distributed mainly in the area near the lakeshore.
Remote Sens. 2020, 12, x FOR PEER REVIEW 14 of 32 Figure 6. Changes in the spatial distribution of the water area in different months.

Analysis of the Interannual Spatial Changes
For a more detailed understanding of the spatial patterns and regional differences in the water surface changes in different areas, the per pixel frequency of water detection in each season from 2001 to 2017 was calculated, and the slope (k) (Figure 7) is shown. The frequency of water detection was calculated as the number of times water was identified in each pixel from a stack of all water/non-water classifications for each season of the year. From 2001 to 2017, strong interannual variation in spatial water frequency patterns across the Dongting Lake was detected from our pixel-level analysis (Figure 7), revealing striking differences between years. The frequency of water detection decreased significantly in all four seasons in the East Dongting Lake. The frequency of South Dongting Lake decreased in the three seasons including spring, summer and autumn. The frequency decrease of the West Dongting Lake was mainly observed in summer and autumn. The frequency reduction of these three lake areas occurred dominantly in the lakeshore area.

Analysis of the Interannual Spatial Changes
For a more detailed understanding of the spatial patterns and regional differences in the water surface changes in different areas, the per pixel frequency of water detection in each season from 2001 to 2017 was calculated, and the slope (k) (Figure 7) is shown. The frequency of water detection was calculated as the number of times water was identified in each pixel from a stack of all water/non-water classifications for each season of the year. From 2001 to 2017, strong interannual variation in spatial water frequency patterns across the Dongting Lake was detected from our pixel-level analysis (Figure 7), revealing striking differences between years. The frequency of water detection decreased significantly in all four seasons in the East Dongting Lake. The frequency of South Dongting Lake decreased in the three seasons including spring, summer and autumn. The frequency decrease of the West Dongting Lake was mainly observed in summer and autumn. The frequency reduction of these three lake areas occurred dominantly in the lakeshore area.
pixel-level analysis (Figure 7), revealing striking differences between years. The frequency of water detection decreased significantly in all four seasons in the East Dongting Lake. The frequency of South Dongting Lake decreased in the three seasons including spring, summer and autumn. The frequency decrease of the West Dongting Lake was mainly observed in summer and autumn. The frequency reduction of these three lake areas occurred dominantly in the lakeshore area.

Seasonal Variations in the Algal Bloom Area
The statistics of the algal bloom outbreak areas in Dongting Lake during different seasons from 2001 to 2017 (Figure 8) indicate that the algal bloom area in Dongting Lake exhibits a normal distribution within a 1-year period and that most algal blooms occur in the summer and fall seasons. This abundance of algal blooms is associated primarily with the high water temperatures in summer and autumn that provide suitable conditions for algal growth. The area of algal blooms in Dongting Lake was the largest in summer, followed by autumn, and the smallest areas were observed in spring and winter. Between

Seasonal Variations in the Algal Bloom Area
The statistics of the algal bloom outbreak areas in Dongting Lake during different seasons from 2001 to 2017 (Figure 8) indicate that the algal bloom area in Dongting Lake exhibits a normal distribution within a 1-year period and that most algal blooms occur in the summer and fall seasons. This abundance of algal blooms is associated primarily with the high water temperatures in summer and autumn that provide suitable conditions for algal growth. The area of algal blooms in Dongting Lake was the largest in summer, followed by autumn, and the smallest areas were observed in spring and winter.  The regions of algal bloom occurrence in Dongting Lake during different seasons from 2001 to 2017 were statistically analyzed, and distribution maps of the algal blooms in different seasons were obtained (Figure 9). Eastern Dongting Lake experienced more frequent algal blooms than the other two regions of the lake, and the area of algal bloom occurrence was concentrated mainly in the easternmost part of the lake. The largest area occurred in summer, followed by autumn and spring, and the smallest algal bloom area occurred in winter. The largest algal bloom areas occurred near the lakeshore in western and southern Dongting Lake, while the center of the lake exhibited the smallest algal bloom area. Furthermore, western and southern Dongting Lake experienced more frequent The regions of algal bloom occurrence in Dongting Lake during different seasons from 2001 to 2017 were statistically analyzed, and distribution maps of the algal blooms in different seasons were obtained (Figure 9). Eastern Dongting Lake experienced more frequent algal blooms than the other two regions of the lake, and the area of algal bloom occurrence was concentrated mainly in the easternmost part of the lake. The largest area occurred in summer, followed by autumn and spring, and the smallest algal bloom area occurred in winter. The largest algal bloom areas occurred near the lakeshore in western and southern Dongting Lake, while the center of the lake exhibited the smallest algal bloom area. Furthermore, western and southern Dongting Lake experienced more frequent algal blooms in autumn than during the other three seasons. For the whole lake, the algal blooms expanded to the center of Dongting Lake after spring and then decreased after summer.

Interannual Variations in the Algal Bloom Area
The frequencies of algal blooms in different areas of Dongting Lake from 2001 to 2017 were similarly analyzed ( Figure 10). Algal blooms occurred in all 17 years and were distributed primarily near the lakeshore. Blooms occurred very frequently in 2001 with an occurrence rate exceeding 50% near the western shore of eastern Dongting Lake, and eastern Dongting Lake showed the largest area with significant bloom coverage. The frequency of algal blooms was low in 2003 and 2004, mainly at 20-30%, except in the northwestern and southwestern parts of eastern Dongting Lake.

Interannual Variations in the Algal Bloom Area
The frequencies of algal blooms in different areas of Dongting Lake from 2001 to 2017 were similarly analyzed ( Figure 10). Algal blooms occurred in all 17 years and were distributed primarily near the lakeshore. Blooms occurred very frequently in 2001 with an occurrence rate exceeding 50% near the western shore of eastern Dongting Lake, and eastern Dongting Lake showed the largest area with significant bloom coverage. The frequency of algal blooms was low in 2003 and 2004, mainly at 20-30%, except in the northwestern and southwestern parts of eastern Dongting Lake.
The western and northeastern sections of eastern Dongting Lake experienced more frequent algal blooms than the other two sections of the lake, where the algal blooms tended to expand from the shore to the center of the lake. The frequency of algal bloom occurrence decreased gradually from the shore to the interior of the lake area. There are 17 major pollution sources around eastern Dongting Lake, with a daily discharge volume of 126,359 m 3 , accounting for 85.55% of the total daily discharge; accordingly, large quantities of pollutants have been discharged into the lake from agriculture, industries, and aquaculture, resulting in nutrient concentrations that are significantly higher in the eastern part of Dongting Lake than in the western and southern parts of the lake [92,93]. The northeastern region of eastern Dongting Lake is simultaneously located near both the city of Yueyang, which shows considerable growth potential with a large population and widespread industrial activities, and the confluence of eastern Dongting Lake and the Yangtze River, where numerous pollutants accumulate and are released from the lake sediments due to disturbances from the Yangtze River. Therefore, significant algal blooms appeared in the northeastern part of eastern Dongting Lake. In contrast, the western region of eastern Dongting Lake is restricted by its geographical location and is not impacted by the estuary, and the amount of water in this area is low. In addition, the terrain of Dongting Lake is high in the northwest and low in the southeast. The water flows from west to east through Chenglingji into the Yangtze River, and thus, water exchange does not occur frequently in the western part of the lake. Moreover, the water variability is small, and pollutants accumulate more easily in this area than in other areas, resulting in water quality deterioration. This nutrient pattern is the primary reason for the obvious west-to-east decrease in the algal bloom occurrence frequency. The western and northeastern sections of eastern Dongting Lake experienced more frequent algal blooms than the other two sections of the lake, where the algal blooms tended to expand from the shore to the center of the lake. The frequency of algal bloom occurrence decreased gradually from the shore to the interior of the lake area. There are 17 major pollution sources around eastern In western and southern Dongting Lake, algal blooms were distributed sporadically along the shore, and the area of algal bloom occurrence in western Dongting Lake was smaller than that in southern Dongting Lake. Although 19 sources of industrial pollution are situated around southern Dongting Lake, the daily pollutant discharge of these sources is 179,237 m 3 , accounting for only 12.14% of the total daily discharge of Dongting Lake. Furthermore, the economy around southern Dongting Lake is poorly developed, and the nonpoint-source pollution therein is lower than that in eastern Dongting Lake. Therefore, algal blooms in southern Dongting Lake encompass a smaller area than those in eastern Dongting Lake. In contrast, the environment around western Dongting Lake is better than that around the other two areas of the lake. There are only six industrial pollution sources, with a daily discharge of 34,191 m 3 , accounting for only 2.31% of the total daily discharge of Dongting Lake. At the same time, the surrounding population is small, and economic development has been slow, so the algal bloom area was smallest in western Dongting Lake.
Algal blooms in Dongting Lake occur mainly between April and September. To more clearly ascertain the interannual changes in algal blooms, we analyzed the average area of algal blooms in Dongting Lake during April to September from 2001 to 2017. The monthly maxima is more representative of the monthly water bloom coverage than the monthly average value due to the influence of wind on the formation of the surface mats [64,65]. Therefore, on the basis of this argument, the annual bloom coverage for the lake were estimated from the mean of monthly maxima, in which these months were from April to September. Then we also calculated the average of the water area extracted from April to September, and finally a time-series of bloom changes was determined using the annual bloom coverage compared with the average of the water area extracted. The results are presented in Figure 11. The changes in the algal bloom area showed decreasing, increasing, decreasing and stabilizing trends over time; accordingly, the changes in the algal bloom area of Dongting Lake were divided into four periods. First, the algal bloom area of Dongting Lake gradually decreased from 15.

Intensity Analysis of Algal Blooms in the Summer and Autumn Seasons
To further evaluate the intensity of algal blooms, the equal-interval grading method was used to divide the algal blooms into three levels according to the degree of algal bloom severity: low risk (−0.004 ≤ FAI ≤ 0), medium risk (0 ≤ FAI ≤ 0.004) and high risk (0.004 ≤ FAI ≤ 0.2).

Intensity Analysis of Algal Blooms in the Summer and Autumn Seasons
To further evaluate the intensity of algal blooms, the equal-interval grading method was used to divide the algal blooms into three levels according to the degree of algal bloom severity: low risk (−0.004 ≤ FAI ≤ 0), medium risk (0 ≤ FAI ≤ 0.004) and high risk (0.004 ≤ FAI ≤ 0.2).
Because algal blooms occur frequently in many lakes throughout China during summer and autumn [94], three different risk levels of algal blooms were analyzed during these two seasons from 2001 to 2017 (Figure 12). Of all the summers during the study period, 2007 was the worst year in terms of the algal bloom abundance. The summer high-risk area in Dongting Lake was the largest in 2007, reaching 366.74 km 2 , after which the area of algal blooms declined. In the autumns between 2001 and 2017, the most serious outbreak of algal blooms in Dongting Lake occurred in 2007. Hence, the autumn high-risk area in Dongting Lake was the largest in 2007, reaching 228.99 km 2 , after which the area of algal blooms in Dongting Lake obviously decreased. A comparison between the algal bloom areas in the summer and autumn seasons of the same year reveals that the bloom area at the same risk level was larger in summer than in autumn for most years; the area in summer was smaller than that in autumn only in 2008. Generally, algal blooms exhibit a larger range of outbreaks in summer than in autumn, but this range is affected by precipitation and other factors, which can result in an increase in the area of algal blooms in autumn. Figure 11. Percentage of the algal bloom area relative to the entire lake area from 2001-2017. (K is the annual change in Y, and Y is the ratio of the bloom area to the area of Dongting Lake).

Intensity Analysis of Algal Blooms in the Summer and Autumn Seasons
To further evaluate the intensity of algal blooms, the equal-interval grading method was used to divide the algal blooms into three levels according to the degree of algal bloom severity: low risk (−0.004 ≤ FAI ≤ 0), medium risk (0 ≤ FAI ≤ 0.004) and high risk (0.004 ≤ FAI ≤ 0.2).
Because algal blooms occur frequently in many lakes throughout China during summer and autumn [94], three different risk levels of algal blooms were analyzed during these two seasons from 2001 to 2017 (Figure 12). Of all the summers during the study period, 2007 was the worst year in terms of the algal bloom abundance. The summer high-risk area in Dongting Lake was the largest in 2007, reaching 366.74 km 2 , after which the area of algal blooms declined. In the autumns between 2001 and 2017, the most serious outbreak of algal blooms in Dongting Lake occurred in 2007. Hence, the autumn high-risk area in Dongting Lake was the largest in 2007, reaching 228.99 km 2 , after which the area of algal blooms in Dongting Lake obviously decreased. A comparison between the algal bloom areas in the summer and autumn seasons of the same year reveals that the bloom area at the same risk level was larger in summer than in autumn for most years; the area in summer was smaller than that in autumn only in 2008. Generally, algal blooms exhibit a larger range of outbreaks in summer than in autumn, but this range is affected by precipitation and other factors, which can result in an increase in the area of algal blooms in autumn.  The area of a cyanobacteria algal bloom outbreak is not a good indicator of the severity of an outbreak because the outbreak severity is also related to the total lake area. Therefore, it is more reasonable to use the ratio of the algal bloom outbreak area to the total water body area at the corresponding time. Once again, we analyzed the high-risk area of algal blooms in Dongting Lake in summer and autumn between 2001 and 2017 ( Figure 13). The pollution level was higher in the summer of 2007 than during the summers of other years, with the high-risk area accounting for 24.56% of the total water area. The pollution level was low in 2003 and 2004, and the algal bloom area was less than 90 km 2 , with the high-risk area accounting for 4.12% and 6.23% of the water area, respectively. The pollution level was higher in autumn 2007 than during the autumns of other years, with the high-risk area accounting for 22.29% of the total water area, while the pollution level was low in 2011, with the high-risk area accounting for 8.29% of the water area. On 1 April 2007, the government shutdown 234 polluting enterprises around Dongting Lake, resulting in a marked decrease in pollutants being discharged into the lake. Consequently, the ratio of the high-risk area to the total lake area in summer decreased from 24.56% in 2007 to 13.06% in 2017; the same ratio in autumn decreased from 22.29% in 2007 to 12.91% in 2017. Remarkably, the rate of decline in the high-risk area in summer was higher than that in autumn from 2007 to 2008. This result is related to the floods that occurred in Dongting Lake in late April 2008, resulting in a much larger water area in summer than in autumn. decrease in pollutants being discharged into the lake. Consequently, the ratio of the high-risk area to the total lake area in summer decreased from 24.56% in 2007 to 13.06% in 2017; the same ratio in autumn decreased from 22.29% in 2007 to 12.91% in 2017. Remarkably, the rate of decline in the high-risk area in summer was higher than that in autumn from 2007 to 2008. This result is related to the floods that occurred in Dongting Lake in late April 2008, resulting in a much larger water area in summer than in autumn.

Environmental Drivers of Interannual Water Body and Algal Bloom Dynamics
In this section, the water area represents the average annual area of Dongting Lake, and the temperature and rainfall are the average values of the annual temperature and annual precipitation from the eight meteorological stations located around Dongting Lake.
To determine the possible linkage of the interannual water body dynamics with local weather, several meteorological parameters of the past 17 years were plotted ( Figure 14). In 2002, and 2012, when the water area reached a maximum, rainfall simultaneously reached a maximum, and the temperature was at a minimum. In 2004In , 2007In , 2011, and 2017, the water areas and rainfall levels were at a minimum, and the temperature was at the maximum in the study period. In addition, we separately calculated the Pearson correlation coefficients between the water area and the temperature and precipitation during this period. There is positive correlation between precipitation and the water area, and it is negatively correlated with temperature. However, from 2001 to 2017, the water area of Dongting Lake decreased gradually, while the temperature and precipitation showed a very slight upward trend. Thus, in some extreme years, such as drought or rain, the water area of Dongting Lake is closely related to the local weather [95]. However, as a

Environmental Drivers of Interannual Water Body and Algal Bloom Dynamics
In this section, the water area represents the average annual area of Dongting Lake, and the temperature and rainfall are the average values of the annual temperature and annual precipitation from the eight meteorological stations located around Dongting Lake.
To determine the possible linkage of the interannual water body dynamics with local weather, several meteorological parameters of the past 17 years were plotted ( Figure 14). In 2002, and 2012, when the water area reached a maximum, rainfall simultaneously reached a maximum, and the temperature was at a minimum. In 2004In , 2007In , 2011, and 2017, the water areas and rainfall levels were at a minimum, and the temperature was at the maximum in the study period. In addition, we separately calculated the Pearson correlation coefficients between the water area and the temperature and precipitation during this period. There is positive correlation between precipitation and the water area, and it is negatively correlated with temperature. However, from 2001 to 2017, the water area of Dongting Lake decreased gradually, while the temperature and precipitation showed a very slight upward trend. Thus, in some extreme years, such as drought or rain, the water area of Dongting Lake is closely related to the local weather [95]. However, as a whole, the shrinkage of the lake area during the study period could not be explained by either temperature or precipitation. The impact of human activities can be inferred to be the most important cause for the reduction in the water area of DongtingLake [91].
Remote Sens. 2020, 12, x FOR PEER REVIEW 21 of 32 whole, the shrinkage of the lake area during the study period could not be explained by either temperature or precipitation. The impact of human activities can be inferred to be the most important cause for the reduction in the water area of DongtingLake [91]. To determine the possible linkage of the interannual algal bloom dynamics with local weather and water quality data (TN, TP), several water quality parameters from in situ field of the past 17 years were plotted ( Figure 15) and principal component analysis (PCA) was used to investigate the combined effects of water quality and local weather on the growth of algal blooms more To determine the possible linkage of the interannual algal bloom dynamics with local weather and water quality data (TN, TP), several water quality parameters from in situ field of the past 17 years were plotted ( Figure 15) and principal component analysis (PCA) was used to investigate the combined effects of water quality and local weather on the growth of algal blooms more quantitatively in 2001-2003, 2003-2009, 2009-2013, and 2013-2017 (Figure 16). In PCA ordination diagram, the correlation is shown as the cosine of the angle between any two axes in the figure, while the influence intensity is the distance one axis projected on another, and the distance can quantify how strong the correlation between the two is (the longer length, the stronger the correlation) [96]. Figure 16a shows that the algal blooms in 2001-2003 were mainly affected by the positive influence of TN and temperature. During this period, the annual bloom area of Dongting Lake decreased by 3.66%. This reduction may be due to the many floods that affected the Dongting Lake region between 1996 and 1999, during which surface runoff carried a large amount of nutrients into Dongting Lake, resulting in eutrophication [92]. Consequently, the 2001 algal bloom outbreak in Dongting Lake was quite serious. Subsequently, as the temperature dropped ( Figure 14) and the nutrients (TN) in the water decreased ( Figure 15) in 2001-2003, the water bloom area of Dongting Lake also decreased. In 2003-2009 (Figure 16b), algal blooms were mainly affected by the positive influence of TN and temperature, and the negative influence of precipitation, while TP had a small influence. During this period, secondary industries in the Dongting Lake surrounding cities, especially the machinery, chemical, and metal processing industries, developed somewhat rapidly, leading to the introduction of additional pollutants, and an increase in the direct discharge of domestic sewage [97]. In addition, the annual growth rate of grain crops, fertilizers and pesticides in Dongting Lake area is 1:4.7:14.6 in 2003-2009, resulting in serious agricultural nonpoint-source pollution [93]. As a result, the water quality of Dongting Lake worsened, and many coastal areas suffered from agricultural nonpoint-source pollution and industrial pollution, leading to algal blooms. The changes of algal blooms from 2009 to 2013 (Figure 16c) were mainly affected by the positive influence of TN and TP. After 2009, the algal bloom outbreak area in Dongting Lake gradually decreased, with the algal bloom area reached 9.91% in 2013 as a result of the numerous measures implemented by the Hunan Provincial Government beginning in 2008. These measures shut down many polluting enterprises in the region and increased the protection and management policies concerning Dongting Lake. As a result, the water quality of Dongting Lake improved significantly [97] ( Figure 15), and the algal bloom area gradually decreased. The influencing factors of algal blooms variation from 2013 to 2017 (Figure 16d) were consistent with those in the previous years and were also mainly affected by the positive influence of TN and TP. The phase was a relatively stable period, during which no significant increase or decrease in TN and TP showed in Dongting Lake, and there was no significant increase or decrease in the algal bloom area.
Remote Sens. 2020, 12, x FOR PEER REVIEW 22 of 32 period, during which no significant increase or decrease in TN and TP showed in Dongting Lake, and there was no significant increase or decrease in the algal bloom area.

Environmental Drivers of Monthly Water Body and Algal Bloom Dynamics
The average monthly water surface area of Dongting Lake from 2001 to 2017 coincided well with the rainfall pattern in different seasons. According to Figure 17, large water surface area

Environmental Drivers of Monthly Water Body and Algal Bloom Dynamics
The average monthly water surface area of Dongting Lake from 2001 to 2017 coincided well with the rainfall pattern in different seasons. According to Figure 17, large water surface area occurred between July and September during the second half of the wet season; the average water surface area of Dongting Lake reached 1325.87 km 2 during the wet period. In addition, small water surface area occurred between November and February during the dry season, with an average water area of 593.46 km 2 . Over the past 17 years, the average water area of Dongting Lake during the wet season was 2.2 times greater than that during the dry season. From 2001 to 2017, the Dongting Lake water surface area changed greatly in August, September, and November, with the surface area fluctuating by more than 320 km 2 . In addition, the Dongting Lake surface area changed less in January and March (fluctuating by only 142 km 2 ) than in other months. Compared with the monthly precipitation climatology for the 17-year study period, it suggests that the significant seasonality in the water surface area might be primarily controlled by precipitation. We also calculated the Pearson correlation coefficients between the monthly water area and precipitation during this period, with a Pearson correlation coefficients of 0.62 (p-values < 0.01), which also indicated that the precipitation is an important cause for the monthly water body dynamics.
precipitation climatology for the 17-year study period, it suggests that the significant seasonality in the water surface area might be primarily controlled by precipitation. We also calculated the Pearson correlation coefficients between the monthly water area and precipitation during this period, with a Pearson correlation coefficients of 0.62 (p-values < 0.01), which also indicated that the precipitation is an important cause for the monthly water body dynamics. To determine the possible linkage of the monthly algal bloom dynamics with local weather and water quality data (TN, TP), principal component analysis (PCA) was used to investigate the combined effects of water quality and local weather on the growth of algal blooms more quantitatively. Figure 18b shows that the monthly algal bloom dynamics in 2001-2017 were mainly affected by the positive influence of temperature or precipitation, while TP and TN denote an insignificant relationship. Therefore, we further studied the relationship between meteorological factors and monthly bloom dynamics. The average precipitation, average temperature and average bloom area during each month from 2001-2017 were calculated. Figure 18a shows that the changes in the algal bloom area in Dongting Lake were consistent with the precipitation and temperature trends, showing a trend of first increasing and then decreasing, which further indicates that the extents of algal blooms in Dongting Lake throughout the year are affected by precipitation and the water temperature. Theoretically, sunlight is limited when precipitation is abundant, and algae in water exhibit strong phototaxis; thus, cloudy and rainy weather is not conducive to the formation of algal blooms or algal floating. Simultaneously, the water area increases with increased precipitation, and the lake water becomes diluted; correspondingly, the concentration of algal blooms is reduced [98,99]. A decrease in the algal bloom area was detected in the remote sensing images. However, the result of this study show that the monthly algal bloom dynamics in Dongting Lake is affected by the positive influence of precipitation. This positive influence may be because increased precipitation increases runoff into the lake, resulting in nonpoint-source pollutants flowing into Dongting Lake, thereby providing favorable conditions for algae growth [92,93]. To determine the possible linkage of the monthly algal bloom dynamics with local weather and water quality data (TN, TP), principal component analysis (PCA) was used to investigate the combined effects of water quality and local weather on the growth of algal blooms more quantitatively. Figure 18b shows that the monthly algal bloom dynamics in 2001-2017 were mainly affected by the positive influence of temperature or precipitation, while TP and TN denote an insignificant relationship. Therefore, we further studied the relationship between meteorological factors and monthly bloom dynamics. The average precipitation, average temperature and average bloom area during each month from 2001-2017 were calculated. Figure 18a shows that the changes in the algal bloom area in Dongting Lake were consistent with the precipitation and temperature trends, showing a trend of first increasing and then decreasing, which further indicates that the extents of algal blooms in Dongting Lake throughout the year are affected by precipitation and the water temperature. Theoretically, sunlight is limited when precipitation is abundant, and algae in water exhibit strong phototaxis; thus, cloudy and rainy weather is not conducive to the formation of algal blooms or algal floating. Simultaneously, the water area increases with increased precipitation, and the lake water becomes diluted; correspondingly, the concentration of algal blooms is reduced [98,99]. A decrease in the algal bloom area was detected in the remote sensing images. However, the result of this study show that the monthly algal bloom dynamics in Dongting Lake is affected by the positive influence of precipitation. This positive influence may be because increased precipitation increases runoff into the lake, resulting in nonpoint-source pollutants flowing into Dongting Lake, thereby providing favorable conditions for algae growth [92,93].

Potential Impact of the Three Gorges Dam on the Water Area Changes
The changes in the water and sediment discharge regimes in the middle and lower reaches of the Yangtze River due to the operation of the Three Gorges Dam (TGD) have inevitably led to

Potential Impact of the Three Gorges Dam on the Water Area Changes
The changes in the water and sediment discharge regimes in the middle and lower reaches of the Yangtze River due to the operation of the Three Gorges Dam (TGD) have inevitably led to changes in the water level and surface area of Dongting Lake [100][101][102]. The volume of water trapped by the TGD often changes the flow regime and leads to an overall reduction in the water level in the downstream river reaches, thereby affecting the water surface areas and water levels of downstream lakes [91,103]. The water level at the Chenglingji hydrological station and the lake water area during the pre-TGD (2001-2003) and post-TGD (2004-2017) periods are shown in Figure 19. The variations in the annual minimum water surface area (WS min ) and annual minimum water level (WL min ) displayed no clear pattern. However, the annual maximum water surface area (WS max ) and annual maximum water level (WL max ) showed decreasing trends. Furthermore, the difference between WS min and WS max and the difference between WL min and WL max in the post-TGD period were lower than those in the pre-TGD period. The annual water surface area and annual water level in the post-TGD period were lower than those in the pre-TGD period. Specifically, the annual mean surface area of Dongting Lake decreased by 13.3% (from 1177 km 2 in the pre-TGD period to 1034 km 2 in the post-TGD period). Chenglingji recorded an annual mean water level of 25.41 m during the pre-TGD period, which decreased by 0.43 m to 24.98 m in the post-TGD period. The fluctuations in the annual average water level and annual average water surface area were small (0.5 m and 100 km 2 ) before TGD operations began. After the completion of the TGD, these fluctuations increased significantly, with the annual average water level and water surface area reaching 2.0 m and 300 km 2 , respectively. In addition, after the impoundment of the TGD, the supplies of water from the Songzi, Hudu, and Ouchi channels and the Xiang, Zi, Yuan, and Li tributaries were significantly reduced. Since the operation of the TGD began, the annual mean runoff in the three channels has decreased from 625 ×  (Table 4). After the impoundment of the TGD, the streamflow decreased sharply in the middle and lower reaches of the Yangtze River, contributing to the shrinkage of Dongting Lake [104]. In addition to being impacted by the abovementioned reduction in the water supply, the area of Dongting Lake has been affected by sediment deposition within the lake. Although the TGD has acted as an interceptor of sediment, that is, the annual mean sediment inflow into Dongting Lake from the three channels and four tributaries decreased markedly due to the operation of the TGD, decreasing from 6835 × 10 4 t (1999-2002) to 1730 × 10 4 t (2003-2017), which has reduced the rate of sedimentation in Dongting Lake; however, siltation still occurrs and this sedimentation phenomenon in the lake area is not projected to stop [105,106]. Therefore, Dongting Lake will remain in a siltation state for at least a short period of time. The slow increase in the elevation of the lake bottom and the increase in siltation in shallow areas will gradually decrease the surface area of

Four Tributaries (FT) Three Channels (TH) Annual Mean Runoff (FT and TH)
In addition to being impacted by the abovementioned reduction in the water supply, the area of Dongting Lake has been affected by sediment deposition within the lake. Although the TGD has acted as an interceptor of sediment, that is, the annual mean sediment inflow into Dongting Lake from the three channels and four tributaries decreased markedly due to the operation of the TGD, decreasing from 6835 × 10 4 t (1999-2002) to 1730 × 10 4 t (2003-2017), which has reduced the rate of sedimentation in Dongting Lake; however, siltation still occurrs and this sedimentation phenomenon in the lake area is not projected to stop [105,106]. Therefore, Dongting Lake will remain in a siltation state for at least a short period of time. The slow increase in the elevation of the lake bottom and the increase in siltation in shallow areas will gradually decrease the surface area of Dongting Lake.

Implications and Prospects
Based on the above analysis, the variations in the water area of Dongting Lake are driven not only by human factors but also by temperature, precipitation and the operation of the TGD. Monitoring the variations in the water area and the factors that influence Dongting Lake can provide guidance for the prevention and control of flood disasters in the lake. In particular, during seasons with abundant precipitation, monitoring the dynamic change of the surface water can provide early warning for the occurrence of flood disasters. Nutrient concentrations are extremely important for the formation of algal blooms, although local meteorological conditions also influence the outbreak of algal blooms [107,108]. When the nutrient concentrations in the lake are relatively high, the algal bloom frequency will certainly increase as the climate warms. Especially in the period from March to May, when the temperature begins to rise, the algal bloom area will increase. During this period, water quality monitoring in regions with high algal bloom frequencies can provide early warning for the emergence of blooms. Nevertheless, the internal sources of contamination in the lake should be removed to improve the water quality and inhibit algal blooms when the water quality is poor.
Spatial and temporal monitoring and analyses of the water surface and algal bloom areas in Dongting Lake using remote sensing data can reflect the pollution emissions from surrounding enterprises and cities and can therefore provide real-time monitoring information for the government to control the water quality of Dongting Lake. However, due to the influences of factors such as clouds and aerosols and the limitations of existing algorithms, the extraction accuracy of the water area and algal bloom area must be improved. The use of the traditional algorithms optimized by the DL-NN to extract the target not only avoids setting a threshold for each image but also improves the extraction accuracy to some extent. However, there is still a need to combine remote sensing monitoring with real-time in situ monitoring in the future. Ground observation stations should be established in regions with high algal bloom frequencies (the hardest-hit areas) to improve the accuracy of water area and algal bloom area monitoring. Such stations would provide the government with real-time monitoring information regarding the water area and water quality variations in Dongting Lake, allowing for early warnings concerning the occurrence of disasters such as algal blooms and floods. This information will also provide guiding significance for controlling Dongting Lake, which will help ensure the safe supply of drinking water for local residents and the balance of aquatic ecosystems in the lake.

Conclusions
Dongting Lake has experienced increased eutrophication and a reduction in the water surface area, increasing the risks to the environment and freshwater supply sectors. Therefore, to improve the lake water quality, it is of great significance to monitor the water area and the development of algal blooms in Dongting Lake and enact measures to ensure the safe supply of drinking water and balance among the ecosystems therein. The long-term dynamics of the water surface and algal bloom areas in Dongting Lake were systematically analyzed and quantified based on 17 years of MODIS observations. We attribute the success of the work to four factors: (1) obtaining 17 consecutive years of study area with few cloudy observations from the MODIS instruments, (2) optimizing the traditional algorithms to extract the areas of the water surface and algal blooms more accurately and easily, (3) applying multiple methods to improve the extraction accuracy, and (4) collecting meteorological data and various verification data that are now widely available. The use of remote sensing technology to quickly and accurately determine large-area, long-term sequence changes in the water surface area and water quality of a lake is a major technological advancement.
Among the most striking results of this study, the water area of Dongting Lake gradually decreased during the study period coincident with a decrease in the bloom area, of which four different stages were In addition, the changes in the water surface and algal bloom areas in different lake sections varied due to the influences of human activities in the surrounding areas. Eastern Dongting Lake changed the most, followed by southern Dongting Lake and western Dongting Lake. Precipitation and temperature are two important factors that affected the changes in the water and bloom areas. Temperature can promote the growth of algal blooms, while precipitation has a short-term inhibitory effect on bloom outbreaks.
Water quality management in Dongting Lake has become a heated topic. The government is currently taking measures to control the deterioration of the water environment. Based on the remote sensing monitoring results, the water quality of Dongting Lake could be further improved by combining the removal of internal contamination sources from the lake and the interception of external contaminants in regions with high algal bloom frequencies. These measures could be implemented to control the water quality in Dongting Lake and ensure the safe supply of drinking water for residents around Dongting Lake. The methodology used in this study can be easily extended to other inland lakes to assess and monitor long-term temporal and spatial variations. This research topic is important for many environmental pollution and governance issues, which have considerable scientific and socioeconomic value.