Coastal Water Clarity in Shenzhen: Assessment of Observations from Sentinel-2

: Shenzhen is a crucial city in the Guangdong–Hong Kong–Macao Greater Bay Area (GBA). With high-intensity land development and rapid population growth, the ocean has become an essential space for expansion, leading to signiﬁcant variations in water quality in the coastal area of Shenzhen. Water clarity (Z sd ) is a key indicator for evaluating water quality. We applied the quasi-analytical algorithm (QAA) to Sentinel-2 data and retrieved the Z sd of the coastal area of Shenzhen. By adjusting the red band for distinguishing water types, we avoided underestimating Z sd for clear water. This study pioneered the production of a 10 m Z sd product for the coastal area of Shenzhen from 2016 to 2021. The results showed that the coastal area of Shenzhen exhibited a spatial distribution pattern with low Z sd in the west and high in the east, with Pearl River Estuary (PRE: 0.41–0.67 m) and Shenzhen Bay (SZB: 0.30–0.58 m) being lower than Dapeng Bay (DPB: 2.7–2.9 m) and Daya Bay (DYB: 2.5–2.9 m). We analyzed the seasonal and interannual variations and driving factors of the four areas, where PRE and SZB showed similar variation patterns, while DPB and DYB showed similar variation patterns. PRE and SZB are important estuaries in southern China, signiﬁcantly affected by anthropogenic activities. DPB and DYB are important marine aquaculture areas, mainly affected by natural factors (wind speed, precipitation, and sea level). The Z sd of the coastal area of Shenzhen, along with the analysis of its results and driving factors, contributes to promoting local water resource protection and providing a reference for formulating relevant governance policies. It also provides a practical method for assessing and monitoring near-shore water quality.


Introduction 1.Background
Shenzhen, located in the South China Plain of China, is an important city in the GBA agglomeration.The western part of Shenzhen is adjacent to the PRE, SZB, and Lingdingyang Estuary (LDYE).The eastern part of Shenzhen is adjacent to the DPB and DYB, with a vast water area and rich marine resources.As a coastal city with scenic beauty, the water quality variations in the coastal area of Shenzhen are attributed to anthropogenic and urbanization activities.Decades of rapid development and urbanization have resulted in an overload of pollutant emissions, accelerating the deterioration of the coastal area of Shenzhen, causing shortages in water resources and affecting local development [1].For example, plenty of marine pollution incidents-such as red tides and oil spills-have been reported in PRE and DYB [2,3].Therefore, sustainable and daily monitoring of the coastal area is increasingly important for the environmental protection of Shenzhen, in order to ensure the life quality of its residents.
Water clarity provides visual perception of and vital information on the quality of aquatic ecosystems.At present, some parameters can quantitatively express Z sd .Previous

