Satellite Observation of the Long-Term Dynamics of Particulate Organic Carbon in the East China Sea Based on a Hybrid Algorithm

: The distribution pattern and flux variation of POC in the continental shelf seas are essential for understanding the carbon cycle in marginal seas. The hydrodynamic environment and com-plicated estuarine processes in the East China Sea result in challenging estimates and substantial spatio-temporal variability in terms of POC concentrations. A hybrid retrieval model based on the mutual combination of the color index algorithm (CI POC ) and the empirical band ratio algorithm was applied in this study to effectively and dynamically monitor the surface POC concentration in the East China Sea in a long-term series for the first time using MODIS/Aqua remote sensing satellite data from 2003 to 2020. A hybrid retrieval model based on the mutual combination of the color index algorithm (CI POC ) and the empirical band ratio algorithm was applied in this study. The MODIS/Aqua remote sensing satellite data from 2003 to 2020 were employed for the first time to dynamically monitor the surface POC concentrations in the East China Sea for a long time series. The results demonstrated that the performance (R 2 = 0.84, RMSE = 156.14 mg/m 3 , MAPE = 43.30%, bias = −64.79 mg/m 3 ) exhibited by this hybrid retrieval algorithm confirms the usability of inversion studies of surface POC in the East China Sea. Different drivers such as river discharge, phytoplankton, wind, and the sea surface current field jointly influence the spatial and temporal distribution of POC concentrations in the East China Sea. This paper also verifies that the hybrid algorithm can be applied to retrieval tasks for POC in different seas with similar optical properties to the waters of the East China Sea. In conclusion, the long-term series East China Sea POC data record, which was established based on MODIS/Aqua, provides supplementary information for in-situ sampling, which will aid the long-term monitoring of POC fluxes in shelf seas. At the same time, it has also improved our understanding of the transport and spatio-temporal variability of POC in the East China Sea, enhancing our comprehension of the impact of POC on environmental changes and carbon cycling in marginal seas.


