Spectral Classification of the Yellow Sea and Implications for Coastal Ocean Color Remote Sensing

Remote sensing reflectance (Rrs) classification of coastal waters is a useful tool to monitor environmental processes and manage marine environmental resources. This study presents classification work for data sets that were collected in the Yellow Sea during six cruises (spring and autumn, 2003; summer and winter, 2006/2007; and spring and autumn, 2007). Specifically, we analyzed classification features of Rrs spectra and obtained spatio-temporal characteristics of reflectance and bio-optical properties in the coastal waters. Yellow Sea waters were classified into the following four typical regions based on their spatial distribution characteristics: middle of the Yellow Sea (MYS), north Yellow Sea (NYS), coastal Shandong (CS), and Jiangsu shoal (JS), and five water type categories consisting of Classes A–E were used to represent water colors from clear to very turbid. Application of this classification scheme to Medium Resolution Imaging Spectrometer (MERIS) imagery revealed seasonal variations in the data, which suggests that the water types have both significant temporal and spatial distributions. In particular, the area of Class E waters in the Jiangsu shoal tended to gradually shrink in summer and expand in winter. The spatio-temporal variability was due to the influence of various environmental factors such as currents, tidal activity, fresh water discharges, monsoon winds, and typhoons.


Introduction
Coastal and shelf zones occupy about 18% of the surface of the globe, and these areas supply about 90% of global fishery catches and account for some 25% of global marine primary production [1].Because of the heterogeneity of water contents, the diversity of inputs, variations in anthropogenic impacts, and the variability of physical forcing factors in shelf seas [2], water quality can vary dramatically among different regions.Poor water quality in coastal and shelf areas can lead to break outs of red tides [3], green tides [4], and depleted fisheries [5].In order to monitor environmental processes and manage coastal and shelf resources, integrated measurements from in-situ and satellite sensors can be used.Specifically, these data can provide valuable spatial and temporal information about the state of coastal waters for decision makers, which can facilitate the implementation of sustainable actions in the fields of ecology, fisheries, tourism, marine transportation, offshore engineering, and so forth [6].
Numerous high frequency processes occur in coastal waters, and the bio-optical properties in such zones often span a wide range of variability [2,7]; this complexity can make the use of regional algorithms difficult.A feasible approach to analyze the spatial and temporal distribution characteristics of complex coastal waters can be achieved based on optical classifications, which aim to group waters with similar optical traits.The definition of optical water types implicitly assumes that different coastal regions can present similar or different optical characteristics in regard to their marine components, and the bio-optical properties can be evaluated by using spectra reflectance data.Based on these fundamental principles, various researchers have devoted large amounts of time over the past several years to the development of optical data sets that can be used in coastal waterways to identify and solve problems.Roesler and Perry [8] compiled measured irradiance spectra for estuarine, fjord, coastal, and oceanic water types and developed an inverse model for the in-situ spectral reflectance.Arnone [9] used remote sensing reflectance (R rs ) spectra of coastal waters and individual constituents of the water to divide the water masses into six types.Lubac and Loisel [10] applied an unsupervised hierarchical cluster analysis to a data set of normalized R rs spectra and identified five spectrally distinct classes, including three obvious classes involved in bloom situations and two water masses dominated by mineral and non-living particles.Vantrepotte et al. [11] obtained four optical classes through the use of normalized in-situ reflectance spectra and employed statistical analyses to emphasize the shape of the reflectance rather than the associated magnitude data.Moore et al. [12] applied a fuzzy classification approach to NOMAD (NASA Bio-Optical Marine Algorithm Data Set, where NASA stands for the National Aeronautics and Space Administration) remote sensing reflectance and identified eight optical water types for which the uncertainties of the Moderate Resolution Imaging Spectroradiometer (MODIS) ocean chlorophyll product have been calculated.Mélin and Vantrepotte [13] applied an unsupervised clustering technique to a global seven-year SeaWiFS data set and derived a set of 16 classes that covered conditions from very turbid to oligotrophic waters.Other bio-optical parameters have also been applied to the classification of coastal waters, such as, diffuse attenuation coefficients for downwelling plane irradiance (K d ), which resulted in the derivation of 9 classes [14], ternary plots of absorption coefficients from detritus, chlorophyll, and colored dissolved organic matter (CDOM), which resulted in the derivation of four classes [15], and concentrations of chlorophyll a (Chl a) from NOMAD, which resulted in the derivation of 12 classes [16].The discrete Forel-Ule scale has also been used [17].
In regard to the application of such techniques in specific coastal regions, such as the Yellow Sea, various classification methods can be difficult to apply because of problems associated with unsuitable or deficient data.Previous studies have flagged large areas of the Yellow Sea close to river outflows as "no data areas" in regard to Medium Resolution Imaging Spectrometer (MERIS) products [17].Furthermore, in the studied areas, the standard cloud mask algorithm used with SeaWiFS data usually fails since turbid water pixels are generally classified as clouds, thus leading to the loss of data [11].In highly turbid coastal waters, it is possible for dynamic ranges to be exceeded, such as when MODIS bands saturate and the true signal is unknown [18].Additionally, there has been a general lack of in-situ data set training and validation studies in the Yellow Sea, which makes it difficult to apply many of the most common classification methods for coastal waters, i.e., this can lead to too few or too many classification types.Thus, new studies that focus on characterizing the optical diversity of coastal water masses from in-situ measurements, satellite applications of coastal and ocean optical data for monitoring coastal water quality, and improvements in ocean color products would be valuable for the Yellow Sea.
The aim of this paper is to characterize and investigate the mechanisms of the spatio-temporal variability of remote sensing reflectance spectra in the Yellow Sea, which contains diverse water types ranging from very turbid to clear waters as well as a wide range of bio-optical properties.In-situ R rs (λ) spectra data sets allowed us to classify these waters into different classes, which differed significantly in regard to their water constituents and characteristic optical properties.The spectra types were applied to MERIS sensor's products in order to assess the spatio-temporal distribution of the classes, and therefore, the potential applicability of using optical class-based approaches to analyze inter-annual changes in coastal turbid water masses.The development of optical water type classification approaches based on R rs (λ) spectra is important to the study of the spatial and temporal dynamics of ecologically and biogeochemically important coastal regions.

Study Area
The Yellow Sea is semi-enclosed shelf sea surrounded by the mainland of China (Liaoning Peninsula, Shandong Peninsula, and Jiangsu province), the Korean Peninsula (119 ˝-125 ˝E, 31 ˝-40 ˝N), the Bohai Sea in the west, and the East China Sea in the south (Figure 1).The Yellow Sea is a shallow basin with a mean depth of 44 m.A trough with a maximum depth of 152 m lies in the center.North and northwest winds in the autumn and winter are strong and typhoons are common; about nine typhoons pass the region every year, and these bring heavy rains in the summer and autumn [19].The Changjiang River, or Yangtze River, is the largest discharge source into the Yellow Sea and East China Sea, and the discharge contains large amounts of suspended sediments; hence, coastal waters near the estuary have extremely high turbidity values.Surface circulation in the Yellow Sea and East China Sea is driven by three main currents including the Kuroshio, the Tsushima Current, and the Yellow Sea Warm Current (YSWC) [20].
Remote Sens. 2016, 8, 321 3 of 24 approaches based on Rrs(λ) spectra is important to the study of the spatial and temporal dynamics of ecologically and biogeochemically important coastal regions.