Purpose of this Study
Numerous previous studies have analyzed water in the coastal area of Shenzhen, especially PRE.For example, in turbid PRE, the research and evaluation of SPM were performed using retrieval algorithms and remote sensing data, including its long-term variations and driving factors [18,19,[33][34][35][36][37][38].Nazeer et al. (2016) proposed a Chla prediction and classification model using remote sensing data for PRE with complex optical characteristics [39,40].Ma et al. (2022) used Machine Learning and MODIS to retrieve the Chla of PRE, and conducted research on long-term variations [41].These studies focused on the variations in the water quality of PRE in terms of spatial distribution and interannual variations; however, limited research exists on the spatiotemporal seasonal variations in water quality near Shenzhen.
Given the urgent local needs, this study investigates the spatiotemporal variations and driving factors of Z sd in the coastal area of Shenzhen.The study aims to achieve the following objectives: (1) To construct a Z sd retrieval model for the coastal area of Shenzhen based on Sentinel-2, and to distinguish clean and turbid water using the red band to avoid the underestimation of clean water.(2) Using the model proposed in this study, to Water 2023, 15, 4102 3 of 18 evaluate the spatiotemporal variations of Z sd for the period of 2016-2021.(3) Analyze the relationship between Z sd and several key natural factors (precipitation, wind speed, and sea level).

Study Area
Shenzhen is an important port located in southeastern China.It is also an important industrial maritime core city in China.Shenzhen is adjacent to PRE and LDYE in the west and includes DYB and DPB in the east, with a vast water area under its jurisdiction.PRE is China's second estuary (runoff) and an important waterway.The PRE, with a total length of approximately 22.6 km and a drainage area of about 4.53 × 10 5 km 2 , is the largest basin in southern China.The PRE mainly comprises the West, North, and East Rivers.These rivers flow into the South China Sea through eight western outlets (see in Figure 1).The water quality of both DYB and DPB is the best in Guangdong province.With an excellent geographical location, they are natural ports, scenic tourist areas, and aquaculture areas.DYB is located between Honghai Bay and DPB in the eastern Guangdong Province, with a total area of 650 km 2 ; DPB is located between the Dapeng and Kowloon Peninsulas of Hong Kong, with a total area of about 335 km

Study Area
Shenzhen is an important port located in southeastern China.It is also an important industrial maritime core city in China.Shenzhen is adjacent to PRE and LDYE in the west and includes DYB and DPB in the east, with a vast water area under its jurisdiction.PRE is China's second estuary (runoff) and an important waterway.The PRE, with a total length of approximately 22.6 km and a drainage area of about 4.53 × 10 5 km 2 , is the largest basin in southern China.The PRE mainly comprises the West, North, and East Rivers.These rivers flow into the South China Sea through eight western outlets (see in Figure 1).The water quality of both DYB and DPB is the best in Guangdong province.With an excellent geographical location, they are natural ports, scenic tourist areas, and aquaculture areas.DYB is located between Honghai Bay and DPB in the eastern Guangdong Province, with a total area of 650 km 2 ; DPB is located between the Dapeng and Kowloon Peninsulas of Hong Kong, with a total area of about 335 km 2 .Precipitation data for Guangzhou (https://swj.gz.gov.cn,accessed on 10 August 2023) and Shenzhen (https://swj.sz.gov.cn,accessed on 10 August 2023) were obtained from the monthly report of the Water Authority.Sea level data were obtained from the China Sea Level Bulletin (https://www.mnr.gov.cn/sj/sjfw/hy/gbgg/zghpmgb,accessed on 29 July 2023).Wind speed was obtained from the National Centers for Environmental Information (NCEI) of the National Oceanic and Atmospheric Administration (NOAA) (https://www.ncei.noaa.gov/data/global-summary-of-the-day/archive/,accessed on 29 June 2023).Maritime transportation information was obtained from Shiptracker.org(https://ship-tracker.org/marinetraffic,accessed on 20 June 2023).

Pre-Processing
We used images with cloud cover of ≤30% as valid data.If cloud cover was ≥30%, such data were excluded.The study selected Sentinel-2 L1C TOA (started in 2016) and Sentinel-2 L2A SR (started in 2019 in China) data from Sentinel-2.First, L1C TOA data for 2016-2018 were corrected using the Sen2Cor 02.11.00 software released by ESA to obtain SR data.Second, SNAP 9.0.0 software released by the ESA was used to resample SR data for 2016-2021, resulting in SR data with a spatial resolution of 10 m.Finally, cloud and mist removal was performed on SR data using the QA60 band (see in Figure 2).2023).Wind speed was obtained from the National Centers for Environmental Information (NCEI) of the National Oceanic and Atmospheric Administration (NOAA) (https://www.ncei.noaa.gov/data/global-summary-of-the-day/archive/,accessed on 29 June 2023).Maritime transportation information was obtained from Shiptracker.org(https://ship-tracker.org/marinetraffic,accessed on 20 June 2023).

Pre-Processing
We used images with cloud cover of ≤30% as valid data.If cloud cover was ≥30%, such data were excluded.The study selected Sentinel-2 L1C TOA (started in 2016) and Sentinel-2 L2A SR (started in 2019 in China) data from Sentinel-2.First, L1C TOA data for 2016-2018 were corrected using the Sen2Cor 02.11.00 software released by ESA to obtain SR data.Second, SNAP 9.0.0 software released by the ESA was used to resample SR data for 2016-2021, resulting in SR data with a spatial resolution of 10 m.Finally, cloud and mist removal was performed on SR data using the QA60 band (see in Figure 2).
To remove residual noise from SR data, we applied a simple correction method for Sentinel-2 SR [42].The following correction formula was used:

Mean Aggregation
Owing to the significant randomness of the spatiotemporal variations of water clarity in the coastal area of Shenzhen analyzed using single scene Zsd, the calculation of the monthly Zsd can be more accurate.If there are multiple valid data in a month, their mean value is calculated as the monthly result, and used as the calculation unit to obtain the To remove residual noise from SR data, we applied a simple correction method for Sentinel-2 SR [42].The following correction formula was used: where R c rs (λ) is the corrected water-leaving reflectance of the Sentinel-2 band centered at (λ) i , R rs (λ) is the original reflectance of the Sentinel-2 band centered at (λ) i , and min (R rs (NIR): R rs (SWIR)) is the minimum positive reflectance value of the NIR and SWIR bands.

Mean Aggregation
Owing to the significant randomness of the spatiotemporal variations of water clarity in the coastal area of Shenzhen analyzed using single scene Z sd , the calculation of the monthly Z sd can be more accurate.If there are multiple valid data in a month, their mean value is calculated as the monthly result, and used as the calculation unit to obtain the quarterly and annual results.Figure 3 is the quantities of Sentinel-2 covering the coastal area of Shenzhen during 2016-2021.

Algorithm to Retrieve Zsd
Lee et al. (2002) found that the classic underwater visibility theory cannot accurately explain the observation of the human eye on the Secchi disk [43], and they proposed an innovative model based on radiation transfer theory to estimate Zsd, which ranges from 0.1 to 30 m.The validation results showed that the estimated results of the model were consistent with those of on-site measurements, with an average absolute difference of 18% (note: the 18% difference includes 10% uncertainty in the on-site in situ measurement of Zsd).
The QAA has good stability and universality without adjusting for any coefficients, and its application is not limited by time, region, or sensors.The QAA has been widely used for marine, nearshore, and inland water.The latest version of QAA, that is, QAA_V6 [44], is widely used [10,16,25,45].Through the relationship between Rrs and inherent optical property, the total absorption and backscattering coefficients of water are retrieved, and then the Zsd is obtained through the total absorption coefficient (See in Figure4 and Table1).The formula is as follows [44]: where ( ) is the minimum diffuse attenuation coefficient at the transparent window of the water within the visible domain (443-665 nm) and  is the remote sensing reflectance corresponding to the wavelength.We used the method described by Qiu et al. (2023) to distinguish between clean and turbid water [25], which improves the usability of the QAA model in retrieving Zsd in the coastal area of Shenzhen and avoids the underestimation of the results of clean water.This provides valuable reference information for improving the use of QAA models to retrieve the Zsd of water with complex optical characteristics.This method provides sustainable good data for the subsequent analysis of the 2.5.Algorithm to Retrieve Z sd Lee et al. (2002) found that the classic underwater visibility theory cannot accurately explain the observation of the human eye on the Secchi disk [43], and they proposed an innovative model based on radiation transfer theory to estimate Z sd , which ranges from 0.1 to 30 m.The validation results showed that the estimated results of the model were consistent with those of on-site measurements, with an average absolute difference of 18% (note: the 18% difference includes 10% uncertainty in the on-site in situ measurement of Z sd ).
The QAA has good stability and universality without adjusting for any coefficients, and its application is not limited by time, region, or sensors.The QAA has been widely used for marine, nearshore, and inland water.The latest version of QAA, that is, QAA_V6 [44], is widely used [10,16,25,45].Through the relationship between R rs and inherent optical property, the total absorption and backscattering coefficients of water are retrieved, and then the Z sd is obtained through the total absorption coefficient (See in Figure 4 and Table 1).The formula is as follows [44]: where min K tr d is the minimum diffuse attenuation coefficient at the transparent window of the water within the visible domain (443-665 nm) and R tr rs is the remote sensing reflectance corresponding to the wavelength.We used the method described by Qiu et al. (2023) to distinguish between clean and turbid water [25], which improves the usability of the QAA model in retrieving Z sd in the coastal area of Shenzhen and avoids the underestimation of the results of clean water.This provides valuable reference information for improving the use of QAA models to retrieve the Z sd of water with complex optical characteristics.This method provides sustainable good data for the subsequent analysis of the spatiotemporal variations of Z sd in the coastal area of Shenzhen.The results obtained can serve as a reference for subsequent studies on water resource management.Step Property Calculation If  (670) < 0.005 ( ) =  ( ) + 10 .Notes: Rrs: above-surface remote-sensing reflectance; rrs: below-surface remote sensing reflectance; u(λ): ratio of backscattering coefficient to the sum of absorption and backscattering coefficients (bb/(a + bb)); λ0: reference wavelength; aw: absorption coefficient of pure water; bbp: backscattering coefficient of suspended particles; bbw: backscattering coefficient of pure seawater; η: spectral power of particlescattering coefficient;  : solar zenith angle.

Evaluation of Zsd
Because we lack synchronous (quasi-synchronous) in situ data in the coastal area of Shenzhen, it proves the applicability and reliability of the QAA algorithm by evaluating the consistency between Landsat8 OLI and Sentinel-2.Because the QAA algorithm has been widely applied to Landsat [30,45,46], we use the result of Landsat8 OLI as a reference.Validation points are randomly selected in PRE, SZB, DPB, and DYB (see in Figure 1).With each validation point as the center, 9 × 9 pixels (Sentinel-2) and 3 × 3 pixels

Evaluation of Zsd
Because we lack synchronous (quasi-synchronous) in situ data in the coastal area of Shenzhen, it proves the applicability and reliability of the QAA algorithm by evaluating the consistency between Landsat8 OLI and Sentinel-2.Because the QAA algorithm has been widely applied to Landsat [30,45,46], we use the result of Landsat8 OLI as a reference.Validation points are randomly selected in PRE, SZB, DPB, and DYB (see in Figure 1).With each validation point as the center, 9 × 9 pixels (Sentinel-2) and 3 × 3 pixels (Landsat8 OLI) are selected.Both Landsat8 OLI and Sentinel-2 flew over on the same day, with times of 20190328, 20190927, 20191123, 20210119, 20210317, 20211205, and 20211227.
Two factors need to be considered while selecting the sampling points.First, if relevant validation points are selected on the edge, they may be affected by optical shallow water or land proximity effect.Second, the spatial resolution of Landsat8 OLI and Sentinel-2 is inconsistent at 30 and 10 m, respectively.If a relevant verification point selected is on the shore, mixed pixels may also affect it.To avoid the relevant influence of such confounding factors, the verification points selected in the study should be more than 1 km away from the shore.
Thus, the evaluation criteria were the coefficient of determination (R 2 ), Pearson correlation coefficient (PCC), and unbiased percentage difference (UPD).
where A is the Z sd of Landsat8 OLI, A' is the Z sd of Sentinel-2, and n is quantity.

Evaluation of Model Applicability
In order to show the universality and stability of the QAA algorithm, the mean values of Z sd at each validation point in 9 × 9 pixels (Sentinel-2) and 3 × 3 pixels (Landsat8 OLI) were calculated separately, and the results showed that R 2 = 0.95, PCC = 0.96, and UPD = 33.6%.The results of Z sd from Landsat8 OLI and Sentinel-2 had a high correlation, further proving that the QAA algorithm has good reliability and applicability in Sentinel-2 (See in Figure 5).

Variations in Water Clarity
Using Sentinel-2 as input, the monthly Zsd products of Shenzhen coastal areas (PRE, SZB, DPB, and DYB) for the period 2016-2021 were produced.The obvious spatiotemporal variations of Zsd were reflected by the monthly climatology products as high in summer and low in spring (Figures 6 and 7).To provide a more intuitive explanation of the variations of Zsd in the coastal area of Shenzhen, the monthly spatial mean values of Zsd

Variations in Water Clarity
Using Sentinel-2 as input, the monthly Z sd products of Shenzhen coastal areas (PRE, SZB, DPB, and DYB) for the period 2016-2021 were produced.The obvious spatiotemporal variations of Z sd were reflected by the monthly climatology products as high in summer and low in spring (Figures 6 and 7).To provide a more intuitive explanation of the variations of Z sd in the coastal area of Shenzhen, the monthly spatial mean values of Z sd for PRE, SZB, DPB, and DYB were calculated separately (Figure 8).The Z sd of DPB and DYB are significantly better than those of PRE and SZB, with DPB having the best Z sd among the four water areas.Owing to the distribution, Z sd in the western coastal areas of Shenzhen (PRE and SZB) was lower than in the east (DPB and DYB).Both PRE and SZB exhibited the pattern of being low in the north and high in the south, low near the inlet and high on the seaward side.PRE and SZB became increasingly similar, and their variation ranges were 0.41~0.67m and 0.39~0.58m, respectively, with relatively small fluctuations.The Z sd of PRE was only slightly higher than that of SZB in February, September, and November.However, PRE and SZB were different from DPB and DYB.The seasonal variations of DPB and DYB were the highest in summer (2.7~2.9 m) and lowest in spring (2.5~2.9 m), with relatively huge fluctuations.DPB and DYB showed a monthly upward trend in Z sd during winter and spring, while DPB showed a monthly upward trend in autumn, and a monthly first downward, and then upward trend in summer.In summer, the Z sd variations in DPB and DYB were opposite, with DYB exhibiting a downward trend and then an upward trend.These spatial distribution patterns are consistent with PRE and SZB, with both showing low values in the north and high values in the south.
Trend variations were analyzed using the Mann-Kendall (M-K) significance test [47,48].Figures 9 and 10 show the statistical results of the interannual variations and mean annual variations of Z sd .Linear regression results are presented in Table 2.The Z sd of DPB (1.26~2.98 m) and DYB (1.41~3.27m) were higher than those of PRE (0.41~0.61 m) and SZB (0.39~0.63 m), and the Z sd of PRE and SZB were very similar.The interannual variations of Z sd for DPB and DYB were significantly higher than those of PRE and SZB.The trends of DPB and DYB were similar, while those of PRE and SZB were also similar.DPB and DYB only reached their highest levels in the autumn of 2019, while the rest reached their highest in the summer.In 2016, 2017, 2019, and 2020, Z sd reached its lowest point in the spring, and in 2018 and 2021 it did so during the autumn.Analyses of the annual mean of Z sd from different waters during 2016-2021 showed that the Z sd variations in PRE and SZB showed a gentle trend, but both showed a slow upward trend.The Z sd variations in DPB and DYB were significant, showing a downward trend from 2016 to 2019 and an upward trend from 2019 to 2021.

Uncertainties in Z sd Estimation
According to previous research, the QAA algorithm has good application capabilities for retrieving Z sd (Z sd : 0.1~30 m) in marine and coastal water [44,49], and it is not limited by region, time, or sensors.However, owing to regional differences, the water quality of PRE, SZB, DPB, and DYB was relatively clear, thus making the application of semianalytical algorithms in the study area a challenge.However, QAA_V6 applies a red band threshold to distinguish turbid and clean water [25,45,46], modifies the QAA_V6 parameters [23,25], demonstrates good universality with significant regional differences, and further improves the ability of QAA to identify regional differences in water.On the one hand, owing to the lack of quasi-synchronous field measurement data in the study area, a correlation was conducted with the results of Landsat8 OLI Z sd .To widely apply the semi-analytical algorithm to regional differences in water, a large amount of synchronous or quasi-synchronous in situ data may improve the application of the algorithm in relevant areas.On the other hand, Z sd is managed by Kd(λ) and R rs (λ), and Kd(λ) is determined by the calculation of a and bb using QAA.Owing to regional differences, the QAA has more accurate results for a and bb by changing the reference wavelength.Because the absorption coefficient of clean water is low, QAA_V6 is more suitable.Owing to the high absorption of turbid water, QAA_V6 generates estimation errors due to the underestimation of a [25].

Driving Forces of Water Clarity
The coastal area of Shenzhen is an important estuary in southeastern China; the driving factors affecting water quality were mainly divided into natural factors (precipitation, wind speed, and sea level) and human factors (navigation and land reclamation) for the purpose of this study.The correlation coefficients between Z sd and three natural factors were counted, respectively.PRE has a negative correlation with precipitation and positive correlations with wind speed and sea level (See in Figure 11a and Table 3).When precipitation increases, the discharge of pollutants and sediments in multiple estuaries of PRE also increases, leading to plenty of pollutants and sediments being discharged into the sea through the estuary, resulting in a decrease in Z sd .As PRE is a typical subtropical estuary with ocean monsoon climate characteristics (southwesterly winds in summer and northeasterly winds in winter), variations in wind speed may cause Z sd at the estuary mouth to be higher than that on the seaward side due to different forcing mechanisms and synergistic effects, making the hydrodynamic conditions complex [50].Sea level variations accelerate water stratification and extend the residence time of river inputs, which may increase Z sd [51,52].On the other hand, PRE is a micro-tidal estuary, with tidal ranges increasing from 0.8-0.9m near the Wanshan Islands to 1.7 m near Humen [50].Tidal mixing is increasingly important for weakening stratification, especially during the spring tide period, as river discharge is relatively low in winter when tides propagate from southeast to northwest in PRE.Therefore, this phenomenon is weaker in winter and stronger in summer, and the seasonal variations in these hydrodynamic conditions will also cause seasonal variations in Z sd in PRE.
Water 2023, 15, x FOR PEER REVIEW 15 of 20 led to groundwater exploitation and sand excavation activities in the PRE, resulting in settlement and variations in sedimentation.The rise in sea level caused by global warming and local land subsidence may weaken sediment resuspension, and push the maximum turbidity of the estuary towards land [52].In addition, with the increasingly changing structure of the coastline, a decrease in nonpoint source loads from the coast will lead to a decrease in sediment supply [53].In previous studies, factors such as soil erosion [54], meteorological and hydrological conditions [55], dam construction [56], and sand mining [57] have all been well explained.PRE and SZB are important marine transportation areas (see in Figure 12), with the daily heavy discharge of pollutants from marine transportation potentially debasing water quality.Moreover, the trend of Zsd variations in PRE in this study is negatively correlated with the results of SPM variations reported by Nukapothula (2019) [18], further proving that the algorithm used in this study has good stability and universality, and the Zsd results obtained are reliable.11b and Table 3).When precipitation increases, the discharge of pollutants and sediments in various estuaries of SZB also increases, leading to plenty of pollutants and sediments being discharged into the sea through the estuary, thus reducing Z sd .As SZB is a semi-enclosed pocket bay with weaker marine dynamics and poor exchange capacity, it may be that SZB has a low correlation with wind speed and sea level.
As PRE and SZB are important estuaries in southern China, the water environment is significantly affected by humans.The low Z sd near the shore of SZB may be influenced by human activities such as land reclamation.In recent decades, rapid urbanization has led to groundwater exploitation and sand excavation activities in the PRE, resulting in settlement and variations in sedimentation.The rise in sea level caused by global warming and local land subsidence may weaken sediment resuspension, and push the maximum turbidity of the estuary towards land [52].In addition, with the increasingly changing structure of the coastline, a decrease in nonpoint source loads from the coast will lead to a decrease in sediment supply [53].In previous studies, factors such as soil erosion [54], meteorological and hydrological conditions [55], dam construction [56], and sand mining [57] have all been well explained.PRE and SZB are important marine transportation areas (see in Figure 12), with the daily heavy discharge of pollutants from marine transportation potentially debasing water quality.Moreover, the trend of Z sd variations in PRE in this study is negatively correlated with the results of SPM variations reported by Nukapothula (2019) [18], further proving that the algorithm used in this study has good stability and universality, and the Z sd results obtained are reliable.DPB and DYB are not only beautiful vacation destinations, but also important aquaculture zones in the local area.The ecological environment is superior, so the Zsd in this area is mainly affected by natural factors.DPB and DYB are positively correlated with precipitation (See in Figure13 and Table3).When the precipitation increases, it dilutes SPM in water and excreta produced by the aquaculture areas, making the water much clearer.They are negatively correlated with wind speed.The increase in wind speed will disturb SPM in water and release pollutants from sediments, decreasing Zsd.They are also negatively correlated with sea level.Due to the small number of rivers entering the sea in DPB and DYB, that is, the variations of sea level reducing the stratification of water, the amount of river input is decreased, which may be another reason for the decrease of Zsd.
To protect the marine environment and ecological balance, we can adjust and improve the maritime functional areas; strengthen land-based pollution control; and improve the quality of the marine environment, thus enhancing comprehensive utilization of the ocean and ecological restoration of polluted areas, as well as strengthening public education, regional communication, and cooperation in comprehensive prevention and control of marine pollution.DPB and DYB are not only beautiful vacation destinations, but also important aquaculture zones in the local area.The ecological environment is superior, so the Z sd in this area is mainly affected by natural factors.DPB and DYB are positively correlated with precipitation (See in Figure 13 and Table 3).When the precipitation increases, it dilutes SPM in water and excreta produced by the aquaculture areas, making the water much clearer.They are negatively correlated with wind speed.The increase in wind speed will disturb SPM in water and release pollutants from sediments, decreasing Z sd .They are also negatively correlated with sea level.Due to the small number of rivers entering the sea in DPB and DYB, that is, the variations of sea level reducing the stratification of water, the amount of river input is decreased, which may be another reason for the decrease of Z sd .
To protect the marine environment and ecological balance, we can adjust and improve the maritime functional areas; strengthen land-based pollution control; and improve the quality of the marine environment, thus enhancing comprehensive utilization of the ocean and ecological restoration of polluted areas, as well as strengthening public education, regional communication, and cooperation in comprehensive prevention and control of marine pollution.

