Water Occurrence in the Two Largest Lakes in China Based on Long-Term Landsat Images: Spatiotemporal Changes, Ecological Impacts, and Inﬂuencing Factors

: Poyang Lake and Dongting Lake are the two largest freshwater lakes in China. They play a critical ecological role in regional sustainable development, which is closely related to water occurrence (WO). This study aimed to characterize the spatiotemporal changes, ecological impacts, and inﬂuencing factors of WO in the two lakes. First, this study utilized the Google Earth Engine (GEE) platform and all available Landsat images to retrieve the time series of WO in Poyang Lake and Dongting Lake on a two-year scale from 1989 to 2020. Then, the study explored the landscape variation intensity, landscape structure, and their ecological impacts on net primary production (NPP) and total nitrogen (TN) based on the landscape transfer rate (LTR), landscape indices, and geostatistical connectivity analysis. Finally, the study investigated the natural and human factors that inﬂuence the changes in WO. The major conclusions include the following: (1) the mean WO showed a signiﬁcantly decreasing trend in the two lakes from 1989 to 2020. According to the LTR, the landscape variation in Dongting Lake was more intense before 2002, while after 2002, the landscape variation in Poyang Lake was more intense. (2) The landscape structure of WO gradually fragmented and diversiﬁed, along with decreasing hydrological connectivity in the two lakes. (3) The increase in NPP in Poyang Lake and the increase in TN in Dongting Lake were closely correlated with the changes in WO. (4) The upstream runoff showed a signiﬁcantly positive correlation with the mean WO in Poyang Lake and the annual precipitation showed a signiﬁcantly positive correlation with the mean WO in Dongting Lake; the inﬂuences of natural factors on WO changes became weaker. This study reinforces the understanding of long-term changes in the water body in Poyang Lake and Dongting Lake. calculated on a two-year scale; therefore, the study calculated the two-year average TN series for linear regression analysis.


Introduction
Inland lakes play an important role in regional sustainable development for a variety of ecosystem services, such as recreational fisheries, outdoor activity spaces, climate regulation, and water resources [1,2]. Therefore, it is especially important to monitor water body changes in inland lakes [3]. However, under the background of climate change and human activities, some inland lakes are facing a series of worsening problems, such as the shrinkage of lake water bodies [4], decreases in water quality [5], and increases in lacustrine

Study Area
Dongting Lake and Poyang Lake are the largest freshwater lakes in China. They are located in the Dongting Lake Basin and Poyang Lake Basin in the Middle Yangtze River Basin (Figure 1). The water from five major tributaries (Xiu Shui, Gan Jiang, Fu He, Xin Jiang, and Rao He) in Poyang Lake Basin flows into Poyang Lake and then flows into the Yangtze River. The water from four major tributaries (Li Shui, Yuan Jiang, Zi Jiang, and Xiang Jiang) in Dongting Lake Basin flows into Dongting Lake and then flows into the Yangtze River. The seasonal variations in precipitation result in high fluctuation in the extent of water bodies in the two lakes. The maximum water body extent usually occurs from June to September, and the minimum water body extent usually occurs from December to February [17].

Study Data
Landsat satellite images were selected to extract water bodies in the study area. On the one hand, Landsat images can provide long-term observations as far back as the 1980s. On the other hand, Landsat images have a relatively high spatial resolution (30 m). In this study, all available Landsat 4, 5, 7, and 8 surface reflectance images from 1989 to 2020 were adopted based on the GEE platform. The total number of images collected was 1976 and 2856 for Dongting Lake and Poyang Lake, respectively. Moreover, the study also used the WO dataset from the Joint Research Center of European Commission (JRC) Global Surface Water Mapping Layers [10] to compare with the WO extracted by our study.
This study explored the ecological impacts of WO changes based on two case studies. One case study is the impacts of WO changes on net primary productivity (NPP) in

Study Data
Landsat satellite images were selected to extract water bodies in the study area. On the one hand, Landsat images can provide long-term observations as far back as the 1980s. On the other hand, Landsat images have a relatively high spatial resolution (30 m). In this study, all available Landsat 4, 5, 7, and 8 surface reflectance images from 1989 to 2020 were adopted based on the GEE platform. The total number of images collected was 1976 and 2856 for Dongting Lake and Poyang Lake, respectively. Moreover, the study also used the WO dataset from the Joint Research Center of European Commission (JRC) Global Surface Water Mapping Layers [10] to compare with the WO extracted by our study.
This study explored the ecological impacts of WO changes based on two case studies. One case study is the impacts of WO changes on net primary productivity (NPP) in Poyang Lake. NPP indicates how much carbon dioxide remains after accounting for carbon dioxide the plants released during respiration, and is regarded as an important ecosystem service Poyang Lake. NPP indicates how much carbon dioxide remains after accounting for carbon dioxide the plants released during respiration, and is regarded as an important ecosystem service of wetland vegetation. Poyang Lake has more vegetation than Dongting Lake; therefore, this study explores the impacts of WO changes on NPP. Annual NPP data from 2000 to 2015 in Poyang Lake were collected from Ye et al. [23]. Another case study is the impacts of WO changes on total nitrogen (TN) in Dongting Lake. TN is an indicator of the nitrogen content in water. Excessive TN is an important driving factor for lake eutrophication. The TN data from 1997 to 2014 in Dongting Lake were collected from Tian et al. [24]. Moreover, to analyze the natural and human factors that influence water occurrence changes, this study collected precipitation, runoff, and land cover datasets. The annual precipitation from 1989 to 2020 was calculated based on the PERSIANN-CDR daily precipitation product provided in the GEE platform [25]. The annual runoff from 1989 to 2020 was collected from the Changjiang Water Resources Commission of the Ministry of Water Resources. The annual runoff comes from 11 hydrological stations (Figure 1). We calculated the sum of the runoff from 5 hydrological stations around Poyang Lake as the total inflow runoff into Poyang Lake and 4 hydrological stations around Dongting Lake as the total inflow runoff into Dongting Lake. The runoff at Shashi Station and Hankou Station was regarded as the upstream runoff in this study for Dongting Lake and Poyang Lake, respectively. The land cover dataset with a resolution of 30 m in 1990, 2000, 2010, and 2020 was collected from the China Land Cover Dataset (CLCD) [26] to reflect the human factor in the study area.

