Sampling Uncertainties of Long-Term REMOTE-Sensing Suspended Sediments Monitoring over China’s Seas: Impacts of Cloud Coverage and Sediment Variations

: Satellite-based ocean color sensors have provided an unprecedentedly large amount of information on ocean, coastal and inland waters at varied spatial and temporal scales. However, observations are often adversely a ﬀ ected by cloud coverage and other poor weather conditions, like sun glint, and this inﬂuences the accuracy associated with long-term monitoring of water quality parameters. This study uses long-term (2013–2017) and high-frequency (eight observations per day) datasets from the Geostationary Ocean Color Imager (GOCI), the ﬁrst geostationary ocean color satellite sensor, to quantify the cloud coverage over China’s seas, the resultant interrupted observations in remote sensing, and their impacts on the retrieval of total suspended sediments (TSS). The monthly mean cloud coverage for the East China Sea (ECS), Bohai Sea (BS) and Yellow Sea (YS) were 62.6%, 67.3% and 69.9%, respectively. Uncertainties regarding the long-term retrieved TSS were a ﬀ ected by a combination of the e ﬀ ects of cloud coverage and TSS variations. The e ﬀ ects of the cloud coverage dominated at the monthly scale, with the mean normalized bias (P bias ) at 14.1% ( ± 2.6%), 7.6% ( ± 2.3%) and 12.2% ( ± 4.3%) for TSS of the ECS, BS and YS, respectively. Cloud coverage-interfering observations with the Terra / Aqua MODIS systems were also estimated, with monthly P bias ranging from 6.5% ( ± 7.4%) to 20% ( ± 13.1%) for TSS products, and resulted in a smaller data range and lower maximum to minimum ratio compared to the eight GOCI observations. Furthermore, with approximately 16.7% monthly variations being missed during the periods, signiﬁcant “missing trends” e ﬀ ects were revealed in monthly TSS variations from Terra / Aqua MODIS. For the entire region and the Bohai Sea, the most appropriate timeframe for sampling ranges from 12:30 to 15:30, while this timeframe was narrowed to from 13:30 to 15:30 for observations in the East China Sea and the Yellow Sea. This research project evaluated the e ﬀ ects of cloud coverage and times for sampling on the remote sensing monitoring of ocean color constituents, which would suggest the most appropriate timeframe for ocean color sensor scans, as well as in situ data collection, and can provide design speciﬁcation guidance for future satellite sensor systems.


Introduction
Effective monitoring of water quality in coastal and open seas, including variations of phytoplankton levels, suspended sediments and colored dissolved organic matter, calls for uninterrupted sampling over appropriate time scales. Field survey-based water quality monitoring, using in-situ sampling, land or shipboard laboratory analysis, are precise but time-and labor-costly. Remote sensing data are favored for the advantages of rich spatial-temporal-spectral information. The benefits of remote sensing include its wider temporal and spatial range (coverage and scale) for water resource monitoring, and relatively low costs compared to the conventional field sampling approach [1,2]. Significant efforts have been made to applying remote sensing technologies to monitor water quality from local to global scales, and as a result, considerable progress has been achieved in our understanding of the behavior of coastal and inland waters [3][4][5][6].
Important ocean color sensors orbiting the earth include the Moderate Resolution Imaging Spectroradiometer (Terra/Aqua MODIS, 1999-present), the US Visible Infrared Imager Radiometer Suite (VIIRS, 2011-present), the Ocean and Land Color Instrument (Sentinel-3A OLCI, 2016-present) and the Korean Geostationary Ocean Color Imager (GOCI, 2010-present). Prepared to launch are several other ocean color sensors, including the Chinese HY-series satellites and NASA PACE mission [7,8]. Ideally, ocean color remote sensing unprecedentedly provides data for water quality properties at various spectral, temporal and spatial scales [9][10][11]. Unfortunately, applications of both in-situ and remote sensing methods are often limited by either low temporal or spatial resolutions, for regions that experience diurnal or semi-diurnal variations [12,13]. For instance, the spatial coverage or sampling frequency of traditional field sampling is often too small to derive meaningful conclusions [14,15]. The reliability of these methods is primarily determined by sampling protocol strategies with respect to sampling times and intervals. In most conditions, valid samplings are often limited: routine field surveys are generally conducted monthly by local authorities, such as the environmental or water quality agency. The most widely used Terra/Aqua MODIS satellites, designed to perform two observations a day, are particularly prone to disruption from cloud coverage and other poor weather conditions, as well as sun glints [16][17][18].
Although ocean color sensors, including the MODIS system, have been widely used for monitoring both short-and long-term water quality, previous studies have shown that observations are often hampered by cloud coverage, thick aerosol or sun glint [19]. As a result, it was found that each day, only 12%, on average, of the Aqua/MODIS ocean products were valid [10,19,20], and for some regions under persistent cloud coverage, monthly averaged valid observation could be less than one. High-latitude and equatorial oceans experienced the lowest levels of cloud-free period, with the valid number of observations being less than 2.4 per month [21]. Some research has been dedicated to improving temporal coverage for a single observation area by merging data from multiple satellite sources. For instance, observations merged from SeaWiFS and MODIS systems can produce valid data for up to 23% of the total ocean area, compared to approximately 18% from the MODIS system alone. Similar improvements can be derived by merging data over a four-day period, increasing the percentage of coverage to 30%-35% [22,23]. Despite efforts to improve the sampling frequency and coverage, the availability of valid observations from all the ocean color sensors remains limited.
These limitations greatly reduce the number of feasible applications of ocean color data, and can result in erroneous estimates of water quality indicators, such as mean Chl-a and suspended sediments levels. For coastal and open sea waters, where highly dynamic suspended sediment has been observed, high-frequency monitoring is required. Observations should preferably be taken every 30 min to 1 h [24], while inappropriate sampling frequencies and/or times could introduce large biases and Remote Sens. 2020, 12,1945 3 of 26 eventually lead to unscientific decisions in water management. Previous studies have indicated that inappropriate sampling time and frequency in water quality monitoring have led to statistical errors larger than 50% [7,25,26]. The facts of missing data from CZCS remote sensing resulted in a root mean-square error (RMSE) of 8% for annual mean chlorophyll concentrations, and a RMSE of 33% when determining the peak chlorophyll values [18]. For the Gulf of Mexico region, Barnes and Hu demonstrated that the difference in monthly mean Chl-a levels between MODIS and SeaWiFS datasets could be larger than 10% if monthly valid observations are fewer than 15 [27]. Despite these challenges and efforts to improve sampling strategies, systematic research on sampling coverage and its impacts on ocean color remote sensing data remains scarce.
For the aforementioned reasons, there is an urgent need to quantify the influence of cloud coverage in remote sensing and evaluate the impact of insufficient valid observations. Accordingly, appropriate sampling times for ocean color applications could be determined, and a benchmark for future satellite sensor planning be set. Pioneering research, using polar-orbit satellites like the MODIS system, showed that simulated datasets [18,25] with limited temporal coverage were inadequate for measuring diurnal variability. Thus, the geostationary-orbit (GEO) ocean color sensors have become prevalent to meet the demands of the high-frequency sampling of water quality. The first GEO ocean color satellite, the Geostationary Ocean Color Imager (GOCI), was launched by South Korea in 2010, and takes hourly measurements during the daytime [16]. GEO satellites, including GOCI, provide high-frequency sampling to capture the dynamic of coastal/inland waters [28]. GEO data have been widely used in ocean color remote sensing applications, for monitoring parameters such as the chlorophyll, suspended sediment and colored dissolved organic matter concentrations [29][30][31]. Most commonly, studies on the remote sensing of water quality involve time-series analysis of, for example, TSS in the regions of interest. Therefore, a fundamental question is how these sampling uncertainties revealed above would affect the estimations of long-term trends in TSS observations. The differences among these time-series TSS data were examined to answer the following questions: do they show the same magnitudes of mean monthly TSS? Do they show the same seasonality (i.e., monthly maximum/monthly mean ratio)? Do they show the same inter-annual changes and long-term trends?
By using the long-term high-frequency GOCI datasets for the China seas, this study aimed, firstly, to quantify cloud coverage over targeted waters at long-term scales, and secondly, to provide a solution to resolve the effects of sampling times on the remote sensing monitoring of total suspended sediment (TSS) in highly dynamic waters. Thirdly, the level of uncertainty in the data produced by the satellite sensors most commonly used for ocean color observations, the Terra/Aqua and MODIS systems, in ocean color remote sensing were assessed, thus providing a benchmark for the development of sampling strategies.