Implications
Zsd is an important evaluation indicator that can intuitively reflect the clarity and quality of water [56,57].Using Sentinel-2 to retrieve Zsd to evaluate water quality can not only help to obtain the spatial distribution of water quality, but also the continuous and dynamic water quality variations.Therefore, from 2016 to 2021, a 10 m dataset was produced using Sentinel-2 data in the coastal area of Shenzhen, and based on this dataset, an analysis of the spatiotemporal variations and driving factors of Zsd was conducted.According to the analysis results, the overall trend of PRE and SZB was flat, but both showed a slow upward trend.The interannual trends of DPB and DYB were significant, showing a downward trend from 2016 to 2019 and an upward trend from 2019 to 2021.The results indicate that the overall water quality in the coastal area of Shenzhen has improved from 2016 to 2021; Zsd variations in PRE were negatively correlated with the results of SPM variations reported by Nukapothula et al. (2019) [18].Therefore, the results further prove the reliability and applicability of the model.

Conclusions
The study analyzed the spatiotemporal variations of Zsd obtained from Sentinel-2 in the coastal areas of Shenzhen (such as PRE, SZB, DPB, and DYB) for the period of 2016-2021.First, this study processed Sentinel-2 SR data for cloud and mist removal and waterleaving reflectance correction.Second, to avoid the underestimation of Zsd in clean water, the red band was used to distinguish between clean and turbid water, which is more suitable for retrieving the Zsd results in the coastal area of Shenzhen.Finally, spatiotemporal variations of the Zsd and driving factors in the coastal area of Shenzhen from 2016 to 2021 were conducted.The overall spatial distribution pattern of the coastal area of Shenzhen is low in the west and high in the east, and low in the north and high in the south.From 2016 to 2021, the overall water quality of the coastal area of Shenzhen has improved, which may be attributed to the good governance measures taken locally.