Study Area
The Yellow Sea is semi-enclosed shelf sea surrounded by the mainland of China (Liaoning Peninsula, Shandong Peninsula, and Jiangsu province), the Korean Peninsula (119°-125° E, 31°-40° N), the Bohai Sea in the west, and the East China Sea in the south (Figure 1).The Yellow Sea is a shallow basin with a mean depth of 44 m.A trough with a maximum depth of 152 m lies in the center.North and northwest winds in the autumn and winter are strong and typhoons are common; about nine typhoons pass the region every year, and these bring heavy rains in the summer and autumn [19].The Changjiang River, or Yangtze River, is the largest discharge source into the Yellow Sea and East China Sea, and the discharge contains large amounts of suspended sediments; hence, coastal waters near the estuary have extremely high turbidity values.Surface circulation in the Yellow Sea and East China Sea is driven by three main currents including the Kuroshio, the Tsushima Current, and the Yellow Sea Warm Current (YSWC) [20].1).During the cruises, all bio-optical parameters and ancillary data were obtained strictly according to NASA's ocean optics protocols for satellite ocean color sensor validation [21].The data sets and adopted   1).During the cruises, all bio-optical parameters and ancillary data were obtained strictly according to NASA's ocean optics protocols for satellite ocean color sensor validation [21].The data sets and adopted measurement methods for the remote sensing reflectance, R rs , diffuse attenuation coefficient for downwelling plane irradiance, K d , total absorption coefficient, a, absorption by CDOM, a g , absorption by particles, a p , absorption by detritus, a d , absorption by phytoplankton pigments, a phy , scattering coefficient, b, backscattering coefficient, b b , total beam attenuation coefficient, c, particle beam attenuation, c p , suspended particle matter (SPM), Chl a and accessory pigment concentrations, Secchi depth and Forel-Ule value have been described systematically by Tang et al. [22], Zhang et al. [23], and Wang et al. [24].This work is briefly reviewed below.
R rs (λ) was measured following the water methods described by Mueller et al. [21] through the use of an ASD FieldSpec 931 spectroradiometer with 512 bands, a spectral resolution of 3 nm at 700 nm, and a sampling interval of 1.4 nm over the spectral range of 350-1050 nm.The radiometric stability of the instrument was tracked in-situ during the cruises of 2003 with the SeaWiFS Quality Monitor (Satlantic, Inc., SQM-II); the instrument was carried back periodically for calibration with a NIST (National Institute of Standards and Technology) reference lamp in the laboratory during the cruises of 2006 and 2007.Radiance spectra were collected 15 times (span of one more than the wave period) each for a reference panel, the water, and the sky.All procedures were repeated three times at each station and then data were examined visually to eliminate abnormal spectra; measurements were repeated if necessary.Averaged data were used in the calculations.R rs (λ) (Figure 2) was derived according to Equation (1): R rs pλq " pL t pλq ´r ˆLsky pλqq{pL p pλq ˆπ{ρ p q where L t (λ) is the total upwelling spectral radiance above the water surface; r ˆLsky (λ) is the direct upwelling radiance reflected on the water surface contributed by the sky; r is calculated from the Fresnel formula; and L p (λ) is the simultaneously observed radiance of the reference panel, which has an accurately calibrated reflectance, ρ p , of approximately 35%.The R rs (λ) measurement standard deviations were between 2% and 4%.
Remote Sens. 2016, 8, 321 4 of 24 measurement methods for the remote sensing reflectance, Rrs, diffuse attenuation coefficient for downwelling plane irradiance, Kd, total absorption coefficient, a, absorption by CDOM, ag, absorption by particles, ap, absorption by detritus, ad, absorption by phytoplankton pigments, aphy, scattering coefficient, b, backscattering coefficient, bb, total beam attenuation coefficient, c, particle beam attenuation, cp, suspended particle matter (SPM), Chl a and accessory pigment concentrations, Secchi depth and Forel-Ule value have been described systematically by Tang et al. [22], Zhang et al. [23], and Wang et al. [24].This work is briefly reviewed below.Rrs(λ) was measured following the water methods described by Mueller et al. [21] through the use of an ASD FieldSpec 931 spectroradiometer with 512 bands, a spectral resolution of 3 nm at 700 nm, and a sampling interval of 1.4 nm over the spectral range of 350-1050 nm.The radiometric stability of the instrument was tracked in-situ during the cruises of 2003 with the SeaWiFS Quality Monitor (Satlantic, Inc., SQM-II); the instrument was carried back periodically for calibration with a NIST (National Institute of Standards and Technology) reference lamp in the laboratory during the cruises of 2006 and 2007.Radiance spectra were collected 15 times (span of one more than the wave period) each for a reference panel, the water, and the sky.All procedures were repeated three times at each station and then data were examined visually to eliminate abnormal spectra; measurements were repeated if necessary.Averaged data were used in the calculations.Rrs(λ) (Figure 2) was derived according to Equation (1): where Lt(λ) is the total upwelling spectral radiance above the water surface; r × Lsky(λ) is the direct upwelling radiance reflected on the water surface contributed by the sky; r is calculated from the Fresnel formula; and Lp(λ) is the simultaneously observed radiance of the reference panel, which has an accurately calibrated reflectance, ρp, of approximately 35%.The Rrs(λ) measurement standard deviations were between 2% and 4%.For SPM measurements, the water samples were filtered through previously weighted 47 mm diameter filters (Whatman GF/F filters, pore size of 0.45 μm).The filters were then dried and reweighed.Chl a concentrations were measured by using a spectrophotometer after the water samples had been filtered through 47 mm diameter filters (membrane filter, pore size of 0.45 μm).Chl a data were based on in vitro fluorescence measurements following dark 24 h extractions in acetone, and quantifications were performed via both high performance liquid chromatography (HPLC) and fluorometric methods.The HPLC methods were applied to separate the different types of pigments [25].The HPLC derived total Chl a (TChl a) values in the data sets represent the sum of divinyl chlorophyll a plus chlorophyllide a.The absorptions of constituents such as a, ap, aphy, ad, and For SPM measurements, the water samples were filtered through previously weighted 47 mm diameter filters (Whatman GF/F filters, pore size of 0.45 µm).The filters were then dried and reweighed.Chl a concentrations were measured by using a spectrophotometer after the water samples had been filtered through 47 mm diameter filters (membrane filter, pore size of 0.45 µm).Chl a data were based on in vitro fluorescence measurements following dark 24 h extractions in acetone, and quantifications were performed via both high performance liquid chromatography (HPLC) and fluorometric methods.The HPLC methods were applied to separate the different types of pigments [25].The HPLC derived total Chl a (TChl a) values in the data sets represent the sum of divinyl chlorophyll a plus chlorophyllide a.The absorptions of constituents such as a, a p , a phy , a d , and a g were measured with the quantitative filter technique (QFT) by using a GBC UV/VIS Cintra20 spectrophotometer [26].The AC9 sensor was used to measure a and c, and it was regularly calibrated with pure (Milli-Q) water, i.e., before, during, and after the Yellow Sea field campaigns.The raw absorption and attenuation spectral values were corrected for temperature and salinity effects following the recommended procedure [27] for the visible and near infrared spectral regions.The proportional method was then used to correct for residual scattering effects on absorption measurements [28].Hydroscat 6 (HS6) was used to measure b b , and in-situ data were calibrated with deionized water and corrected by the sigma method to eliminate the underestimates of scattering in turbid waters [29].The AC9 and HS6 instruments were bonded together and used to take profiles of the water column from the ship.Transparency (Secchi depth) and Forel-Ule water color were measured by a Secchi disc and the Forel-Ule color scale, respectively.K d (λ) was determined by a Satlantic profiling radiometer system (Satlantic, Inc.).The downwelling irradiance profiles were used to calculate the K d (λ) by applying an exponential fit over the depth range from the subsurface (Z = 0 ´) to the penetration depth Z 90 (λ), which was defined as the depth at which the downwelling irradiance decreased to e ´1 of its value at the surface [30].