Study Areas and GOCI Data
GOCI is the first geostationary ocean color satellite sensor, onboard the South Korean Communication, Ocean and Meteorological Satellite (COMS). It has successfully collected data covering geographical regions between 110 • E and 150 • E, and 2 • N and 50 • N, with a 500-m spatial resolution, since its launch in June 2010 [16,32]. GOCI was designed to measure six visible light bands (at 412, 443, 490, 555, 660 and 680 nm) and two near infrared bands (at 745 and 865 nm) with hourly observations during the daytime. GOCI aims to monitor the marine environment in real time with a rapid data-acquisition rate [33]. More than two-thirds of the coastal and sea waters within the Chinese maritime territory are imaged hourly by GOCI, which greatly increases the sampling frequency of China's seas. Figure 1 shows a red-green-blue (RGB) composite image from GOCI acquired on May 13, 2017, covering the China seas, including the Bohai Sea (BS), the Yellow Sea (YS) and the East China Sea (ECS), when there was massive cloud coverage across the region. The Bohai Sea is a semi-enclosed sea located Figure 1 shows a red-green-blue (RGB) composite image from GOCI acquired on May 13, 2017, covering the China seas, including the Bohai Sea (BS), the Yellow Sea (YS) and the East China Sea (ECS), when there was massive cloud coverage across the region. The Bohai Sea is a semi-enclosed sea located in the northeastern region of China, with an area of approximately 7700 km 2 and an average water depth of 18 m, extending 550 km from south to north and 300 km from west to east [34]. The Bohai Sea is dominated by an irregular semi-diurnal tide. The constituent of surface sediment is mainly silt and clayey silt, and the spatial variations of suspended sediment in the Bohai Sea was dominated by the river input, together with coastal resuspension caused by hydrodynamics [35,36]. The Yellow Sea flows to the Bohai Sea through the narrow Bohai straight in the east. The Yellow Sea has an area of approximately 380,000 km 2 with average depth of about 44 m [37]. In recent decades, abundant sediments were exported to the Bohai and Yellow Sea, caused by human activities together with the accumulative effects of winds, currents and the resuspension of particles [38]. Previous research showed that wind and tide are the principal agents for sediment resuspension [39,40]. The East China Sea is connected to the Yellow Sea at the north, with a continental shelf of approximately 640 km and an average water depth of 72 m [41]. This wide and shallow shelf receives a great amount of sediment from the Yangtze river, one of the largest rivers in the world. The annual sediment load from the Yangtze River is about 480 million tons [42]. GOCI obtains eight hourly measurements during daytime from 08:30 to 15:30 local time and provides unparalleled high-frequency data for close monitoring of the highly dynamic aquatic environment [43]. The GOCI L1B data products were downloaded from the website (http://kosc.kiost. ac.kr/eng/). High-frequency observations from 2013 to 2017, 17,120 images in total, were obtained to analyze the cloud coverage, interrupting and reducing frequency of valid samplings, and to analyze the associated uncertainty levels. Note that although GOCI data were evenly selected from different months, the averaged valid observations varied monthly, as listed in Table 1 for each sub-region. The valid pixels in each image were first determined using "Cloud or Ice" flags provided along with the GOCI products. Then the valid observations (%) for each month were calculated for each pixel as the ratio of the valid number of data to total images of the month. In general, more valid observations were available in summer than in winter for the East China Sea, but opposite trends were observed for the Bohai Sea and the Yellow Sea, with more valid data in winter.