Implications
Z sd is an important evaluation indicator that can intuitively reflect the clarity and quality of water [56,57].Using Sentinel-2 to retrieve Z sd to evaluate water quality can not only help to obtain the spatial distribution of water quality, but also the continuous and dynamic water quality variations.Therefore, from 2016 to 2021, a 10 m dataset was produced using Sentinel-2 data in the coastal area of Shenzhen, and based on this dataset, an analysis of the spatiotemporal variations and driving factors of Z sd was conducted.According to the analysis results, the overall trend of PRE and SZB was flat, but both showed a slow upward trend.The interannual trends of DPB and DYB were significant, showing a downward trend from 2016 to 2019 and an upward trend from 2019 to 2021.The results indicate that the overall water quality in the coastal area of Shenzhen has improved from 2016 to 2021; Z sd variations in PRE were negatively correlated with the results of SPM variations reported by Nukapothula et al. (2019) [18].Therefore, the results further prove the reliability and applicability of the model.

Conclusions
The study analyzed the spatiotemporal variations of Z sd obtained from Sentinel-2 in the coastal areas of Shenzhen (such as PRE, SZB, DPB, and DYB) for the period of 2016-2021.First, this study processed Sentinel-2 SR data for cloud and mist removal and water-leaving reflectance correction.Second, to avoid the underestimation of Z sd in clean water, the red band was used to distinguish between clean and turbid water, which is more suitable for retrieving the Z sd results in the coastal area of Shenzhen.Finally, spatiotemporal variations of the Z sd and driving factors in the coastal area of Shenzhen from 2016 to 2021 were conducted.The overall spatial distribution pattern of the coastal area of Shenzhen is low in the west and high in the east, and low in the north and high in the south.From 2016 to 2021, the overall water quality of the coastal area of Shenzhen has improved, which may be attributed to the good governance measures taken locally.
Although the 10 m Z sd dataset produced in this study can provide valuable reference information for achieving sustainable utilization of marine resources, it has certain limi-tations in terms of accuracy.First, we lack the Z sd satellite-ground synchronous in situ data from the study area to evaluate the accuracy of the QAA algorithm and prove its universality.Second, the time series of this study is relatively short, and the data sources and evaluation parameters are limited.This hinders both an in-depth analysis of spatiotemporal variations in water quality and the evaluation of correlations among various water quality parameters.To enhance the reliability of our results, future research will involve multi-source remote sensing data and an analysis of multiple evaluation parameters to assess spatiotemporal variations.Additionally, we will collect more driving factors, such as marine discharge volume and distribution of marine outfalls.By conducting a detailed analysis of such data, we can gain an accurate understanding of the factors contributing to the variations in water quality in the study area, in order to undertake effective measures for the protection and management of water resources.