MERIS Image L2R Remote Sensing Reflectance Data
The remote sensing reflectance of MERIS Level 2 (L2) products can be directly acquired from the ESA website (http://mer-merci-uk.eo.esa.int/merci/), which were retrieved by the neural network (NN) algorithm as the Case II water algorithm for atmospheric correction [31].The software of BEAM was used for visualization of the MERIS L2 products, and Interface Description Language (IDL) was programed to batch process, which included radiometric correction, geometric correction, band synthesis, project changes, and MERIS imagery reflectance classification.
with pure (Milli-Q) water, i.e., before, during, and after the Yellow Sea field campaigns.The raw absorption and attenuation spectral values were corrected for temperature and salinity effects following the recommended procedure [27] for the visible and near infrared spectral regions.The proportional method was then used to correct for residual scattering effects on absorption measurements [28].Hydroscat 6 (HS6) was used to measure bb, and in-situ data were calibrated with deionized water and corrected by the sigma method to eliminate the underestimates of scattering in turbid waters [29].The AC9 and HS6 instruments were bonded together and used to take profiles of the water column from the ship.Transparency (Secchi depth) and Forel-Ule water color were measured by a Secchi disc and the Forel-Ule color scale, respectively.Kd(λ) was determined by a Satlantic profiling radiometer system (Satlantic, Inc.).The downwelling irradiance profiles were used to calculate the Kd(λ) by applying an exponential fit over the depth range from the subsurface (Z = 0 − ) to the penetration depth Z90(λ), which was defined as the depth at which the downwelling irradiance decreased to e −1 of its value at the surface [30].

MERIS Image L2R Remote Sensing Reflectance Data
The remote sensing reflectance of MERIS Level 2 (L2) products can be directly acquired from the ESA website (http://mer-merci-uk.eo.esa.int/merci/), which were retrieved by the neural network (NN) algorithm as the Case II water algorithm for atmospheric correction [31].The software of BEAM was used for visualization of the MERIS L2 products, and Interface Description Language (IDL) was programed to batch process, which included radiometric correction, geometric correction, band synthesis, project changes, and MERIS imagery reflectance classification.
MERIS imagery L2R data were available for the study area from 30 April 2002 to 8 April 2012, and data were collected about 1/2-1 times per day.Figure 3 (left) shows the pixels spectra from one cloudless MERIS 2P product image (16 April 2003) corresponding to the data from all field stations sampled during the spring cruise in 2003 without use of a synchronous time window.The results showed that several pixel spectra had negative values as result of over correcting for the atmosphere in very turbid water regions [17].Figure 3 (right) shows the strictly matched-up validation results [32] based on 3 × 3 pixel boxes, 0.15 satellite coefficient of variation thresholds and less than ± 3 h between the 2003 cruise's in-situ sets and MERIS data.The data were spread around the 1:1 line, with Rrs overestimated in the blue bands and underestimated in the red bands.The median of absolute percentage difference ranged from 62% to 15% in the visible bands of 412 and 665 nm, respectively.The results also demonstrate that the shape and value of MERIS pixel spectra were similar to the field spectra.MERIS imagery L2R data were available for the study area from 30 April 2002 to 8 April 2012, and data were collected about 1/2-1 times per day.Figure 3 (left) shows the pixels spectra from one cloudless MERIS 2P product image (16 April 2003) corresponding to the data from all field stations sampled during the spring cruise in 2003 without use of a synchronous time window.The results showed that several pixel spectra had negative values as result of over correcting for the atmosphere in very turbid water regions [17].Figure 3 (right) shows the strictly matched-up validation results [32] based on 3 ˆ3 pixel boxes, 0.15 satellite coefficient of variation thresholds and less than ˘3 h between the 2003 cruise's in-situ sets and MERIS data.The data were spread around the 1:1 line, with R rs overestimated in the blue bands and underestimated in the red bands.The median of absolute percentage difference ranged from 62% to 15% in the visible bands of 412 and 665 nm, respectively.The results also demonstrate that the shape and value of MERIS pixel spectra were similar to the field spectra.

Method for the Classification of In-Situ R rs Spectra
The Yellow Sea (YS) contains typical turbid class II waters, and the peak of R rs spectra is to a great extent related to the dominant components and concentrations in the water.This was especially obvious in waters dominated by suspended sediments; see the appearance of the red shift with increasing concentrations (Figure 2 red arrow).Maximum-minimum values for supervised classification (referred to in shortened form as max-classification) [33] of R rs spectra were adopted for the whole in-situ data set in order to divide the data into homogeneous optical groups.One key advantage of the max-classification method is that it is sensitive to the shape of the reflectance rather than the magnitude data, which is more relevant to the dominant constituents and their relative concentrations; it is particularly suitable for waters dominated by suspended sediments.The max-classification algorithm was used along with the max-min values for all R rs data sets between 400-580 nm to define five water types corresponding to Classes A-E as shown below (Figure 4).

‚
Class A: Phytoplankton dominant waters.The spectral peak was at 400 nm.
‚ Class B: Phytoplankton and CDOM co-dominant waters.The peak band was at approximately 490 nm, and the minimum one was near 580 nm.

‚
Class C: Mixed waters with no dominant components.A flat peak was present at 500-550 nm, and the minimum was near 400 nm.

‚
Class D: Total suspended matter dominant waters.The peak was around 565 nm.
‚ Class E: Very high suspended sediment concentrations in water.The peak exceeded 570 nm.

Method for the Classification of In-Situ Rrs Spectra
The Yellow Sea (YS) contains typical turbid class II waters, and the peak of Rrs spectra is to a great extent related to the dominant components and concentrations in the water.This was especially obvious in waters dominated by suspended sediments; see the appearance of the red shift with increasing concentrations (Figure 2 red arrow).Maximum-minimum values for supervised classification (referred to in shortened form as max-classification) [33] of Rrs spectra were adopted for the whole in-situ data set in order to divide the data into homogeneous optical groups.One key advantage of the max-classification method is that it is sensitive to the shape of the reflectance rather than the magnitude data, which is more relevant to the dominant constituents and their relative concentrations; it is particularly suitable for waters dominated by suspended sediments.The maxclassification algorithm was used along with the max-min values for all Rrs data sets between 400-580 nm to define five water types corresponding to Classes A-E as shown below (Figure 4).

In-Situ R rs Classification and Spectra Features
Class A (N = 27): phytoplankton dominant water or oligotrophic water.The shape of spectra peaked at wavelengths of 400 nm (or less than 412 nm) and values declined with increasing wavelengths; flat or near zero values were found over 600 nm (Figure 5).In general, this type of water was closely related to the ocean Case I water type, which exists farther away from the shore and river discharges; hence, this water is only slightly influenced by coastal matter.The water color always appeared as blue, Forel-Ule values ranged from 4 to 7, and the transparency was very high; the median value for the Secchi depth was 12 m (Table 1).This type of water contained relatively low Chl a concentrations, but they dominated the change of the spectral shape; low concentrations of CDOM and SPM were also present.This type water was not observed during the two 2003 cruises, which were restricted in terms of the sampling region and time.In general, these data were scarce in the other data sets (Figure 2).

In-Situ Rrs Classification and Spectra Features
Class A (N = 27): phytoplankton dominant water or oligotrophic water.The shape of spectra peaked at wavelengths of 400 nm (or less than 412 nm) and values declined with increasing wavelengths; flat or near zero values were found over 600 nm (Figure 5).In general, this type of water was closely related to the ocean Case I water type, which exists farther away from the shore and river discharges; hence, this water is only slightly influenced by coastal matter.The water color always appeared as blue, Forel-Ule values ranged from 4 to 7, and the transparency was very high; the median value for the Secchi depth was 12 m (Table 1).This type of water contained relatively low Chl a concentrations, but they dominated the change of the spectral shape; low concentrations of CDOM and SPM were also present.This type water was not observed during the two 2003 cruises, which were restricted in terms of the sampling region and time.In general, these data were scarce in the other data sets (Figure 2).Class B (N = 167): phytoplankton and CDOM co-dominant water.The shape of spectra peaked at around 490 nm and values declined with increasing wavelengths; values were flat over 600 nm, and most spectra had absorption peaks at 443 nm (Figure 5).The change of spectral shape in this type of water appeared to be dominated by changes in the Chl a and CDOM absorptions.In general, this type water was relatively far away from the shore and river discharges, and water depths were relative deep.The water color always appeared as green-blue, Forel-Ule values ranged from 5 to 11, and the transparency was relatively high; the median value for the Secchi depth was 9 m (Table 1).This type of water was common in clear coastal zones.
Class C (N = 185): mixed water with no dominant components.The shape of spectra showed a flat peak over 500-550 nm (Figure 5), and decreases or increases in sediment concentrations would necessitate new classifications as Class B or Class D, respectively.The SPM concentration was relative low, but SPM contributed to the reflectance spectra notably via backscattering processes.The water color always appeared as green, Forel-Ule values ranged from 6 to 14, and the transparency was moderate; the median value for the Secchi depth was 6 m (Table 1).This type of water displayed high spatio-temporal variability.
Class D (N = 211): suspended particle matter dominant water.The shape of spectra showed a peak around 565 nm (Figure 5), which was remarkably governed by backscattering SPM processes.The water color always appeared as bright yellow, Forel-Ule values ranged from 7 to 18, and the transparency was relatively low; the median value for the Secchi depth was 2.75 m (Table 1).This class of water was influenced by coastal suspended sediments that were either re-suspended from the seafloor in shallow waters or discharged by rivers.The SPM concentration was usually about 10 mg/L in magnitude.Class B (N = 167): phytoplankton and CDOM co-dominant water.The shape of spectra peaked at around 490 nm and values declined with increasing wavelengths; values were flat over 600 nm, and most spectra had absorption peaks at 443 nm (Figure 5).The change of spectral shape in this type of water appeared to be dominated by changes in the Chl a and CDOM absorptions.In general, this type water was relatively far away from the shore and river discharges, and water depths were relative deep.The water color always appeared as green-blue, Forel-Ule values ranged from 5 to 11, and the transparency was relatively high; the median value for the Secchi depth was 9 m (Table 1).This type of water was common in clear coastal zones.
Class C (N = 185): mixed water with no dominant components.The shape of spectra showed a flat peak over 500-550 nm (Figure 5), and decreases or increases in sediment concentrations would necessitate new classifications as Class B or Class D, respectively.The SPM concentration was relative low, but SPM contributed to the reflectance spectra notably via backscattering processes.The water color always appeared as green, Forel-Ule values ranged from 6 to 14, and the transparency was moderate; the median value for the Secchi depth was 6 m (Table 1).This type of water displayed high spatio-temporal variability.
Class D (N = 211): suspended particle matter dominant water.The shape of spectra showed a peak around 565 nm (Figure 5), which was remarkably governed by backscattering SPM processes.The water color always appeared as bright yellow, Forel-Ule values ranged from 7 to 18, and the transparency was relatively low; the median value for the Secchi depth was 2.75 m (Table 1).This class of water was influenced by coastal suspended sediments that were either re-suspended from the seafloor in shallow waters or discharged by rivers.The SPM concentration was usually about 10 mg/L in magnitude.Class E (N = 28): very high suspended sediment concentrations in water.The shape of spectra showed a peak at wavelengths exceeding 570 nm and a line appeared at 400-565 nm; major reflectance changes occurred at 560-700 nm or longer wavelengths (Figure 5), and these characteristics were always obtained in shallow water regions or near the mouths of rivers.The water color always appeared as yellow-brown, Forel-Ule values ranged from 18 to 21, and the transparency was very low; the median value for the Secchi depth was 0.4 m (Table 1).The SPM concentrations were usually near or above 100 mg/L in magnitude.

Distribution Features of In-Situ Reflectance and Bio-Optical Properties
The remote sensing reflectance data from all field stations were classified by the max-classification method.As can be seen, obvious differences in the spatial distribution characteristics were observed for four typical regions in the Yellow Sea (Figure 6), including coastal Shandong (CS, the water area in the inner 30 m isobaths near Shangdong Peninsula, north of 35°N), the north Yellow Sea (NYS, north of 37°N excluding the part consisting of coastal Shangdong), in the middle of the Yellow Sea (MYS, located approximately in the trough of the Yellow Sea above the 60 m isobaths), and Jiangsu shoal (JS, the water area in the inner 20 m isobaths near Jiangsu province between 32°N-34.5°N)[19,34]; the water in these areas ranged from clear to very turbid.Moreover, the suspended particle  Class E (N = 28): very high suspended sediment concentrations in water.The shape of spectra showed a peak at wavelengths exceeding 570 nm and a line appeared at 400-565 nm; major reflectance changes occurred at 560-700 nm or longer wavelengths (Figure 5), and these characteristics were always obtained in shallow water regions or near the mouths of rivers.The water color always appeared as yellow-brown, Forel-Ule values ranged from 18 to 21, and the transparency was very low; the median value for the Secchi depth was 0.4 m (Table 1).The SPM concentrations were usually near or above 100 mg/L in magnitude.

Distribution Features of In-Situ Reflectance and Bio-Optical Properties
The remote sensing reflectance data from all field stations were classified by the max-classification method.As can be seen, obvious differences in the spatial distribution characteristics were observed for four typical regions in the Yellow Sea (Figure 6), including coastal Shandong (CS, the water area in the inner 30 m isobaths near Shangdong Peninsula, north of 35°N), the north Yellow Sea (NYS, north of 37°N excluding the part consisting of coastal Shangdong), in the middle of the Yellow Sea (MYS, located approximately in the trough of the Yellow Sea above the 60 m isobaths), and Jiangsu shoal (JS, the water area in the inner 20 m isobaths near Jiangsu province between 32°N-34.5°N)[19,34]; the water in these areas ranged from clear to very turbid.Moreover, the suspended particle  Class E (N = 28): very high suspended sediment concentrations in water.The shape of spectra showed a peak at wavelengths exceeding 570 nm and a line appeared at 400-565 nm; major reflectance changes occurred at 560-700 nm or longer wavelengths (Figure 5), and these characteristics were always obtained in shallow water regions or near the mouths of rivers.The water color always appeared as yellow-brown, Forel-Ule values ranged from 18 to 21, and the transparency was very low; the median value for the Secchi depth was 0.4 m (Table 1).The SPM concentrations were usually near or above 100 mg/L in magnitude.

Distribution Features of In-Situ Reflectance and Bio-Optical Properties
The remote sensing reflectance data from all field stations were classified by the max-classification method.As can be seen, obvious differences in the spatial distribution characteristics were observed for four typical regions in the Yellow Sea (Figure 6), including coastal Shandong (CS, the water area in the inner 30 m isobaths near Shangdong Peninsula, north of 35°N), the north Yellow Sea (NYS, north of 37°N excluding the part consisting of coastal Shangdong), in the middle of the Yellow Sea (MYS, located approximately in the trough of the Yellow Sea above the 60 m isobaths), and Jiangsu shoal (JS, the water area in the inner 20 m isobaths near Jiangsu province between 32°N-34.5°N)[19,34]; the water in these areas ranged from clear to very turbid.Moreover, the suspended particle  Class E (N = 28): very high suspended sediment concentrations in water.The shape of spectra showed a peak at wavelengths exceeding 570 nm and a line appeared at 400-565 nm; major reflectance changes occurred at 560-700 nm or longer wavelengths (Figure 5), and these characteristics were always obtained in shallow water regions or near the mouths of rivers.The water color always appeared as yellow-brown, Forel-Ule values ranged from 18 to 21, and the transparency was very low; the median value for the Secchi depth was 0.4 m (Table 1).The SPM concentrations were usually near or above 100 mg/L in magnitude.

Distribution Features of In-Situ Reflectance and Bio-Optical Properties
The remote sensing reflectance data from all field stations were classified by the max-classification method.As can be seen, obvious differences in the spatial distribution characteristics were observed for four typical regions in the Yellow Sea (Figure 6), including coastal Shandong (CS, the water area in the inner 30 m isobaths near Shangdong Peninsula, north of 35°N), the north Yellow Sea (NYS, north of 37°N excluding the part consisting of coastal Shangdong), in the middle of the Yellow Sea (MYS, located approximately in the trough of the Yellow Sea above the 60 m isobaths), and Jiangsu shoal (JS, the water area in the inner 20 m isobaths near Jiangsu province between 32°N-34.5°N)[19,34]; the water in these areas ranged from clear to very turbid.Moreover, the suspended particle  Class E (N = 28): very high suspended sediment concentrations in water.The shape of spectra showed a peak at wavelengths exceeding 570 nm and a line appeared at 400-565 nm; major reflectance changes occurred at 560-700 nm or longer wavelengths (Figure 5), and these characteristics were always obtained in shallow water regions or near the mouths of rivers.The water color always appeared as yellow-brown, Forel-Ule values ranged from 18 to 21, and the transparency was very low; the median value for the Secchi depth was 0.4 m (Table 1).The SPM concentrations were usually near or above 100 mg/L in magnitude.

Distribution Features of In-Situ Reflectance and Bio-Optical Properties
The remote sensing reflectance data from all field stations were classified by the max-classification method.As can be seen, obvious differences in the spatial distribution characteristics were observed for four typical regions in the Yellow Sea (Figure 6), including coastal Shandong (CS, the water area in the inner 30 m isobaths near Shangdong Peninsula, north of 35°N), the north Yellow Sea (NYS, north of 37°N excluding the part consisting of coastal Shangdong), in the middle of the Yellow Sea (MYS, located approximately in the trough of the Yellow Sea above the 60 m isobaths), and Jiangsu shoal (JS, the water area in the inner 20 m isobaths near Jiangsu province between 32°N-34.5°N)[19,34]; the water in these areas ranged from clear to very turbid.Moreover, the suspended particle matter concentrations ranged from low to very high and optical properties of the water bodies and the associated causes were different from each other.Representative water types Class A and Class B occurred in the middle of the Yellow Sea, Class B and Class C occurred in the north Yellow Sea, Class C and Class D occurred in coastal Shandong, and Class D and Class E occurred in the Jiangsu shoal.There were a few instances of Class D waters near the mouth of the Yalu River, and rarely, Class A waters were observed in the center of the north Yellow Sea.Data trended from Class D to Class C waters away from the shoreline in 30 m isobaths off coastal Shandong.The external driving factors included the warm current of the YSWC, coastal currents along the coast, tidal resuspension in shallow waters, Changjiang fresh water discharges with suspended sediments, tidal activity, Class E (N = 28): very high suspended sediment concentrations in water.The shape of spectra showed a peak at wavelengths exceeding 570 nm and a line appeared at 400-565 nm; major reflectance changes occurred at 560-700 nm or longer wavelengths (Figure 5), and these characteristics were always obtained in shallow water regions or near the mouths of rivers.The water color always appeared as yellow-brown, Forel-Ule values ranged from 18 to 21, and the transparency was very low; the median value for the Secchi depth was 0.4 m (Table 1).The SPM concentrations were usually near or above 100 mg/L in magnitude.

Distribution Features of In-Situ Reflectance and Bio-Optical Properties
The remote sensing reflectance data from all field stations were classified by the max-classification method.As can be seen, obvious differences in the spatial distribution characteristics were observed for four typical regions in the Yellow Sea (Figure 6), including coastal Shandong (CS, the water area in the inner 30 m isobaths near Shangdong Peninsula, north of 35 ˝N), the north Yellow Sea (NYS, north of 37 ˝N excluding the part consisting of coastal Shangdong), in the middle of the Yellow Sea (MYS, located approximately in the trough of the Yellow Sea above the 60 m isobaths), and Jiangsu shoal (JS, the water area in the inner 20 m isobaths near Jiangsu province between 32 ˝N-34.5 ˝N) [19,34]; the water in these areas ranged from clear to very turbid.Moreover, the suspended particle matter concentrations ranged from low to very high and optical properties of the water bodies and the associated causes were different from each other.Representative water types Class A and Class B occurred in the middle of the Yellow Sea, Class B and Class C occurred in the north Yellow Sea, Class C and Class D occurred in coastal Shandong, and Class D and Class E occurred in the Jiangsu shoal.There were a few instances of Class D waters near the mouth of the Yalu River, and rarely, Class A waters were observed in the center of the north Yellow Sea.Data trended from Class D to Class C waters away from the shoreline in 30 m isobaths off coastal Shandong.The external driving factors included the warm current of the YSWC, coastal currents along the coast, tidal resuspension in shallow waters, Changjiang fresh water discharges with suspended sediments, tidal activity, monsoon winds, and typhoons.Overall, the optical properties and values for each region were variable, but the leading water types were relatively stable in one or two adjoining water bodies.monsoon winds, and typhoons.Overall, the optical properties and values for each region were variable, but the leading water types were relatively stable in one or two adjoining water bodies.The Rrs is an apparent optical property parameter that depends on inherent optical properties (such as absorption and scattering) and matter concentration within water [7].Each region in this study displayed unique optical properties that were influenced by different external factors (Figure 7 and Table 2), such as, summer which is influenced by southwest monsoon and fresh rain water discharged by Changjiang River, winter which is influenced by north and northwest winds, and spring and autumn which are the transition seasons [19,34].The distribution of diffuse attenuation coefficient, total absorption and backscattering (Kd(490), a(488) and bb(488)) are directly related to reflectance (Rrs(490)).ad(400) and ag(400) are related to matter concentration within water [7].The hypothesis for making contour lines was that the field station data were representative of the region and the bio-optical properties of surrounding water were fairly uniformly distributed.A contour map could be generated from dispersed site data through applying an interpolation method in a geographic statistical model.The interpolation method used in this study was ordinary Kriging interpolation, which can be performed by using the "Geostatistical Analyst" function module in ArcMap software.The statistical results are provided in Table 2.
(1) In the middle of the Yellow Sea, the average values of Rrs were very small, i.e., lower than 0.0070 sr −1 (Figure 7a,b, Table 2), and values peaked at 400 nm or 490 nm in accordance with the appearance of Class A and Class B waters; most of the Class B spectra appeared as a feature of absorption at 443 nm (Figure 6, MYS).Optical property values were the lowest in the Yellow Sea, and a(488), c(488), and bb(488) were respectively lower than 0.175 m −1 , 1.14 m −1 , and 0.005 m −1 in winter (Figure 7e,f, Table 2); moreover, ad(400) and ag(400) were respectively lower than 0.1 m −1 and 0.17 m −1 in spring, 0.084 m −1 and 0.354 m −1 in summer (Figure 7g,d), and 0.045 m −1 and 0.15 The R rs is an apparent optical property parameter that depends on inherent optical properties (such as absorption and scattering) and matter concentration within water [7].Each region in this study displayed unique optical properties that were influenced by different external factors (Figure 7 and Table 2), such as, summer which is influenced by southwest monsoon and fresh rain water discharged by Changjiang River, winter which is influenced by north and northwest winds, and spring and autumn which are the transition seasons [19,34].The distribution of diffuse attenuation coefficient, total absorption and backscattering (K d (490), a(488) and b b (488)) are directly related to reflectance (R rs (490)).a d (400) and a g (400) are related to matter concentration within water [7].The hypothesis for making contour lines was that the field station data were representative of the region and the bio-optical properties of surrounding water were fairly uniformly distributed.A contour map could be generated from dispersed site data through applying an interpolation method in a geographic statistical model.The interpolation method used in this study was ordinary Kriging interpolation, which can be performed by using the "Geostatistical Analyst" function module in ArcMap software.The statistical results are provided in Table 2. (1) In the middle of the Yellow Sea, the average values of R rs were very small, i.e., lower than 0.0070 sr ´1 (Figure 7a,b, Table 2), and values peaked at 400 nm or 490 nm in accordance with the appearance of Class A and Class B waters; most of the Class B spectra appeared as a feature of absorption at 443 nm (Figure 6, MYS).Optical property values were the lowest in the Yellow Sea, and a(488), c(488), and b b (488) were respectively lower than 0.175 m ´1, 1.14 m ´1, and 0.005 m ´1 in winter (Figure 7e,f, Table 2); moreover, a d (400) and a g (400) were respectively lower than 0.1 m ´1 and 0.17 m ´1 in spring, 0.084 m ´1 and 0.354 m ´1 in summer (Figure 7g,d), and 0.045 m ´1 and 0.15 m ´1 in autumn.The K d (490) was lower than 0.27 m ´1 in summer.The water optical properties were mainly influenced by the warm current of the YSWC, and clear waters or Class A areas were little influenced by turbid tides from Changjiang River discharges [20,34].(2) In the north Yellow Sea, which extends from the middle of the Yellow Sea, the average values of R rs were relatively small and lower than 0.01 sr ´1 (Figure 7a,b, Table 2).Reflectance peaks appeared near 490 nm and flat peaks appeared at 500-565 nm; these were affected by non-pigmented particles and CDOM (Figure 6, NYS), similar to the features of Class B and Class C waters spectra.The north Yellow Sea optical property values were relatively low near the middle of the Yellow Sea.The average values of a d (400), a phy (400), and a g (400) were 0.084 m ´1, 0.055 m ´1, and 0.15 m ´1 in autumn, respectively, which were larger than the values in the middle of the Yellow Sea.The relative percentage of a phy (400) was lower than a g and a d , and this could suggest that R rs optical properties were affected by non-pigmented particles and CDOM.K d (490) was 0.072-0.46m ´1 in spring and 0.109-0.507m ´1 in autumn ((Figure 7c), and the values were distributed as concentric circles with minimum values in the center.The water optical properties were influenced by the warm current branch of the YSWC that extends through the MYS region and the NYS region [20,34].Four stations near the Yalu River mouth were influenced by river discharges as well as by coastal and bottom erosion processes; here, the reflectance peaked at 565 nm and values were lower than 0.025 sr ´1.(3) In coastal Shandong, the average values of R rs were relatively high, i.e., 0.01-0.02sr ´1 (Figure 7a,b, Table 2), and flat peaks appeared at 500-560 nm as well as a peak at 565 nm in accordance with the appearance of Class C and Class D waters (Figure 6, CS).According to the statistical data (Table 1), the suspended sediment and CDOM concentrations in coastal water bodies were relatively high compared to the north and middle regions of the Yellow Sea.The values of a d (400), a phy (400), and a g (400) were relative larger than the values in the north Yellow Sea, and the absorption coefficients showed contours parallel with the coastline in spring with relatively higher values in the inner 30 m isobaths than the outside ones.Additionally, the absorption coefficient isoclines were nearly perpendicular with the coastline in summer, and the values decreased from high values along the Qingdao coast to low values in the middle of the Yellow Sea (Figure 7g,d The water optical properties were mainly influenced by coastal currents and re-suspended sediments from shallow waters [19,34].(4) In Jiangsu shoal, the values of R rs were higher than those in the other regions of the Yellow Sea; specifically, the values were larger than 0.015 sr ´1 (Figure 7a,b, Table 2).A peak at 565 nm and a peak at wavelengths larger than 570 nm were observed in accordance with the appearance of Class D and Class E waters (Figure 6, JS).The spectra of Class D in terms of the amount were different between coastal Shandong and Jiangsu shoal; the maximum values of R rs in Jiangsu shoal were more than 0.025 sr ´1, and these were seriously influenced by suspended sediments [34,35].a g (400) values were larger than 0.3 m ´1 in the inner 10 m isobaths, and the contours were approximately parallel to the shoreline when influenced by tides and currents in spring and larger than 0.25 m ´1 in summer when the maximum center of values shifted to north of the shoal (Figure 7d). a d (400) and a phy (400) were, respectively, larger than 3.5 m ´1 and 0.2 m ´1 in spring and larger than 0.5 m ´1 and 0.17 m ´1 in summer (Figure 7g).sediments from shallow waters that were disturbed by the actions of coastal currents, mooninduced tides, and monsoons.

Classification Variation in Different Regions and Seasons
On the basis of an interpolation method implemented with a geographic statistical model, dispersed site data were also used to generate continuous surfaces.Again, ordinary Kriging interpolation was performed using the "Geostatistical Analyst" function module of ArcMap software.The data for the interpolation were derived from four cruises in 2006-2007, and these data were selected to represent the four seasons of the year as follows: July (summer), 2006, January (winter), April (spring), October (autumn), 2007.The spatio-temporal characteristic of the Yellow Sea were then analyzed with the interpolated classification results for the R rs spectra (Figure 8).winter and in the north Yellow Sea at all seasons except winter.In the north Yellow Sea, the largest area of Class B waters occurred in summer.In coastal Shandong, the inshore region was mainly dominated by Class D waters, which adjoined the coastline all year round, and the offshore region contained Class C waters in spring and winter; here, Class B waters appeared occasionally in spring.
In Jiangsu shoal, Class E waters were present all year round and the area experienced tremendous seasonal variations.Because Class A's optical properties were closely related to ocean clear water properties, its origin was presumably the Tsushima Current in the west Pacific, which is carried by the Yellow Sea warm current, and Class B waters evolved there from Class A. The area of Class A or Class B waters in the middle of the Yellow Sea demonstrates the strength of the Yellow Sea warm current to a certain extent.When the area of Class A or Class B grew gradually, the water turned clear and the water quality was good.Class E's optical properties were closely related to high concentrations of suspended particle matter (Table 1).Class D was distributed along the coastline or at the periphery of Class E waters.The areas and shapes of Class E and Class D were strongly influenced by discharges of the Changjiang River, which were turbid and carried nutrients such as nitrogen and phosphorus.The area and shape of Class E and Class D waters near the Changjiang River were closely related to the quantity of terrigenous materials discharged by the Changjiang River.In comparison to Class A and Class B waters, the water quality of Class E and Class D waters was more turbid, the transparency was very poor, and the euphotic layer was shallow; these features may have serious impacts on fisheries and tourism in the coastal zones.Apparent optical properties (AOPs) depend on the medium (inherent optical properties), and such data typically display enough regular features and stability to be useful descriptors of water bodies [7].The distribution of R rs classification types, which were representative of the optical properties of regional waters, changed temporally and spatially in the Yellow Sea.In the middle of the Yellow Sea, Class A waters were not observed in winter, and these waters were the most widely distributed in autumn.Class B waters were the dominant type in the middle of the Yellow Sea in winter and in the north Yellow Sea at all seasons except winter.In the north Yellow Sea, the largest area of Class B waters occurred in summer.In coastal Shandong, the inshore region was mainly dominated by Class D waters, which adjoined the coastline all year round, and the offshore region contained Class C waters in spring and winter; here, Class B waters appeared occasionally in spring.In Jiangsu shoal, Class E waters were present all year round and the area experienced tremendous seasonal variations.
Because Class A's optical properties were closely related to ocean clear water properties, its origin was presumably the Tsushima Current in the west Pacific, which is carried by the Yellow Sea warm current, and Class B waters evolved there from Class A. The area of Class A or Class B waters in the middle of the Yellow Sea demonstrates the strength of the Yellow Sea warm current to a certain extent.When the area of Class A or Class B grew gradually, the water turned clear and the water quality was good.Class E's optical properties were closely related to high concentrations of suspended particle matter (Table 1).Class D was distributed along the coastline or at the periphery of Class E waters.The areas and shapes of Class E and Class D were strongly influenced by discharges of the Changjiang River, which were turbid and carried nutrients such as nitrogen and phosphorus.The area and shape of Class E and Class D waters near the Changjiang River were closely related to the quantity of terrigenous materials discharged by the Changjiang River.In comparison to Class A and Class B waters, the water quality of Class E and Class D waters was more turbid, the transparency was very poor, and the euphotic layer was shallow; these features may have serious impacts on fisheries and tourism in the coastal zones.

Operational Classification Tree
In general, classification of the waters with satellite imagery data proved to be a useful tool for providing timely spatial and temporal information about the state of coastal waters.Based on the characteristic peaks for the five types of waters (Figure 4), band wave centers for MERIS data (Figure 2), and the quantity of matched-up pixel spectra (Figure 3), a classification tree (Figure 9) was designed by optimizing maximum-minimum values with a supervised classification method.The classification rules in Figure 4 were partly modified as discussed below. (1) Setting threshold value to preserve the valid water data.To avoid the loss of turbid water data and the removal of other water data, all reflectance pixels in the study area took part in the classification and threshold values were used.Up to 90% of the total radiance received at the sensor for most cloud-free scenes was affected by the intervening atmosphere [36], and the field spectra and MERIS matched-up site pixel spectra in Figures 2 and 3 show that the B1 (413 nm) band and B2 (443 nm) band reflectances of coastal waters were less than 0.05 sr ´1; the peak of reflectance at the B5 (560 nm) band was greater than 0. Because of overestimates due to atmosphere scattering at short wavelengths, B2 was chosen as the reference band and a 0.05 sr ´1 threshold value was used.Additionally, positives for B5 that included most of the water spectra peaks between 413-709 nm were examined.Through the threshold values and peak bands, we were able to completely discard the pixels for cloud, glint, white caps, and other non-relevant objects.(2) Partition in water types by band peaks and band ratios.Waters were classified according to the shape of spectra, and for this, thresholds were used to the classification in overcoming the shifts of peaks and the influence of quantity changes in the adjacent bands caused by atmosphere correction.It is reasonable to categorize pixels to the corresponding water types based on spectra shapes by comparing band peaks with adjacent band values within 5% or 10% in practice, such as R rs (490)/R rs (510) < 1.05 or R rs (490)/R rs (560) < 1.1 for Class C when the peak is R rs (490), and R rs (560)/R rs (510) < 1.1 for Class C when the peak is R rs (560).Class E spectra peaks displayed an obvious red-shift phenomenon near 570 nm (Figures 2 and 6), and the inflexion point near 600 nm was not changed very much in comparison to that for Class D. B5 and B6 band data from MERIS were obtained at 560 nm and 620 nm, respectively, and these data were useful for categorizing the pixels in which the peaks appeared; this was particularly true for over 50% of the B5 data for Class E in this paper.

Operational Classification Tree
In general, classification of the waters with satellite imagery data proved to be a useful tool for providing timely spatial and temporal information about the state of coastal waters.Based on the characteristic peaks for the five types of waters (Figure 4), band wave centers for MERIS data (Figure 2), and the quantity of matched-up pixel spectra (Figure 3), a classification tree (Figure 9) was designed by optimizing maximum-minimum values with a supervised classification method.The classification rules in Figure 4 were partly modified as discussed below.
(1) Setting threshold value to preserve the valid water data.To avoid the loss of turbid water data and the removal of other water data, all reflectance pixels in the study area took part in the classification and threshold values were used.Up to 90% of the total radiance received at the sensor for most cloud-free scenes was affected by the intervening atmosphere [36], and the field spectra and MERIS matched-up site pixel spectra in Figures 2 and 3 show that the B1 (413 nm) band and B2 (443 nm) band reflectances of coastal waters were less than 0.05 sr −1 ; the peak of reflectance at the B5 (560 nm) band was greater than 0. Because of overestimates due to atmosphere scattering at short wavelengths, B2 was chosen as the reference band and a 0.05 sr −1 threshold value was used.Additionally, positives for B5 that included most of the water spectra peaks between 413-709 nm were examined.Through the threshold values and peak bands, we were able to completely discard the pixels for cloud, glint, white caps, and other non-relevant objects.(2) Partition in water types by band peaks and band ratios.Waters were classified according to the shape of spectra, and for this, thresholds were used to assist the classification in overcoming the shifts of peaks and the influence of quantity changes in the adjacent bands caused by atmosphere correction.It is reasonable to categorize pixels to the corresponding water types based on spectra shapes by comparing band peaks with adjacent band values within 5% or 10% in practice, such as Rrs(490)/Rrs(510) < 1.05 or Rrs(490)/ Rrs(560) < 1.1 for Class C when the peak is Rrs(490), and Rrs(560)/ Rrs(510) < 1.1 for Class C when the peak is Rrs(560).Class E spectra peaks displayed an obvious red-shift phenomenon near 570 nm (Figure 2 and 6), and the inflexion point near 600 nm was not changed very much in comparison to that for Class D. B5 and B6 band data from MERIS were obtained at 560 nm and 620 nm, respectively, and these data were useful for categorizing the pixels in which the peaks appeared; this was particularly true for over 50% of the B5 data for Class E in this paper.

Imagery Classification and Frequency Distribution
The imageries of MERIS products were applied to the classification tree (Figure 9) to categorize water masses into the five optical types.From the cloudless imagery example (MER_RR_PRBCM_20070428_021351_Data) shown in Figure 10, which was almost the same as the interpolation result obtained with in-situ point data (Figure 8, spring), we found that the classification tree performed well and that the water types described the regional optical properties very well.The accumulative water types of a site in one season can show the spatial distribution of classes throughout the region; moreover, these data can be used to analyze seasonal or annual dynamic changes over time when compiled as imagery time series.

Imagery Classification and Frequency Distribution
The imageries of MERIS products were applied to the classification tree (Figure 9) to categorize water masses into the five optical types.From the cloudless imagery example (MER_RR_PRBCM_ 20070428_021351_Data) shown in Figure 10, which was almost the same as the interpolation result obtained with in-situ point data (Figure 8, spring), we found that the classification tree performed well and that the water types described the regional optical properties very well.The accumulative water types of a site in one season can show the spatial distribution of classes throughout the region; moreover, these data can be used to analyze seasonal or annual dynamic changes over time when compiled as imagery time series.Figure 11 shows that the Class A-E seasonal accumulative frequency distributions were generally in accordance with the interpolation map shown in Figure 8 for the four seasons.As to seasonal distributions over the entire MERIS archive (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012), the area of Class E at Jiangsu shoal changed tremendously, and it was at its largest in winter; in summer, the area of Class E water in this region decreased to about 30% to 60% of the wintertime area size (Figure 12).This seasonality was present during every year analyzed.The Changjiang River discharged large amounts of fresh water and suspended sediments in August, and the lowest levels of discharge were observed in winter.However, other environmental parameters (e.g., wind stress) also played an important role and contributed to the magnitude of the Class E turbid plume extension.In contrast, Class A waters were rare or non-existent in winter and reached their maximum extent in summer.Class B waters were relatively rare in winter, and they were practically non-existent in the north Yellow Sea too.Class C waters appeared at the lowest frequency in summer and showed obvious overlaps with Class A waters, which may have been related to large fresh water discharges from the Changjiang River, typhoons passing through the region, and clear ocean water carried by the warm current.Southeast and south Pacific monsoons in the summer were strong and brought heavy rains to the area.Additionally, north and northwest winds in the autumn and winter were strong and stirred up the sediments in shallow waters.Figure 11 shows that the Class A-E seasonal accumulative frequency distributions were generally in accordance with the interpolation map shown in Figure 8 for the four seasons.As to seasonal distributions over the entire MERIS archive (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012), the area of Class E at Jiangsu shoal changed tremendously, and it was at its largest in winter; in summer, the area of Class E water in this region decreased to about 30% to 60% of the wintertime area size (Figure 12).This seasonality was present during every year analyzed.The Changjiang River discharged large amounts of fresh water and suspended sediments in August, and the lowest levels of discharge were observed in winter.However, other environmental parameters (e.g., wind stress) also played an important role and contributed to the magnitude of the Class E turbid plume extension.In contrast, Class A waters were rare or non-existent in winter and reached their maximum extent in summer.Class B waters were relatively rare in winter, and they were practically non-existent in the north Yellow Sea too.Class C waters appeared at the lowest frequency in summer and showed obvious overlaps with Class A waters, which may have been related to large fresh water discharges from the Changjiang River, typhoons passing through the region, and clear ocean water carried by the warm current.Southeast and south Pacific monsoons in the summer were strong and brought heavy rains to the area.Additionally, north and northwest winds in the autumn and winter were strong and stirred up the sediments in shallow waters.suspended sediments and organic matter [11].These features pose a great challenge for field-based optical measurements, which influence the validation results for satellite bio-optical products.The accuracy of measurements for the absorption coefficients a, a p , a phy , and a d were seriously affected by QFT measurements and absorption in waters containing organic particulate matter discharged from terrigenous sources.Previous research results have indicated that the incomplete removal of absorption from organic debris particles would lead to overestimates for the chlorophyll absorption coefficient [40], and the correction error for the Transmissivity method (T method) can be very large in coastal regions with high concentrations of inorganic particulate matter [41].Hence, only reasonable data that are not affected by the absorption of the residual particles can be used.Data from other field optical equipment such as AC9/ACS, BB9, VSF3, HS6, Lisst100X, and MPRO are usually outside of the measurement range, and the correction methods may be invalid too.Although, it may be possible to modify the measuring principle of AC9/ACS for shorter light paths in the measuring tube, i.e., use of 25 cm for ocean water to 10 cm for turbid water (Class D-E).Additionally, optical path reduction modules are available to reduce the optical paths by 50, 80, and 90% [42], and iterative scattering correction schemes for in-situ AC9/ACS measurements can be applied [43], albeit, there will still exist regions that are out of the dynamic measuring range of the device in very turbid waters (the R rs spectral peak above 620 nm as band 6 of MERIS).Thus, it should be emphasized that most current in-situ optical instruments, measurement protocols, and correction methods need to be checked to ensure that they are performing in an appropriate range and optimized for Class E and Class D water regions for validation of satellite products.
(4) Inversion algorithm.The spectra were relatively not sensitive to changes in the concentrations of Chl a and CDOM under high concentrations of suspended sediments [34,35], but at lower concentrations of suspended sediments, the data may have been influenced by the concentrations of Chl a and CDOM based upon the absorption characteristics.Previous research used empirical orthogonal function (EOF) analysis to study similar turbid waters and the results indicated that 74% [10] or nearly two-thirds [44] of the total variance of R rs can be partly explained by particulate backscattering (b bp ), whereas particulate and dissolved absorption only explained 15% [10] or 30% [44] of the ocean color variability based on different data sets.The associated classification results were shown to be a good tool for an inversion algorithm in previous research too [10,45].Because water bodies are optically complex in Yellow Sea, there needs to be more data to derive algorithms for concentrations of Chl a, CDOM and suspended sediments with R rs classifications.
(5) Eutrophic water.Class B-D water types are widespread in global coastal Case II waters, and most coastal water classifications have included them in their analysis procedures; these waters can be categorized even further as subdivisions such as eutrophic waters and red tides [10,32], for which the spectra peaks are about 560 nm, and the amplitudes of spectra are as low as 0.015 sr ´1.Waters in this study were categorized as eutrophic Class C waters according to imagery data showing flat spectra before the peak band; red tide spectra were rare in the data sets from the six cruises analyzed in this study.The peaks near 565 nm, for which the peak values were more than 0.015 sr ´1, were categorized into Class D waters, and appropriate values for classes of eutrophic waters and red tides need further analysis in the future.

Test of the Max-Classification on Other Turbid Waters and Comparisons with Other Operational Classification Methods
The five water types of the Yellow Sea distinguished here were based on the peaks and shapes of reflectance spectra, which were easy to apply to satellite products, and the classification results reflect different absorption and scattering properties of the water and concentrations of various water quality parameters.The maximum-minimum value supervised method may be too simple or too general to use to describe the detailed features of other complex waters, but it was found to be convenient for processing imagery data and suitable for representing the water types of the Yellow Sea during all four seasons.Moreover, the results were robust when the method was applied to turbid waters, and hence, it could be suitable for other turbid coastal waters such as in the Mississippi region and the English Channel (Figure 13).The max-classification results for the English Channel Class C-E waters were very similar to Classes 1-3 of Vantrepotte et al. [11] which were based on the average reflectance spectra derived for Ward's hierarchical clustering.The Loire and Dordogne estuaries in France contain Class D-E waters as confirmed by field measurements taken by Doxaran et al. [46] who developed an algorithm to estimate the concentrations of biogeochemical constituents.The five water types of the Yellow Sea, the Mississippi region and the English Channel were also consistent with the most frequently selected classes of maximum membership (averaged 5.2 dominant class) in the work of Mélin and Vantrepotte [13], which was based on a global seven-year SeaWiFS data set and the distribution of 16 classes in global coastal ocean water.France contain Class D-E waters as confirmed by field measurements taken by Doxaran et al. [46] who developed an algorithm to estimate the concentrations of biogeochemical constituents.The five water types of the Yellow Sea, the Mississippi region and the English Channel were also consistent with the most frequently selected classes of maximum membership (averaged 5.2 dominant class) in the work of Mélin and Vantrepotte [13], which was based on a global seven-year SeaWiFS data set and the distribution of 16 classes in global coastal ocean water.The partitions for the five classes of Yellow Sea water relevant to in-situ hyperspectral data were robust for MERIS multispectral data.There is an abundance of classification methods that can be applied to such data from Case II waters, such as k-means or fuzzy c-means [10,47], the ISODATA (Iterative Self-Organizing Data Analysis Technique) clustering method [13], the spectral angle method (SAM), vector machines, NN [48], and the Forel-Ule method [16].However, very few methods are available that can be used operationally with satellite sensors in the Yellow Sea because of the lack of in-situ data sets supporting their use and the complexity of coastal waters.Wernand et al. [16] classified the Yangtze (Changjiang) River with the Forel-Ule MERIS (FUME) algorithm to FU 19, and a major area close to the river outflow and Jiangsu shoal in MERIS imagery was flagged as a "no data area" that could not cover all water types because of the loss of turbid area data such as Class E data for which the FU values were larger than FU 19.The losses of turbid area data were also restricted by the signal saturation of MODIS data [18] and a mask in SeaWiFS [38].The fuzzy c-means method [10,47] needs to characterize optically distinct water classes a priori and to parameterize algorithms for each class, and it could not directly classify spectra to Class A-E because of the changeable shape of spectra in each class; however, it always clusters the normalized spectra together under the same feature, so it could be used to subdivide the max-classification result in some situations.
In some areas, the water classes appeared as continuous transitions, and in other areas, the closely related classes overlapped (Figure 10).These data reflect the optical complexity of coastal waters and high temporal variability in coastal processes, which makes it difficult to assess regional water types with global water classification methods.The findings show that the bio-optical properties of the Yellow Sea are typically dominated by the influence of suspended inorganic matter, and results for these properties may be a useful reference for global research on turbid water bodies or regions influenced by river discharge and sea bottom matter resuspension effects.

Conclusions
Based on spectral classifications with Rrs data in Yellow Sea waters, the temporal and spatial distributions of remote sensing reflectance were obtained along with diffuse attenuation coefficients, The partitions for the five classes of Yellow Sea water relevant to in-situ hyperspectral data were robust for MERIS multispectral data.There is an abundance of classification methods that can be applied to such data from Case II waters, such as k-means or fuzzy c-means [10,47], the ISODATA (Iterative Self-Organizing Data Analysis Technique) clustering method [13], the spectral angle method (SAM), vector machines, NN [48], and the Forel-Ule method [16].However, very few methods are available that can be used operationally with satellite sensors in the Yellow Sea because of the lack of in-situ data sets supporting their use and the complexity of coastal waters.Wernand et al. [16] classified the Yangtze (Changjiang) River with the Forel-Ule MERIS (FUME) algorithm to FU 19, and a major area close to the river outflow and Jiangsu shoal in MERIS imagery was flagged as a "no data area" that could not cover all water types because of the loss of turbid area data such as Class E data for which the FU values were larger than FU 19.The losses of turbid area data were also restricted by the signal saturation of MODIS data [18] and a mask in SeaWiFS [38].The fuzzy c-means method [10,47] needs to characterize optically distinct water classes a priori and to parameterize algorithms for each class, and it could not directly classify spectra to Class A-E because of the changeable shape of spectra in each class; however, it always clusters the normalized spectra together under the same feature, so it could be used to subdivide the max-classification result in some situations.
In some areas, the water classes appeared as continuous transitions, and in other areas, the closely related classes overlapped (Figure 10).These data reflect the optical complexity of coastal waters and high temporal variability in coastal processes, which makes it difficult to assess regional water types with global water classification methods.The findings show that the bio-optical properties of the Yellow Sea are typically dominated by the influence of suspended inorganic matter, and results for these properties may be a useful reference for global research on turbid water bodies or regions influenced by river discharge and sea bottom matter resuspension effects.

Conclusions
Based on spectral classifications with R rs data in Yellow Sea waters, the temporal and spatial distributions of remote sensing reflectance were obtained along with diffuse attenuation coefficients, absorption coefficients, and scattering coefficients.This study focused on seasonal changes in the spectral classifications in the Yellow Sea, and the main causes of these changes were investigated.The values of R rs spectra were similar in spring and autumn, and values were highest in winter and lowest in summer.The R rs of Jiangsu shoal was constantly at a high value all year round, and relatively low values were detected in the middle of the Yellow Sea and the central region of the north Yellow Sea.Influenced by the northward Yellow Sea warm current, the optical properties of the north Yellow Sea showed seasonal changes that were larger than those in other areas.
Class D spectra values in coastal Shandong were less than those in out-layers of Class E waters characterized by low SPM concentrations.Scattering, absorption, and attenuation of water were strong in winter and weak in summer.(3) The five spectra types in the Yellow Sea were mainly categorized by the R rs spectral peak features, which were easy to obtain from satellite products and analyze with an operational classification tree.The imagery classification frequency distribution maps showed a tendency whereby the area of Class E waters at Jiangsu shoal gradually shrank in summer and expanded in winter; this change was tremendous and the coverage during the summer was about 30% to 60% of the wintertime area size.
Optical properties are very complex in coastal waters, and the development of a suitable optical classification approach that can be used operationally with satellite sensors needs more in-situ data set collection and validation work.Five types of water classes were identified in this work, and the approach used needs to be tested with other sensors.Several mechanisms can affect ocean color variability in complex coastal environments, and influencing factors include currents, tidal activity, fresh water discharges, monsoon winds, and typhoons.

Figure 1 .
Figure 1.Locations of the in-situ sampling sites in the Yellow Sea (left); and coastal water depths from 0 to 120 m, which were referenced from the GEBCO_2014 Grid, www.gebco.net(right).

2. 2 .
Data Acquisition and Pre-Processing 2.2.1.In-Situ Data Sets This work uses Yellow Sea in-situ data sets that contain data from 618 stations.These data were acquired from six cruises including: (i) two Cal/Val cruises in spring (22 March-23 April) and autumn (3-27 September) of 2003; (ii) two coastal background survey cruises in summer (15 July-6 August) and winter (2 January-4 February) of 2006/2007; and (iii) two coastal background survey cruises in spring (4 April-4 May) and autumn (14 October-4 November) of 2007 (symbols in Figure

Figure 1 .
Figure 1.Locations of the in-situ sampling sites in the Yellow Sea (left); and coastal water depths from 0 to 120 m, which were referenced from the GEBCO_2014 Grid, www.gebco.net(right).

2. 2 .
Data Acquisition and Pre-Processing 2.2.1.In-Situ Data Sets This work uses Yellow Sea in-situ data sets that contain data from 618 stations.These data were acquired from six cruises including: (i) two Cal/Val cruises in spring (22 March-23 April) and autumn (3-27 September) of 2003; (ii) two coastal background survey cruises in summer (15 July-6 August) and winter (2 January-4 February) of 2006/2007; and (iii) two coastal background survey cruises in spring (4 April-4 May) and autumn (14 October-4 November) of 2007 (symbols in Figure

Figure 2 .
Figure 2. Field measured remote sensing reflectance spectra and the visible-near infrared bands (B 1-15) of the Medium Resolution Imaging Spectrometer (MERIS) sensor.Red arrow covers the peak of spectra.

Figure 2 .
Figure 2. Field measured remote sensing reflectance spectra and the visible-near infrared bands (B 1-15) of the Medium Resolution Imaging Spectrometer (MERIS) sensor.Red arrow covers the peak of spectra.

Figure 3 .
Figure 3. (a) Pixel reflectance spectra of one cloudless imagery MERIS 2P products (16 April 2003) corresponding to the station analyzed during the spring cruise in 2003.(b) Strictly match-up scatter plot of the MERIS Rrs(λ) versus in-situ measurements during the spring and autumn cruises in 2003; the solid line is the 1:1 line, and dashed lines are the 1:2 and 2:1 lines.

Figure 3 .
Figure 3. (a) Pixel reflectance spectra of one cloudless imagery MERIS 2P products (16 April 2003) corresponding to the station analyzed during the spring cruise in 2003.(b) Strictly match-up scatter plot of the MERIS R rs (λ) versus in-situ measurements during the spring and autumn cruises in 2003; the solid line is the 1:1 line, and dashed lines are the 1:2 and 2:1 lines.

•
Class A: Phytoplankton dominant waters.The spectral peak was at 400 nm.• Class B: Phytoplankton and CDOM co-dominant waters.The peak band was at approximately 490 nm, and the minimum one was near 580 nm.• Class C: Mixed waters with no dominant components.A flat peak was present at 500-550 nm, and the minimum was near 400 nm.• Class D: Total suspended matter dominant waters.The peak was around 565 nm.• Class E: Very high suspended sediment concentrations in water.The peak exceeded 570 nm.

Figure 4 .
Figure 4.The flowchart for the maximum-minimum values for supervised classification (referred to in shortened form as the max-classification) of Rrs spectra.Min(Rrs) was used only to assist with the confirmation of each class; the minimums of Classes A and B were at 580 nm and Classes C and D were near 400 nm.

Figure 4 .
Figure 4.The flowchart for the maximum-minimum values for supervised classification (referred to in shortened form as the max-classification) of R rs spectra.Min(R rs ) was used only to assist with the confirmation of each class; the minimums of Classes A and B were at 580 nm and Classes C and D were near 400 nm.
matter concentrations ranged from low to very high and optical properties of the water bodies and the associated causes were different from each other.Representative water types Class A and Class B occurred in the middle of the Yellow Sea, Class B and Class C occurred in the north Yellow Sea, Class C and Class D occurred in coastal Shandong, and Class D and Class E occurred in the Jiangsu shoal.There were a few instances of Class D waters near the mouth of the Yalu River, and rarely, Class A waters were observed in the center of the north Yellow Sea.Data trended from Class D to Class C waters away from the shoreline in 30 m isobaths off coastal Shandong.The external driving factors included the warm current of the YSWC, coastal currents along the coast, tidal resuspension in shallow waters, Changjiang fresh water discharges with suspended sediments, tidal activity, matter concentrations ranged from low to very high and optical properties of the water bodies and the associated causes were different from each other.Representative water types Class A and Class B occurred in the middle of the Yellow Sea, Class B and Class C occurred in the north Yellow Sea, Class C and Class D occurred in coastal Shandong, and Class D and Class E occurred in the Jiangsu shoal.There were a few instances of Class D waters near the mouth of the Yalu River, and rarely, Class A waters were observed in the center of the north Yellow Sea.Data trended from Class D to Class C waters away from the shoreline in 30 m isobaths off coastal Shandong.The external driving factors included the warm current of the YSWC, coastal currents along the coast, tidal resuspension in shallow waters, Changjiang fresh water discharges with suspended sediments, tidal activity, matter concentrations ranged from low to very high and optical properties of the water bodies and the associated causes were different from each other.Representative water types Class A and Class B occurred in the middle of the Yellow Sea, Class B and Class C occurred in the north Yellow Sea, Class C and Class D occurred in coastal Shandong, and Class D and Class E occurred in the Jiangsu shoal.There were a few instances of Class D waters near the mouth of the Yalu River, and rarely, Class A waters were observed in the center of the north Yellow Sea.Data trended from Class D to Class C waters away from the shoreline in 30 m isobaths off coastal Shandong.The external driving factors included the warm current of the YSWC, coastal currents along the coast, tidal resuspension in shallow waters, Changjiang fresh water discharges with suspended sediments, tidal activity, matter concentrations ranged from low to very high and optical properties of the water bodies and the associated causes were different from each other.Representative water types Class A and Class B occurred in the middle of the Yellow Sea, Class B and Class C occurred in the north Yellow Sea, Class C and Class D occurred in coastal Shandong, and Class D and Class E occurred in the Jiangsu shoal.There were a few instances of Class D waters near the mouth of the Yalu River, and rarely, Class A waters were observed in the center of the north Yellow Sea.Data trended from Class D to Class C waters away from the shoreline in 30 m isobaths off coastal Shandong.The external driving factors included the warm current of the YSWC, coastal currents along the coast, tidal resuspension in shallow waters, Changjiang fresh water discharges with suspended sediments, tidal activity,

Figure 6 .
Figure 6.In-situ remote sensing reflectance (Rrs) spectra classification for distributed features, with 10 m, 20 m, 30 m and 60 m isobaths.Four regions with typical bio-optical properties are shown; these regions include: coastal Shandong (CS, inner 30 m isobaths, north of 35°N), the north Yellow Sea (NYS, north of 37°N, where the part containing CS is excluded), the middle of the Yellow Sea (MYS, above 60 m isobaths), and Jiangsu shoal (JS, inner 20 m isobaths).

Figure 6 .
Figure 6.In-situ remote sensing reflectance (R rs ) spectra classification for distributed features, with 10 m, 20 m, 30 m and 60 m isobaths.Four regions with typical bio-optical properties are shown; these regions include: coastal Shandong (CS, inner 30 m isobaths, north of 35 ˝N), the north Yellow Sea (NYS, north of 37 ˝N, where the part containing CS is excluded), the middle of the Yellow Sea (MYS, above 60 m isobaths), and Jiangsu shoal (JS, inner 20 m isobaths).
K d (490) was larger than 2.8 m ´1 in spring and 0.67 m ´1 in summer.a(488), c(488), and b b (488) were respectively larger than 6.85 m ´1, 30 m ´1, and 0.11 m ´1 in the inner 10 m isobaths, and values gradually decreased in the middle of the Yellow Sea; their contours basically paralleled the coast.The contour lines expanded in the northeast direction as a result of the northeastern flow of the Changjiang River discharges.In winter, Jiangsu shoal became a high absorption water body, and the backscattering in Jiangsu shoal was obviously stronger than that in the other areas (Figure7e,f).The water optical properties were influenced by discharges from the Changjiang River and re-suspended sediments from shallow waters that were disturbed by the actions of coastal currents, moon-induced tides, and monsoons.Remote Sens. 2016, 8, 321 11 of 24

Figure 8 .
Figure 8. Yellow Sea water class types diagrams for spring, summer, autumn, and winter.

Figure 8 .
Figure 8. Yellow Sea water class types diagrams for spring, summer, autumn, and winter.

Figure 9 .
Figure 9.The classification tree used to categorize water masses into the five optical types (A-E) for MERIS data; the max is the wavelength where the peak of MERIS B1-B9 lies in.

Figure 9 .
Figure 9.The classification tree used to categorize water masses into the five optical types (A-E) for MERIS data; the max is the wavelength where the peak of MERIS B1-B9 lies in.

Figure 10 .
Figure 10.(a) The color composite shows the 709, 560, and 413 nm channels as red, green, and blue, respectively.(b) Cloudless imagery classified into Class A (blue), Class B (green), Class C (light green), Class D (yellow), and Class E (brown); cloud, glint, and so forth are shown in black.Data for MERIS imagery (MER_RR_PRBCM_20070428_021351_Data) are from 28 April 2007 (spring).

Figure 10 .
Figure 10.(a) The color composite shows the 709, 560, and 413 nm channels as red, green, and blue, respectively.(b) Cloudless imagery classified into Class A (blue), Class B (green), Class C (light green), Class D (yellow), and Class E (brown); cloud, glint, and so forth are shown in black.Data for MERIS imagery (MER_RR_PRBCM_20070428_021351_Data) are from 28 April 2007 (spring).

Figure 11 .
Figure 11.Seasonal accumulative frequency distzributions of water types (Classes E-A, up to down) for the four seasons (summer (1 June-31 August) and winter (1 December-28 February) in 2006/2007, and spring (1 March-30 May) and autumn (1 September-30 November) in 2007) corresponding to four survey cruises in 2006-2007.The color bar at the bottom represents the number of days over the season.

Figure 11 .
Figure 11.Seasonal accumulative frequency distzributions of water types (Classes E-A, up to down) for the four seasons (summer (1 June-31 August) and winter (1 December-28 February) in 2006/2007, and spring (1 March-30 May) and autumn (1 September-30 November) in 2007) corresponding to four survey cruises in 2006-2007.The color bar at the bottom represents the number of days over the season.

Figure 13 .
Figure 13.Cloudless imagery classified into Class A (blue), Class B (green), Class C (light green), Class D (yellow), and Class E (brown); cloud, glint, and so forth are shown in black: (a) the Mississippi region (27 April 2007); and (b) the English Channel (8 April 2007).

Figure 13 .
Figure 13.Cloudless imagery classified into Class A (blue), Class B (green), Class C (light green), Class D (yellow), and Class E (brown); cloud, glint, and so forth are shown in black: (a) the Mississippi region (27 April 2007); and (b) the English Channel (8 April 2007).
spectra could be divided into five categories, from clear to very turbid waters, which were termed as Class A to Class E, and these classes showed seasonal and spatial variations in four regions, namely, MYS, NYS, CS, and JS.Class E waters were mainly distributed in the estuary of the Changjiang River and in Jiangsu shoal.Class D waters were distributed along the coast of the Shandong Peninsula and along the outer edge of Class E waters in Jiangsu shoal.Class C waters were distributed in the north Yellow Sea and along the southern coast.Class B waters were mainly distributed in the central area of the Yellow Sea and north Yellow Sea.Class A waters appeared in the north Yellow Sea central region during the summer and in the middle of the Yellow Sea except in winter.Regarding the origins of the water classes, Class E waters in Jiangsu shoal were mainly due to tidal-induced bottom sediment re-suspension and terrigenous inputs.Class D waters along the coast of the Shandong Peninsula were mainly caused by turbid waters created by alongshore currents.Class C waters appeared when concentrations of suspended sediments in former Class D type waters decreased or when concentrations of suspended sediments in former Class B type waters increased.Class B waters that appeared in the middle of the Yellow Sea and in the central region of the north Yellow Sea were mainly influenced by the Yellow Sea warm current.Class A waters existed in all four seasons except for a very small area in winter; compared with Class B, these waters were less affected by suspended particle matter and chlorophyll absorption dominated the optical property characteristics.(2)

Table 1 .
Statistics for the in-situ data sets for Classes A-E.Min, Max, Avg., Std., and N represent the minimum, maximum, average, standard deviation, and total numbers of stations, respectively.

Table 1 .
Statistics for the in-situ data sets for Classes A-E.Min, Max, Avg., Std., and N represent the minimum, maximum, average, standard deviation, and total numbers of stations, respectively.

Table 1 .
Statistics for the in-situ data sets for Classes A-E.Min, Max, Avg., Std., and N represent the minimum, maximum, average, standard deviation, and total numbers of stations, respectively.

Table 1 .
Statistics for the in-situ data sets for Classes A-E.Min, Max, Avg., Std., and N represent the minimum, maximum, average, standard deviation, and total numbers of stations, respectively.

Table 1 .
Statistics for the in-situ data sets for Classes A-E.Min, Max, Avg., Std., and N represent the minimum, maximum, average, standard deviation, and total numbers of stations, respectively.

Table 1 .
Statistics for the in-situ data sets for Classes A-E.Min, Max, Avg., Std., and N represent the minimum, maximum, average, standard deviation, and total numbers of stations, respectively.

Table 2 .
Statistics for the in-situ data set values for coastal Shandong (CS), north Yellow Sea(NYS), the middle Yellow Sea (MYS), and Jiangsu shoal (JS) during spring, summer, autumn, and winter.Minimum and maximum values are linked with dashes, as in Min-Max; averages are shown in bold, with standard deviations in parentheses; N is the number of samples.