Satellite Remote Sensing of Water Quality Variation in a Semi-Enclosed Bay (Yueqing Bay) under Strong Anthropogenic Impact

: The semi-enclosed bays impacted by heavy anthropogenic activities have weak water exchange and puriﬁcation capacities. Most of the sea bays have suffered severe eutrophication, water quality deterioration, ecosystem degradation and other problems. Although many countries and local governments have carried out corresponding environmental protection actions, the evaluation of their effectiveness still requires monitoring technology and data support for long-term water environment change. In this study, we take Yueqing Bay, the fourth largest bay in China, as a case to study the satellite-based water quality monitoring and variation analysis. We established a nutrient retrieval model for Yueqing Bay to produce a long-term series of nutrient concentration products in Yueqing Bay from 2013 to 2020, based on Landsat remote sensing images and long-term observation data, combined with support vector machine learning and water temperature and satellite spectra as input parameters, and then we analyzed its spatiotemporal variations and driving factors. In general, nutrient concentrations in the western part of the bay were higher than those in the eastern part. Levels of dissolved inorganic nitrogen (DIN) were lower in summer than in spring and winter, and reactive phosphate (PO 4 -P) levels were lower in summer and higher in autumn. of in nutrient concentrations (the variation rates: a DIN > 0, a PO 4 − P > 0), which directly or indirectly caused high nutrient concentrations in some areas of the bay (e.g., Southwest Shoal at the mouth of Yueqing Bay). After 2015, the various ecological remediation policies adopted by cities around Yueqing Bay have mitigated, to some extent, the increasing nutrient concentration trends (the variation rates: a DIN < 0, a PO 4 − P < 0), but not signiﬁcantly ( P > 0.1). The environmental restoration of Yueqing Bay also requires continuous and long-term ecological protection and restoration work to be effective. This research can provide a reference for ecological environment monitoring and remote sensing data application for similar semi-enclosed bays, and support the sustainable development of the bay. three main coastal cities, Yueqing, and Yuhuan. We investigated the spatiotemporal variational patterns in nutrient concentration in the waters from 2013 to 2020 as well as the main driving factors. These results indicate that the intensity of coastal anthropogenic activity


Introduction
With the rapid development of the coastal economy and continuous urbanization, over 60% of the world's population live within 200 km of a coastline [1], and more than 40% of China's population are centered in coastal areas [2]. As such, the impact of anthropogenic activity on the regional environment has also increased, natural resources are exploited on an ever-larger scale, and natural coastal areas and marine ecosystems are under pressure from intensive economic development and coastal zone projects [3][4][5]. At present,~41% of sea areas worldwide, especially estuaries and bays, have been seriously disturbed by anthropogenic activities [3]. As a result of ever-increasing urbanization, continuous mariculture scale-up, and rapid development of the ocean economy, land-derived pollutants (e.g., nitrogen and phosphorus) are continuously entering into waters, resulting in water quality deterioration and eutrophication in local areas. In particular, as typical ecologically sensitive areas heavily disturbed by land-sea interactions and anthropogenic activities, semi-enclosed and enclosed bays exhibit weaker water exchange and purification capacities, and thus are highly vulnerable environments [6][7][8]. As reported by the Ministry of Ecology and Environment of the People's Republic of China, 9.4% of nearshore seawaters in China are worse than grade IV in quality, with particularly high levels of dissolved inorganic nitrogen (DIN) and reactive phosphate (PO 4 -P) [9]. Therefore, to monitor variations in water quality of nearshore seas and bays, it is essential to conduct dynamic monitoring of nitrogen and phosphorus elements in water and investigate their driving factors, and then implement corresponding measures and policies to control and improve water quality status.
The accumulation of observation data provides an important basis for regional water quality monitoring. Based on measured DIN, PO 4 -P, and chlorophyll data from Chesapeake Bay, Harding et al. reported that DIN concentration doubled from 1945 to 1980, leading to an increase in chlorophyll concentration and serious ecosystem damage, but DIN concentration and eutrophication decreased somewhat from 1981 to 2012 due to the control of freshwater influx and DIN concentration [10]. Based on annual and monthly field measurement data of dissolved inorganic carbon (DIC), dissolved organic carbon (DOC), chlorophyll, and nutrients, Atsushi et al. investigated long-term trends in water quality status in Tokyo Bay, and showed that total input of land-derived pollutants has decreased to some extent due to gradual improvement in coastal wastewater treatment measures in recent years [11,12]. Based on field measurement data at monitoring stations, hydrological characteristics, and status of land-derived sewage, Cao et al. investigated nearshore marine pollution in Hangzhou Bay and found that serious nutrient pollution was primarily affected by topography, hydrological environment, and industrial development in the bay area [13]. The addition of land-derived nutrients, eutrophication, and other issues have led to water quality deterioration in bays heavily affected by anthropogenic activities, prompting governments to take corresponding measures, including strict control of the total amount of pollutants discharged into the sea, strengthening the pollutants treatment, and restricting the human activities in the bay [14,15], as well as protecting the natural shoreline, restoring the ecological environment of mudflats and wetlands, promoting the ecological functions of the bay, etc. [15][16][17]. To evaluate remediation performance and assess and modify existing management measures, continuous monitoring data are still needed, while the large-scale and long-term simultaneous observation of water quality data in major bays is generally insufficient. Satellite remote sensing technology has the advantages of long time series, large-scale simultaneous observation, and historical backtracking, which can overcome the shortcomings of in situ surveys and enable the observation of spatiotemporal variations in water quality.
Water quality parameters with optical activity, such as chlorophyll, colored dissolved organic matter, total suspended matters and turbidity, can be retrieved by the satellite monitoring spectrum [18][19][20][21][22][23][24] and are widely used for the satellite-based monitoring of water quality [25][26][27][28]. However, nutrients are not directly linked to the spectral signal and the variation mechanism is extremely complex, making nutrient retrieval a difficult Remote Sens. 2022, 14, 550 3 of 26 task [29]. Although some studies had attempted to retrieve nutrients by optically active variables [30][31][32], these algorithms were extremely regionally applicable and depended heavily on field measurement data, making them difficult to apply to other regions. Nevertheless, remote sensing is still an effective tool for water quality monitoring, with satellites such as the Landsat series and MODIS having been applied quite successfully to this aspect [33][34][35][36]. Therefore, the establishment of nutrient retrieval models based on field measured data can help to clarify the spatiotemporal change patterns of water quality in a region and explore the influencing factors, thereby providing theoretical and scientific support for water and ecological environment governance decisions.
Yueqing Bay is a typical semi-enclosed tidal bay and the fourth largest in area in China. Its tidal range, tidal influx, water exchange capacity, and other hydrological properties are representative of semi-enclosed tidal bays, to some extent. Moreover, Yueqing Bay has complex water exchange processes, including exchange between incoming offshore waters and internal bay water and between bay mouth waters and diluted waters from the Ou River. From the outlet to the bay head, internal waters vary considerably in water exchange capacity, as it takes less than 5 days for 90% of the waters near the outlet to be replaced by offshore waters, compared to 40 days for 90% of the waters in the whole bay to be replaced [37,38]. In addition, more than 30 streams crisscross along Yueqing Bay, and the bay mouth is located on the north side of the Ou River estuary, one of the eight major river systems in Zhejiang, China. Hence, given these natural geographic conditions, Yueqing Bay is greatly affected by land area change. In recent decades, cities along Yueqing Bay have witnessed rapid economic development, dense population growth, and surging demand for urban land for industries and agriculture. Human demand for ocean space development has increased significantly, and since the 1950s, reclaimed tidal flats along Yueqing Bay have increased to 165.1 km 2 , accounting for 38% of total tidal flat area and 18% of total sea area in the bay [39]. Along with accelerated urbanization, intensified anthropogenic activities, and increased domestic sewage and industrial and agricultural wastewater discharge, total input of land-derived pollutants into the bay has increased, further affecting water quality. Under highly intensive anthropogenic activities, the environment along Yueqing Bay has experienced rapid and dramatic deterioration. Since 2006, local governments have undertaken various measures for marine environment pollution remediation and controlled the use of phosphorus and land-derived pollutants. Based on nutrient data observed in May and August from 2000 to 2014, Guo et al. studied interannual variations in nutrient salt concentrations in the fishing waters of Yueqing Bay and found that concentrations increased continuously from 2002 to 2006 due to high river runoff and industrial and agricultural wastewater and domestic sewage, but that the rising trend slowed slightly after 2006 due to the implementation of national policies and measures to control pollutant discharge [40]. Given the intense anthropogenic activities, in-depth study of variations in environmental status and the main driving factors of Yueqing Bay could offer effective reference for environmental remediation of this and other semi-enclosed bays worldwide and may be of importance to the rapidly developing global bay economy.
Based on water quality data of Yueqing Bay collected from 2014 to 2020 (provided by Wenzhou Environmental Monitoring Center) and Landsat-8 remote sensing data, we built a nutrient concentration retrieval model for Yueqing Bay using support vector machine (SVM) learning with water temperature, red-light waveband, near-infrared waveband, and short infrared waveband as input parameters, and then, spatiotemporal variations in nutrient concentration in the bay retrieved from remote sensing data (2013-2020) were studied along with influencing factors, such as hydrodynamic and environmental information, river basinrelated management policies, and anthropogenic activities. This study can provide an effective reference for the remote sensing dynamic monitoring of water quality and provide technical support for the management of environmental protection and remediation of similar bays. Situated in the subtropical monsoon climate zone, Yueqing Bay is subject to various environmental factors, such as topography and marine climate, resulting in a land-bay microclimate with good heat, moisture, and temperature, resulting in a mean annual air temperature of 17-17.5 • C and mean annual rainfall of 1191.7-1506.8 mm [44]. Yueqing Bay contains 54,527.9 acres of tidal flats, and over 20 species of commercial fish, 58 species of shellfish, and 60 species of crustaceans, and is thus a natural farm of shellfish breeding. Furthermore, with a mean tidal range of more than 4 m (maximum of 8.53 m) [41], Yueqing Bay has considerable tidal energy, which is useful for tidal power plants, e.g., the Jiangxia pilot tidal power plant has an installed capacity of up to 550,000 kW [39,42].
Our study area included Yueqing Bay and its three main coastal cities, namely, Yueqing, Wenling, and Yuhuan. We investigated the spatiotemporal variational patterns in nutrient concentration in the waters from 2013 to 2020 as well as the main driving factors.

