Retrieval of Secchi Disk Depth in Turbid Lakes from GOCI Based on a New Semi-analytical Algorithm

: The accurate remote estimation of the Secchi disk depth (Z SD ) in turbid waters is essential in the monitoring the ecological environment of lakes. Using the field measured Z SD and the remote sensing reflectance (Rrs(λ)) data, a new semi-analytical algorithm (denoted as Z SDZ ) for retrieving Z SD was developed from Rrs(λ), and it was applied to Geostationary Ocean Color Imager (GOCI) images in extremely turbid waters. Our results are as follows: (1) the Z SDZ performs well in estimating Z SD in turbid water bodies (0.15 m < Z SD < 2.5 m). By validating with the field measured data that were collected in four turbid inland lakes, the determination coefficient (R 2 ) is determined to be 0.89, with a mean absolute square percentage error (MAPE) of 22.39%, and root mean square error (RMSE) of 0.24 m. (2) The Z SDZ improved the retrieval accuracy of Z SD in turbid waters and outperformed the existing semi-analytical schemes. (3) The developed algorithm and GOCI data are in order to map the hourly variation of Z SD in turbid inland waters, the GOCI-derived results reveal a significant spatiotemporal variation in our study region, which are significantly driven by wind forcing. This study can provide a new approach for estimating water transparency in turbid waters, offering important support for the management of inland waters.