Figure 1 .
Figure 1.The coastal area of Shenzhen.

2. 2 .
Data Sentinel-2A/2B was launched in 2015 and 2017 by the European Space Agency (ESA), with a revisit period of 10 d per sensor.When the two sensors complement each other, the revisit period can be shortened to 5 d, and in mid-latitude areas, it is reduced to 2-3 d.Sentinel-2 has high resolutions (including temporal and spatial).During 2016-2018, Sentinel-2 only included Level-1C Top of Atmosphere (L1C TOA), which obtained Level-2A Surface Reflectance (L2A SR) by Sen2Cor.Since 2019, Level-2 can be downloaded from the ESA.

Figure 1 .
Figure 1.The coastal area of Shenzhen.

2. 2 .
Data Sentinel-2A/2B was launched in 2015 and 2017 by the European Space Agency (ESA), with a revisit period of 10 d per sensor.When the two sensors complement each other, the revisit period can be shortened to 5 d, and in mid-latitude areas, it is reduced to 2-3 d.Sentinel-2 has high resolutions (including temporal and spatial).During 2016-2018, Sentinel-2 only included Level-1C Top of Atmosphere (L1C TOA), which obtained Level-2A Surface Reflectance (L2A SR) by Sen2Cor.Since 2019, Level-2 can be downloaded from the ESA.
15, x FOR PEER REVIEW 5 of 20 quarterly and annual results.Figure3 is the quantities of Sentinel-2 covering the coastal area of Shenzhen during 2016-2021.