Satellite Remote Sensing Data
Land satellite Landsat-8 Operational Land Imager (OLI) high-spatial-resolution (30 m) multi-spectral remote sensing data (released by United States Geological Survey (USGS)) were used for analysis. Since its launch in 2013, data from the OLI has been widely applied in long-term water quality monitoring for inland rivers, lakes, nearshore seas, and bays. The Landsat-8 OLI remote sensing data were downloaded from the USGS website (https: //earthexplorer.usgs.gov/, accessed on 30 August 2021), and remote sensing images of the study area (strip no.: 118040) during 2013-2020 with sky cover less than 30% were selected. In total, 63 Landsat Level-1 remote sensing scenes met the above criteria (Table 1), including more than 10 remote sensing images per season and more than five remote sensing images per year covering most months, seasons, and years. To obtain the remote sensing reflectance product, we performed radiometric calibration and atmospheric correction of the data using ENVI v5.3 (Exelis Visual Information Solutions (Boulder, CO, USA)), with atmospheric correction conducted using the FLAASH algorithm in ENVI v5.3 [45]. Table 1. Valid Landsat-8 remote sensing images from January 2013 to December 2020 (green cells and numbers represent presence of data and quantity of valid data, respectively, while white cells represent absence of data).

Field Measurement Data
For this study, environmental water quality survey data of Yueqing Bay were provided by Wenzhou Environmental Monitoring Center. During 2014-2020, monthly voyage observations were performed, and in situ sampling and observation and data acquisition were conducted at nine fixed stations and several non-fixed stations. Monitoring parameters included DIN (nitrite-nitrogen, nitrate-nitrogen, and ammonia-nitrogen), PO 4 -P, water temperature, chlorophyll, dissolved oxygen, salinity, and pH. Both sampling and testing were performed in accordance with the Code of Practice for Marine Monitoring Technology [46]. During 2014-2020, the in situ sampling period and statistical sampling size were relatively stable and more than 45 pairs of observed DIN and PO 4 -P data per month were retained, representing primary characteristics of nutrients in Yueqing Bay.

Matching Satellite and Field Measurement Data
In this study, satellite data were matched with field measurement data. Spatially, 1 × 1 pixels (30 m × 30 m) closest to a survey station in latitude and longitude were matched; temporally, satellite data differing from field measurement data in observation time by less than ±3 days were matched, which was an optimal choice to obtain enough representative data in the shortest possible time interval. We could not obtain enough matchup data for the satellite algorithm development if we set up a shorter time interval, Remote Sens. 2022, 14, 550 6 of 26 due to the long revisiting period of 16 days for Landsat and the limited number of field measurements. To extract reflectance and other information, the "ExtractMultiValuesTo-Points" feature of ArcGIS v10.7 (Environmental Systems Research Institute, Inc. (Redlands, CA, USA)) was used to generate a vector file based on latitude and longitude information of a data observation station, with corresponding remote sensing reflectance information then extracted.
Using the above matching approach for data during 2014-2020, 80 matched DIN data pairs and 73 matched PO 4 -P data pairs were generated. We randomly sampled 65% of the data as a training dataset for algorithm construction, while the remaining 35% of data served as a validation dataset for algorithm accuracy assessment. In general, the matched data encompassed and reflected spatiotemporal variability of the field measurement data, and thus could be used to construct and assess the nutrient retrieval algorithm for waters of Yueqing Bay over a long time series.

SVM Algorithm and Statistical Parameters of Model Assessment
Machine learning-based regression algorithms, e.g., random forest, artificial neural network, and SVM algorithms, have been widely applied in water quality research [47][48][49]. In this study, based on matched data pairs, we used the SVM algorithm to construct a nutrient retrieval algorithm for waters of Yueqing Bay. This was mainly due to the smaller data size, as the random forest and artificial neural network algorithms are prone to overfitting. In contrast, the SVM algorithm attempts to map training data points to high-dimensional eigenspace based on kernel function, resulting in many decision-making boundaries (hyperplanes) in eigenspace that can separate different classes of data points. In addition, the SVM algorithm does not involve probability measures or the law of large numbers; moreover, the purpose of SVM optimization is to minimize structuring risk rather than empirical risk, thereby reducing the risk of overfitting [50]. Regarding prediction, the training process generated an SVM model with multiple input parameters for automatic optimization. Based on the trained SVM model, we predicted dependent variables (e.g., DIN and PO 4 -P) by inputting new training parameters (e.g., remote sensing reflectance).
The fitting effect during algorithm construction was expressed by coefficients of determination (R 2 ) and significance level (P) resulting from statistical analysis. In addition, to assess the accuracy of the model-predicted results quantitatively, the following statistical parameters were used: (1) bias, i.e., general bias information of the predicted result; (2) Mean Absolute Error (MAE), i.e., mean value of absolute error between predicted and true values; and (3) Root Mean Squared Error (RMSE), i.e., square root of the ratio of squared bias between observed and true values to number of observations. Results were statistically analyzed using Microsoft Excel 2016 (Microsoft (Redmond, Washington, DC, USA)): where DI N(PO 4 − P) i RS and DI N(PO 4 − P) i in-situ are the predicted and observed results, respectively, and N is the number of data pairs, i ∈ [1, N].

Hydrodynamic Model for Yueqing Bay
The Finite Volume Coastal Ocean Model (FVCOM) was used to simulate threedimensional (3D) hydrodynamics of Yueqing Bay. FVCOM is a 3D model for nearshores, estuaries, and oceans based on the finite volume method [51,52]. Unstructured triangular meshes were used in the horizontal direction and the σ coordinate system was used in the vertical direction to better resolve the complex shorelines and topography of Yueqing Bay.
This model was previously validated and applied by our research team in the Changjiang River estuary [53].
The model consisted of 14,197  , and background flow field data from HYCOM, among which the predicted tidal levels had a temporal resolution of 1 h and the background flow field data had a temporal resolution of 1 month. On the river bounds, climatological mean discharges of the Yangtze, Qiantang, Tiao, Yong, Jiao, Ou, Feiyun, and Ao rivers, and three major mountain creek-based riverlets, flowing into Yueqing Bay were superposed. Discharge data for the Yangtze and Qiantang rivers were sourced from Bulletin of River Sediments in China for 2002-2020, while discharge data for other rivers were estimated from the TRMM 3B43 climatological mean precipitation product for 1998-2019 using the river discharge estimation formula of Liu et al. (2020) [54]. For the sea-air upper bounds, wind data were sourced from the European Centre for Medium-Range Weather Forecasts (ECMWF), with a spatial resolution of 0.125 • × 0.125 • and a temporal resolution of 6 h.

Characterization of Intensity of Anthropogenic Activities in Cities along Yueqing Bay
Based on the administrative boundaries of Yueqing, Wenling, and Yuhuan along Yueqing Bay, point-of-interest (POI) data of the various cities in 2020 were acquired using the "POI retrieval" feature in Baidu Maps (http://lbsyun.baidu.com/index.php?title= webapi accessed on 30 August 2021), which are a type of vector point data containing latitude and longitude coordinates, facility name, and facility type [55]. Each vector point in space represents one facility space in real life and corresponds to commercial and human socioeconomical activities, such as traveling and dining, and governmental services. The POI data are applicable for analysis of anthropogenic activities, economic development status, and distribution from a spatial perspective [56,57].
Since many categories of facilities are closely associated with each other in functional properties, POI data are often manifested in aggregated form in geographic space, i.e., hot spots [58]. Usually, "density" is used to characterize aggregation intensity of POI data in a local area quantitatively, while kernel density estimation is a common and effective density estimation method [59] for nonparametric analysis of geospatial distribution, reflecting the effect of a geographic phenomenon attenuating with increasing spatial distance. Given an arbitrary point in space, kernel density is the sum of weights of all sample points within a circle with a given radius whose center is at the point. Such a weight is inversely proportional to the distance between the sample point and point [60,61]. Here, kernel density was calculated using the following formula: where f (x) represents the kernel density function of point x; h is the bandwidth, i.e., radius of the circle; n is the number of sample points within the circle; K is the spatial weight function; and (x − x 1 ) is the distance between the sample point and point x.
Given the above understanding, intensities of anthropogenic activities in cities along the Yueqing Bay were spatially analyzed based on POI data by means of kernel density estimation.

Algorithm Construction and Validation
Nutrients are non-optically active substances, and their concentrations have no significant spectral features [62]. However, nutrients are to some extent related to the physical properties of water, such as water temperature, and optically active substances in water, such as suspended solids and chlorophyll-a [63][64][65], which are associated with spectra [66][67][68][69][70]. In this study, waveband information was used as a variable input to retrieve nutrient concentrations in the waters of Yueqing Bay.
Correlation analysis was carried out for temperature and waveband information. As shown in Figure 2, DIN was highly correlated with sea surface temperature (SST) (R = −0.84), as well as the red-light (Band 4), near-infrared (Band 5) (R > 0.1), and short infrared wavebands (Bands 6 and 7) (R < −0.1); PO 4 -P was highly correlated with SST (R = −0.70), as well as the red-light (Band 4) and near-infrared wavebands (Band 5) (R > 0.15). Therefore, SST and the red-light, near-infrared, and two short infrared wavebands were selected as input parameters of the DIN retrieval model, and SST, red-light waveband, and near-infrared waveband were selected as input parameters of the PO 4 -P retrieval model. Given the above understanding, intensities of anthropogenic activities in cities along the Yueqing Bay were spatially analyzed based on POI data by means of kernel density estimation.

Algorithm Construction and Validation
Nutrients are non-optically active substances, and their concentrations have no significant spectral features [62]. However, nutrients are to some extent related to the physical properties of water, such as water temperature, and optically active substances in water, such as suspended solids and chlorophyll-a [63][64][65], which are associated with spectra [66][67][68][69][70]. In this study, waveband information was used as a variable input to retrieve nutrient concentrations in the waters of Yueqing Bay.
Correlation analysis was carried out for temperature and waveband information. As shown in Figure 2, DIN was highly correlated with sea surface temperature (SST) (R = −0.84), as well as the red-light (Band 4), near-infrared (Band 5) (R > 0.1), and short infrared wavebands (Bands 6 and 7) (R < −0.1); PO4-P was highly correlated with SST (R = −0.70), as well as the red-light (Band 4) and near-infrared wavebands (Band 5) (R > 0.15). Therefore, SST and the red-light, near-infrared, and two short infrared wavebands were selected as input parameters of the DIN retrieval model, and SST, red-light waveband, and nearinfrared waveband were selected as input parameters of the PO4-P retrieval model. Based on the data arrays acquired by matching, a nutrient retrieval model for Yueqing Bay was obtained by SVM. As shown in Figure 3a  Based on the data arrays acquired by matching, a nutrient retrieval model for Yueqing Bay was obtained by SVM. As shown in Figure 3a

Spatial Distribution of Nutrient Concentrations in Yueqing Bay
From the perspective of the spatial distribution of nutrient concentration based on remote sensing retrieval, in spring, there were two high DIN regions in the narrowest area of the middle part of the bay and Southwest Shoal at the bay mouth, while DIN concentrations were more uniform in most other areas. In summer, high DIN values remained in the narrowest area of the middle part of the bay and Southwest Shoal, while at the bay head, DIN concentrations on the west side of Ximen Island were slightly higher than those on the east side. In autumn, high DIN concentrations were mainly observed in the stream on the east side of Ximen Island at the bay head, Southwest Shoal at the bay mouth, and the sea area on the north side of Damen Island, while low concentrations were mainly observed in the narrowest area of the middle part of the bay and coastal area to the southwest side of Yuhuan Island. In winter, high DIN concentrations were observed in the stream on the east side of Ximen Island at the bay head, Southwest Shoal at the bay mouth, and sea areas on the north side of Damen and Xiaomen islands, while low concentrations were observed in the narrowest area of the middle part of the bay through the sea area in the middle of the bay mouth (Figure 4a-d). In spring, there was a high PO4-P concentration zone in the Southwest Shoal thr the sea area in the middle of the Yueqing Bay mouth, while low concentrations wer served in the sea areas near Maoyan Island and the southwest side of Yuhuan Islan summer, high PO4-P concentrations were observed in the Southwest Shoal, while concentrations were observed at the bay head and sea area between Yuhuan Island Damen Island. In autumn, a high PO4-P concentration zone was observed in the m part of the bay through the sea area on the north side of Damen Island, while low co trations were mainly observed in the coastal area of the Southwest Shoal. In winter, PO4-P concentrations were observed in the Southwest Shoal and sea area on the sides of Damen and Xiaomen islands, while low concentrations were observed i stream on the west side of Ximen Island at the bay head and in the narrowest area o middle part of the bay (Figure 5a-d).
Based on field measurement data from 2014 to 2020, the spatial distribution of n ent concentrations in Yueqing Bay was obtained by Kriging interpolation (Figures  and 5e-h), and a comparison was made with the remote sensing retrieval results (Fi 4a-d and 5a-d). In general, the spatial distribution of the satellite-derived nutrient centration was consistent with that of the field measured data; high-value zones we uated in the streams on both sides of Ximen Island at the bay head, narrowest area o middle part of the bay, and Southwest Shoal at the bay mouth, while low-value z were situated in the sea area on the west side of Yuhuan Island. Due to the limited pling data and coverage, there were no data in the nearshore and the distribution of In spring, there was a high PO 4 -P concentration zone in the Southwest Shoal through the sea area in the middle of the Yueqing Bay mouth, while low concentrations were observed in the sea areas near Maoyan Island and the southwest side of Yuhuan Island. In summer, high PO 4 -P concentrations were observed in the Southwest Shoal, while low concentrations were observed at the bay head and sea area between Yuhuan Island and Damen Island. In autumn, a high PO 4 -P concentration zone was observed in the middle part of the bay through the sea area on the north side of Damen Island, while low concentrations were mainly observed in the coastal area of the Southwest Shoal. In winter, high PO 4 -P concentrations were observed in the Southwest Shoal and sea area on the north sides of Damen and Xiaomen islands, while low concentrations were observed in the stream on the west side of Ximen Island at the bay head and in the narrowest area of the middle part of the bay (Figure 5a-d).
Based on field measurement data from 2014 to 2020, the spatial distribution of nutrient concentrations in Yueqing Bay was obtained by Kriging interpolation (Figures 4e-h and 5e-h), and a comparison was made with the remote sensing retrieval results (Figures 4a-d and 5a-d). In general, the spatial distribution of the satellite-derived nutrient concentration was consistent with that of the field measured data; high-value zones were situated in the streams on both sides of Ximen Island at the bay head, narrowest area of the middle part of the bay, and Southwest Shoal at the bay mouth, while low-value zones were situated in the sea area on the west side of Yuhuan Island. Due to the limited sampling data and coverage, there were no data in the nearshore and the distribution of interpolation results might be unreasonable in the nearshore. As shown in Figure 5f, there was an obvious mutation in the concentration of PO 4 -P from the middle of the bay to the bay mouth, which might be the artificial error by interpolation. The spatial distribution of satellite-derived nutrients was more complete in spatial coverage. The distribution of nutrient concentrations reported by Zeng et al. [71] in Yu Bay from 2006 to 2007 was basically consistent with our satellite retrieval results, th the high nutrient concentrations were mainly concentrated in the Southwest Shoa the narrowest part of the middle part of the bay, and the low concentrations were m concentrated in the west side of Yuhuan Island. Through model simulation, Li et a found that the nutrient concentrations in the western part of the bay were higher those in the eastern part. However, Li et al. [72] reported that high nutrient concentra were observed in the bay head and low nutrient concentrations were observed in th mouth, and there were no obvious high nutrient concentrations in the narrowest ar the middle part of the bay, which was slightly different from the results of this study main reason would be the different input parameters of the nutrient concentratio trieval model. The water quality model constructed by Li et al. [72] mainly considere influence of water flow and diffusion, while we used the temperature and spectral mation as the main parameters to retrieve the nutrient concentrations. At the same on both sides, the narrowest part in the middle part of the bay might be affected by power plants, the Huaneng Yuhuan Power Plant (put into operation in 2007) an Zheneng Yueqing Power Plant (put into operation in 2008). Although the influence r of warm water drainage of the power plant was limited, it may still have a certain im on the satellite-derived nutrient concentrations in the narrowest part of the middle p the bay. Therefore, more data are needed for refined verification and analysis. The distribution of nutrient concentrations reported by Zeng et al. [71] in Yueqing Bay from 2006 to 2007 was basically consistent with our satellite retrieval results, that is, the high nutrient concentrations were mainly concentrated in the Southwest Shoal and the narrowest part of the middle part of the bay, and the low concentrations were mainly concentrated in the west side of Yuhuan Island. Through model simulation, Li et al. [72] found that the nutrient concentrations in the western part of the bay were higher than those in the eastern part. However, Li et al. [72] reported that high nutrient concentrations were observed in the bay head and low nutrient concentrations were observed in the bay mouth, and there were no obvious high nutrient concentrations in the narrowest area of the middle part of the bay, which was slightly different from the results of this study. The main reason would be the different input parameters of the nutrient concentration retrieval model. The water quality model constructed by Li et al. [72] mainly considered the influence of water flow and diffusion, while we used the temperature and spectral information as the main parameters to retrieve the nutrient concentrations. At the same time, on both sides, the narrowest part in the middle part of the bay might be affected by two power plants, the Huaneng Yuhuan Power Plant (put into operation in 2007) and the Zheneng Yueqing Power Plant (put into operation in 2008). Although the influence range of warm water drainage of the power plant was limited, it may still have a certain impact on the satellite-derived nutrient concentrations in the narrowest part of the middle part of the bay. Therefore, more data are needed for refined verification and analysis.

Seasonal Variations in Nutrient Concentrations in Yueqing Bay
Based on analysis of the remote sensing retrieval results (Table 2)  Based on the field measurement data ( Figure 6) in Yueqing Bay, DIN concentrations were lower in summer (mean: 0.482 mg/L) and higher in spring and winter (means: 0.73 mg/L and 0.66 mg/L, respectively); PO 4 -P concentrations were lower in summer (mean: 0.043 mg/L) and higher in autumn and winter (means: 0.047 mg/L and 0.048 mg/L, respectively). Nutrient concentrations varied greatly, with the DIN concentration ranging from 0.2 to 1.0 mg/L and PO 4 -P concentration ranging from 0.01 to 0.08 mg/L.  Based on the field measurement data ( Figure 6) in Yueqing Bay, DIN concentrations were lower in summer (mean: 0.482 mg/L) and higher in spring and winter (means: 0.73 mg/L and 0.66 mg/L, respectively); PO4-P concentrations were lower in summer (mean: 0.043 mg/L) and higher in autumn and winter (means: 0.047 mg/L and 0.048 mg/L, respectively). Nutrient concentrations varied greatly, with the DIN concentration ranging from 0.2 to 1.0 mg/L and PO4-P concentration ranging from 0.01 to 0.08 mg/L.

Long Time-Series Variations in Nutrient Concentrations in Yueqing Bay
As shown in Figure 7a, the DIN concentration in Yueqing Bay exhibited a notable rising trend (P = 0.0262, t = 2.5341) during 2013-2015 and a non-significant decreasing trend during 2015-2020 (P = 0.1233, t = −1.5729). As shown in Figure 7b, the PO 4 -P concentration in Yueqing Bay exhibited a non-significant increasing trend (P = 0.3196, t = 1.0383) during 2013-2015 and a non-significant decreasing trend (P = 0.8906, t = −0.1383) during 2015-2020. Comparing field measurement data with the remote sensing retrieval results ( Figure  7), the long time-series variation trends in remote sensing retrieval generally agreed with the observed results, but due to the acquisition time of satellite remote sensing images and other objective factors, the field measurement data of high DIN concentration could not be matched during the data matching process, resulting in the DIN retrieval results being slightly lower than the field measurement data in some periods. At the same time, abnormally low nutrient concentrations were found in Yueqing Bay in the summer of 2017. To further elucidate the water quality trends in local regions of Yueqing Bay under strong anthropogenic impact, three 500 m × 500 m regions located at coastal stream estuaries were selected (Figure 8a), i.e., East Stream of the bay head on the east side of Ximen Island (#1), West Stream of the bay head on west side of Ximen Island and at Qing River estuary (#2), and Southwest Shoal at the bay mouth (#3), where mean, maximum, and minimum nutrient concentrations were extracted. As shown in Figure 8b-g, during 2013-2015, the increasing nutrient concentrations in these three regions were all higher than the mean nutrient concentration rate in Yueqing Bay (see slopes in Table 3). Among these results, the PO4-P concentrations in the East Stream (#1) and West Stream of the bay head (#2) increased at the fastest rate (Figure 8c,e), aPO4-P(#1) = 1.98 × 10 −5 , aPO4-P(#2) = 2.03 × 10 −5 ). After 2015, nutrient Comparing field measurement data with the remote sensing retrieval results (Figure 7), the long time-series variation trends in remote sensing retrieval generally agreed with the observed results, but due to the acquisition time of satellite remote sensing images and other objective factors, the field measurement data of high DIN concentration could not be matched during the data matching process, resulting in the DIN retrieval results being slightly lower than the field measurement data in some periods. At the same time, abnormally low nutrient concentrations were found in Yueqing Bay in the summer of 2017.
To further elucidate the water quality trends in local regions of Yueqing Bay under strong anthropogenic impact, three 500 m × 500 m regions located at coastal stream estuaries were selected (Figure 8a), i.e., East Stream of the bay head on the east side of Ximen Island (#1), West Stream of the bay head on west side of Ximen Island and at Qing River estuary (#2), and Southwest Shoal at the bay mouth (#3), where mean, maximum, and minimum nutrient concentrations were extracted. As shown in Figure 8b-g, during 2013-2015, the increasing nutrient concentrations in these three regions were all higher than the mean nutrient concentration rate in Yueqing Bay (see slopes in Table 3). Among these results, the PO 4 -P concentrations in the East Stream (#1) and West Stream of the bay head (#2) increased at the fastest rate (Figure 8c,e), a PO4-P(#1) = 1.98 × 10 −5 , a PO4-P(#2) = 2.03 × 10 −5 ). After 2015, nutrient levels in the waters of the three regions decreased to variable extents. Notably, the nutrient concentrations in the East Stream (#1) decreased at a slower rate than Yueqing Bay overall (Figure 8b,e), especially DIN (a DIN(#1) = −1.30 × 10 −5 ), while the decreasing nutrient concentrations in the West Stream (#2) and Southwest Shoal (#3) were higher than the mean decreasing rate in Yueqing Bay (Figure 8d-g), especially PO 4 -P in Southwest Shoal (#3) (a PO4-P(#3) = −2.90 × 10 −6 ). Table 3 shows the statistical parameters of the time-series regression of nutrient concentrations in the waters of various regions and Yueqing Bay.  Table 3 shows the statistical parameters of the time-series regression of nutrient concentrations in the waters of various regions and Yueqing Bay.

Influences of Rivers, Hydrodynamics, and Biological Effects
Yueqing Bay is the largest semi-enclosed bay in Zhejiang and is influenced by the offshore Fujian-Zhejiang Coastal Current and Taiwan Warm Current [73,74]. Although there is no large river along the bay providing land-derived input, Yueqing Bay is on the north side of the Ou River estuary and more than 30 streams are distributed in the bay area [42]. Based on HYCOM background flow field, runoff, wind farm, and tidal data, the climatological flow field of nearshore sea areas in the middle and southern parts of Zhejiang and Yueqing Bay in summer and winter was acquired by FVCOM numerical simulation.
Based on variations in flow fields in the nearshore sea areas in summer ( Figure 9a) and winter (Figure 9c), the high-temperature low-nutrient Taiwan Warm Current flowing from south to north was the dominant water mass in summer [74], while in winter, the Taiwan Warm Current weakened significantly and retreated towards the south, and the low-temperature, nutrient-rich Fujian-Zhejiang Coastal Current flowing from north to south gradually became the dominant water mass, with a flow field intensity much higher than that of the Taiwan Warm Current [73].
Waters in Yueqing Bay are affected by streams inside the bay and water from the Ou River ( Figure 9). These waters are mainly manifested by bay head water flowing past the bay center and turning to the west side of Yuhuan Island, then flowing out; by waters in the Southwest Shoal flowing to the east and converging with water from the bay head; and by a small amount of water from the Ou River flowing past the north side of Damen Island and Xiaomen Island, draining likewise into the water stream from the bay head, and finally flowing out of the channel between Luxi Island and Damen Island. As seen in Figures 4a-d and 5a-d, the spatial distribution and transport direction of nutrient concentrations in Yueqing Bay remained generally consistent with the above results, indicating that river transport has a great impact on the spatial distribution of nutrients in the bay. In addition, Yueqing Bay is affected by multiple influencing factors, including inflowing streams, the Ou River, tides, and flow fields, and thus the bay mouth zone is hydrodynamically active, with frequent resuspension activity [75]. Sediment resuspension results in the redistribution of nutrients at the sediment-water interface, which can have an obvious impact on migration, transformation, and storage of nitrogen and phosphorus [76,77]. Thus, the complex water environment of Yueqing Bay is subjected to many factors.
In addition, observed chlorophyll concentration was used to characterize biomass and photosynthesis intensity and their relationship to nutrient concentration. The results show that, in summer (June to August), chlorophyll concentration increased significantly, resulting in high primary productivity, while in autumn and winter, chlorophyll concentration started to decline significantly (Figure 10a), resulting in lower primary productivity; correspondingly, in summer, nutrient concentration in Yueqing Bay decreased markedly until autumn and winter when the nutrient concentration began to increase again. Based on multi-linear curve fitting, we fit the climatological mean values of DIN, PO 4 -P, and chlorophyll concentration during 2014-2020. As shown in Figure 10b, a significant inverse correlation existed between chlorophyll concentration and DIN and between chlorophyll concentration and PO 4 -P (R 2 = 0.7192, P < 0.005), i.e., DIN and PO 4 -P concentrations decreased with increasing chlorophyll. Thus, seasonal variation in photosynthesis intensity has a certain impact on nutrient concentration in Yueqing Bay. In spring and summer, plankton photosynthesis was intense, and high levels of DIN and PO 4 -P were consumed, while in autumn and winter, plankton photosynthesis was weakened, thereby contributing to the accumulation of DIN and PO 4 -P.
Notes: Regression equation y = ax + b, where a is slope and b is intercept. Values in bold indicate significance P < 0.05.

Influences of Rivers, Hydrodynamics, and Biological Effects
Yueqing Bay is the largest semi-enclosed bay in Zhejiang and is influenced by the offshore Fujian-Zhejiang Coastal Current and Taiwan Warm Current [73,74]. Although there is no large river along the bay providing land-derived input, Yueqing Bay is on the north side of the Ou River estuary and more than 30 streams are distributed in the bay area [42]. Based on HYCOM background flow field, runoff, wind farm, and tidal data, the climatological flow field of nearshore sea areas in the middle and southern parts of Zhejiang and Yueqing Bay in summer and winter was acquired by FVCOM numerical simulation.
Based on variations in flow fields in the nearshore sea areas in summer ( Figure 9a) and winter (Figure 9c), the high-temperature low-nutrient Taiwan Warm Current flowing from south to north was the dominant water mass in summer [74], while in winter, the Taiwan Warm Current weakened significantly and retreated towards the south, and the low-temperature, nutrient-rich Fujian-Zhejiang Coastal Current flowing from north to south gradually became the dominant water mass, with a flow field intensity much higher than that of the Taiwan Warm Current [73]. (Black arrow represents flow field direction, length of black arrow represents flow field intensity; green, yellow, red, and purple arrows represent directions of water from bay head, Southwest Shoal, Ou River, and final outflowing stream, respectively).
uting to the accumulation of DIN and PO4-P.
According to the available data analysis, from the perspective of natural factors, physical effects, such as the Fujian-Zhejiang Coastal Current and Taiwan Warm Current, combined with biological effects, such as seasonal variation in the intensity of plankton photosynthesis in the bay have an important impact on the seasonal variation of nutrient concentration in Yueqing Bay, but more parameters are needed for detailed comprehensive analysis. Figure 10. Inter-monthly variations (a) and correlation analysis (b) of nutrient and chlorophyll concentrations in Yueqing Bay during 2013-2020. (In (b), abscissa represents sum of DIN concentration multiplied by a constant and PO 4 -P concentration multiplied by another constant after multi-linear curve fitting, which was used to characterize nutrient level in water. Ordinate represents chlorophyll concentration in water. Formula represents a relational formula for fitting nutrient and chlorophyll concentrations. As DIN and PO 4 -P concentrations are multiplied by negative constants, nutrient concentrations are inversely correlated to chlorophyll concentration).
According to the available data analysis, from the perspective of natural factors, physical effects, such as the Fujian-Zhejiang Coastal Current and Taiwan Warm Current, combined with biological effects, such as seasonal variation in the intensity of plankton photosynthesis in the bay have an important impact on the seasonal variation of nutrient concentration in Yueqing Bay, but more parameters are needed for detailed comprehensive analysis.

Influences of Intensity of Anthropogenic Activities in River Basin
Changes in the intensity of anthropogenic activities inside the river basin can have an immediate impact on nutrient level in land-derived water. To investigate the impact of anthropogenic activities and stream transport on the spatial distribution of nutrient concentration in the bay, the intensity of anthropogenic activities along Yueqing Bay and spatial distribution status were analyzed by kernel density using POI data. As shown in Figure 11a, six regions along Yueqing Bay had strong anthropogenic activities, namely, Yuecheng Town (#1), Hongqiao Town (#2), and Dajing Town (#3) in Yueqing City, Chengxi Subdistrict (#4) in Wenling, and Chumen Town (#5) and Yucheng Subdistrict (#6) in Yuhuan. Based on water system distribution, elevation, and other geographic information of the bay (Figure 11b,c), streams flowing past Yuecheng Town (#1), Hongqiao Town (#2), and Dajing Town (#3) eventually drain into Yueqing Bay. Dajing Town (#3) is nearest to Yueqing Bay, but anthropogenic intensity in this region was weaker than that in the other regions (low kernel density); streams passing the Chengxi Subdistrict (#4) eventually flow into Yueqing Bay, but they are farther away. Streams flowing past Chumen Town (#5) are situated in Xuanmen Bay, but water flow is blocked by dams inside Xuanmen Bay and thus water is not continuously transported to Yueqing Bay; streams flowing past Yucheng Subdistrict (#6) eventually drain into the sea area on the east side of Yuhuan Island without having an immediate impact on the waters in Yueqing Bay. Based on the spatial distribution of nutrient concentration in Yueqing Bay (Figure 11b), high-concentration zones were situated in shoals southwest of the bay mouth (corresponding to Yuecheng Town (#1) and Hongqiao Town (#2)), streams on both sides of Ximen Island at the bay head (corresponding to Dajing Town (#3) and Chengxi Subdistrict (#4)), and coastal sea areas beyond the Xuanmen Bay dams (corresponding to Chumen Town (#5)), whereas low-concentration zones were situated in the sea area on the west side of Yuhuan Island (corresponding to the Yucheng Subdistrict (#6)). These results indicate that the intensity of coastal anthropogenic activity and transport of inflowing streams are the main driving factors contributing to the high or low concentrations of nutrients in waters in local regions inside the bay.

Influences of Urban Management in River Basin and Changes in Bay Governance Policy
Since its development in 1980, Yueqing Bay has also experienced many problems similar to the world-famous bay areas (e.g., Tokyo Bay), such as increased pollutant discharge and ecosystem damage [11,12]. As shown in Figure 12, compared to 1980, nutrient concentrations in Yueqing Bay have increased about 5-7-fold. From 1988 to 2005, cities along Yueqing Bay have undertaken significant land reclamation in the bay to address urban land shortages, totaling 61.81 km 2 . As a continuation of the coastal reclamation plan of Xuanmen Bay in Yuhuan, Xuanmen phase II further strengthened the reclamation of land from the sea on the basis of the construction of Xuanmen dam in Xuanmen phase I, and the reclamation area reached 37.34 km 2 . The implementation of Xuanmen phase I reclamation project directly blocked the exchange of water in Xuanmen Port and waters inside and outside the bay, while the large-scale reclamation of land from the sea in Xuanmen phase II further resulted in a remarkable decline in tidal influx in the bay. During 1988-2005, tidal influx from the bay center to the bay head decreased by 10%-20% and the water exchange rate declined significantly. Nutrients in the bay failed to be transported to the open sea and thus accumulated inside the bay, resulting in a marked nutrient concentration increase during this period [72]. After 2005, sea reclamation, urban areas, and anthropogenic activities increased, with an accompanying increase in the total discharge of land-derived pollutants and their transportation along inflowing streams to Yueqing Bay, resulting in further enhancement of nutrient concentrations in Yueqing Bay [40,43]. By 2018, the Yueqing Bay River basin contained a population of about 3 million and small enterprises had expanded, with industrial wastewater discharge in the bay increasing accordingly (up to 3 million tons each year). In addition, the aquatic farming industry in the bay has continued to develop. Since 2011, annual gross output of aquatic products in the three cities along Yueqing Bay has exceeded 80 million tons [43,78]; urban domestic sewage, coastal sewage outfalls, and pollutants from aquatic farming have drained continuously into Yueqing Bay.

Influences of Urban Management in River Basin and Changes in Bay Governance Policy
Since its development in 1980, Yueqing Bay has also experienced many problems similar to the world-famous bay areas (e.g., Tokyo Bay), such as increased pollutant discharge and ecosystem damage [11,12]. As shown in Figure 12 [72]. After 2005, sea reclamation, urban areas, and anthropogenic activities increased, with an accompanying increase in the total discharge of land-derived pollutants and their transportation along inflowing streams to Yueqing Bay, resulting in further enhancement of nutrient concentrations in Yueqing Bay [40,43]. By 2018, the Yueqing Bay River basin contained a population of about 3 million and small enterprises had expanded, with industrial wastewater discharge in the bay increasing accordingly (up to 3 million tons each year). In addition, the aquatic farming industry in the bay has continued to develop. Since 2011, annual gross output of aquatic products in the three cities along Yueqing Bay has exceeded 80 million tons [43,78]; urban domestic sewage, coastal sewage outfalls, and pollutants from aquatic farming have drained continuously into Yueqing Bay. As seen in Figure 13, gross output values of the primary sectors and industrial in cities along Yueqing Bay have shown a growing tendency. During 2011-2020, gross output values of both sectors increased by about RMB 90 billion, while gross output of aquatic products generally increased by ~16 million tons over the seven years, decreased somewhat after 2017, and decreased to 86.5 million tons by 2020. Most inflowing streams along Yueqing Bay run past rural towns. Many are classified as grade IV-V surface water due to increasing input from industrial wastewater, agricultural irrigation run-off, domestic sewage, and aquatic farming pollutants. Gao et al. [79] reported that most inflowing streams in Yueqing Bay have a total nitrogen and phosphorus content of 1.5-6.0 mg/L and 0.06-0.2 mg/L, respectively. At present, the output value of various industries continues to grow, the total population in the basin continues to increase, and rural towns continue to develop. Although the influence of the current pollution control plan, which highlighted the control of total nitrogen and total phosphorus emissions from rivers entering the sea and direct discharge sea pollution sources, achieved results and controlled the nutrient concentration in surface water to a certain extent, land-based pollution is still one of the main sources of pollutants in Yueqing Bay at this stage. As seen in Figure 13, gross output values of the primary sectors and industrial in cities along Yueqing Bay have shown a growing tendency. During 2011-2020, gross output values of both sectors increased by about RMB 90 billion, while gross output of aquatic products generally increased by~16 million tons over the seven years, decreased somewhat after 2017, and decreased to 86.5 million tons by 2020. Most inflowing streams along Yueqing Bay run past rural towns. Many are classified as grade IV-V surface water due to increasing input from industrial wastewater, agricultural irrigation run-off, domestic sewage, and aquatic farming pollutants. Gao et al. [79] reported that most inflowing streams in Yueqing Bay have a total nitrogen and phosphorus content of 1.5-6.0 mg/L and 0.06-0.2 mg/L, respectively. At present, the output value of various industries continues to grow, the total population in the basin continues to increase, and rural towns continue to develop. Although the influence of the current pollution control plan, which highlighted the control of total nitrogen and total phosphorus emissions from rivers entering the sea and direct discharge sea pollution sources, achieved results and controlled the nutrient concentration in surface water to a certain extent, land-based pollution is still one of the main sources of pollutants in Yueqing Bay at this stage. Time-series variations in nutrient concentrations in Yueqing Bay during 2013-2020 ( Figure 5) showed abnormally low nutrient concentrations in the summer of 2017, which was mainly due to prolonged drought, with monthly mean rainfall lower than perennial monthly mean rainfall; indeed, only 2 days in August showed precipitation >5 mm and monthly total rainfall was only 51.2 mm in Yueqing City [82]. The drought led to a significant decline in water levels and river discharge and a lower total content of nutrient trans- Time-series variations in nutrient concentrations in Yueqing Bay during 2013-2020 ( Figure 5) showed abnormally low nutrient concentrations in the summer of 2017, which was mainly due to prolonged drought, with monthly mean rainfall lower than perennial monthly mean rainfall; indeed, only 2 days in August showed precipitation >5 mm and monthly total rainfall was only 51.2 mm in Yueqing City [82]. The drought led to a significant decline in water levels and river discharge and a lower total content of nutrient transported by streams into Yueqing Bay, resulting in abnormally low nutrient concentrations in the bay. This climatic event showed that industrial wastewater, agricultural irrigation water, domestic sewage, and aquatic farming pollutants in cities along Yueqing Bay are the main sources of nutrient pollution, and inflowing streams are the main mode of pollutant transport.
In May 2014, the "Five Water Treatment Actions" program (i.e., wastewater treatment, flood control, stagnant water drainage, water supply assurance, and water conservation) was implemented in Zhejiang for the purpose of water control, and mandatory water quality indicators have been used to control pollutant discharge and increase sustainable development [83]. In the context of this program, various ecological governance and remediation policies were proposed to address issues such as pollutant discharge control and environmental protection in cities along Yueqing Bay and to improve water quality in the bay by human intervention. In 2016, the Yueqing municipal government proposed policies to rationalize marine resource exploitation and utilization, control entry of polluting enterprises, and develop the ecological sector to control land-derived pollutants and pollutant discharge entering Yueqing Bay [84]. In 2021, the Yuhuan municipal government also proposed actions to govern and monitor pollutants discharged into the sea to improve water quality status in the bay by ecosystem remediation and pollutant control [85]. In addition to relevant pollution control and ecological remediation policies, after 2017, the development of the agricultural, forestry, husbandry, fishing, and aquatic product sectors along Yueqing Bay increased at somewhat slower rates, which had a positive impact on total input of land-derived pollutants and pollution control.
As seen in Figure 7, after 2015, the DIN and PO 4 -P concentrations in Yueqing Bay decreased slightly but non-significantly (P > 0.1), also with no further growth noted, indicating that the ecological governance and remediation policies proposed by the cities along Yueqing Bay have had a positive impact on water quality status in Yueqing Bay.
As seen in Figure 8, for stream estuaries corresponding to regions along Yueqing Bay with intense anthropogenic activity (i.e., Yuecheng Town (#1), Hongqiao Town (#2), Dajing Town (#3), and Chengxi Subdistrict (#4) in Figure 12), nutrient concentrations changed distinctly before and after 2015. During 2013-2015, nutrient concentrations in all three regions increased at higher rates than the mean increasing rate in Yueqing Bay. In the East Stream (#1 in Figure 8) and West Stream of the bay head (#2 in Figure 8), the PO 4 -P concentration increased at the fastest rate (Figure 8c,e and Table 3). After 2015, nutrient levels in the three regions declined to variable extents. In the East Stream (#1), nutrient concentrations in the water decreased at a slower rate, in particular, DIN decreased at a slightly lower rate than that in Yueqing Bay, while in the West Stream (#2 in Figure 8) and Southwest Shoal (#3 in Figure 8), nutrient concentrations decreased at a higher rate than that in Yueqing Bay, especially in the Southwest Shoal (#3 in Figure 8), where PO 4 -P decreased markedly. On the west side of the bay, remediation polies to control PO 4 -P were initiated relatively earlier by Yueqing city [40], resulting in significantly lower PO 4 -P concentrations on the west side of Yueqing Bay. After 2015, a series of ecological governance and remediation policies were issued and implemented by local city governments along Yueqing Bay. These measures slowed the rise in nutrient concentrations in the bay and mitigated the water quality deterioration trend, especially in Yueqing City, where PO 4 -P concentrations in inflowing streams were effectively controlled, resulting in a reduction in concentrations in the sea area on the west side of the bay. However, despite the mitigation effects, nutrient concentrations in Yueqing Bay remain high. For example, DIN concentrations only reach grade III seawater criteria (0.30 mg/L < DIN ≤ 0.40 mg/L) in summer, with levels higher than the upper limit specified for grade IV seawater (0.40 mg/L < DIN ≤ 0.50 mg/L) in the other seasons; furthermore, PO 4 -P concentrations exceed grade IV seawater criteria (0.30 mg/L < PO 4 -P ≤ 0.045 mg/L) all year [86].
The restoration of bay ecosystems is a long process; thus, continuous long-term ecological remediation efforts are essential for ecosystem restoration and regeneration. As an important semi-enclosed bay in the Kanto Area of Japan, Tokyo Bay is surrounded by several large cities, including Tokyo, Yokohama, and Kawasaki. During its development, coastal engineering, rapid urbanization, and other anthropogenic activities resulted in significant sea reclamation, pollutant discharge, ecosystem damage, and water quality deterioration in the bay [11,12]. In response, the Japanese government launched the "Tokyo Bay Restoration Plan" in 2001 to address these ecological problems, including the specific policies on wastewater discharge control, river dredging, wetland ecosystem restoration, and monitoring system establishment, etc. After years of governance and remediation, total pollutants (e.g., nitrogen and phosphorus) in Tokyo Bay have decreased, and water quality, while still poor, has not deteriorated further [87]. That is, although ecological regeneration policies have shown some success, the environment inside Tokyo Bay remains poor due to its early and high intensity development and deterioration. In comparison, Yueqing Bay developed much later, and ecological remediation has been implemented gradually over the past few years. Thus, from the perspective of the developmental and remediation history of Tokyo Bay, long-term ecological remediation efforts are still needed to improve the environment in Yueqing Bay and to eliminate or reduce the negative impacts of coast engineering and other anthropogenic activities.

Conclusions
This paper focused on water quality monitoring and variation analysis of a semienclosed bay (Yueqing Bay) in China, which is under strong anthropogenic impact. By constructing a regional remote sensing retrieval algorithm, the long time series remote sensing monitoring of water quality changes in Yueqing Bay is realized. On this basis, the impact of relevant management policies and anthropogenic activities in the basin on water quality variations in the bay can be understood, which can provide an effective reference basis for ecological environmental protection and restoration work in Yueqing Bay and other similar bays.
Based on long time-series Landsat-8 satellite remote sensing data, field measurement data, and supplementary statistical data, we constructed a model to retrieve nutrient concentrations in Yueqing Bay using SVM, and further investigated the spatiotemporal variations and main influencing mechanisms of nutrient concentrations in the bay during 2013-2020.
The nutrient retrieval results show that DIN concentration varied seasonally within 0.211-0.687 mg/L, with the lowest level in summer (averaging 0.3940 mg/L); PO 4 -P concentration varied seasonally within 0.016-0.056 mg/L, with the lowest level in summer (averaging 0.0414 mg/L) and highest level in autumn (averaging 0.0508 mg/L). From a spatial distribution perspective, nutrient concentrations in Yueqing Bay were generally higher in the west and lower in the east, with high-concentration zones observed in the narrowest area of the middle part of the bay, Southwest Shoal at the bay mouth, and sea area on the north side of Damen and Xiaomen islands, and low-concentration zones observed in coastal waters on the west side of Yuhuan Island. During 2013-2015, the DIN and PO 4 -P concentrations in Yueqing Bay showed an increasing tendency (DIN, P = 0.0262, t = 2.5341; PO 4 -P, P = 0.3196, t = 1.0383), whereas, during 2015-2020, the DIN and PO 4 -P concentrations showed a decreasing though non-significant tendency (P > 0.1).
The observed and remote-sensing-retrieved data were highly consistent in terms of seasonal variations, long time-series variations, and the spatial distribution of nutrient concentrations in Yueqing Bay. The field measurement data were mainly sourced from in situ surveys by marine research vessels, with high labor and material costs, limited observation stations, small spatial coverage, and low observation frequency. Based on its spatiotemporal continuity, field measurement data are less effective for investigations into the complex issues in Yueqing Bay, such as large-scale pollution prevention, dynamic monitoring of environmental pollution, long time series recording, multiple driving factors, and environmental remediation and restoration. In contrast, remote sensing technology not only shortens the data acquisition period but also improves the timeliness and spatial coverage of marine data. Hence, in comparison to field measurement data, the satellite remote sensing imagery of bays can be used to acquire environmental data with larger coverage and more continuous time series to investigate variation trends and driving factors.
Based on numerical simulation, kernel density estimation, and related literature data, seasonal variations in nutrient concentration in Yueqing Bay were mainly controlled by physical effects, including the Taiwan Warm Current and Fujian-Zhejiang Coastal Current, and biological effects, such as seasonal variations in the intensity of plankton photosynthesis in the bay. Furthermore, spatial distribution characteristics of nutrient concentration were mainly affected by synergy between the intensity of anthropogenic activities and transport by rivers along Yueqing Bay. Compared to 1980, nutrient concentrations in Yueqing Bay have increased 5-7-fold, mainly due to significantly enhanced anthropogenic activities from 1998 to 2015; in particular, during 1988-2005, large-area reclamation projects (e.g., Xuanmen phase II project) have resulted in a remarkable decline in tidal influx in the bay and a significant decrease in the water exchange rate, with nutrients thus accumulating inside the bay, rather than being transported to the open sea, resulting in a significant increase in nutrient concentration during this period. After 2005, the urban area along Yueqing Bay increased significantly, resulting in an increase in anthropogenic activities, discharge of land-derived pollutants, and transportation of pollutants along inflowing streams into Yueqing Bay, and thus caused the second rapid increase in nutrient concentrations in the bay waters [40,43]. Since 2011, the overall GDP of the three cities along Yueqing Bay has increased continuously, e.g., from 2011 to 2020, the GDP of the primary and secondary sectors increased by~RMB 90 billion and annual aquatic production exceeded 80 million tons [43,76]. Due to this rapid development, industrial wastewater, agricultural irrigation run-off, domestic sewage, and aquatic farming pollutants increased, resulting in the gradual rise in nutrient concentrations in inflowing streams and the bay, especially the estuaries. After 2015, the main cities along Yueqing Bay implemented various restoration policies, and after 2017, the rising trend in gross aquatic product output slowed, leading to a decline (though non-significant) in the nutrient concentrations in waters inside the bay and improvement in water quality deterioration, with an overall positive impact on ecological remediation in Yueqing Bay.
Compared with Tokyo Bay, Yueqing Bay developed much later, and ecological remediation has been implemented gradually over the past few years. Thus, from the perspective of the developmental and remediation history of Tokyo Bay, long-term ecological remediation efforts are still needed to improve the environment in Yueqing Bay and to eliminate or reduce the negative impacts of coast engineering and other anthropogenic activities.
This paper demonstrates the satellite remote sensing techniques to monitor water quality variation in a bay under strong anthropogenic impact. An updated satellite algorithm is still needed in future work with more data collection and a deep mechanism understanding of the nutrient variation. The ±3 day time interval we used in this study to guarantee enough matched data for regional algorithm development may induce errors. Moreover, satellite remote sensing only relates to the surface layer, and the detection of profile information should also be considered, although the water column was mixing in the coast and bay much more than in the open ocean under the highly dynamic hydrological condition. In general, satellite remote sensing data with a high spatial resolution have a strong ability to monitor the water quality of bays and gulfs, but they still need to be updated and combined with more field observations to provide better data support for environmental management and in-depth mechanism understanding.