Research Framework
The research framework of this study is shown in Figure 2. It is divided into three sections. The first section is about the data described in the contents above.  The second section describes how to generate the WO dataset for further applications: (1) perform image preprocessing; (2) extract water bodies from the preprocessed Landsat images based on the water detection rule and conduct the accuracy assessment, and (3) calculate and classify WO on a two-year scale. Then, we obtain the WO dataset in Poyang Lake and Dongting Lake from 1989 to 2020. The last section describes how we will apply the obtained WO dataset, which includes: (1) analysis of the temporal changes in WO from three aspects, including variation intensity, landscape pattern, and hydrological connectivity; (2) exploration of the ecological impacts of WO changes; and (3) investigation of the influencing factors on WO changes. In Section 2.4, we will provide details about the methods we used.

Image Preprocessing
For each image, the study removed the clouds and cloud shadows according to the Landsat image quality band. Then, the preprocessed Landsat image was used to extract water bodies.

Extraction of Water Bodies
The water detection rule is a common method for identifying water bodies from remote sensing images for convenience and high accuracy [13,27]. In this study, we also used the water detection rule proposed by Zou et al. [12] to retrieve the water bodies from 1989 to 2020. The water detection rule combines the modified normalized difference water index (MNDWI), normalized difference vegetation index (NDVI), and enhanced vegetation index (EVI). If a pixel meets the following criteria: EVI < 0.1 and MNDWI > EVI or MNDWI > NDVI, it is regarded as a water body. The water detection rule has been widely used in other studies [12,28] due to its simplicity and high accuracy. Then, the accuracy of the extracted water body was assessed by 2000 random sample points and visual interpretation. Moreover, we compared the results with the JRC Global Surface Water Mapping Layers (V 1.4).

Calculation and Classification of Water Occurrence
According to extracted water bodies, WO was calculated as: where ε represents the corresponding cell value (1 for water and 0 for nonwater) of the ith lake water body map, and N indicates the total valid observations (the clouds and cloud shadows were masked according to the Landsat quality band) for the corresponding cell during a certain period. WO ranges from 0 to 100%. Here, given that the number of images within a year is small, we calculated WO every two years, namely, 1989-1990, 1991-1992 ..., and 2019-2020. In a previous study, breakpoints of 5% and 75% were used to extract ephemeral water bodies (0 < WO < 5%), seasonal water bodies (5% ≤ WO ≤ 75%), and permanent water bodies (WO > 75%) in the contiguous United States [11]. In this study, to make the classification finer, the study adopted breakpoints of 5%, 25%, 50%, 75%, and 95% to divide the study area into six types of landscape (Table 1). Here, 95% was used as the threshold to distinguish seasonal and permanent water bodies. This study suggests that 95% is more suitable for Poyang Lake and Dongting Lake because the two lakes show high seasonal variation in the water bodies. Moreover, seasonal water bodies in this study (5% ≤ WO ≤ 95%) were further divided into four types based on roughly equal intervals. To quantify the variation intensity among landscape types, an indicator called the landscape transfer rate (LTR) was proposed in this study. The LTR is defined as follows: where A t,i indicates the total area of type i at time t, and TA t→t 1 ,i→j stands for the transition area from type i to type j (j = i) during the time phase from the time t to time t 1 . i and j are the IDs of the landscape types (Table 1). LTR t,i→j indicates the variation intensity from type i to type j during the time phase from the time t to time t 1 . According to LTR t,i→j , the study can obtain multitype average LTR at time t by averaging all LTR t,i→j to reflect the landscape variation intensity of the whole study area.

Landscape Pattern Index
In this study, the patch density (PD), largest patch index (LPI), contagion index (CON-TAG), patch cohesion index (COHESION), and Shannon diversity index (SHDI) were selected to reveal the landscape pattern features of Poyang Lake and Dongting Lake. Among them, PD, LPI, CONTAG, SHDI, and COHESION can reflect fragmentation, dominance, aggregation, connectivity, and heterogeneity of landscape patterns, respectively [29]. The calculation of landscape indices was based on Fragstats4.2 software.