Figure 3 .
Figure 3.The quantities of Sentinel-2 covering the coastal area of Shenzhen during 2016-2021.

Figure 3 .
Figure 3.The quantities of Sentinel-2 covering the coastal area of Shenzhen during 2016-2021.
Water 2023, 15, 4102 6 of 18 Water 2023, 15, x FOR PEER REVIEW 6 of 20 spatiotemporal variations of Zsd in the coastal area of Shenzhen.The results obtained can serve as a reference for subsequent studies on water resource management.
8(λ) b b (λ) Notes: R rs : above-surface remote-sensing reflectance; r rs : below-surface remote sensing reflectance; u(λ): ratio of backscattering coefficient to the sum of absorption and backscattering coefficients (b b /(a + b b )); λ 0 : reference wavelength; a w : absorption coefficient of pure water; b bp : backscattering coefficient of suspended particles; b bw : backscattering coefficient of pure seawater; η: spectral power of particle-scattering coefficient; θ s : solar zenith angle.

Figure 8 .Figure 9 .
Figure 8. Seasonal variations of Z sd in PRE, SZB, DPB, and DYB.(a) Seasonal variations of Z sd in PRE and SZB.(b) Seasonal variations of Z sd in DPB and DYB.Water 2023, 15, x FOR PEER REVIEW 13 of 20