Introduction
Particulate organic carbon (POC) is a crucial component of organic matter in the water column [1] and an essential portion of the ocean carbon cycle, participating in various biogeochemical processes and influencing organic and inorganic carbon cycles [2][3][4][5]. Although shelf seas account for less than 10% of the global ocean, they also play a crucial role in the global carbon cycle [6]. As the border between the East Asian continent and the Pacific Ocean, the East China Sea shelf is an area with a very active carbon cycle [7] and one of the most productive marginal seas in the world [8,9]. On the one hand, obtaining the distribution pattern of POC in the East China Sea can help us understand the potential size of POC as a form of carbon storage. On the other hand, it can help to clarify the spatial and temporal characteristics and drivers of organic carbon in the East China Sea and deepen our understanding of the role played by particulate organic carbon within the carbon cycle. Therefore, studying POC in the East China Sea is vital to understand the global ocean carbon cycle.
Unfortunately, the spatial and temporal variability of POC content in the East China Sea has not been fully characterized due to sparse conventional oceanographic sampling. In contrast, satellite remote sensing is labor-saving, has extensive coverage, and can be observed for long periods, providing a unique method for monitoring marine POC and a beneficial alternative tool for exploring the long-term dynamics of POC [10]. Since the first POC remote sensing retrieval algorithm by Stramski [11], many empirical and semi-analytical algorithms have been utilized to estimate ocean surface POC concentrations [6,[12][13][14][15][16][17][18][19][20]. Although these algorithms provide reasonable estimates of POC in waters at regional or global scales, they perform poorly when applied to complex optical properties in waters or coastal waters with numerous POC sources [17,19]. In addition, IOP observations are not available for many coastal waters, which hinders the development, validation, and application of IOP-based algorithms [21].
Large-scale mapping of POC concentrations in coastal waters is challenging due to the lack of stable estimation inversion models for complex waters with multiple POC origins. Le [21] created an ocean color index method (CIPOC) for the retrieval of POC concentrations in the global open ocean and coastal waters based on absorption coefficients of waters (three-band difference) with robust application to minimize redundant noise effects and enhance algorithm performance. Nevertheless, the algorithm still has significant uncertainties, particularly in turbid and complex coastal water environments. To this end, a hybrid retrieval model is proposed to accurately estimate POC concentrations in turbid waters using a simple water classification method based on the different optical properties of waters in this study. This study investigates the applicability and transferability of the hybrid inversion model for retrieving sea surface POC in coastal waters with complex optical properties and explores the mechanism of POC dynamics. In this paper, we collected measured POC data, matched remote sensing images in the East China Sea, and applied the constructed hybrid retrieval algorithm to the East China Sea. The spatial distribution and annual variation of POC in the East China Sea from 2003 to 2020 are revealed for the first time in a long-term series.
The results not only validate the fantastic performance and applicability of the hybrid algorithm for POC inversion in coastal waters but also demonstrate the potential of MODIS/Aqua remote sensing data for the long-term series monitoring of POC concentrations in the surface layer of the East China Sea. The achievement also provides a basis for further applications in other coastal waters. At the same time, the main drivers affecting the distribution pattern of POC in the sea were investigated, and our understanding of the future POC concentration and flux variation and the carbon cycle in the East China Sea was strengthened.

Study Area and Field Measurements
The East China Sea, located east of mainland China, is one of the most productive marginal seas on Earth [9] and receives a substantial supply of terrigenous material imported from large rivers such as the Changjiang and Qiantang [8,[22][23][24]. In addition, the East China Sea is a highly variable marine ecosystem where different water masses introduced or discharged by the Kuroshio Current, Taiwan warm current, Northern Jiangsu current, and Zhejiang-Fujian coastal current intersect and interact, resulting in complex optical properties and significant seasonal variations in water quality parameters in this region [25]. In this paper, field surveys of four cruises in the East China Sea were conducted from May to June 2014, July-August 2017, and May-September 2019 (as illustrated in Figure 1A,B), respectively, and a total of 120 sample points were collected for the development and verification of the inversion model. Another different bio-optical dataset (N = 33) was collected in the Chaoshan sea area (as demonstrated in Figure 1C) from June to July 2021 to verify the migration ability of the POC retrieval model under different water optical conditions.

Satellite Data
MODIS/Aqua's Rrs products were downloaded from the NASA website (http://oceancolor.gsfc.nasa.gov/ (accessed on 10 September 2021). Based on the SeaWiFS Data Analysis System software (SeaDAS, version 8.0, NASA Ocean Biology Processing Group, Greenbelt, Maryland, USA), the remote sensing data were reprocessed utilizing the latest calibration coefficients and algorithms (Version 2018.0, NASA Ocean Biology Processing Group, Greenbelt, Maryland, USA) and extracted Rrs of position and time matching the in-situ POC. Due to the preciousness of the measured data, a time gap of ±5 days was used to match satellite data. The median value of a 3 × 3 window centered on each sampling point was substituted for the pixel value of the corresponding sampling point to filter the satellite sensors and remote sensing algorithm's noise. Analysis and comparisons between in-situ POC and satellite-derived POC were only made when the number of valid pixels (i.e., after discarding the pixels that failed quality control flags) in the box was >4 and the coefficient of variation (CV) among pixel values was <0.4 [26]. A total of 79 satellite remote sensing reflectance datasets matched with in-situ measured samples were filtered for MODIS/Aqua (as shown in Figure 1a,b), which were employed to establish the POC time series for 18 years between 2003 and 2020. Another 15 MODIS/Aqua remote sensing reflectance datasets matched with in-situ measured data (as shown in Figure 1c) were observed near the Chaoshan sea area to prove the transferability of the POC remote sensing algorithm.

Classifying Waters and Building a Hybrid POC Algorithm
Hu [27] proposed an algorithm based on the difference in reflectance (i.e., color index algorithm) between the blue, green, and red bands to invert the surface chlorophyll-a concentration in nutrient-poor oceans. Considering the advantages of the CI algorithm in estimating chlorophyll-a and the inheritable correlation between chlorophyll-a and POC, Le [28] improved the algorithm based on the measured POC data collected on SeaBASS and the dataset matched by MODIS/Aqua, which retrieved the POC concentration of the global ocean. Aiming at the high uncertainty of remote sensing reflectance at 412 nm and 443 nm, a three-band difference algorithm (Equation 1) of 488 nm, 547 nm, and 678 nm was developed. The result demonstrates that the algorithm has an advantage over the traditional empirical band-ratio algorithm for inversion in open waters. For a more detailed description of the algorithm, please read [28]. The algorithm is expressed as follows: where CI POC is the color index of POC, and Rrs(x) is the central wavelength at different bands in MODIS/Aqua. Despite the satisfactory performance of the CIPOC algorithm in inverting POC in oceanic waters, the algorithm still has significant uncertainties when targeting particularly complex coastal waters. The optical properties of the East China Sea exhibit significant differences from the estuary of the Changjiang River to the open ocean. As shown in Figure 1a,b, two distinct optical characteristics of the remote sensing reflectance spectra of the Changjiang River estuary and the East China Sea waters exist. The different optical properties of the waters were classified, and different forms of inversion models were constructed to meet the accuracy requirements. The Changjiang River releases enormous amounts of sediment and organic matter into the East China Sea. Under the dual influence of high sediment concentration and CDOM concentration, the Rrs spectra of the Changjiang estuary waters generally show a reflection peak at 547 nm. In contrast, the Rrs in the open waters of the East China Sea gradually decrease with increasing wavelengths. As demonstrated in Figure 2, using the unique patterns of these optical properties, this paper provides a simple classification of the different types of remote sensing reflectance into two major categories, as follows: A separate remote sensing inversion algorithm was developed based on the different Rrs types mentioned above in this study. The water column of type I (as illustrated in Figure 2a) exhibits typical open water characteristics, and therefore the CIPOC algorithm is employed. The continuous import of river plumes has contributed to the resuspension of sediments in the Changjiang estuary, creating a highly turbid sea area. The inability of phytoplankton-related particles to dominate the optical properties in this sea area will prevent the CIPOC algorithm, which was established based on the absorption coefficient of waters, from accurately estimating the POC concentration. For this reason, a band-ratio algorithm needs to be constructed to monitor the POC concentration in turbid water bodies such as the Changjiang estuary based on the empirical relationship between the measured data and the remote sensing reflectance (as shown in Figure 2b). In the process of algorithm construction, this study verified various band-ratio algorithms and found that Rrs(645)/Rrs(547) was the most suitable for measuring POC. A study by Le [29] method is preferred over the existing chlorophyll ocean color algorithm. Therefore, a hybrid inversion model was constructed by combining the CIPOC algorithm and the empirical band-ratio algorithm for application in the East China Sea to verify the algorithm's performance for estimating the retrieved POC concentrations in the shelf sea in this paper.

Environmental Datasets
This paper compares satellite-retrieved POC data with several environmental forcing variables, including the river discharge of the Changjiang River and chlorophyll-a concentration data. For major rivers in China, the Ministry of Water Resources of the People's Republic of China (MWRC) conducts river discharge measurements at the corresponding hydrological stations of the rivers. It compiles a bulletin to be published on its official website (http://www.mwr.gov.cn/sj/#tjgb (accessed on 10 September 2021)) [10]. The monthly mean river discharge data from the most downstream non-tidal hydrographic station, Datong Station, published by the MWRC for 2003-2020 was used in this study to investigate how the Changjiang River's flushing water influences the POC dynamics in the East China Sea. Considering that the study area is located near the mouth of the Changjiang River, ocean color satellite data usually overestimate chlorophyll-a concentrations in turbid coastal waters. Cui [30] evaluated satellite observations and in-situ measurements in the turbid coastal waters of the Bohai Sea from 2003 to 2007. The results revealed that although satellite observation overestimates field measurement results to some extent, the two observations were consistent in terms of spatial distribution and time variation. He [31] compared the measured chlorophyll-a from 2006 to 2007 cruises in the Bohai Sea, Yellow Sea, and East China Sea with that retrieved by MODIS satellites. The experimental results demonstrated that the in-situ measurements were strongly correlated with the satellite-derived results (R = 0.72, SD = 0.40, p < 0.0001), and the chlorophylla products derived from the satellite inversion were calibrated based on this relationship. For this purpose, this paper calibrates the L3 chlorophyll-a (chlor_a) product from MODIS/Aqua processed by NASA's OPBG group based on this method to eliminate the overestimation of chlorophyll-a by turbid waters. The time series of localized monthly chlorophyll-a concentrations in the East China Sea from 2003 to 2020 were estimated using the calibrated chlorophyll-a products to analyze the synergistic relationship between phytoplankton and POC. The chlorophyll-a data used subsequently in this paper were also calibrated using this method.

Accuracy Evaluation of POC Retrieval Algorithm
For the determination of the best regression fit effect, this study quantified the superiority of fit by examining the explained coefficient of determination (R 2 ) (defined in Equation (4)) in this paper. The root mean square error (RMSE) (defined in Equation (5)) and the mean absolute percentage error (MAPE) (defined in Equation (6)) were calculated to assess the predictive performance of the algorithm. The significance test (p) was utilized to indicate errors occurring in the algorithmic model. Bias (defined in Equation (7)) was applied to measure the accuracy of prediction results. The related equations are expressed as follows: where POC measured is the in-situ POC concentration, POC retrieved is the POC concentration estimated by the inverse model, POC mean is the average of the measured POC concentrations, and N is the number of sample points. In this paper, all data were sorted according to the concentration of the POC, and the extraction interval was determined based on the number of sample points. Finally, training and testing samples were extracted from the sample dataset. The newly composed training and testing samples had the same interval (arithmetic sequence) to ensure that the POC concentration gradient had a uniform distribution to meet the requirement of constructing a high-precision inversion model. In this way, the sample dataset was divided into two parts: approximately 70% of the sample points were designated to calibrate the algorithm's coefficients, and the remaining 30% were designated to evaluate the performance of the POC retrieve model. All statistical analyses were undertaken using R.4.12 software.

Performance Evaluation of Hybrid Retrieval Models
The 79 MODIS/Aqua satellite remote sensing reflectances matched to the measured sample points were classified using the methodology in Section 2.3, and 70% of the classification results were employed for derivation model parameters and performance evaluation. The sample points in the type I dataset (N = 32) were applied to the CIPOC algorithm (Equation (1)), and the sample points in the type II dataset (N = 23) were applied to the band-ratio algorithm. These two datasets were linearly related to the measured log10(POC) data to determine the coefficient of the best-fit regression, respectively (as shown in Equations (8) and (9): The correlation coefficients R 2 of the CIPOC algorithm and the band-ratio algorithm are 0.71 and 0.75, respectively, according to Figure 3a,c, and the significance test p is less than 0.001 for both, indicating a favorable linear relationship between the measured

Comparison with Different Forms of Algorithms
Accuracy is critical to ocean retrieval because time-series analyses call for the most robust products in order to stabilize temporal variation. However, current POC algorithms for marginal seas often have significant uncertainties due to the optical complexity of coastal waters. We evaluated the inversion performance of multiple forms of POC algorithms in the East China Sea based on Stramski [32].
Among them, the CAT1 algorithms are referred to as BR-PF (band ratio-power function) and have the form shown in Table 1, where X is usually the ratio of the two bands. The current NASA standard POC algorithms fall into this category. For example, a global standard POC concentration product has been produced using this algorithm based on the MODIS bands [19]. The CAT2 algorithm is called MBR-OCx (maximum band ratio-OCx), where x indicates the number of bands involved in the equation of the algorithm. This algorithm has also been used to manufacture global ocean color (OC) chlorophyll-a products [33,34]. For this reason, the CAT2 algorithm was migrated to the East China Sea for POC inversion in this paper. The CAT3 algorithm is a color index algorithm (CI) developed to improve chlorophyll algorithms in waters with low chlorophyll-a based on the concept of three-band reflectance difference [27] and is now also used in global POC estimation [28]. The CAT4, CAT5, and CAT6 algorithms are conceptually similar to the normalized difference vegetation index (NDVI), initially developed for terrestrial vegetation, and the algorithm formulation was later rewritten and used by Son [17] for POC estimation in the Gulf of Mexico. The CAT7 algorithms are based on the multiple linear regression (MLR) approach. Le [21] and Chen [35] applied their similar MLR algorithm in POC studies of the Louisiana shelf and the Ross Sea, respectively, and obtained excellent inversion results.    Table 1 mainly defines the inversion models for the seven categories mentioned above, and many alternative candidate formulas are further derived in each category. In the comparative analysis for the purpose of evaluating all algorithms, the best coefficients were fitted using the same training dataset based on ordinary least squares. The accuracy of each algorithm in estimating the POC was quantified using the same testing dataset using the statistical measures in Section 2.5. For brevity, here, we limit the presentation of these algorithms to the case of MODIS bands. Since the MODIS sensor lacks a band centered at or near 510 nm, a virtual band of 510 v was constructed based on the reflectivity of 488 nm versus 547 nm, drawing on the method of Stramski [32]. For more details on the description, please refer to the literature by Stramski [32]. As can be seen from Table 1, the POC concentrations retrieved by the locally optimized inversion models are widely underestimated, with high uncertainties. It is worth noting that the CIPOC-based CAT3 algorithm performs the worst in terms of statistical metrics (slope = 0.36, R 2 = 0.34, RMSE = 374.26 mg/m 3 ), which also demonstrates that the CIPOC algorithm is not sensitive enough to accurately estimate POC in particularly turbid and complex waters. Although the hybrid inversion model features an overestimation phenomenon (slope = 1.14, bias = −64.79 mg/m 3 ), it still demonstrates better performance (R 2 = 0.84, RMSE = 156.14 mg/m 3 , MAPE = 43.30%) and can be used for the retrieval of POC concentrations in the East China Sea shelf area and to capture the dynamic variation of POC in coastal waters.

Spatial and Temporal Variation of POC
After verifying the POC hybrid remote sensing retrieval model, Equations 8 and 9 in Section 3.1 were applied to the MODIS/Aqua remote sensing images of the East China Sea area from 2003 to 2020. The calculated annual average concentration of POC in this sea area based on the East China Sea area (25-32°N, 120.5-126°E) is plotted in Figure 1. The seasonal average POC concentration values for spring (March to May), summer (June to August), autumn (September to November), and winter (December to February) were also calculated using MODIS/Aqua satellite remote sensing data concerning POC concentration values from 2003 to 2020 (as shown in Figure 4). Figure 5 presents the long-term series distribution of the average POC concentration in the sea area near the sampling site from 2003 to 2020. Each scene image is synthesized from all available remote sensing data for that year.

Driving Factors of POC Dynamics on the East China Sea
The effects of different drivers such as phytoplankton, TSM, tidal processes, precipitation, temperature, river discharge, surface sea currents, and wind speed on the distribution pattern of POC concentrations in seawater have been the subject of previous studies [3,21,36,37]. Given that the distribution of POC concentrations in the East China Sea shows significant spatial and temporal variations, this paper attempts to identify the causes of this phenomenon. In contrast, the components of POC in the East China Sea are extremely intricate. The spatiotemporal variability of sources such as high primary productivity, riverine flux, sediment resuspension, and DOC transformation in different waters of the East China Sea exacerbates the conclusion that the origins of POC are inconsistent [38], so the synergistic effects of different driving factors for POC in the East China Sea need to be considered. The monthly average POC values of the region were counted according to the drawn sea area range (25-32°N,120.5-126°E) and combined with the monthly average river discharge from the Datong hydrological station on the Changjiang River mainstream for correlation analysis, which is exhibited in Figure 6. The river discharge of the Changjiang River is significantly negatively correlated with POC concentration (r = −0.60, p < 0.001). During the spring and autumn, the river discharge remains stable. The dilution of water brought by the Changjiang River will further dilute the original POC content in the East China Sea and contribute to a decrease in the POC concentration; thus, the river discharge in spring and autumn has a strong negative correlation with POC concentration (spring r = −0.34, p < 0.05; Autumn r = −0.48, p < 0.05). Moreover, extreme weather events such as heavy rainfall weather and increased frequency of flooding during the summer period lead to a sharp increase in river discharge and affect the flux and composition of the terrestrial material transported from the river to the sea [39]; thus, releasing more terrestrial organic matter into the marginal sea means that the river discharge exhibits a degree of positive correlation with POC (r = 0.41, p < 0.001). Despite the annual rainy season (June-August), which is the peak discharge period of the Changjiang River and the main period for the import of organic matter from terrestrial sources [40], substantial amounts of soil are washed into the river by precipitation [41]. As a result, the POC decreases with the increase in suspended matter due to the adsorption of inorganic minerals [10]. The annual rainy season is also the season with the lowest POC concentration. As precipitation gradually decreases in winter and the river discharge diminishes, there is also a decrease in terrestrial organic matter entering the ocean from rivers. Figure 6 also illustrates no statistical correlation between the two in this period. Due to it being one of the most productive marginal seas in the world [9], high primary productivity and phytoplankton blooms are frequently observed in the East China Sea. The occurrence of such events also leads to the export of significant POC. In recent years, spring phytoplankton blooms have been a recurring event in the East China Sea [42]. Phytoplankton in the ocean's surface layer utilizes photosynthesis to convert dissolved CO2 into its cellular components and POC. For this purpose, the effect of phytoplankton bloom events on POC distribution in the East China Sea requires investigation. As an essential proxy for phytoplankton biomass in the ocean, chlorophyll-a can be identified by satellite-derived data to characterize the area of phytoplankton blooms [31,43]. This study applied remote sensing images with low cloudiness in spring and a chloro-phyll-a threshold of >10 mg/m 3 to evaluate whether a phytoplankton bloom event occurred in this sea area [31]. Figure 7 plots the spatial and temporal distribution of chlorophyll-a and POC concentrations on April 8 and 19, 2018. Interestingly, the chlorophyll-a concentration in the red boxed area is abnormally high and greater than 10 mg/m 3 , thus determining that a phytoplankton bloom event occurred in the East China Sea at this time (as shown in Figure 7b). In contrast, high POC concentrations can also be observed at the corresponding location (as demonstrated Figure 7c). After about ten days, the chlorophylla concentration at the corresponding location returned to normal, and no POC concentration abnormality was observed (as represented in Figure 7d,e). This phenomenon demonstrates that a short phytoplankton bloom can produce abnormally high POC concentrations in the region of the East China Sea. With the depletion of nutrients, the phytoplankton will also disappear rapidly, contributing to the gradual disappearance of the POC concentrations initially present in the region. In this study, long-term series of the dynamics of chlorophyll-a and POC were also determined using satellite archived data for the last two decades. Considering that there are also geographical differences in the East China Sea, the original sea area was divided into two regions for analysis: the Changjiang estuary region (30-32°N,120.5-126°E) and the East China Sea region (25-30°N,120.5-126°E). The results in Figure 8 show that the proliferation of marine phytoplankton in winter plays an important regulatory role regarding POC ((a) r = 0.93, p < 0.001; (b) r = 0.94, p < 0.001). This conclusion is consistent with the research of [21,36], who demonstrated the POC concentration in seawater increases with higher chlorophyll-a concentrations, that there is no correlation between winter river discharge and POC, and that POC dynamics during this period is dominated by phytoplankton (chlorophyll-a). The same statistical results exist for chlorophyll-a and POC in summer ((a) r = 0.75, p < 0.001; (b) r = 0.85, p < 0.001), but both are regulated by the import of nutrients from high-intensity rivers. The Changjiang River carries enormous amounts of sediment into the East China Sea, making the optical properties of the waters in this area exceptionally complex [9], interfering with the correlation between chlorophyll-a and POC and affecting the dynamics of POC concentrations. Although phytoplankton bloom events frequently occur during spring and autumn, Figure 8 reveals that the correlation between POC concentrations and phytoplankton is insignificant, which also indicates that sudden changes in chlorophyll-a concentrations over a short period are not capable of causing a trend shift, while other processes may also regulate the spatial and temporal distribution of POC. Wind and sea surface currents can also influence the migration of seawater, thereby indirectly affecting the spatial distribution of POC. For example, extreme weather events such as tropical cyclones or typhoons cause POC migration driven by solid winds that trigger sea surface currents. To further investigate the dynamic variation of POC before and after the transit of a typhoon, this study selected Typhoon TAPAH (No. 1917), which transited in the East China Sea waters on 21-22 September 2019, with a maximum wind speed of 33 m/s and a wind circle radius of 500 km during the transit (http://typhoon.nmc.cn/web.html (accessed on 24 June 2022), as shown in Figure 9a). Figure 9b demonstrates the composite map of POC concentration distribution for seven days before the typhoon's transit. The POC in the coastal waters of the Changjiang estuary was still at a high concentration level, while the POC concentration in the open waters of the East China Sea was low. Figure 9c indicates a composite map of the POC concentration distribution for seven days after the passage of Typhoon TAPAH, where the POC concentration in the coastal waters of the Changjiang estuary showed a significant increase and the POC concentration in the open waters of the East China Sea increased significantly (lower right area of Figure 9b). In this study, two sub-regions were mapped to highlight the effect of typhoons on POC transportation. There is a distinct downward trend in the POC exhibited in Zone 1 compared to the pre-transit period, and higher concentrations of POC in Zone 2 spread from coastal waters to more distant waters (>122°E). Such extreme weather events induce water mixing, enhance land runoff, and promote sediment resuspension, causing an increase in sea surface POC concentrations and diffuse migration to more distant open waters [44]. In contrast, when no extreme weather events occur, wind and sea surface currents also carry out these hydrodynamic processes, but the processes are relatively slow; however, they affect the dynamic distribution of POC to some extent. Other studies have revealed that ocean currents in different regions of the East China Sea impact the distribution of POC. For example, the sea surface currents in the northeastern region flow roughly southeastward in spring and autumn, accelerating the rapid transport of the Changjiang River plume, which results in higher POC concentrations in this sea area. The prevalence of the northwest monsoon in winter enhances Ekman pumping, and the surface currents also reach their maximum at this time, thus driving convective mixing in the upper ocean and increasing the depth of the mixed layer, thus transporting bottom nutrients to the surface [36,45,46]. These phenomena also promote the rapid growth of phytoplankton and increase biomass and POC concentrations. In summer, the south wind reverses the direction of local surface currents, and the Changjiang river's dilution water is restricted to the northeastern part of the East China Sea, which also prevents the transfer of terrestrial substances to more distant waters, resulting in low POC concentrations in summer. Similar effects in terms of sea surface currents on POC concentrations can also vary seasonally. Winds also interact with the distribution of POC in the East China Sea by driving surface currents, as the prevailing northwesterly winds in winter accelerate the movement of the Jiangsu coastal current to the southeast. In contrast, during summer, prevailing southeasterly winds contribute to the northwestward flow of the Taiwan warm current [47]. These current movements contribute to higher POC concentrations in the East China Sea in winter than in summer, since the Jiangsu coastal current brings high POC concentrations while the Taiwan warm current does not [3].
As mentioned above, the seasonal variation in terms of POC concentrations in the East China Sea is mainly influenced by river discharge and phytoplankton. At the same time, the relationship between POC and chlorophyll-a is also influenced by the Changjiang river's dilution water, resulting in significant territorial differences. The import of terrestrial material and phytoplankton reproduction affect the spatial distribution of POC in the East China Sea, where winds and surface currents also accelerate the migration and transformation of POC in seawater, contributing to spatial variability in terms of POC in the East China Sea.

Performance and Applicability of Algorithms in Other Marine Areas
Satellite remote sensing is inexpensive, consumes extremely little labor, and is an excellent alternative tool for exploring long-time series POC variability [10]. However, obtaining accurate POC predictions in nearshore waters remains challenging due to the complex optical properties of coastal waters and the many different sources of water column components that produce a range of biochemical contaminants of POC [21]. Recognizing the difficulty of retrieving POC in coastal waters, this paper successfully estimates POC in the coastal and open waters of the East China Sea using a newly constructed hybrid retrieval algorithm. Since the CI algorithm effectively exploits the effect of variations in absorption and backscattering coefficients on remote sensing reflectance based on the water column, it can improve uncertainties caused by atmospheric corrections, inversion errors, and sensor noise. It also provided a satisfactory performance in the retrieval application in the open waters of the East China Sea (R 2 = 0.71, RMSE = 45.97 mg/m 3 , MAPE = 52.21%, bias = 11.54 mg/m 3 ).
From the remote sensing reflectance spectra matched to the sampling points in Figure  1c, the optical properties of water bodies in the tidal waters remain consistent with the spectral distribution characteristics of water bodies in type I. For this reason, the algorithm form of type I in the hybrid algorithm (i.e., Equation (8)) was applied to the Chaoshan sea area to examine its potential applicability. The reasonably favorable statistical parameters (slope = 0.67, R 2 = 0.81, RMSE = 23.39 mg/m 3 , MAPE = 19.58%, bias = −14.29 mg/m 3 ) are exhibited in Figure 10, and the POC concentration distribution in the Chaoshan sea area retrieved from MODIS/Aqua remote sensing images using this hybrid algorithm is presented in Figure 11. The results prove that the hybrid inversion algorithm can accurately capture the dynamic variation of POC in the Chaoshan sea area (Figure 11). At the same time, the inverse POC values also have good consistency compared with the in-situ data, with similar spatio-temporal distribution patterns to East China Sea waters. The reliability shown by the model in the Chaoshan Sea implies that the hybrid algorithm used in this study can be extended to other sea areas with similar water types. However, it is essential to be cautious before validating the migration of the model. Empirically, this algorithm works stably only in similar environments, for which it was designed. When the model is transferred and implemented in waters with typically complex optical properties, the performance often fails to meet the requirements. The model form needs to be dynamically adjusted for a local optimization to meet the application requirements based on the sampled measured sample points.

Assessment of Uncertainty
It has been suggested that the backscatter coefficient may be strongly correlated with POC and can be a proxy for POC concentration levels in the open ocean [11,19,28,32,48]. However, Hunter [49] concludes that organic particles may contribute enormously to the absorption and slightly less to the backscattering coefficients. This inference indicates that absorption may be a better optical proxy for retrieving POC than the backscattering coefficient and that the total absorption coefficient of seawater significantly causes the variation in remote sensing reflectance. Therefore, some connections can be established by remote sensing reflectance for POC retrieval. However, the performance of the hybrid algorithm based on remote sensing reflectance may be influenced by many factors during the remote sensing application.
The Rrs variables utilized in the hybrid algorithm are estimated from satellite inversions, where the algorithm established between Rrs and POC does not relay the measured relationships between these variables. The Rrs variables are also subject to processes such as specific algorithms for atmospheric correction and the processing of raw top-of-atmosphere radiation signal data, which can be incorrectly calculated and introduce uncertainty into the inversion POC [28]. In extremely turbid waters with complex optical properties, strong backscattering influenced by high concentrations of non-algal suspended matter can significantly increase water-leaving reflectance in the green and red wavelength bands. Although the increase in water-leaving reflectance in the blue band is relatively small compared to the longer band, the strong absorption caused by extensive detritus and chromophoric dissolved organic matter (CDOM) reduces the blue-green band reflectance [31,50].
These same factors can make Rrs potentially erroneous in the inversion process. As shown in Figure 3. (a), it can be observed that some sample points in the lower right cluster are more densely distributed in the area with CIPOC values of 0-0.001 and log10(POC) values of 1.8-2.4. After data analysis, the positions of sampling points that deviate from the regression line are mainly distributed near the mouth of the Changjiang River, resulting in higher POC concentrations in the sampled waters. In this case, the signal of suspended sediment and CDOM will reach the MODIS sensor through the reflection and scattered light from the waters when the satellite is conducting the waters signal observation, which causes the reflectance error. The sensor fails to promptly capture the rapid variation of POC in the waters over a short period, resulting in a significant difference between the POC retrieved by the model constructed based on the remote sensing reflectance and the actual measured value. These errors also bring some uncertainties to the algorithm in terms of calibration and validation.
All the factors described above may vary with river discharge. Firstly, river discharge affects phytoplankton biomass by influencing freshwater export fluxes and water column stratification [26]. The Changjiang river dilution water also carries abundant nutrients and accelerates the dynamic distribution of phytoplankton in water column mixing regulating regions, thereby affecting POC concentrations. Secondly, rivers also release enormous amounts of silt to resuspend sediments in shallow waters [51] or deepen the surface mixing layer, which reduces the availability of light to the phytoplankton community as a whole and impacts its photosynthesis, indirectly decreasing the production of POC. Finally, high concentrations of suspended sediment and abundant organic matter exported by rivers can significantly affect the relationship between Rrs and POC, contributing to errors in the estimation of POC. In addition, higher sea surface temperatures may reduce phytoplankton processes, such as enhanced phytoplankton respiration (decreased net primary productivity) and increased zooplankton grazing. Variations in phytoplankton community structure associated with changes in nutrient, temperature, and grazing regimes [26] can indirectly affect POC estimates. Simultaneous matching validation was performed using MODIS/Aqua remote sensing imagery with a time window of ±5 days and in-situ data. Although the modeling and validation results are more satisfactory, the lack of timeliness means this method cannot quickly capture temporal and spatial variations in POC over a short period, limiting the accuracy of the inversion results. In field measurement data, inconsistent sample acquisition methods and measurement steps may introduce errors. Factors such as limited samples, a small number of sample points, and an uneven spatio-temporal distribution may also affect the dynamic analysis of POC and yield incorrect conclusions.
In summary, the results of this study are encouraging, even though these factors can introduce more uncertainties in the retrieved POC. The performance exhibited by the hybrid model constructed based on the CIPOC algorithm and the empirical band-ratio algorithm demonstrates the potential applicability of the remote sensing inversion of the POC of the Chinese shelf seas, especially for applications in complex coastal waters.

Conclusions
The CIPOC algorithm based on the absorption coefficient of waters does not perform well in extremely complex water environments. Still, it can be adopted in the remote sensing and retrieval of POC in coastal and open waters when combined with an empirical band-ratio algorithm. The new hybrid algorithm provides acceptable verification results. The estimated POC concentration is consistent with the in-situ POC at the height of the matching points, for which the spatial and temporal distribution of POC along the coastal and open waters of the East China Sea from 2003 to 2020 was successfully retrieved. Factors such as river discharge and phytoplankton largely determine the dynamics of POC concentrations in the East China Sea. At the same time, wind and sea surface currents also influence the transport and transformation of POC.
Future improvements in the remote sensing of complex waters by new sensors could improve accuracy in the monitoring of POC fluxes in different regions of the shelf seas by using Sentinel-3/OLCI data with higher sensitivity and better resolution in coordination with MODIS/Aqua data. This study also provides a new method for monitoring coastal POC distribution dynamics in the long term and on a large scale. It is believed that with an improvement in the quality and quantity of field measurement data and continuous research into the remote sensing inversion mechanism, the regional algorithm of inversion for POC will continue to develop, improve, and better serve ocean monitoring research. It can also provide valuable clues for other marine studies in similar water environments, assist in estimations of the carbon budget of nearshore and open waters, and aid in formulating related policies such as peak carbon dioxide emissions and carbon neutrality.