Introduction
Water transparency, which has been widely and regularly measured with Secchi disk in both limnology and oceanography, is key information for assessing the water quality and trophic state. Secchi disk depth (ZSD, m) is the depth at which the disk, lowered from the surface, is no longer visible to the observer, and is a measure of transparency or vertical visibility in the water body. The values of ZSD depend on the intensity of light, water molecules, and optical components (e.g., chlorophyll-a and suspended matter), and they play a key role in understanding aquatic environment variations and biogeochemical processes [1]. Traditional methods of estimating ZSD largely depend on cruise surveys, which are laborious and time-consuming. The large coverage characteristics of satellite data makes it a perfect tool for describing the spatial and temporal variability of ZSD [2].
Substantial efforts, which are based on remote sensing methods, have been made to obtain ZSD for describing the spatiotemporal variability of water properties in marine, coastal, and inland waters [3][4][5][6][7]. They can be mainly divided into two strategies: empirical and semi-analytical approaches. (1) empirical methods: Binding et al. [8] used remote-sensing reflectance at 550 nm to estimate ZSD in Great Lakes based on multi-source satellite data. Shi et al. [5] also estimated ZSD in Lake Taihu by constructing a linear regression model between ZSD and MODIS-Aqua reflectance at 645 nm. However, the empirical methods have the characteristic of location-dependency and may not be portable to other waters. (2) semi-analytical methods: Preisendorfer et al. proposed a semi-analytical approach for estimating inherent optical properties (IOPs) [9] based on classical underwater visibility theory; this approach has been applied to estimate ZSD [7]; however, the results from this approach did not agree well with field observations and human experiences, and a new semi-analytical model (denoted as ZSDV6) based on the radiative transfer theory was proposed by Lee et al. [10] to retrieve ZSD. Shang et al. subsequently verified the encouraging performance of the new method [11].
This new model was developed and validated with a large data covering different IOPs conditions, including marine, coastal, and inland waters (ZSD range of 0.1-30 m), and received excellent results (~18% average absolute difference, R 2 = 0.96) [12]. By this method, the inherent optical properties (IOPs) were first estimated while using the quasi-analytical algorithm (QAA), and then positing the reference wavelength at 550 or 670 nm further derived the ZSD. However, several obvious limitations of the original QAA become evident for applications in inland turbid waters [6,[13][14][15][16]. Firstly, the estimation method of the total absorption coefficient at the reference wavelength (550 or 670 nm) does not work in inland turbid waters with lower absorption, but high scattering coefficients; secondly, the original coefficients of the IOPs were calibrated while using the data from marine and coastal waters and inland waters. The ZSD in turbid water is much lower than in oceanic water or coastal waters, and so the ZSD model must be reevaluated for these turbid conditions. Moreover, vertical heterogeneous distribution of water constituents and complex optical properties in turbid lakes challenge the validation of parameters in the semi-analytical model of ZSD that was proposed by Lee et al. [10]. Therefore, an improved semi-analytical algorithm is needed to retrieve ZSD in turbid waters.
The Geostationary Ocean Color Imager (GOCI), an ocean color satellite imager that has been set in geostationary orbit, provides eight images during the daytime from 8:30 to 16:30 local time. With its eight bands from visible to near-infrared spectral region (412 ~ 865 nm) and one hourly short temporal resolutions, GOCI has been confirmed as an ideal satellite sensor with high signal-to-noise ratios (SNR) for mapping suspended particulate matter, colored dissolved organic matter (CDOM), and harmful algal blooms [13,[17][18][19]. Although the application of GOCI data for ZSD estimation has seldom been reported in inland turbid waters, it could be of great significance for the observation of ZSD in inland turbid lakes.
The ZSD of Lake Taihu and Lake Hongze were estimated based on in-situ data and GOCI images to improve the algorithm of retrieving ZSD in inland turbid water. The main objectives of this study were: (1) to develop a new semi-analytical algorithm of ZSD for inland turbid waters, (2) compare the performance of the new scheme with the existing semi-analytical algorithm, and (3) to obtain the spatiotemporal dynamic characteristics of ZSD from GOCI observations.

Study Area
Four turbid waters-Lake Dongting (28°30′-30°20′N, 111°40′-113°10′E), Lake Taihu (30°55′-31°32′N, 118°52′-120°36′E), Lake Hongze (33°06′-33°40′N, 118°10′-118°52′E), and Lake Erhai (25°61′-25°98′N, 100°08′-100°31′E)-are observed in our study. As the second, third, and fourth largest freshwater lakes in China, Lake Dongting, Lake Taihu, and Lake Hongze have water areas of ~2794 km 2 , ~2338 km 2 , and ~1577 km 2 , respectively ( Figure 1). Lake Erhai is the second largest freshwater lake in Yunnan Province. Lake Dongting has a mean water depth of 1.9m and a maximum depth of ~5.5m. Lake Taihu has a mean water depth of 1.9m and a maximum depth of ~2.6m. Lake Hongze has a mean water depth of 1.8 m and a maximum depth of ~4.5m. Lake Erhai has a mean water depth of 10m and maximum depth of ~ 20m. Over the past few decades, because of environmental overload from catchments, the aquatic ecosystems of lakes have faced water quality deterioration and reduced biodiversity. Blue-green algae outbreaks in these lakes seriously affect aquatic biodiversity, tourism, and the health of nearby residents. Lake Taihu is characterized by hypereutrophic and high turbidity due to high levels of primary production and high concentrations of suspended sediments; Lake Hongze is characterized by mesotrophic and high turbidity due to a large amount of sediment and terrestrial debris being carried from the Huai River and human activities [20]. Similar to the above two lakes, Lake Dongting has also experienced eutrophication in recent years due to anthropogenic activities and environmental change [21]. Figure 1. The location of Lake Hongze (B), Lake Dongting (C), Lake Taihu (D), and Lake Erhai (E) in China (A). The symbols represent the sampling date in year-month forma.

In-Situ Water Quality Data and Spectra Data Collection
From 2014 to 2018, the remote-sensing reflectance ( (λ)) and the underwater downward irradiance (E (λ, z)) were measured in-situ by taking 232 samples (87 samples in Lake Hongze, 58 in Lake Taihu, 40 in Lake Dongting, and 47 in Lake Erhai) in ten cruise campaigns ( Figure 1). The (λ) data were field-measured while using an Analytical Spectral Devices (ASD) Inc. FieldSpec Pro. Full details of the measurement process of (λ) are provided by previous studies [22,23]. At each sampling site, the total radiance ( ), the sky-viewing radiance ( ), and the radiance reflected by a standard gray panel ( ) were collected. The (λ) can be computed with the following equation: is the reflectance of the grey panel. The diffuse attenuation coefficient ( ) was field collected while using an underwater spectroradiometer (TriOS Mess-und Datentechnik GmbH, Rastede, Germany). With a spectral channels range of 320 to 950 nm, the spectral sampling interval of spectroradiometer is 3.3 nm, which was subsequently interpolated to a new resolution of 1 nm in indoor data processing. The spectra of underwater downward irradiance ( , ) at different sampling depths ( = 0.2, 0.4,…….,1.2 m) were measured by the spectroradiometer [21].
Secchi disk depth (ZSD) measurements and water samples were simultaneously collected at each sampling site. The Secchi disk is lowered into a body of water by an observer until it disappears from view; the depth at which it is no longer visible is recorded as the ZSD value. One skilled person carried out all of the measurements in order to decrease artificial error. At the same time, in situ water samples were collected using Niskin bottles at a depth of 50 cm under the surface and were kept in containers at −20 °C for laboratory analysis. Whthin the next five days, the water samples were filtered and measured in the laboratory. The measured parameters of water samples in laboratory analysis included chlorophyll-a (Chla), total suspended matter (TSM), and organic/ inorganic suspended matter (OSM/ISM) concentrations, which were determined using the method that was described by Chen et al. [25]. Note that Lake Erhai has not field measured ( ) data used in this study. Finally, the total datasets of 232 samples were divided into two groups, one of which consisted of 137 samples, approximately 58% of the total number of samples, was used to calibrate the model; the remaining group was used to validate the accuracy of the model.

Satellite Data and Preprocessing
GOCI images, which were captured in Lake Hongze and Lake Taihu from 2014 to 2018, were downloaded from the Korea Ocean Satellite Center (KOSC). The 6S method (Second Simulation of the Satellite Signal in the Solar Spectrum, 6S) was used to correct the GOCI Level-1B data, which are the same as the latest research [19,22]. After image data pre-processing, these data were used to validate the performance of 6s model and estimate the ZSD in Lake Hongze and Lake Taihu. During the cruise campaigns that were conducted from 2014 to 2018, the measurement time of a total of 125 points (69 points in Lake Hongze and 56 in Lake Taihu) was close to GOCI imaging time (± 1 h). These sampling points were selected for further validation of the accuracy of 6S atmospheric correction.

Wind Speed Data Collection
Hourly wind data from 2015-2018 in Dongshan meteorological station (31°4′N, 120°26′E) near Lake Taiyhu and Xuyi meteorological station (33°42′N, 118°32′E) near Lake Hongze were obtained from the China Meteorological Data Sharing Service System (http://cdc.c,a.gov.cn/home.do). These data were recorded per 1 h, which could be matched with the GOCI imaging time.

Data analyses and Accuracy Assessment
Correlation analyses were performed with OriginPro 2018 software. Three evaluation indicators, including the coefficient of determination (R 2 ), the mean absolute square percentage error (MAPE), and root mean square error (RMSE), were selected to characterize the performance of the model: where the and are the estimated and measured values, respectively; n is the number of samples.

ZSDZ Algorithm
A new semi-analytical algorithm (denoted as ZSDZ) that is suitable for estimating ZSD in inland turbid waters was developed based on the original ZSDV6 by Lee et al. [10] and modified based on the in-situ data. Table 1 lists the derivation flowchart of ZSDZ algorithm and the corresponding expressions. Two parts (part I and part II) are contained in the new algorithm. Table 1. Steps of ZSDZ algorithm (a new semi-analytical algorithm).
Step Property Expression Approach Semi-analytical 2.6.1. Part I In step 1, ( ) can be calculated as [26]: In step 2, (λ) can be calculad as [26]: where g0 =0.084 and g1 = 0.17 in ZSDZ algorithm, which were suggested for high scattering waters by [15]. In step 3, (λ ) can be obtained based on the assumption that the (λ ) is dominant at λ .
In original ZSDV6 algorithm that was proposed by Lee et al. [10], (λ ) is the absorption of dissolved components and detritus at original reference wavelength (550 nm or 670 nm). However, poor accuracy can be obtained by deriving △ (550) or △ (670), and further leading to poorer performance of (λ ) in turbid waterbodies. Moreover, an ambiguous extrapolation index of the aerosol model usually leads to the poor estimation of from satellite images at short wavelengths in highly turbid inland waters [27]. These defects will propagate to the next steps and affect the accuracy of the model. One way to overcome this shortcoming is to first shift the reference wavelength to longer bands [6,28]. There are two longer NIR (745 nm and 845 nm) that meet the above conditions given the band setting of the GOCI image. However, poor performance of atmospheric correction existed at 865nm. Therefore, in the modified algorithm of ZSDZ, the reference wavelength of (λ ) will be shifted to 745 nm, and the (λ ) is considered to be equal to the absorption coefficient of pure water at 745 nm, (745) . Samplings of 745 nm of GOCI are determined in this step, other sensors with longer NIR band, such as the 740 nm of Sentinel-2, can also be selected as reference wavelength when applying the ZSDZ to other study areas.
In step 4, (λ ) can be calculated as: In the step 5, (λ ) can be calculated as: In ZSDV6 algorithm, the original values of parameters (m0~3) were determined while using Hydrolight simulations on oligotrophic waters and Case-1 model of Morel et al. [29] for optical properties, where the input inherent optical properties (IOPs) were kept vertically constant [30]. Therefore, parameters (m0~3) should be retuned for turbid waters. Parameters (m0~3) were retuned using our calibration data (137 samples) in order to meet the application of turbid waters, similar to Lee et al. [30]. In detail, we fixed the same values for m1, m2, and m3, but varied m0, and the optimal value of m0 was further confirmed based on nonlinear best fit. The values of m1~m3 were obtained in the same way. The final values of the four model parameters (m0~3) were 0.0124, 3.16, 0.52, and 10.8, respectively. After retuning, Figure 2a presents the match up points between the field measured and the retrieved (745). All of the points were closely distributed along the 1:1 line with low MAPE (10.24%) and RMSE (0.81 m -1 ), which suggests that the retuned values of parameters (m0~3) are acceptable for estimation of (745). In step 6, ( ) can be calculated based on line regression between ( ) and ( ), as: ( ) = * ( ) + (10) where b and c are the coefficient and intercept of the line regression equation, respectively, which can be determined by the field measured data.

Part II
In part II, the ZSD is derived based on part I: where represents the minimum diffuse attenuation coefficient value in the transparent spectral window within the visible domain (443-665 nm). For GOCI images that were captured in Lake Hongze and Lake Taihu, the value of (555) was used to represent in this study area. At the same time, a good line regression between (555) and (745) was found (R 2 = 0.89). After that, the values of b and c in Equations (2)-(10) are 0.99 and -1.96, respectively. The is the remotesensing reflectance (sr−1) of the maximum transparency at the transparent window of the water body, and was taken as the value that corresponds to the wavelength with minimum [10].

Noise-Equivalent ZSD
The signal to noise ratio (SNR) determines whether a sensor can be used for parameter retrieval of water quality and what level of spatial or temporal averaging is required to reach the desired performance, and the noise-equivalent reflectance ( ) was used to measure the performance of the sensors [31]: where 0 is the extraterrestrial solar irradiance and is the sun zenith angle. The earth-sun distance in Astronomical Units, , is approximately set to 1 AU, as the variability of due to changes in is much larger. is the noise-equivalent radiance: where is the signal-to-noise ratio at the reference radiance, . 0, , and values for different bands are given in [32]. The noise-equivalent ZSD can be estimated from while using our developed algorithm. Table 2 summarized the basic statistics of the field measured water quality parameters in Lake Hongze (HZ), Lake Taihu (TH), and Lake Dongting (DT). The parameters from these lakes showed a coefficient of variation (CV) between 17.36% and 142.64%, which indicated high dynamic ranges and substantial variability amongst the samples in the study areas. Although differences in the concentration of water parameters are present in these three lakes, they are all characterized by high turbidity. Table 2. Descriptive statistics from Lake Hongze (HZ), Lake Taihu (TH), Lake Dongting, and Lake Erhai (EH).

Parameters
Statistics Aver-average, SD-standard deviation, CV-coefficient of variation In Lake Hongze, the maximum values of Chla, TSM, and ZSD were more than 15 times larger than the minimum values. The Chla, TSM, and ZSD were 1.39-149.32 μg/L, 7.18-193.33 mg/L, and 0.1-0.8 m, respectively. In Lake Hongze, the share of OSM in TSM is smaller than that of ISM, which demonstrated that ISM dominated the TSM in Lake Hongze. In Lake Taihu, the Chla, TSM, and ZSD were 9.23-301.93 μg/L, 8.66-96.47 mg/L, and 0.14-1.1 mg/L, respectively. The maximum values of Chla, TSM and ZSD were 32, 11, and 7 fold larger than the minimum values. The ratio of OSM to TSM (OSM / TSM > 0.5) indicates that suspended particulates in Lake Taihu are mainly dominated by organic suspended particulates. In Lake Dongting, the Chla, TSM and ZSD were 2.79-52.08 μg/L, 3.75-200.53 mg/L and 0.15-1.05 mg/L, respectively. In Lake Erhai, the Chla, TSM and ZSD were 7.27-34.34 μg/L, 1.96-7.5 mg/L, and 0.15-4 mg/L, respectively. Overall, Lake Hongze, Lake Taihu, Lake Dongting, and Lake Erhai showed high levels of suspended matter and Chla, but low ZSD, which showed huge concentration differences with other inland waterbodies [4,[33][34][35], indicating that Lake Hongze and Lake Taihu have the characteristics of highly turbid waters.
The average Rrs(λ) spectra in four lakes (HZ, DT, TH, and EH) had distinct characteristics, as shown in Figure 3a. Note that Lake Erhai has not field measured ( ) data used in this study. When compared with Lake Hongze and Lake Erhai, Lake Taihu and Lake Dongting have a higher spectral reflection curve due to the high intensity reflections. In addition, there is an obvious reflection peak and absorption trough appearing at 550 nm and 680 nm, respectively. Similar curve variations of ( ) were observed in Lake Taihu, Lake Dongting and Lake Hongze, all of which had a decreasing tendency from 400 to 570 nm and an increasing tendency from 700 to 750 nm Figure 3b. At the same time, the diffusion attenuation coefficient at 555 nm (K (555)) presents a minimum value within the visible domain (443-665 nm) rather than other bands in GOCI images.

Algorithm Validation and Noise-equivalent ZSD
Although the model presented a significant improvement in the robustness to turbid water, comprehensive testing using in situ data that were obtained under various water conditions is still critical for algorithm validation and subsequent application to these extremely turbid waters. The remaining dataset (95 samples) was used to validate the semi-analytical algorithm of ZSD. Note that only 48 validated samples have field measured ( ) data. The in-situ measured and corresponding retrieved values of (555) and ZSD were distributed well along the 1:1 line (Figure 4). The MAPE and RMSE of match-up points between the measured and retrieved (555) were 11.31% and 0.63 m -1 , respectively (Figure 4a). This suggests that the modified steps in the new algorithm (step 3, 5, and 6) performed well with satisfactory accuracy and had great potential to further estimate the ZSD in turbid lakes. The MAPE and RMSE of match-up points between the measured and retrieved ZSD were 22.39% and 0.24 m, respectively (Figure 4b), which demonstrates that the performance of the semi-analytical algorithm of ZSD had great potential for estimating ZSD in satellite images.
Remote sensing of water clarity relies not only on validated algorithms for ZSD, but also on the quality of remotely sensed imagery. Some previous studies [31,36] showed significant differences of imagery quality in the difference sensors between open sea and inland turbid waters, raising the question of whether GOCI imagery are suitable for the algorithm that we proposed for mapping ZSD. Finally, the noise-equivalent ZSD of 0.084 m was calculated based on our developed algorithm. When considering its well noise-equivalent ZSD, GOCI imagery can be used for ZSD estimation while using our developed algorithm in extremely turbid inland waters, such as Lake Taihu and Lake Hongze.

Atmospheric Correction Assessment by Synchronized Images
The performance of the 6S model mentioned in Section (2.3) and the ZSD deriving algorithm were evaluated by comparing the GOCI-derived and the synchronized field measured Rrs at 443, 560, 665, 674, 709, and 754 nm ( Figure 5), while the field measured reflectance obtained by ASD measurement were integrated while using the GOCI spectral response functions (SRF) to simulate the Rrs at the ground level. There are 42 synchronized in-situ samples in Lake Hongze and 30 samples in Lake Taihu.
For both turbid lakes, the MAPE values at 490, 555, 660, 680, and 745 nm were less than 30%, and poor performance of atmospheric correction at 412 nm was found (MAPE > 30%), which might be related to the characteristic that the spectrum is prone to saturation at the short blue wavelengths. [37]. At the same time, poor performance existed at 865 nm with MAPE > 32% due to the high absorbing properties in NIR wavelengths for turbid water [38]. Overall, the GOCI-derived and field obtained Rrs(λ) were close to the 1:1 line. Given the superior performance of the 6S algorithm, it is reasonable to use it for atmospheric correction of GOCI data.
For the performance of ZSDZ algorithm, the in situ measured ZSD and the corresponding retrieved values were close to the 1:1 line (Figure 5c). The MAPE and RMSE for the validation of the developed models were 29.07% and 0.12 m, respectively. This suggests that our proposed retrieval models that were based on GOCI data performed well.   Figure 5. The field measured reflectance corresponded to the Geostationary Ocean Color Imager (GOCI) bands for Lake Taihu (A) and Lake Hongze (B). The performance of new algorithm based on synchronized images (C).

Mapping ZSD from GOCI Based on Developed Algorithm
The GOCI Level-1B data, captured on, 4 August 2018 in Lake Taihu and 5 November 2016 in Lake Hongze, were used to estimate ZSD based on our developed algorithm. The newly developed semi-analytical algorithm was applied to GOCI atmosphere corrected images in order to derive the hourly scale of ZSD in Lake Taihu and Lake Hongze ( Figure 6). The ZSD retrieval results in the eastern part of Lake Taihu are not shown due to the impact of submerged vegetation and bottom reflectance. The hourly scale maps show the consecutive spatiotemporal dynamic characteristics of ZSD in both Lake Taihu and Hongze with low ZSD values within one day. Lake Taihu and Lake Hongze are shallow inland turbid waters, which are susceptible to environmental and anthropogenic factors, and they have significantly different characteristics from other inland waterbodies.
Significant spatiotemporal variability in ZSD was present from 8:30 to 15:30 in Lake Taihu. The spatiotemporal distribution of ZSD was significantly regulated by particulate suspended matter in Lake Taihu as a result of wind speed and wind direction [5,39]. The distribution of low ZSD is mostly concentrated in the southwest areas of Taihu Lake, whereas high values were observed in central open zones and lake bays (such as Meiliang Bay, as shown in Figure 1). Sheltered conditions and large amounts of submerged aquatic vegetation, which maintain relatively high transparency due to the water filtration, wind-induced wave attenuation, and the inhibition of sediment resupension by aquatic vegetation characterize these areas [5,39]. In conclusion, wind conditions and submerged aquatic vegetation can both regulate the distribution of ZSD in Lake Taihu. In addition, a key finding was that the ZSD value of Lake Taihu shows a distinct increasing trend from 8:30 in the morning to 15:30 in the afternoon, which is probably due to the continued weakening of wind speed from morning to afternoon. Figure 6. Hourly variations of GOCI-derived ZSD based on ZSDZ model in Lake Taihu (A) and Lake Hongze (B). The retrieval results of ZSD in the eastern part of Lake Taihu are not shown due to the dense distribution of aquatic vegetation (see Figure 1).
A different distribution pattern of ZSD was observed in Lake Hongze than in Lake Taihu, which is, the transparency value showed a significant downward trend from 8:30 am to 15:30 pm. Temporally, the ZSD in the eastern and northern part of Lake Hongze go through a process from high to low, whereas the ZSD in the southwestern Lake remains low from 8:30 to 15:30. Spatially, the distribution of high ZSD is concentrated in the north and west of Lake Hongze, whereas the low values were observed in the southern region, especially in the upstream riverine region (Figure 6b). The headwater of Huai River accounts for more than 70% of the annual mean incoming runoff of Lake Hongze; the vigorous movement of water is one of main factors for the low transparency in these areas [19].

The Relationship between In-Situ ZSD and Water Constituent Concentrations
In the present study, the field measured ZSD values were low, with 0.15-0.8 m in Lake Hongze and 0.15 m-1.0 m in Lake Taihu, indicating that extreme turbidity characterizes these two lakes. These findings were similar to those found in other turbid inland lakes [40,41], as well as previous studies in Lake Taihu [42]. However, the ZSD values in Lake Hongze and Lake Taihu are obviously lower than those for oceanic and coastal waters [10,11], and slightly turbid waters [4,6]. For example, the mean value of the field measured ZSD in Lake Hongze is six times lower than that for Lake Qiandao [4].
Suspended particulate matter, phytoplankton, submerged vegetation, and the surface turbulence of the lake water affect ZSD [5,43]. Among these, the concentrations of main optical components (such as Chl-a and TSM) and their relative contributions to ZSD vary for different waters [4,5,8]. For example, in open marine waters, the optical properties are mainly dominated by phytoplankton and its attachments substantially contribute to the variability in ZSD. In contrast, particulate sediment and debris may have a significant effect on water clarity in turbid lakes [5].
Robust positive correlation coefficients between ZSD and TSM were found in both Lake Hongze and Lake Taihu (Figure 7a), indicating that ZSD has a tight relationship with TSM; more than 60% of the ZSD dynamic can be explained by the TSM in these two lakes. Relative higher correlation coefficients between ZSD and ISM were also found in Lake Hongze (R 2 = 0.60, p < 0.01) and Lake Taihu (R 2 = 0.59, p < 0.01), whereas low values between ZSD and OSM were in both lakes. Note that the coefficient in Lake Hongze (R 2 = 0.31, p < 0.01) is higher than that in Lake Taihu (R 2 = 0.08, p < 0.01), indicating that there are different water compositions between these two lakes. These results conclude that inorganic suspended matter (ISM), rather than organic suspended matter (OSM), plays an important role in water quality parameters that affects transparency. In addition, the relationships between field measured ZSD and field measured Chla in both lakes were poor (Figure 7c,d). The ratio of OSM/TSM is commonly used to indicate the relative quantity of ISM and OSM in water bodies. The mean values of OSM/TSM in Lake Hongze and Lake Taihu are 0.14 and 0.46, respectively, which suggests that inorganic suspended matter is dominant in Lake Hongze and Lake Taihu. This is different from other inland lakes (like Lake Chaohu and Lake Dianchi) that are dominated by organic suspended matter [44,45]. The low ZSD values are characterized by high turbidity and high TSM. High TSM caused by sediment resuspension significantly decreased the water clarity [46]. Moreover, sediment resuspension can be induced by meteorological and hydrological events, such as rainstorms, strong winds, and human activities, such as sailing and trawling [47]. The intense water currents and sand mining activities in Lake Hongze cause sediment resuspension and reduce water transparency, especially in the estuarine part of the lake, according to the study by Lei et al. [20].

Comparison with the Exist ZSD Algorithms for GOCI
An original and several modified algorithms for estimating ZSD in slightly turbid or turbid waters have been proposed [6,11,12]. These algorithms were selected for comparison with our developed model in this study (Table 3 and Figure 8). Model 4 (ZSDV6) was the original semi-analytical algorithm of ZSD, and this model was developed and validated with large data covering inland, marine, and coastal waters [10]. Model 1 (denoted as ZSDM14) and 2 (denoted as ZSDR17) were developed by parameterizing and calibrating for application in highly turbid coastal and slightly turbid waters [6,48].
The results show the diverse performance of the different algorithms. Model 3 and 4 performed better than model 1 and 2, and the performance of Model 4 was best with relatively low MAPE (MAPE = 22.39%) and RMSE (RMSE = 0.24 m), but high R2 (R 2 = 0.89). Model 3 performed well for low ZSD, but it had an overestimation effect for high ZSD, which was consistent with a previous finding in slightly turbid waters [6]. Among the existing models, model 1 and 2 significantly underestimated ZSD, which suggested that the original ZSD algorithm that was proposed by Lee et al. [10] was more suitable than others in turbid inland waters. The results indicate that, although the three existing models can be used to estimate ZSD in turbid inland waters, their performance was highly impacted by the local parameters and parameterization of the inherent optical properties, the new proposed model showed the best estimation effect.  Figure 8. Comparison of the measured and derived ZSD of the four Models. ZSDM14, ZSDR17, and ZSDV6 were provided by previous studies, ZSDZ is provided by this study.

Model Parameterization
A suitable model is needed to estimate the total absorption and backscattering coefficients and further retrieve the water transparency (ZSD) in turbid lakes. The values of ZSD were underestimated or overestimated based on existing algorithms ( Figure 6), which was probably due to these algorithms introducing larger uncertainty into the derived IOPs and further leading to the poorer performance of ZSD in turbid inland waters. Lee et al. [12] noticed the underestimated IOPs using original QAA for turbid waters and concluded that the larger absorption effect of particulates at short wavelengths results in the weak accuracy of the model. Although several modified versions have been proposed for slightly turbid and coastal waters [6,48], the uncertainty analysis of semi-analytical algorithms for deriving ZSD in extremely turbid lakes is seldom reported.
The uncertainty analysis that was carried out by previous studies [6,28] suggested that the reference wavelength shifted to longer wavelengths in turbid waterbodies that have low absorption but high scattering coefficients. Röttgers et al. [49] also found that the near infrared (NIR) absorption of particulates was approximately equal to the absorption of water in turbid waters. Liu et al. [50] noticed that selecting the appropriate NIR band as the reference wavelengths depends on varying types of turbid waters. In this study, the absorption in reference wavelength was assumed at 745 nm (step 3 in Table 1), which took into consideration that the absorption of pure water should dominate in this band. However, numerical difference was found between a(λ ) and a (λ ) in step 3, which would propagate to the estimation of IOPs in step 3-6, and then to the estimation of ZSD in step 7.
In step 2, g0 and g1 are both empirically constant and the optimized values can improve the estimation of the IOPs value. Ogashawara et al. [51] indicated that the derived u(λ) based on fixing the g0 and g1 values can result in an underestimation of a(λ) for different waters. The g0 = 0.089 and g1 = 0.1245 are suggested in the three versions of QAA (v4-6) for Case II waters [10,28]. Lee et al. [52] also demonstrated that g0 = 0.084 and g1 = 0.17 for higher scattering coastal waters and the same values are also applied in this study when considering the water properties of Lake Hongze and Lake Taihu. For (λ ) retrieval in step 5, the values of coefficients (m0~3) were also retuned while using the field measured datasets and the model brought the estimated K (745) close to the field measured values (Figure 3a), indicating the necessity of parameter readjustment in our study.

Measurement Uncertainties
The uncertainty from field measured and satellite-derived Rrs(λ) also has a negative impact on the estimation of the IOPs and ZSD. The measurement of spectral data followed the above-water measurement method in this study. It is worth noting that the water-leaving radiance was obtained by correcting the total radiance of the water surface ( ) while using the coefficient of skylight reflectance ( ), which depended on the water surface conditions that are caused by wind conditions [21]. Zheng et al. [21] suggested that = 0.022 for calm conditions, 0.025 for < 5 m/s wind speed, and 0.026-0.028 for 10 m/s wind speed. According to field measured speed data that we collected, low wind speeds are dominant in Lake Hongze and Lake Taihu, and = 0.022 was applied in this study. However, the variation of would result in a high variation of mean coefficient variance value (CV = ~50%) of Rrs at the spectral range of 400 to 800 nm [15], which definitely has a negative impact on the performance of ZSD models. Moreover, although a reliable method was used to reduce sun glint and ensure the accurate measurement of reflectance, it is still possible for glint to remain. In addition, the inherent noise that is generated by the sensor is also an influencing factor.
In this study, the atmospheric correction was based on the 6S method to obtain the GOCIderived Rrs. However, slight overcorrection of Rrs(745) and Rrs(555) were also found, which, as key input parameters, may introduce uncertainty into the derived IOPs and further lead to the poorer performance of ZSD. The sensibility of using improved semi-analytical model of ZSD for the input parameters (Rrs(745) and Rrs(555)) were analyzed, the result indicated that this model is highly sensitive to the atmospheric corrected GOCI-derived Rrs (Figure 9). Therefore, the retrieved IOPs rely on the accuracy of corrected Rrs directly, indicating that an excellent atmospheric correction model with high accuracy is necessary.

Limitation of Model in Retrieving ZSD
For the ZSDZ algorithm, one of the key steps (step 3 in Table 1) is to shift the reference wavelength to longer NIR band in order to satisfy an assumed condition that the (λ ) is dominant in NIR wavelength. Additionally, the relatively high (λ ) value in extremely turbid waters in NIR wavelength provides a higher signal-to-noise ratio to ensure the algorithm development. Conversely, for some slightly turbid waters or clear waters, a relatively low (λ ) magnitude at longer NIR wavelength implies a low SNR, which might introduce large uncertainty to estimate (λ ). Although we do not have enough datasets to validate whether the developed algorithm can effectively estimate ZSD in slightly turbid or clear waters, caution should be exercised in applying the newly developed algorithm to the above waters. One empirical step (step 5 in Table 1) was constructed based on field measured data that we collected, which is location dependent and may not be suitable to other waters. Therefore, when our algorithm is applied to other water bodies, it might be necessary to reconstruct similar empirical equations based on the optical properties of water.
The water quality and optical parameters of Lake Hongze and Lake Taihu have high variability, as presented in Table 2 and Figure 3. The Rrs, for example, has high variability of coefficient variance at 555 nm (CV = 54%), which is probably due to larger contribution of backscattering from particulates in extremely turbid water. Additionally, high variability of Chla (CV = 142.64%) in Lake Hongze and the ISM (CV = 90.66%) in Lake Taihu was observed. Obviously, high spatial dynamic of bio-optical properties in this study area limit the accuracy of ZSD model. Therefore, for waterbodies with a high spatial dynamic of bio-optical properties, the retrieval of IOPs and the estimation of ZSD may be compromised based on ZSDZ algorithm.

The Response of GOCI-Derived ZSD to Wind Speed
Lake Hongze and Lake Taihu are characterized by shallow and turbid water, with an average depth of no more than 2 m, the wind forcing, river sediment discharge, and the short-term rainfall often lead to strong water mixing and suspension of particles [20,46]. This process results in an increase in TSM concentration, which then leads to weak light transmission and low water transparency. In addition, low wind speeds in enclosed areas often promote the growth of submerged aquatic vegetation with the functions of water filtration, water purification, and sediment resuspension reduction [5,47]. Figure 10 presents the hourly variation of GOCI-derived ZSD and the corresponding wind speed. Strong relationships between GOCI-derived ZSD and wind speed are found in both Lake Hongze (R 2 = 0.57, p < 0.01) and Lake Taihu (R 2 = 0.63, p < 0.01), and the ZSD decreases with a reduction in wind speed, indicating that wind greatly affects the ZSD. To further understand the connection between ZSD and wind speed, wind speed was divided into three levels: low (0-3 m/s), medium (3-4 m/s), and high (> 4 m/s), according to the intensity of the wind [53]. When the wind speed was low, the ZSD varied from 0.08 m to 0.9 m with a mean value of 0.29 ± 0.21 m, and 84% of the variation of ZSD can be explained by wind speed. However, when the wind condition was medium or high, the ZSD ranged from 0.05 to 0.24 m with a mean value of 0.13 ± 0.05 m, and a weak relationship is observed between the wind speed and ZSD (Figure 10b). The results indicate that low wind speeds are dominant in Lake Taihu and Lake Hongze, which is good for the growth of submerged aquatic vegetation. However, an increase in ZSD is observed in both Lake Hongze and Lake Taihu, which might be related to anthropogenic activities [20,42].

Conclusion
An improved semi-analytical scheme (ZSDZ) was developed to estimate water clarity in turbid waters. When compared with original ZSDV6, the new method has three main extensions. Firstly, the reference wavelength for estimating (λ ) was relocated at 745 nm to reduce the adverse contribution from the optical components in extremely turbid waters. Second, four coefficients (m0 -m3) were modified based on field measured datasets to better estimate (λ ). Third, target (λ) was derived based on the good empirical relationship between (λ) and (λ ). The validation on a comprehensive field measured ZSD dataset (69 samples) demonstrated that ZSDZ improved the retrieval accuracy of ZSD in turbid waters and outperformed the three existing semi-analytical schemes. The hourly variation patterns of ZSD in Lake Hongze and Lake Taihu were also discovered after applying the ZSDZ to GOCI data. A distinct difference was presented in the spatiotemporal distribution pattern of ZSD between the two lakes. The transparency value showed a significant downward trend from 8:30 am to 15:30 pm in Lake Hongze; conversely, a positive change in transparency was observed during the same time in Lake Taihu. Moreover, it has been found that wind is an important regulating factor in the spatiotemporal variation of ZSD. Our findings suggest that our proposed algorithm that is based on GOCI imagery can be suitable for capturing and analyzing the hourly variation and the dynamic characteristics of transparency of turbid waterbodies in a good manner, providing an important source of information on monitoring and mitigating the adverse effects on ecosystems.