Figure 10 .
Figure 10.Mean annual Z sd of PRE, SZB, DPB, DYB.(a) Mean annual of Z sd of PRE and SZB.(b) Mean annual of Z sd of DPB and DYB.

Figure 11 .
Figure 11.The correlation between the Zsd of PRE and SZB and wind speed, precipitation, and sea level.(a) The correlation between the Zsd of PRE and wind speed, precipitation, and sea level.(b) The correlation between the Zsd of SZB and wind speed, precipitation, and sea level.

Figure 11 .
Figure 11.The correlation between the Z sd of PRE and SZB and wind speed, precipitation, and sea level.(a) The correlation between the Zsd of PRE and wind speed, precipitation, and sea level.(b) The correlation between the Zsd of SZB and wind speed, precipitation, and sea level.

Water 2023 , 20 Figure 12 .
Figure 12.The route of marine traffic over the coastal areas of Shenzhen.

Figure 12 .
Figure 12.The route of marine traffic over the coastal areas of Shenzhen.

Figure 13 .
Figure 13.The correlation between the Zsd of DPB and DYB and wind speed, precipitation, and sea level.(a) The correlation between the Zsd of DPB and wind speed, precipitation, and sea level.(b) The correlation between the Zsd of DYB and wind speed, precipitation, and sea level.

Figure 13 .
Figure 13.The correlation between the Z sd of DPB and DYB and wind speed, precipitation, and sea level.(a) The correlation between the Zsd of DPB and wind speed, precipitation, and sea level.(b) The correlation between the Zsd of DYB and wind speed, precipitation, and sea level. 2 To construct a Zsd retrieval model for the coastal area of Shenzhen based on Sentinel-2, and to distinguish clean and turbid water using the red band to avoid the underestimation of clean water.(2)Using the model proposed in this study, to evaluate the spatiotemporal variations of Zsd for the period of 2016-2021.(3)Analyze the relationship between Zsd and several key natural factors (precipitation, wind speed, and sea level).

Table 2 .
Results of regression analysis performed on an interannual scale of Z sd .

Table 2 .
Results of regression analysis performed on an interannual scale of Zsd.

Table 3 .
Results of correlation analyses performed on interannual Z sd and wind speed, precipitation, and sea level over the PRE, SZB, DPB, and DYB.

Table 3 .
Results of correlation analyses performed on interannual Zsd and wind speed, precipitation, and sea level over the PRE, SZB, DPB, and DYB.