Geostatistical Connectivity
This study adopted the geostatistical connectivity analysis method to investigate the hydrological connectivity of Dongting Lake and Poyang Lake. The geostatistical connectivity analysis method has been used in wetland hydrological connectivity studies [30,31].
The calculation of geostatistical connectivity in this study is based on a binary image. Here, the binary image can be obtained based on the WO image and threshold segmentation. If a pixel in the WO image meets the criterion of pixel value (WO) > threshold, it is regarded as a wet pixel. In contrast, if a pixel in the WO image meets the criterion of pixel value (WO) ≤ threshold, it is regarded as a dry pixel. The study set two thresholds to investigate the hydrological connectivity under the different WO situations. One was to set the threshold as 75%, and the other was to set the threshold as 50%. Finally, the study can obtain a wet-dry binary image.
Based on the wet-dry binary image, the connectivity function (CF) was calculated. CF is the core point of geostatistical connectivity analysis. It quantifies the probability that wet cells are connected along a given direction at a given distance [30]. A simple illustrated example is presented in Figure 3. The CF value was calculated as the fraction of wet cell pairs at the given distance that were connected. More details can be found in the paper [31,32]. illustrated example is presented in Figure 3. The CF value was calculated as the fraction of wet cell pairs at the given distance that were connected. More details can be found in the paper [31,32]. This study calculated the CF in the x-axis and y-axis directions at different distances using 8-way connectivity, and then averaged all the CF values. The final average CF value was used to reflect the hydrological connectivity in Poyang Lake and Dongting Lake.

Other Methods
Linear regression and the Pearson correlation coefficient were used to analyze the trend of time series data and investigate the correlation among different time series data. Table 2 shows the accuracy of extracted water bodies according to the water detection rule. The overall accuracy and kappa coefficient are 96.10% and 0.922, respectively, which indicates that the extracted water bodies had high accuracy, and can be used for further analysis. The study also compared the water body results with the existing JRC WO layer. We calculated the WO from 1989 to 2020 in the two lakes and compared the results with the JRC WO ( Figure 4). First, in terms of Poyang Lake, the percentage of WO < 75% according to this study is more than that according to the JRC dataset. Moreover, in terms of This study calculated the CF in the x-axis and y-axis directions at different distances using 8-way connectivity, and then averaged all the CF values. The final average CF value was used to reflect the hydrological connectivity in Poyang Lake and Dongting Lake.

Other Methods
Linear regression and the Pearson correlation coefficient were used to analyze the trend of time series data and investigate the correlation among different time series data. Table 2 shows the accuracy of extracted water bodies according to the water detection rule. The overall accuracy and kappa coefficient are 96.10% and 0.922, respectively, which indicates that the extracted water bodies had high accuracy, and can be used for further analysis. The study also compared the water body results with the existing JRC WO layer. We calculated the WO from 1989 to 2020 in the two lakes and compared the results with the JRC WO ( Figure 4). First, in terms of Poyang Lake, the percentage of WO < 75% according to this study is more than that according to the JRC dataset. Moreover, in terms of Dongting Lake, the percentage of WO < 50% according to this study is greater than that according to the JRC dataset. In general, due to the difference in the water extraction method, WO also shows little difference. Figure 5 shows the spatial pattern of water occurrence in Poyang Lake and Dongting Lake from 1989 to 2020. From the figures, it can be found that higher WO (WO > 75%) in Poyang Lake occurred in rivers and small dish-shaped lakes. The same phenomena can be observed in Dongting Lake. Moreover, some sublakes that are separated from the main lake of Poyang Lake also showed higher WO, and these sublakes are located in eastern and southern Poyang Lake. Moreover. According to Figure 5, Poyang Lake had a higher WO than Dongting Lake in all the periods, except for the period from 1993 to 1994. Moreover, according to Figure 5, the study found that mean WO had a significantly decreasing trend in Poyang Lake and Dongting Lake, and in particular, the decreasing rates in Dongting Lake were faster. It is noted that after 2010, the mean WO showed an increasing trend for Poyang Lake and Dongting Lake. Dongting Lake, the percentage of WO < 50% according to this study is greater than that according to the JRC dataset. In general, due to the difference in the water extraction method, WO also shows little difference. represent Poyang Lake and Dongting Lake, respectively. Here, the water occurrence of this study was calculated for the whole period (1989 to 2020). Figure 5 shows the spatial pattern of water occurrence in Poyang Lake and Dongting Lake from 1989 to 2020. From the figures, it can be found that higher WO (WO > 75%) in Poyang Lake occurred in rivers and small dish-shaped lakes. The same phenomena can be observed in Dongting Lake. Moreover, some sublakes that are separated from the main lake of Poyang Lake also showed higher WO, and these sublakes are located in eastern and southern Poyang Lake. Moreover. According to Figure 5, Poyang Lake had a higher WO than Dongting Lake in all the periods, except for the period from 1993 to 1994. Moreover, according to Figure 5, the study found that mean WO had a significantly decreasing trend in Poyang Lake and Dongting Lake, and in particular, the decreasing rates in Dongting Lake were faster. It is noted that after 2010, the mean WO showed an increasing trend for Poyang Lake and Dongting Lake. . Spatial pattern of and temporal changes in water occurrence in Poyang Lake and Dongting Lake from 1989 to 2020. (a,b) indicate the spatial patterns of Poyang Lake and Dongting Lake, respectively. (c) indicates the temporal changes in the mean WO in Poyang Lake and Dongting Lake from 1989 to 2020. Figure 6 shows the time series of LTR in Poyang Lake and Dongting Lake from 1989 to 2020. Before 2002, the LTR in Dongting Lake was higher than that in Poyang Lake during most periods. After 2002, the LTR in Poyang Lake was higher than that in Dongting Lake during most periods. In general, the average LTR in Poyang Lake and Dongting Lake Figure 5. Spatial pattern of and temporal changes in water occurrence in Poyang Lake and Dongting Lake from 1989 to 2020. (a,b) indicate the spatial patterns of Poyang Lake and Dongting Lake, respectively. (c) indicates the temporal changes in the mean WO in Poyang Lake and Dongting Lake from 1989 to 2020. Figure 6 shows the time series of LTR in Poyang Lake and Dongting Lake from 1989 to 2020. Before 2002, the LTR in Dongting Lake was higher than that in Poyang Lake during most periods. After 2002, the LTR in Poyang Lake was higher than that in Dongting Lake during most periods. In general, the average LTR in Poyang Lake and Dongting Lake were 8.07% and 8.18%, respectively, which indicates that the landscape variation in Dongting Lake was more intense. In addition, the study found that the LTR in Poyang Lake showed a slightly increasing trend. In contrast, the LTR in Dongting Lake showed a significantly decreasing trend. Figure 5. Spatial pattern of and temporal changes in water occurrence in Poyang Lake an Lake from 1989 to 2020. (a,b) indicate the spatial patterns of Poyang Lake and Dongti spectively. (c) indicates the temporal changes in the mean WO in Poyang Lake and Do from 1989 to 2020. Figure 6 shows the time series of LTR in Poyang Lake and Dongting Lake to 2020. Before 2002, the LTR in Dongting Lake was higher than that in Poyang ing most periods. After 2002, the LTR in Poyang Lake was higher than that in Lake during most periods. In general, the average LTR in Poyang Lake and Don were 8.07% and 8.18%, respectively, which indicates that the landscape variatio ting Lake was more intense. In addition, the study found that the LTR in Po showed a slightly increasing trend. In contrast, the LTR in Dongting Lake sho nificantly decreasing trend.   Figure 7 shows the time series of landscape indices in Poyang Lake and Dongting Lake from 1989 to 2020. First, in terms of PD, it had an increasing trend in both Poyang Lake and Dongting Lake. However, only the PD in Dongting Lake showed a significantly increasing trend, with a faster increasing rate. Then, from the view of the LPI, Dongting Lake showed only a slightly decreasing trend, while Poyang Lake showed a significantly decreasing trend. The decreasing rate of LPI in Poyang Lake was approximately 7.15 times faster than that in Dongting Lake. A decreasing trend can also be found for CONTAG and COHESION in the two lakes. Likewise, the decreasing rate in Poyang Lake was faster than that in Dongting Lake. Moreover, the SHDI showed no trend in Dongting Lake; however, it had an increasing trend in Poyang Lake. According to the correlation matrix (Figure 7), CONTAG and SHDI had the highest negative correlation in the two lakes. CONTAG and COHESION had the highest positive correlation in the two lakes.