Data Processing
GOCI Level-1 datasets containing calibrated-at-sensor radiance measurements were batch-processed using a module of the GOCI Data Processing System (GDPS) software package version 1.4.1 (released in December 2015) to generate spectral remote sensing reflectance (RRS) images. Atmospheric correction of the GOCI datasets was performed using an atmospheric correction algorithm similar to that used with the MODIS system. Details of the atmospheric correction scheme used with the GDPS software have already been published [44,45].
After atmospheric correction, two TSS algorithms were applied: the YOC algorithm and the Case-2 algorithm, to obtain the required TSS data products. The YOC algorithm [46] is an algorithm developed for the region and validated with extensive in situ measurements collected in the areas covered by GOCI, including the East China Sea and the Yellow Sea. For extremely turbid water (e.g., the estuarine part of the Yangtze River and the coastal area of the Yellow Sea, where TSS > 100 g/m 3 ), the YOC algorithm does not perform well. For this reason, the new Case-2 algorithm for turbid water was developed using in situ datasets obtained from the western coastal area of the Korean Peninsula [46,47]. Validation of the standard GOCI data products was carried out for both low and highly turbid waters, including the Bohai Sea, the Yellow Sea and the East China Sea [43,[48][49][50][51]. These studies showed that the GOCI ocean color data products are in close agreement with data from the MODIS and MERIS systems and are effective for monitoring water quality in these regions.
Cloud mask products derived using GDPS was used to analyze the cloud coverage in the target regions at both short-and long-term temporal scales. The algorithm for cloud masking is a threshold-based method using atmospheric corrected images. The sum of multiply scattered aerosol reflectance (ρ a ) and interactively scattered reflectance between aerosols and molecular (ρ ra ), was used as [ρ a (865nm) + ρ ra (865nm)], to determine whether a pixel is cloud or not. The default value of the threshold is 0.028 Sr −1 , and thus pixels with [ρ a (865nm) + ρ ra (865nm)] larger than 0.028 were defined as cloud pixels. To account for band registration errors or cloud movements, neighboring pixels were also masked [52].

Cloud Coverage and Total Suspended Sediment Uncertainties Analysis
(1) Cloud statistics: Monthly cloud coverage percentage (CCP) was first calculated for each pixel from five years of GOCI cloud mask data using the ratio of cloud coverage to the total numbers of observations. The spatial and temporal trends of monthly CCP were then obtained. (2) Estimation of TSS uncertainties: The mean normalized bias (P bias ) of TSS observations due to insufficient sampling rates was estimated using the statistical average of the absolute (unsigned) percentile differences between single observations at each hour h and the daily mean TSS from Remote Sens. 2020, 12,1945 6 of 26 eight observations. Note that P bias were only calculated for pixels that satisfy eight valid coverages per day. The value of P bias was calculated through the following equation: where N is the number of valid observations for each pixel from 2013 to 2017, and TSS(ave)i is the bias for the ith day observation: Similar indexes have also been used in previous studies to estimate the sampling bias of different satellite observations [18,25]. During the statistical evaluation of each day's results, pixels that returned fewer than 8 valid observations were neglected. The monthly and annual mean biases for TSS observations were then obtained by averaging the hourly bias over the month and the year.

Sensitivity Analysis of TSS Sampling Uncertainties
Two factors, cloud coverage and TSS variations, were assumed to influence the sampling bias of TSS in highly dynamic waters, and were quantified by coefficient of variation (CV = Standard Deviation/Mean) of each pixel across the entire time-series. An hourly or monthly CV grid was created by stacking all the TSS images per hour or month, and by calculating statistics for the entire image array with a CV value for each pixel. For a specific hour or month, the CV could be calculated as follows: where TSS i is the TSS of the i th observation, TSS is the average TSS for the study period, and N is the total number of observations. Larger CV values indicate that the data are more dispersed and TSS variability is greater over time, whereas smaller CV values mean that the data are more consistent and the TSS are more stable.
To analyze the effect of one factor within the influence of multiple factors, composite effects should be removed during the correlation analysis. The Partial Correlation analysis is an effective method for investigating the relationships between TSS bias and influencing factors, including cloud coverage and TSS variations. In addition, the P value (probability associated with significance) was used to test the significance of the correlations, where a P value < 0.05 is regarded as a significant correlation. The specific formula of partial correlation coefficient is as follows: R xy,z = r xy − r xz r yz where r xy is the Pearson correlation coefficient of variables x and y, x i is the TSS bias value of the i th observation, and y i is the cloud coverage or the CV of the i th observation. x and y are the mean values of the corresponding variables. R xy,z is the Partial correlation coefficient between x and y, with the effects of z removed.

Cloud Coverage
The spatial distribution and the temporal trend of monthly cloud coverage percentage (CCP) are shown in Figures 2 and 3. The studied areas demonstrated significant intra-annual variabilities in CCP, both spatially ( Figure 2) and temporally ( Figure 3). Averaged monthly CCP of the whole study area peaked in March at approximately 77.1%, decreased to the lowest level in July at approximately 52.3%, and then increased again to approximately 72.2% in November. This is consistent with previous results obtained by Xiao et al. [20], based on merged imagery of daily observations from the SeaWiFS, MODIS, MERIS and VIIRS systems. Obvious monthly CCP patterns were identified for sub-regions of the BS, the YS and the ECS area. The BS and YS areas had similar temporal patterns, with peak CCP in July and lower CCP in December, while an opposite pattern was shown in the ECS region, with low CCP in July and high CCP in December.

Cloud Coverage
The spatial distribution and the temporal trend of monthly cloud coverage percentage (CCP) are shown in Figures 2 and 3. The studied areas demonstrated significant intra-annual variabilities in CCP, both spatially ( Figure 2) and temporally ( Figure 3). Averaged monthly CCP of the whole study area peaked in March at approximately 77.1%, decreased to the lowest level in July at approximately 52.3%, and then increased again to approximately 72.2% in November. This is consistent with previous results obtained by Xiao et al. [20], based on merged imagery of daily observations from the SeaWiFS, MODIS, MERIS and VIIRS systems. Obvious monthly CCP patterns were identified for sub-regions of the BS, the YS and the ECS area. The BS and YS areas had similar temporal patterns, with peak CCP in July and lower CCP in December, while an opposite pattern was shown in the ECS region, with low CCP in July and high CCP in December.  The extremely high CCP for the China seas indicates that large portions of the observations from remote-sensing satellites are blinded by clouds. To substantiate this point, the lowest monthly CCP averages for the ECS, BS and YS were found to be 32.5%, 37.3% and 47.9%, equivalently reducing the number of days of cloud-free observations to 20.1, 18.6 and 16, respectively. In the months of the highest CCP, there were only 5, 4 and 5.7 cloud-free days for the ECS, BH and YS area, respectively. The average CCP (standard deviation) values were 62.6% (14.1%), 67.3% (15.5%) and 69.9% (8.6%) for the ECS, BS YS areas, respectively, and it was 65.7% (7.7%) for the whole area. Therefore, it is imperative to understand the effects of the cloud coverage on the remote sensing of TSS in these highly dynamic waters, and to assess the uncertainties in the observations from data from the Terra/Aqua MODI satellite, the most commonly used satellite sensors.  The extremely high CCP for the China seas indicates that large portions of the observations from remote-sensing satellites are blinded by clouds. To substantiate this point, the lowest monthly CCP averages for the ECS, BS and YS were found to be 32.5%, 37.3% and 47.9%, equivalently reducing the number of days of cloud-free observations to 20.1, 18.6 and 16, respectively. In the months of the highest CCP, there were only 5, 4 and 5.7 cloud-free days for the ECS, BH and YS area, respectively. The average CCP (standard deviation) values were 62.6% (14.1%), 67.3% (15.5%) and 69.9% (8.6%) for the ECS, BS YS areas, respectively, and it was 65.7% (7.7%) for the whole area. Therefore, it is imperative to understand the effects of the cloud coverage on the remote sensing of TSS in these highly dynamic waters, and to assess the uncertainties in the observations from data from the Terra/Aqua MODI satellite, the most commonly used satellite sensors.

Biases of TSS Observations
The spatial distribution map of hourly TSS biases was calculated and is presented in Figure 4, while the statistical results representing the temporal trends for the different regions are illustrated in Figure 5. TSS bias is the same as Pbias in the following text of the paper. Notable spatial and temporal variations were revealed throughout the Chinese seas. The average (by standard deviation) biases for hourly TSS estimates were 12.8% (2.7%), 7.0% (1.7%), 10.2% (2.3%) and 11.5% (2.5%), in the ECS, BS, YS, and the whole study area, respectively. The results revealed that the maximum TSS bias was as much as 18.3%, 9.8%, 14.5% and 16.5% for the aforementioned sea regions. The largest TSS bias was found in the ECS, while the smallest TSS bias occurred in the BS. The diurnal variation of TSS bias in the whole area shows a decreasing trend from morning, at around 08:30, to the early afternoon, around 13:30 local time. During this period the bias decreased from 13.3% to 8.9%. It then increased to 16.5% at 15:30. Similar diurnal trends were also observed in the ECS, BS and YS.

Biases of TSS Observations
The spatial distribution map of hourly TSS biases was calculated and is presented in Figure 4, while the statistical results representing the temporal trends for the different regions are illustrated in Figure 5. TSS bias is the same as P bias in the following text of the paper. Notable spatial and temporal variations were revealed throughout the Chinese seas. The average (by standard deviation) biases for hourly TSS estimates were 12.8% (2.7%), 7.0% (1.7%), 10.2% (2.3%) and 11.5% (2.5%), in the ECS, BS, YS, and the whole study area, respectively. The results revealed that the maximum TSS bias was as much as 18.3%, 9.8%, 14.5% and 16.5% for the aforementioned sea regions. The largest TSS bias was found in the ECS, while the smallest TSS bias occurred in the BS. The diurnal variation of TSS bias in the whole area shows a decreasing trend from morning, at around 08:30, to the early afternoon, around 13:30 local time. During this period the bias decreased from 13.3% to 8.9%. It then increased to 16.5% at 15:30. Similar diurnal trends were also observed in the ECS, BS and YS.     The spatial distributions and temporal trends of the monthly bias for TSS are shown in Figure 6 and Figure 7. The pixels that did not satisfy the criterion of eight valid observations a day were labeled as "Nodata" (black color in Figure 6) during the statistical analysis. In general, the average monthly bias (±standard deviation) for the TSS was 14.1% (±2.6%), 7.6% (±2.3%) and 12.2% (±4.3%) for the ECS, BS and YS areas respectively, and was 12.8% (±2.7%) for the whole study area.
The maximum monthly bias for TSS was 18.6% and 19.9%, for the ECS and YS in October, and so too was the maximum monthly bias for the whole study area, at 18.2% in October. The maximum monthly bias was 10.7% for the BS region in August. Among the three areas, the YS area had the largest while the BS area had the smallest monthly bias for TSS, and the BS area also showed a smaller overall monthly biases than the other two area. The spatial distributions and temporal trends of the monthly bias for TSS are shown in Figures 6 and 7. The pixels that did not satisfy the criterion of eight valid observations a day were labeled as "Nodata" (black color in Figure 6) during the statistical analysis. In general, the average monthly bias (±standard deviation) for the TSS was 14.1% (±2.6%), 7.6% (±2.3%) and 12.2% (±4.3%) for the ECS, BS and YS areas respectively, and was 12.8% (±2.7%) for the whole study area.     The changes in monthly bias of TTS, in the ECS and YS, to some extent mirrored that of the whole studied area, with relatively lower biases from May to August, and relatively higher ones from October to December. In comparison, the monthly TSS bias for the BS region displayed almost a reversed pattern, with relatively higher values from June to September, and relatively lower values from October to December and from January to April. The maximum monthly bias for TSS was 18.6% and 19.9%, for the ECS and YS in October, and so too was the maximum monthly bias for the whole study area, at 18.2% in October. The maximum monthly bias was 10.7% for the BS region in August. Among the three areas, the YS area had the largest while the BS area had the smallest monthly bias for TSS, and the BS area also showed a smaller overall monthly biases than the other two area.
The changes in monthly bias of TTS, in the ECS and YS, to some extent mirrored that of the whole studied area, with relatively lower biases from May to August, and relatively higher ones from October to December. In comparison, the monthly TSS bias for the BS region displayed almost a reversed pattern, with relatively higher values from June to September, and relatively lower values from October to December and from January to April.
The spatio-temporal variations in monthly TSS biases arise in a manner agreeable with that of the cloud coverage percentage (Figures 2, 3, 6 and 7). Lower bias values coincided, both temporally and spatially, with higher cloud coverage percentage, and vice versa. For areas with extremely high cloud coverage, such as the majority of the ECS in November and December and the coastal areas of the YS from January to July (Figure 2), the bias for TSS cannot be calculated due to a lack of valid observations ("Nodata" area in Figure 6). Consistently, the monthly CCP (the probability of cloud coverage) for these areas in the corresponding months ( Figure 2) was higher than 90%. The changes of monthly TSS bias for the YS area, however, could not be entirely attributed to the cloud coverage over this area, because the changes of monthly biases for TTS did not correlate with the changes of the monthly CCP all year around (Figures 3 and 5).

Factors Affecting TSS Sampling Uncertainties
Two factors were assumed to potentially influence the uncertainties of retrieved TSS values. The first one is the natural variation of TSS itself, as results would be more susceptible to sampling times and frequencies in waters with highly dynamic TSS. In waters with stable TSS, different sampling times and frequencies make no difference to the TSS estimates. In this study, variations of TSS were indicated using the coefficient of variation (CV) to normalize variations across the whole study region, including both highly turbid and relatively clear waters. The second factor is the cloud coverage, which affects the number of valid observations, and thus affects the variations of obtained TSS values over both the short and long term.
The hourly and monthly trends of TSS CVs across the study areas are presented in Figures 8-11. Significant spatial variations were revealed in the regions of concern at both hourly and monthly scales. At the hourly scale, larger CVs were calculated for offshore waters along the East China Sea (Figure 8), and the coastal zones of the Bohai and Yellow seas indicated higher temporal variations of TSS levels in these areas; they are dominated by tidal currents that control the sediment transport and resuspension [51,53]. Regions with higher TSS bias also displayed larger CVs, which indicates that most of the hourly TSS bias may be due to the TSS variations ( Figure 9). Overall, the temporal trends in the TSS CVs in each subregion were similar to those of the TSS bias, which decreased from morning, around 08:30, to around 13:30 local time, and then increased until 15:30, with the trough found at around 13:30. The magnitudes of the TSS CVs were also consistent with the magnitudes of the TSS bias for the subregions. More specifically, the East China Sea generated the largest TSS bias and CV values, followed by the Yellow Sea. The lowest TSS bias and CV levels were found in the Bohai Sea.
TSS values over both the short and long term.
The hourly and monthly trends of TSS CVs across the study areas are presented in Figures 8-11. Significant spatial variations were revealed in the regions of concern at both hourly and monthly scales. At the hourly scale, larger CVs were calculated for offshore waters along the East China Sea (Figures 8), and the coastal zones of the Bohai and Yellow seas indicated higher temporal variations of TSS levels in these areas; they are dominated by tidal currents that control the sediment transport and resuspension [51,53]. Regions with higher TSS bias also displayed larger CVs, which indicates that most of the hourly TSS bias may be due to the TSS variations (Figures 9). Overall, the temporal trends in the TSS CVs in each subregion were similar to those of the TSS bias, which decreased from morning, around 08:30, to around 13:30 local time, and then increased until 15:30, with the trough found at around 13:30. The magnitudes of the TSS CVs were also consistent with the magnitudes of the TSS bias for the subregions. More specifically, the East China Sea generated the largest TSS bias and CV values, followed by the Yellow Sea. The lowest TSS bias and CV levels were found in the Bohai Sea.  The monthly TSS CV values showed notable spatial and temporal variations (Figures 10 and 11). Since CV was determined as the ratio of standard deviation to mean value, areas that experienced high dynamic TSS variations would have higher CV values. Larger CV values were observed along the shorelines of the Yellow Sea and the East China Sea during the summer and autumn seasons, induced by river inflow in the rainy season, whereas the winter season displayed lower CVs. There are evident plumes where the CV differs from surrounding areas, between the coastal regions and the open oceans, which may have been caused by sediment resuspension driven by the combined effects of river discharges, winds, currents and waves in these zones [54]. Previous studies have demonstrated that suspended matter levels increased significantly from the coastal zones to the open oceans, which were blocked by the warm Taiwan and Kuroshio currents [55,56].  The results above suggests that both cloud coverage and TSS variations contributed to the TSS sampling biases, spatially and temporally, which indicate that areas with higher cloud coverage or TSS variation are prone to having greater sampling biases. For instance, the TSS biases of the Yellow Sea and East China sea in November, December and January were higher than in other months, which were consistent with the temporal trends of cloud coverages in these areas. But for the Bohai Sea, higher TSS biases were revealed in the summer, and lower ones revealed in the winter, and these  The results above suggests that both cloud coverage and TSS variations contributed to the TSS sampling biases, spatially and temporally, which indicate that areas with higher cloud coverage or TSS variation are prone to having greater sampling biases. For instance, the TSS biases of the Yellow Sea and East China sea in November, December and January were higher than in other months, which were consistent with the temporal trends of cloud coverages in these areas. But for the Bohai Sea, higher TSS biases were revealed in the summer, and lower ones revealed in the winter, and these The monthly TSS CV values showed notable spatial and temporal variations (Figures 10 and 11). Since CV was determined as the ratio of standard deviation to mean value, areas that experienced high dynamic TSS variations would have higher CV values. Larger CV values were observed along the shorelines of the Yellow Sea and the East China Sea during the summer and autumn seasons, induced by river inflow in the rainy season, whereas the winter season displayed lower CVs. There are evident plumes where the CV differs from surrounding areas, between the coastal regions and the open oceans, which may have been caused by sediment resuspension driven by the combined effects of river discharges, winds, currents and waves in these zones [54]. Previous studies have demonstrated that suspended matter levels increased significantly from the coastal zones to the open oceans, which were blocked by the warm Taiwan and Kuroshio currents [55,56].
The monthly biases and CVs for TSS did not support the hypothesis that the sampling bias was significantly affected by the TSS variations on a monthly scale for the entire study area. Areas with larger biases did not necessarily coincide with higher CV values, and the changing temporal patterns of the two parameters were not consistent either. These results demonstrate that the monthly observations of TSS levels are not solely determined by the TSS variations. Coastal regions with higher CVs were, however, prone to have larger TSS biases.
The results above suggests that both cloud coverage and TSS variations contributed to the TSS sampling biases, spatially and temporally, which indicate that areas with higher cloud coverage or TSS variation are prone to having greater sampling biases. For instance, the TSS biases of the Yellow Sea and East China sea in November, December and January were higher than in other months, which were consistent with the temporal trends of cloud coverages in these areas. But for the Bohai Sea, higher TSS biases were revealed in the summer, and lower ones revealed in the winter, and these coincident with cloud coverage trends in the area. Thus, the sensitivity of the sampling biases to the two factors should be further analyzed.
Maps of the partial correlation of TSS sampling bias with cloud coverage and with TSS CV were obtained at monthly scales, each month using 60 matchups. A correlation greater than 0.6 or 0.8 is generally described as a high or strong correlation, whereas a correlation less than 0.6 is generally described as a weak correlation. Figure 12 shows the correlations between the hourly TSS sampling biases and the cloud coverage ( Figure 12a) and the TSS CVs (Figure 12b), together with the P value maps (Figure 12c,d). In general, the correlation between the TSS sampling biases and variations was revealed to be higher than that for the relationship between sampling biases and the cloud coverage, at the hourly scale. The results are consistent with the analysis in Sections 3.2 and 4.1, which indicated that the spatial and temporal patterns of hourly scale TSS sampling biases are more closely correlated with the TSS variations than with cloud coverage. Maps of the partial correlation of TSS sampling bias with cloud coverage and with TSS CV were obtained at monthly scales, each month using 60 matchups. A correlation greater than 0.6 or 0.8 is generally described as a high or strong correlation, whereas a correlation less than 0.6 is generally described as a weak correlation. Figure 12 shows the correlations between the hourly TSS sampling biases and the cloud coverage ( Figure 12a) and the TSS CVs (Figure 12b), together with the P value maps (Figure 12c,d). In general, the correlation between the TSS sampling biases and variations was revealed to be higher than that for the relationship between sampling biases and the cloud coverage, at the hourly scale. The results are consistent with the analysis in Sections 3.2 and 4.1, which indicated that the spatial and temporal patterns of hourly scale TSS sampling biases are more closely correlated with the TSS variations than with cloud coverage. A strong correlation between sampling biases and TSS CVs was found in the Bohai Sea, Yellow Sea, and the northeastern section of the East China Sea, while the rest of the ECS showed a lower correlation, of less than 0.6. A total of 61.6% of the total study area was significantly affected (P < 0.05) by the TSS variations, while 47.6% was significantly affected by the cloud coverage. Table 2 presents A strong correlation between sampling biases and TSS CVs was found in the Bohai Sea, Yellow Sea, and the northeastern section of the East China Sea, while the rest of the ECS showed a lower correlation, of less than 0.6. A total of 61.6% of the total study area was significantly affected (P < 0.05) by the TSS variations, while 47.6% was significantly affected by the cloud coverage. Table 2 presents the area statistics for correlation, and P values for each subregion. The results clearly showed that the TSS sampling bias is more sensitive to TSS variation than to cloud coverage, for hourly scale observations. For instance, observations for approximately 87% of the Bohai Sea and 86% of the Yellow Sea were significantly affected by the TSS variations. On the contrary, 48.4% and 46.9% of these two subregions were significantly affected by cloud coverage. About half of the area of the East China sea was significantly affected by the TSS variations.  Figure 13 presents maps of correlation for P values between the monthly bias and the cloud coverage (Figure 13a,c), and between monthly bias and the TSS CV (Figure 13b,d). Both the cloud coverage and the TSS variations obviously generated uncertainties in obtaining TSS. For partial areas of the Bohai Gulf, the central Yellow Sea, and the southeastern area of the East China Sea, amounting to 17.1% of the studied areas, sampling bias strongly correlated with cloud coverage. In total, approximately 51.4% of the entire studied region had correlation coefficients higher than 0.6. Areas where TSS variations correlated more closely with sampling biases were found along the coastal regions and in the southeastern area of the East China Sea.  Figure 13 presents maps of correlation for P values between the monthly bias and the cloud coverage (Figure 13a,c), and between monthly bias and the TSS CV (Figure 13b,d). Both the cloud coverage and the TSS variations obviously generated uncertainties in obtaining TSS. For partial areas of the Bohai Gulf, the central Yellow Sea, and the southeastern area of the East China Sea, amounting to 17.1% of the studied areas, sampling bias strongly correlated with cloud coverage. In total, approximately 51.4% of the entire studied region had correlation coefficients higher than 0.6. Areas where TSS variations correlated more closely with sampling biases were found along the coastal regions and in the southeastern area of the East China Sea.   Secondly, the impact of the cloud coverage on sampling biases overshadows that of the TSS variations for monthly observations. At the P < 0.05 level, 17.6%, 10.9% and 29.5% of the TSS observations in the Bohai Sea, the East China Sea and the Yellow Sea were significantly affected by cloud coverage, respectively, while 4.4%, 6.9% and 4.5% of these seas, in respective order, were significantly affected by TSS variations (Table 3). It is worth noting that more than 50% of the observations displayed correlation coefficients larger than 0.4 for either cloud coverage or TSS variations. For example, 55% and 60% of the TSS observations in the East China Sea showed P values greater than 0.4, for cloud coverage and TSS variations, respectively.

Impacts of Sampling Strategy on Long-Term TSS Trends Monitoring
Time-series of remote sensing data are commonly used for long-term water quality estimations in the regions of interest. Therefore, the impacts of cloud coverage and sampling frequency on sampling uncertainties associated with long-term trends in TSS observations should be resolved. The differences among the TSS time-series retrieved from different satellite sources were examined to determine whether they are consistent regarding the magnitude, seasonality, inter-annual changes and long-term trends, for monthly TSS. To address this, monthly and annual TSS statistics from Terra/MODIS ("simulated" using GOCI 11:30 AM images), Aqua/MODIS (GOCI 13:30 PM images) and Terra/Aqua MODIS (combined observations using GOCI 11:30 and 13:30 images) observations were analyzed and compared to the eight observations of GOCI as references. Figure 14 provides the climatology monthly maps of GOCI-observed TSS from 2013 to 2017. The TSS concentrations are higher along the coasts, especially in the Bohai Sea and the Yangtze River estuary, while being lower in the middle of the Yellow Sea and the southeastern area of the East China Sea. From January to March, and October to December, significantly high TSS values (>~10 g/m 3 ) are observed in most regions of the study area, excluding the middle of the Yellow Sea and the southeastern area of the East China Sea. In general, TSS values are lowest in summer months (July to August), and increase in autumn and winter. These seasonal patterns are consistent with findings in Yuan et al. (2008), whereby the most turbid plumes were observed during winter, and decreased in the following spring and through summer [57]. Figure 15 presents the changes of monthly mean TSS from 2013 to 2017, retrieved from GOCI, Terra/MODIS, Aqua/MODIS and Terra/Aqua, for the entire study regions (ALL) and each subregion (ECS, YS, BS). In addition, the first derivative, the instant slope of the tangent line of the TSS-time curve, is presented. In this study, we defined "missed trends" to quantitatively compare the accuracy of different sensors in capturing long-term monthly variations of TSS. First, the baselines of monthly TSS trends were determined using GOCI observations, and the trend between each two months was calculated as the first derivative. A positive derivative signifies an increasing trend, and a negative derivative signifies a decreasing trend, of TSS between each month. Then, the monthly TSS trends from Terra/Aqua MODIS were also calculated and compared to the GOCI results. If the trends observed from GOCI and Terra/Aqua MODIS were opposing, then the numbers of missed trends were added. This process was checked for all months, and the total numbers of missed trends were obtained. Despite the analogous seasonal pattern, significant difference existed in these time-series of monthly TSS. estuary, while being lower in the middle of the Yellow Sea and the southeastern area of the East China Sea. From January to March, and October to December, significantly high TSS values (>~10 g/m 3 ) are observed in most regions of the study area, excluding the middle of the Yellow Sea and the southeastern area of the East China Sea. In general, TSS values are lowest in summer months (July to August), and increase in autumn and winter. These seasonal patterns are consistent with findings in Yuan et al. (2008), whereby the most turbid plumes were observed during winter, and decreased in the following spring and through summer [57].  Specifically, for the whole region, ECS, YS and BS, the ranges of monthly mean TSS obtained from eight GOCI observations were larger than those obtained from one Terra or Aqua/MODIS observation, as well as those from two Terra/Aqua combined observations. For the entire study region, GOCI observations showed a range of 4.47 g/m 3 (with the maximum and minimum TSS of 1.35 and 5.52 g/m 3 ), whereas the ranges were 2.96, 2.84 and 2.65 g/m 3 for Terra, Aqua and Terra/Aqua MODIS, respectively, as shown in Table 4. Similar differences were also observed in subregions of the ECS, YS and BS. This inconsistency indicates that single, or a limited number of, observations would miss information on variations of TSS during long-term monitoring, as substantiated by the smaller maximum to minimum ratios in comparison with the eight GOCI observations, as listed in Table 4. For extremely dynamic waters, as in the ECS [58], the maximum to minimum ratio obtained from GOCI was 16.57, almost five times larger than that for the single daily Terra or Aqua/MODIS observations. The effects of missing information, due to insufficient sampling frequency, on annual TSS variations were also observed, as shown in Figures 16 and 17. Larger variations indicated by standard deviation (STD) and CV were revealed from GOCI data, despite similar temporal trends and spatial distributions, compared to Terra/Aqua MODIS data. The higher CV from GOCI observations means more detailed TSS variations were observed from high-frequency remote sensing data. For instance, with eight observations a day, GOCI is able to capture TSS variations within each day, from 8:30 to 15:30, which Tera/Aqua MODIS would have missed. Moreover, the single or combined Terra and Aqua/MODIS seemed to overestimate annual trends of TSS for the entire areas and each subregion.  Figure 15 presents the changes of monthly mean TSS from 2013 to 2017, retrieved from GOCI, Terra/MODIS, Aqua/MODIS and Terra/Aqua, for the entire study regions (ALL) and each subregion (ECS, YS, BS). In addition, the first derivative, the instant slope of the tangent line of the TSS-time curve, is presented. In this study, we defined "missed trends" to quantitatively compare the accuracy of different sensors in capturing long-term monthly variations of TSS. First, the baselines of monthly TSS trends were determined using GOCI observations, and the trend between each two months was calculated as the first derivative. A positive derivative signifies an increasing trend, and a negative derivative signifies a decreasing trend, of TSS between each month. Then, the monthly TSS trends from Terra/Aqua MODIS were also calculated and compared to the GOCI results. If the trends observed from GOCI and Terra/Aqua MODIS were opposing, then the numbers of missed trends were added. This process was checked for all months, and the total numbers of missed trends were obtained. Despite the analogous seasonal pattern, significant difference existed in these time-series of monthly TSS. Table 4. Monthly statistics of TSS (g/m3) variations from GOCI, "simulated" Terra, Aqua and Terra + Aqua observations.   standard deviation (STD) and CV were revealed from GOCI data, despite similar temporal trends and spatial distributions, compared to Terra/Aqua MODIS data. The higher CV from GOCI observations means more detailed TSS variations were observed from high-frequency remote sensing data. For instance, with eight observations a day, GOCI is able to capture TSS variations within each day, from 8:30 to 15:30, which Tera/Aqua MODIS would have missed. Moreover, the single or combined Terra and Aqua/MODIS seemed to overestimate annual trends of TSS for the entire areas and each subregion.  Another concern when limited numbers of remote-sensing observations were applied for longterm TSS monitoring was the uncaptured tendencies of changes in the time-series, from Terra, Aqua or Terra/Aqua MODIS, at both monthly and annual scales (Figures 16 and 17). The tendency of changes between successive months or years could be calculated by the first derivative or slope, as shown in Figure 17c-h. Differences between derivatives, from Terra, Aqua or Terra/Aqua MODIS time-series, with regard to TSS and the reference GOCI time-series suggested that the tendencies of changes might be omitted for the Terra/Aqua MODIS systems. Since reference GOCI datasets were Another concern when limited numbers of remote-sensing observations were applied for long-term TSS monitoring was the uncaptured tendencies of changes in the time-series, from Terra, Aqua or Remote Sens. 2020, 12, 1945 19 of 26 Terra/Aqua MODIS, at both monthly and annual scales (Figures 16 and 17). The tendency of changes between successive months or years could be calculated by the first derivative or slope, as shown in Figure 17c-h. Differences between derivatives, from Terra, Aqua or Terra/Aqua MODIS time-series, with regard to TSS and the reference GOCI time-series suggested that the tendencies of changes might be omitted for the Terra/Aqua MODIS systems. Since reference GOCI datasets were free of the impacts of different data source, complications resulted from different calibrations, atmospheric corrections, sensitivity drifts or means of data processing. This is evidenced by Figure 15, the "missed trends" for each region of interest, as comparisons of different time-series to the GOCI reference. In general, approximately 10 monthly tendencies on average were missed in the Terra, Aqua or Terra/Aqua MODIS observations, with the maximum number of 13, and as a consequence, 16.7% of the monthly variations were not captured from 2013 to 2017. Similar results were found at the annual scale, for instance, while an increasing trend from 2013 to 2014 was revealed for the ECS in GOCI observations, a stable or slight decreasing trend was obtained in Terra/Aqua MODIS dataset. Therefore, caution should be taken when limited numbers of observations are used for the monitoring of long-term trends in the water quality parameter, a common approach in ocean color remote-sensing.

Advantages of High-Frequency Observations Compared to Conventional Terra/Aqua MODIS
Ocean color sensors, such as the Terra/Aqua MODIS system, have successfully been used to monitor water quality in inland, coastal and open ocean waters. Previous studies have demonstrated that the observations are often heavily hampered by clouds, and as a result, the average valid MODIS image retrieval rate over the global oceans was only 12% per day [10,19,20]. While some research effort has been dedicated to improving the temporal coverage by merging multiple data sources, systematic assessments of the impact of cloud coverage on MODIS observations, and the knock-on effect on water quality monitoring, remains rare, and there is an outstanding requirement for more of such studies.
This study further revealed that cloud coverage for China's seas was extremely high, for both the Terra/MODIS and Aqua/MODIS systems, with percentages of 67% and 63.2%, respectively, as listed in Table 5. Cloud coverage levels from the Terra/MODIS system were generally higher than afternoon observations from the Aqua/MODIS system, for each sub region as well as for the entire study area. The results are also consistent with previous studies on valid ocean observation rates from the Terra/Aqua MODIS systems [21]. These showed that there were more valid Aqua observations than Terra observations, due to lower cloud coverage in the afternoons. Monthly averaged data products are one of the key resources used for building time-series datasets for the analysis of ocean color parameters, including suspended sediments, chlorophyll and euphotic depth. One issue with using remote sensing data to generate high-level monthly and annually averaged datasets is that they may not reflect actual trends if there are too few valid observations [25]. The mean relative bias in monthly mean observations from the Terra/Aqua MODIS system was affected by cloud coverage. These were compared to the GOCI datasets that included eight observations per day. For the entire study area, the relative biases ranged from 9.5% to 20% for the Terra/MODIS system. For the Aqua/MODIS system, biases from 6.6% to 14.6% were revealed, with the maximum bias found in October. Larger biases were displayed in the Terra MODIS observations than the Aqua MODIS observations, and these may have been caused by more cloud coverage, as shown above.
Cloud coverage on remote sensing images would greatly reduce valid observations, and thus introduce sampling uncertainties in short-or long-term water quality monitoring. As revealed in Section 3.2, the average biases of monthly TSS estimates were 14.1%, 7.6% and 12.2% in the ECS, BS and YS, respectively. Such results were comparable with previous studies, such as a global root mean-square error (RMSE) of 8% observed for annual mean chlorophyll concentrations, due to missing data from CZCS [18]. Barnes and Hu demonstrated that the difference between monthly mean Chl-a products for MODIS and SeaWiFS could be larger than 10% if monthly valid observations are less than 15 in the Gulf of Mexico region [27], which means the monthly sampling frequency should be larger than 50%. In a similar study conducted by Xiao Yanfang et al., a bias greater than 15% could be induced in the Yellow and East China Sea if the monthly data are averaged from 3 days, and 30% bias if the monthly data are from 1 day only [20]. However, one drawback of previous studies is that both the analysis and results were obtained using daily merged products, which means diurnal variations were missed.
Another concern when lacking effective observations is that it reduces the quality of time-series ocean color products. For instance, the monthly and annual TSS products form Terra/MODIS or Aqua/MODIS observations mostly seem to be missing the maximum and minimum values, compared to the GOCI products. The entire study region showed a range of 4.47 g/m 3 when using GOCI, with maximum and minimum TSS of 1.35 and 5.52 g/m 3 , whereas the ranges were 2.96, 2.84 and 2.65 g/m 3 from Terra, Aqua and Terra/Aqua MODIS, respectively, as shown in Figure 15 and Table 4. The benefits of high-frequency GOCI data include the improved ability to capture water quality variations more quickly, in cases a few hours before or after a large cloud/rainstorm has passed over the area and sediment is transported out into these waters. Therefore, stationary satellite remote sensing is preferred for water quality monitoring, especially in highly dynamic water regions.

Implications for Future in situ and Remote Sensing Sampling Strategies
Previous studies suggested that for observations of highly dynamic waters, sampling intervals of 30 min to 1 h should preferably be taken [24], and inappropriate sampling frequencies and times can lead to large biases, and could eventually lead to unscientific decisions being made during water management. A previous study has indicated that inappropriate sampling time and frequency could lead to statistical errors larger than 50% during water quality monitoring [7,25,26]. The results presented here demonstrate that while ocean color sensors, including polar-orbit satellite sensors, such as the Terra/Aqua MODIS system, are regularly used for water quality monitoring, the impacts of cloud coverage on valid observation rates and sampling biases should be taken into consideration. This is more critical when evaluating regions where the water quality undergoes more complex variations, which may be obscured by larger sampling biases. It is important to develop an optimal sampling strategy that incorporates both in situ and remote sensing observations, in order to minimize the effects of cloud coverage, and to obtain more valid observations or matchups for system modeling and validation of the remote sensing products.
Since geostationary ocean color satellite coverage is currently not available for most regions, and will not be for the foreseeable future, polar-orbit satellites are still the primary data source for ocean color research and applications. More polar-orbit ocean color satellite missions are currently scheduled to be launched by space agencies around the world. For instance, three ocean color missions-the Haiyang (HY) series, including the HY-1D and the HY-1E/F-are planned for launch in 2019 and 2021, respectively, by the Chinese National Space Administration (CNSA). These will contribute to surveys of the environment in China's offshore waters, islands and coastal zones, as well as performing research on the global oceanic environment. Other missions include the PACE system from NASA in the USA, and the EnMAP from the DLR, Germany. The results of this study could, therefore, serve as a guide to avoiding cloud coverage and maximizing valid observations. The number of hours corresponding to minimum TSS sampling biases were calculated and are shown in Figure 18 and Table 6. Notable variance was found at different spatial locations and in different months. For the entire region and the Bohai Sea, the most appropriate sampling time window ranges from 12:30 to 15:30, while this time window narrowed from 13:30 to 15:30 for observations in the East China Sea and the Yellow Sea.
Previous studies suggested that for observations of highly dynamic waters, sampling intervals of 30 minutes to 1 hour should preferably be taken [24], and inappropriate sampling frequencies and times can lead to large biases, and could eventually lead to unscientific decisions being made during water management. A previous study has indicated that inappropriate sampling time and frequency could lead to statistical errors larger than 50% during water quality monitoring [7,25,26]. The results presented here demonstrate that while ocean color sensors, including polar-orbit satellite sensors, such as the Terra/Aqua MODIS system, are regularly used for water quality monitoring, the impacts of cloud coverage on valid observation rates and sampling biases should be taken into consideration. This is more critical when evaluating regions where the water quality undergoes more complex variations, which may be obscured by larger sampling biases. It is important to develop an optimal sampling strategy that incorporates both in situ and remote sensing observations, in order to minimize the effects of cloud coverage, and to obtain more valid observations or matchups for system modeling and validation of the remote sensing products.   Table 6. Hours corresponding to minimum TSS sampling bias and cloud coverage for each region. The sampling strategy or temporal coverage of satellite sensors is designed based on a compromise between orbital characteristics, the Field of View (FOV) or maximum scan angle [28], and cloud coverage information. The information set out in this article shows that information regarding the profiles and significance of water quality parameters, such as TSS, can also be applied to the design of future ocean color missions. This information can be used for informing scheduling studies that guide the revisit time windows, as well as for field data collection strategies, remote sensing modeling, and data product validation.

Hours of Minimum TSS Bias Errors Hours of Minimum Cloud Coverage
Although the results of this study are region-and sensor-dependent, it is feasible to adapt the approach to other areas in other time zones, as well as other water quality parameters, including chlorophyll, CDOM, etc. A new generation of remote sensing data is now becoming available, especially from optical remote sensors on geostationary satellites, which would provide a much higher temporal resolution-typically an image once or more per hour during daylight, compared to once per day, for instance, for the PACE system from NASA in the USA and the EnMAP from the DLR, Germany, as well as the GOCI II from Korea. These new ocean color sensors with higher temporal resolutions open up opportunities for dramatically improving the data availability, for resolving fast processes such as tidal or diurnal variability of phytoplankton or suspended particulate matter. Moreover, the results of this study will also benefit the design and application of polar-orbit ocean color satellite missions, such as China HY-1D and the HY-1E/F, planned for launch in 2019 and 2021, and Sentinel 3 OLCI, etc., which are designed to determine the time window of specific regions, or to evaluate the long-term uncertainties of ocean color products caused by insufficient sampling frequency or cloud coverage.

Conclusions
The benefits of ocean color remote sensing in scientific research and quantitative applications lie in their ability to provide high-resolution spatial and temporal observations of the earth's surface. This study, however, reveals that insufficient observations caused by cloud coverage have significant impacts on both the short-and long-term monitoring of ocean color parameters, especially for highly dynamic waters in coastal regions. The effects of the cloud coverage dominated at the monthly scale, with the mean normalized bias (P bias ) for TSS measurements of 14.1% (±2.6%), 7.6% (±2.3%) and 12.2% (±4.3%) in the ECS, BS and YS, respectively. Although satellite sensors, including the Terra/Aqua MODIS system, have been widely used for ocean color research and applications, uncertainties due to insufficient observations were evidently revealed for monthly and annual scale products. Evaluations of the effects of insufficient sampling observations made by the Terra/Aqua MODIS systems showed that P bias ranged from 6.5% (±7.4%) to 20% (±13.1%) with the monthly scale TSS products. Much information would be missed, such as maximum values and variation ranges, when limited remote sensing data are adopted, and significant "missing trends" effects were revealed in monthly TSS variations, which means that approximately 16.7% of monthly variations were not captured by Terra/Aqua MODIS during the periods. For the entire region and the Bohai Sea, the most appropriate sampling time window ranges from 12:30 to 15:30, while this time window narrowed from 13:30 to 15:30 for observations in the East China Sea and the Yellow Sea. Caution should, therefore, be taken when using remote sensing datasets to obtain higher level processed data products, like monthly averaged values, particularly if there are an inadequate number of valid observations. The results of this research could be applied to analyze the uncertainties in ocean color remote sensing, and to determine the appropriate sampling times for collecting ocean color information, as well as providing a benchmark for future satellite sensor design and mission planning.

Conflicts of Interest:
The authors declare no conflict of interest.