Spatial and Temporal Changes in Water Occurrence
Hydrological connectivity is especially important for Poyang Lake and Dongting Lake. Geostatistical connectivity analysis was used to investigate the hydrological connectivity, and the results are shown in Figure 8. The CF value showed a significantly decreasing trend over past years in the two lakes. For Poyang Lake, prior to 2006, the decrease rate of the CF value was relatively low. However, from 2006, the decrease rate of the CF increased abruptly, especially when the study regarded the cells with WO > 75% as wet cells, and the CF value decreased to the minimum value (0.142) during the period from 2013 to 2014. Compared with Poyang Lake, Dongting Lake faced a slightly better situation. In Dongting Lake, the CF value also showed a decreasing trend. The decrease rate of the CF was lower than that in Poyang Lake. Moreover, by averaging the CF value during different periods, the study found that, when WO increased from 50% to 75%, the mean CF decreased by 21.29% and 11.12% in Poyang Lake and Dongting Lake, respectively. decreasing trend. The decreasing rate of LPI in Poyang Lake was approximately 7.15 times faster than that in Dongting Lake. A decreasing trend can also be found for CONTAG and COHESION in the two lakes. Likewise, the decreasing rate in Poyang Lake was faster than that in Dongting Lake. Moreover, the SHDI showed no trend in Dongting Lake; however, it had an increasing trend in Poyang Lake. According to the correlation matrix (Figure 7), CONTAG and SHDI had the highest negative correlation in the two lakes. CONTAG and COHESION had the highest positive correlation in the two lakes. Hydrological connectivity is especially important for Poyang Lake and Dongting Lake. Geostatistical connectivity analysis was used to investigate the hydrological connectivity, and the results are shown in Figure 8. The CF value showed a significantly decreasing trend over past years in the two lakes. For Poyang Lake, prior to 2006, the decrease rate of the CF value was relatively low. However, from 2006, the decrease rate of the CF increased abruptly, especially when the study regarded the cells with WO > 75% as wet cells, and the CF value decreased to the minimum value (0.142) during the period from 2013 to 2014. Compared with Poyang Lake, Dongting Lake faced a slightly better situation. In Dongting Lake, the CF value also showed a decreasing trend. The decrease rate of the CF was lower than that in Poyang Lake. Moreover, by averaging the CF value during different periods, the study found that, when WO increased from 50% to 75%, the mean CF decreased by 21.29% and 11.12% in Poyang Lake and Dongting Lake, respectively.

Influences of WO Changes on NPP in Poyang Lake
Seasonal water inundation results in changes in NPP during different periods. In this study, we explore the influences of water occurrence changes on NPP. The results are shown in Figure 9. From 2000 to 2015, NPP in Poyang Lake experienced an increasing trend, with an increasing rate of 10.03 gC/m 2 . By linear regression analysis, it can be found that NPP showed a significantly negative relationship with mean WO (R 2 = 0.763, p < 0.05). Moreover, from the perspective of the water area of different types, NEWB, SWBME, SWBE, and SWBP exhibited a positive correlation with NPP. However, only SWBME showed a significantly positive relationship (R 2 = 0.700, p < 0.05). In contrast, SWBMP and PWB presented a significantly negative correlation with NPP. In general, NPP is closely

Influences of WO Changes on NPP in Poyang Lake
Seasonal water inundation results in changes in NPP during different periods. In this study, we explore the influences of water occurrence changes on NPP. The results are shown in Figure 9. From 2000 to 2015, NPP in Poyang Lake experienced an increasing trend, with an increasing rate of 10.03 gC/m 2 . By linear regression analysis, it can be found that NPP showed a significantly negative relationship with mean WO (R 2 = 0.763, p < 0.05).
Moreover, from the perspective of the water area of different types, NEWB, SWBME, SWBE, and SWBP exhibited a positive correlation with NPP. However, only SWBME showed a significantly positive relationship (R 2 = 0.700, p < 0.05). In contrast, SWBMP and PWB presented a significantly negative correlation with NPP. In general, NPP is closely related to water occurrence changes.

Influences of WO Changes on TN in Dongting Lake
Here, this study also explores the influences of water occurrence changes on TN. The results are shown in Figure 10. From 1997 to 2014, TN in Dongting Lake remained relatively stable. However, from 2008 to 2010, TN experienced an abrupt increase from 1.22 mg/L to 1.93 mg/L. From 2010 to 2014, TN also had a slightly increasing trend. Linear regression indicates that the mean WO showed a negative relationship with TN. However, the negative relationship was not significant (p > 0.05). From the perspective of water areas of different types, only NEWB showed a significantly positive relationship with TN. Although the R 2 between SWBE and TN was relatively higher (0.4120), the negative relationship was still not significant (p > 0.05). Moreover, the study found that NP and SHDI showed a positive relationship with TN, while LPI, CONTAG, and COHESION showed a negative relationship with TN. However, none of the landscape indices showed a significant relationship with TN. It is worth noting that NP had a higher R 2 (R 2 = 0.2511) than the other four landscape indices. Moreover, the CF value also showed a negative relationship with TN.

Influences of WO Changes on TN in Dongting Lake
Here, this study also explores the influences of water occurrence changes on TN. The results are shown in Figure 10. From 1997 to 2014, TN in Dongting Lake remained relatively stable. However, from 2008 to 2010, TN experienced an abrupt increase from 1.22 mg/L to 1.93 mg/L. From 2010 to 2014, TN also had a slightly increasing trend. Linear regression indicates that the mean WO showed a negative relationship with TN. However, the negative relationship was not significant (p > 0.05). From the perspective of water areas of different types, only NEWB showed a significantly positive relationship with TN. Although the R 2 between SWBE and TN was relatively higher (0.4120), the negative relationship was still not significant (p > 0.05). Moreover, the study found that NP and SHDI showed a positive relationship with TN, while LPI, CONTAG, and COHESION showed a negative relationship with TN. However, none of the landscape indices showed a significant relationship with TN. It is worth noting that NP had a higher R 2 (R 2 = 0.2511) than the other four landscape indices. Moreover, the CF value also showed a negative relationship with TN.

Influences of Precipitation and Runoff on WO Changes
Precipitation plays an important role in controlling the seasonal variations in lake water extents in Poyang Lake and Dongting Lake [19,33]. Usually, more precipitation would result in the expansion of lake water, which increases the mean water occurrence of the lake. Therefore, the study also observed that annual precipitation had a positive correlation with mean WO, although the correlation was not significant (p > 0.05) for Poyang Lake (Figure 11b,c). Poyang Lake and Dongting Lake are facing the fact that mean WO has decreased over the past 39 years. However, the study only found that annual precipitation showed a decreasing trend in Dongting Lake, which may lead to a decrease in mean WO in Dongting Lake. In Poyang Lake, the annual precipitation even showed a slightly increasing trend. However, the mean WO in Poyang Lake has decreased in recent years, which suggests that other factors influence the WO changes in Poyang Lake.

Influences of Precipitation and Runoff on WO Changes
Precipitation plays an important role in controlling the seasonal variations in lake water extents in Poyang Lake and Dongting Lake [19,33]. Usually, more precipitation would result in the expansion of lake water, which increases the mean water occurrence of the lake. Therefore, the study also observed that annual precipitation had a positive correlation with mean WO, although the correlation was not significant (p > 0.05) for Poyang Lake (Figure 11b,c). Poyang Lake and Dongting Lake are facing the fact that mean WO has decreased over the past 39 years. However, the study only found that annual precipitation showed a decreasing trend in Dongting Lake, which may lead to a decrease in mean WO in Dongting Lake. In Poyang Lake, the annual precipitation even showed a slightly increasing trend. However, the mean WO in Poyang Lake has decreased in recent years, which suggests that other factors influence the WO changes in Poyang Lake.
Another natural factor that influences the WO changes is the runoff into the lake. Theoretically, more inflow runoff will enlarge the water extent, thus increasing the mean WO. The inflow runoff changes in Poyang Lake and Dongting Lake are shown in Figure 12. The inflow runoff in Poyang Lake and Dongting Lake showed a decreasing trend, which possibly resulted in a decrease in mean WO in the study area to a certain degree. It is noted that the decreasing trend was not statistically significant (p > 0.05) for these two lakes. Moreover, as expected, inflow runoff showed a positive correlation with mean WO.
In particular, inflow runoff showed a stronger correlation with mean WO in Dongting Lake (R 2 = 0.201) than in Poyang Lake (R 2 = 0.122), although the correlation was not statistically significant (p > 0.05). Another natural factor that influences the WO changes is the runoff into the lake. Theoretically, more inflow runoff will enlarge the water extent, thus increasing the mean WO. The inflow runoff changes in Poyang Lake and Dongting Lake are shown in Figure  12. The inflow runoff in Poyang Lake and Dongting Lake showed a decreasing trend, which possibly resulted in a decrease in mean WO in the study area to a certain degree. It is noted that the decreasing trend was not statistically significant (p > 0.05) for these two lakes. Moreover, as expected, inflow runoff showed a positive correlation with mean WO. In particular, inflow runoff showed a stronger correlation with mean WO in Dongting Lake (R 2 = 0.201) than in Poyang Lake (R 2 = 0.122), although the correlation was not statistically significant (p > 0.05). In addition, as a river-connected lake, the water occurrence changes are naturally affected by the upstream runoff from the Yangtze River. As shown in Figure 13, the upstream runoff of Poyang Lake showed a decreasing trend. At the same time, upstream runoff showed a significantly positive correlation with mean WO in Poyang Lake (R 2 =  Another natural factor that influences the WO changes is the runoff into the lake. Theoretically, more inflow runoff will enlarge the water extent, thus increasing the mean WO. The inflow runoff changes in Poyang Lake and Dongting Lake are shown in Figure  12. The inflow runoff in Poyang Lake and Dongting Lake showed a decreasing trend, which possibly resulted in a decrease in mean WO in the study area to a certain degree. It is noted that the decreasing trend was not statistically significant (p > 0.05) for these two lakes. Moreover, as expected, inflow runoff showed a positive correlation with mean WO. In particular, inflow runoff showed a stronger correlation with mean WO in Dongting Lake (R 2 = 0.201) than in Poyang Lake (R 2 = 0.122), although the correlation was not statistically significant (p > 0.05). In addition, as a river-connected lake, the water occurrence changes are naturally affected by the upstream runoff from the Yangtze River. As shown in Figure 13, the upstream runoff of Poyang Lake showed a decreasing trend. At the same time, upstream runoff showed a significantly positive correlation with mean WO in Poyang Lake (R 2 = 0.331, p < 0.05). Therefore, the upstream runoff decrease plays a critical role in the mean In addition, as a river-connected lake, the water occurrence changes are naturally affected by the upstream runoff from the Yangtze River. As shown in Figure 13, the upstream runoff of Poyang Lake showed a decreasing trend. At the same time, upstream runoff showed a significantly positive correlation with mean WO in Poyang Lake (R 2 = 0.331, p < 0.05). Therefore, the upstream runoff decrease plays a critical role in the mean WO of Poyang Lake. In Dongting Lake, upstream runoff almost retained a stable trend. The positive correlation between upstream runoff and mean WO was also quite weak (R 2 = 0.087, p > 0.05). WO of Poyang Lake. In Dongting Lake, upstream runoff almost retained a stable trend. The positive correlation between upstream runoff and mean WO was also quite weak (R 2 = 0.087, p > 0.05). Moreover, we combined the moving window approach and Pearson correlation coefficient to investigate the correlation between the mean WO and these three natural factors at different periods. The window size is 10 periods (20 years), and the sliding time step is 1 period. For example, according to the mean WO and precipitation from 1989-1990 to 2007-2008, the study calculated the Pearson correlation coefficient between the mean WO and precipitation. Then, according to mean WO and precipitation from the period of 1991-1992 to the period of 2009-2010, the study could obtain another Pearson correlation coefficient between the mean WO and precipitation. Finally, the study could obtain a Pearson correlation coefficient record and investigate the correlation between mean WO and precipitation at different periods. Figure 14 shows the results of a multiscale correlation between mean water occurrence and precipitation, inflow runoff, and upstream runoff in Poyang Lake and Dongting Lake. In Poyang Lake, all the correlation coefficients for mean WO and precipitation, inflow runoff, and upstream runoff showed a decreasing trend. In particular, the decrease in the correlation coefficient for the mean WO and inflow was statistically significant (p < 0.05) and was the fastest. For Dongting Lake, the correlation coefficient for mean WO and precipitation and upstream runoff showed similar decreasing change characteristics with Poyang Lake. However, the correlation coefficient for mean WO and inflow runoff showed a slightly increasing trend, which was different from Poyang Lake. In general, the decreases in the correlation coefficient suggest that the influences of natural factors on WO changes became weaker. Moreover, we combined the moving window approach and Pearson correlation coefficient to investigate the correlation between the mean WO and these three natural factors at different periods. The window size is 10 periods (20 years), and the sliding time step is 1 period. For example, according to the mean WO and precipitation from 1989-1990 to 2007-2008, the study calculated the Pearson correlation coefficient between the mean WO and precipitation. Then, according to mean WO and precipitation from the period of 1991-1992 to the period of 2009-2010, the study could obtain another Pearson correlation coefficient between the mean WO and precipitation. Finally, the study could obtain a Pearson correlation coefficient record and investigate the correlation between mean WO and precipitation at different periods. Figure 14 shows the results of a multiscale correlation between mean water occurrence and precipitation, inflow runoff, and upstream runoff in Poyang Lake and Dongting Lake. In Poyang Lake, all the correlation coefficients for mean WO and precipitation, inflow runoff, and upstream runoff showed a decreasing trend. In particular, the decrease in the correlation coefficient for the mean WO and inflow was statistically significant (p < 0.05) and was the fastest. For Dongting Lake, the correlation coefficient for mean WO and precipitation and upstream runoff showed similar decreasing change characteristics with Poyang Lake. However, the correlation coefficient for mean WO and inflow runoff showed a slightly increasing trend, which was different from Poyang Lake. In general, the decreases in the correlation coefficient suggest that the influences of natural factors on WO changes became weaker. WO of Poyang Lake. In Dongting Lake, upstream runoff almost retained a stable trend. The positive correlation between upstream runoff and mean WO was also quite weak (R 2 = 0.087, p > 0.05). Moreover, we combined the moving window approach and Pearson correlation coefficient to investigate the correlation between the mean WO and these three natural factors at different periods. The window size is 10 periods (20 years), and the sliding time step is 1 period. For example, according to the mean WO and precipitation from 1989-1990 to 2007-2008, the study calculated the Pearson correlation coefficient between the mean WO and precipitation. Then, according to mean WO and precipitation from the period of 1991-1992 to the period of 2009-2010, the study could obtain another Pearson correlation coefficient between the mean WO and precipitation. Finally, the study could obtain a Pearson correlation coefficient record and investigate the correlation between mean WO and precipitation at different periods. Figure 14 shows the results of a multiscale correlation between mean water occurrence and precipitation, inflow runoff, and upstream runoff in Poyang Lake and Dongting Lake. In Poyang Lake, all the correlation coefficients for mean WO and precipitation, inflow runoff, and upstream runoff showed a decreasing trend. In particular, the decrease in the correlation coefficient for the mean WO and inflow was statistically significant (p < 0.05) and was the fastest. For Dongting Lake, the correlation coefficient for mean WO and precipitation and upstream runoff showed similar decreasing change characteristics with Poyang Lake. However, the correlation coefficient for mean WO and inflow runoff showed a slightly increasing trend, which was different from Poyang Lake. In general, the decreases in the correlation coefficient suggest that the influences of natural factors on WO changes became weaker.

Influences of the Encroachment of Cropland on WO Changes
An important human factor is the encroachment of cropland. As shown in Figure 15, from 2000 to 2010, more than 500 km 2 of water bodies were converted into cropland in Poyang Lake and Dongting Lake. Therefore, the increase in cropland area would cause a decrease in the WO of the lake. From Figure 4c, the study also observed that the mean WO started to decrease from 2000. Due to the increasing attention on environmental protection in Poyang Lake and Dongting Lake after 2010, part of the cropland was restored to lake water bodies, which would contribute to the increase in mean WO in the two lakes, which is consistent with the fact that the study also observed the phenomena in Figure 4c.
Remote Sens. 2022, 14, x FOR PEER REVIEW 16 of 20 Figure 14. Multiscale correlation between mean water occurrence and precipitation, inflow runoff, and upstream runoff in Poyang Lake and Dongting Lake. (a-c) indicate precipitation, inflow runoff, and upstream runoff, respectively.

Influences of the Encroachment of Cropland on WO Changes
An important human factor is the encroachment of cropland. As shown in Figure 15, from 2000 to 2010, more than 500 km 2 of water bodies were converted into cropland in Poyang Lake and Dongting Lake. Therefore, the increase in cropland area would cause a decrease in the WO of the lake. From Figure 4c, the study also observed that the mean WO started to decrease from 2000. Due to the increasing attention on environmental protection in Poyang Lake and Dongting Lake after 2010, part of the cropland was restored to lake water bodies, which would contribute to the increase in mean WO in the two lakes, which is consistent with the fact that the study also observed the phenomena in Figure 4c. Figure 15. The area transition between cropland and water body in Poyang Lake and Dongting Lake from 1990 to 2020. (a,b) indicate Poyang Lake and Dongting Lake, respectively.

Similarities and Differences in Water Occurrence Changes and the Influencing Factors
This study showed that there were many similarities in the spatial and temporal changes in the WO of Poyang Lake and Dongting Lake. For example, both experienced a decreasing trend in the mean WO ( Figure 4) and hydrological connectivity (Figure 8). Both Poyang Lake and Dongting Lake were influenced by the decreasing inflow runoff and the increasing encroachment of cropland. These similarities in these external factors caused the changes in WO of the two lakes to show similarities over the past years. The differences between the changes in WO in the two lakes can be attributed to some intrinsic factors, such as lake morphology and topographical characteristics [17]. Of course, external factors also cannot be ignored. In fact, given that the influences of three natural influencing factors (precipitation, inflow runoff, and upstream runoff) on WO changes became weaker (Figure 14), there is a reasonable inference that human activities are playing an increasingly important role in WO changes. Our results indicate that the encroachment of cropland results in a decrease in water bodies. Consequently, the mean WO also decrease. Except for the encroachment of cropland, there are a wide variety of human activities that influence the WO. For example, the construction of water conservancy projects. The Three Gorges Dam (TGD), located in the upper reaches of the Yangtze River, is the world's largest hydroelectric plant. In 2003, it began storing water. Previous studies have demonstrated the profound influences of the TGD on ecohydrology in the middle and lower reaches of the Yangtze River [34][35][36]. For example, the TGD changed the seasonal inundation pattern of Poyang Lake [35].

Ecological Impact of Water Occurrence Changes
NPP is an important ecosystem service provided by Poyang Lake. Ye et al. [22] observed that annual NPP shows a significant negative correlation with water inundation duration in Poyang Lake, and the correlation intensity is far greater than that of the other three influencing factors (air temperature, solar radiation, and precipitation) that control Figure 15. The area transition between cropland and water body in Poyang Lake and Dongting Lake from 1990 to 2020. (a,b) indicate Poyang Lake and Dongting Lake, respectively.

Similarities and Differences in Water Occurrence Changes and the Influencing Factors
This study showed that there were many similarities in the spatial and temporal changes in the WO of Poyang Lake and Dongting Lake. For example, both experienced a decreasing trend in the mean WO ( Figure 4) and hydrological connectivity (Figure 8). Both Poyang Lake and Dongting Lake were influenced by the decreasing inflow runoff and the increasing encroachment of cropland. These similarities in these external factors caused the changes in WO of the two lakes to show similarities over the past years. The differences between the changes in WO in the two lakes can be attributed to some intrinsic factors, such as lake morphology and topographical characteristics [17]. Of course, external factors also cannot be ignored. In fact, given that the influences of three natural influencing factors (precipitation, inflow runoff, and upstream runoff) on WO changes became weaker ( Figure 14), there is a reasonable inference that human activities are playing an increasingly important role in WO changes. Our results indicate that the encroachment of cropland results in a decrease in water bodies. Consequently, the mean WO also decrease. Except for the encroachment of cropland, there are a wide variety of human activities that influence the WO. For example, the construction of water conservancy projects. The Three Gorges Dam (TGD), located in the upper reaches of the Yangtze River, is the world's largest hydroelectric plant. In 2003, it began storing water. Previous studies have demonstrated the profound influences of the TGD on ecohydrology in the middle and lower reaches of the Yangtze River [34][35][36]. For example, the TGD changed the seasonal inundation pattern of Poyang Lake [35].

Ecological Impact of Water Occurrence Changes
NPP is an important ecosystem service provided by Poyang Lake. Ye et al. [22] observed that annual NPP shows a significant negative correlation with water inundation duration in Poyang Lake, and the correlation intensity is far greater than that of the other three influencing factors (air temperature, solar radiation, and precipitation) that control NPP. In this study, we observed that there was a significantly negative correlation between the mean WO and NPP in Poyang Lake. NPP stored by vegetation in Poyang Lake is influenced by hydrological conditions, especially by the water level. When the water level increases, vegetation is flooded by water, and the mean WO increases. In contrast, NPP stored by vegetation decreases. When the water level declines, vegetation appears again, and the mean WO decreases. In contrast, NPP stored by vegetation increases. Therefore, there is no doubt that mean WO showed a significantly negative correlation with NPP in Poyang Lake.
A previous study revealed that nutrient input and hydrological conditions mutually affect water quality in Dongting Lake [36]. Both the components and landscape structure of WO can reflect the hydrological conditions in Dongting Lake, thus impacting the changes in TN. However, most of the components and landscape structure indices of WO showed no significant correlation with TN. A possible reason is that the abrupt increase in TN occurred in 2008 when an abrupt change point of TN in Dongting Lake was detected according to the Mann-Kendall test [37]. From 1997 to 2008, precipitation in Dongting Lake showed a remarkable decreasing trend. Notably, from 2005 to 2006, precipitation reached the historical minimum value. Moreover, the Dongting Lake basin experienced rapid economic development from 2000, resulting in discharging more wastewater. For example, from 2000 to 2010, many lake water bodies became cropland (Figure 15), which would lead to an increase in nitrogen fertilizer [23]. Therefore, comprehensive factors resulted in gradual nitrogen accumulation, and an abrupt increase occurred in 2008, which reduced the correlation between TN and WO. In addition, some studies have demonstrated that high hydrological connectivity is expected to improve the water quality of seasonal lakes [18,38]. Our study found that the CF value (WO > 75%) showed a higher correlation with TN than the CF value (WO > 50%) (Figure 10), suggesting the importance of the hydrological connectivity of the core region of lakes. The regions with WO > 75% represent the core lake region of Dongting Lake, which mainly includes the river channel ( Figure 5). Therefore, the connectivity of the river channel is especially important for the lake.

Uncertainty and Prospects
In this study, only Landsat images were used to extract lake water bodies and calculate the WO. Theoretically, the temporal resolution of Landsat images is 16 days. However, due to the cloud contamination, the total number of clear observations based on Landsat images is spatially and temporally uneven. Therefore, the WO calculated by Landsat images may be uncertain because it cannot accurately capture the dynamic change in lake water.
Spatiotemporal fusion for remote sensing data can generate new image series with a finer spatiotemporal resolution by fusing (temporally dense) coarse resolution images and (temporally sparse) fine resolution images [39]. A previous study used an 8-day series of Landsat resolution (30 m) by blending Landsat and MODIS images to investigate the largescale floodplain inundation dynamics, and the results demonstrated that spatiotemporal fusion showed a high potential [40]. Currently, there are a variety of remote sensing data over the world, such as Sentinel-2, GF-1, and SPOT. Spatiotemporal fusion of these images to generate new image series with finer spatiotemporal resolution will be of great importance for studying the dynamic changes in lake water in support of SDG target 6.6.

Conclusions
In this study, we used all available Landsat images from 1989 to 2020 to extract water bodies and calculate WO in Poyang Lake and Dongting Lake based on the GEE cloud platform. Then, the study area was reclassified into six types (NEWB, SWBME, SWBE, SWBP, SWBMP, and PWB) according to WO breakpoints (5%, 25%, 50%, 75%, and 95%, respectively) for further analysis. The study characterized the spatial and temporal variations in WO from two aspects including variation intensity and landscape structure, and investigated the influences of WO on NPP in Poyang Lake and the influences of WO on TN in Dongting Lake. Moreover, the study analyzed the natural and human factors that influence WO changes. The main conclusions are as follows: