Diurnal Variation of Light Absorption in the Yellow River Estuary

Considering the influence of river discharge and strong winds, the diurnal variability of ocean optical absorption properties in the Yellow River Estuary (YRE) is quantified, using in-situ measurements. The study finds that terrestrial sources due to the Yellow River discharge can cause high diurnal variation of water absorption because of the movement of river plume in the YRE, but such an influence diminishes far away from the Yellow River plume. The diurnal variability of water absorption, affected by strong winds, is found to be strengthened with a rapid increase of particles and colored dissolved organic matter (CDOM) arising from re-suspended sediment induced by wave forcing. The diurnal variability of particle absorption is controlled by non-algal particle absorption in the YRE, and the ratio of non-algal particle absorption (aNAP) and total particle absorption for most wavelengths is more than 0.56. The diurnal variation of spectral slope of non-algal particle absorption (SNAP) is found to vary within a narrow range, although large variability in the aNAP spectrum is observed. The CDOM is correlated negatively with salinity, and such negative correlation becomes weaker with the decreasing influence of riverine input. The spectral slope of CDOM absorption (Sg) may reflect the formation and constituents of CDOM with weak relationship to its concentration, and its relationship with the absorption of CDOM at 440 nm may be associated with the source of CDOM. The value of Sg, which is affected by re-suspended bottom sediment, is much lower than that derived from CDOM affected by Yellow River runoff. Disregarding the absorption of pure water, the diurnal variability of total water absorption stems principally from changes in non-algal particle matter rather than CDOM and Chl-a. By the observations of hourly GOCI (Geostationary Ocean Color Imager) data, the major diurnal variations of remote sensing reflectance at 680 nm are observed in near-coastal waters and the estuary of the Yellow River, which are mainly influenced by the flow discharge of Yellow River and strong winds. Finally, the seasonal differences of diurnal variations of water absorption caused by strong winds and river discharge are determined.


Introduction
Colored dissolved organic matter (CDOM), phytoplankton, and non-algal particles (NAP) are major light absorbing constituents in the ocean.These constituents determine the optical properties of natural waters and directly affect both the availability and spectral quality of light in the water column [1,2].Through their effect on the submarine light field, these optically active constituents have a close relationship with biological activity.The changes in these constituents complicate remote sensing of Chl-a (chlorophyll-a concentration) and hence the estimation of primary product from ocean color imagery [3][4][5].Therefore, a quantitative description of the dynamics and variability of water optical properties is especially important for better interpretation of satellite ocean color data and characterization of biogeochemical processes [6].
The optical properties for Case 1 waters, such as open oceanic waters, are relatively stable over large spatial and temporal scales due to the limited terrestrial and human influence.As opposed to Case 1 waters, the CDOM, phytoplankton, and non-algal particles for Case 2 waters are more abundant and variable due to the intricate hydrodynamic and ecological environment.Especially in the estuary waters, in addition to autogenic constituents which are composed of phytoplankton and its degradation products, large amounts of allogenic substances are transported from land by rivers or by wind, and substances are re-suspended from the sea bottom and are eroded from shorelines [7].Moreover, non-algal particles and CDOM in the coastal regions and estuaries contribute significantly to the optical properties compared to phytoplankton, and the constituent concentrations may be uncorrelated with one another.
Although the optical properties of ocean color constituents have been studied in some coastal regions and estuaries, the diurnal variability of the optical properties of the natural waters is less well documented.Some studies have reported diurnal variability of particle backscattering, phytoplankton absorption [8,9], and ocean optical properties during a coastal algal bloom [6], based on in-situ observations.
The YRE receives a large influx of river runoff, mainly coming from the Yellow River which is the largest river that discharges into the Bohai Sea.The discharge of the Yellow River is a major source of freshwaters and supplies a dominant amount of terrestrial matters to the YRE.Therefore, the optical properties of the water body are complex and variable due to the mix of freshwater from the Yellow River discharge and saline sea water.Moreover, the YRE area is characterized by high suspended particle matter (SPM) concentration and complicated hydrodynamics, which poses a challenge for ocean color remote sensing [10][11][12][13][14].However, the water optical properties in the YRE and its adjacent areas have been rarely reported [15], and the diurnal absorption variations of constituents and their mechanisms of influence on the environment have not been studied.
The objective of the present study therefore was to investigate the diurnal variations of light absorption by ocean constituents in the YRE, using in-situ continuous observation data.To understand the diurnal dynamics of CDOM and particle absorption, the influence of Yellow River discharge and strong winds on the diurnal absorption variability of these constituents was analyzed.Results of this study would improve our understanding of the YRE's constituent abundance, spatial and temporal variability, and optical characteristics, and could provide some reference values for other large estuaries as well.

Study Area
The YRE is located in the southwestern Bohai Sea, which is a semi-enclosed body of large, turbid and shallow water with an average depth of ~18 m in China [16].Some regions of the YRE are very shallow, especially the coastal regions, but in general the turbidity of the waters limits the bottom effect in many instances and therefore does not contaminate ocean color imagery [17].The Yellow River has an annual discharge of 420 × 10 8 m 3 and an annual sediment load estimated at 10 9 tons, nearly accounting for more than half of the total river discharge to Bohai Sea [18].The discharge of the Yellow River changes with season, and reaches a maximum in August and a minimum in April [19].The tides are mainly semi-diurnal in the YRE, and the spatial variability of tidal amplitudes has been given by Hu et al. and Hainbucher et al. [20,21].The basic hydrological features of the YRE, influenced by river inputs and dominated by runoff, can be summarized as low volume, heavy sands and weak tide [22].This region is also characterized by complex and capricious optical properties due to the transport of huge amounts of inorganic sediments and organic matters by the Yellow River carrying to the sea.

Measurement of Suspended Particle Matter Concentration
The suspended particle matter concentration (SPM; units: mg/L) of the surface water was measured by a standard gravimetric method [23][24][25].Water samples were collected from the top several centimeters of water surface using sampling bottles.The water samples were filtered through Whatman GF/F filters (Diameter 47 mm, pore size 0.45 mm).The filter was rinsed with 25-50 mL of pure water to remove as much salts as possible [26].The weight of the filter paper had been measured in the laboratory before the cruise, and the filter paper was dried under infrared light and stored for on-shore analysis [17].After taking to the laboratory, the filter was dried at 40 • C for 48 h and re-weighted at room temperature [27].

Measurement of Chlorophyll-a Concentration
For the estimation of total chlorophyll-a concentration (chl-a; units: µg/L), seawater samples (500 mL) were filtered through Whatman GF/F filter of 0.7 µm pore size under low vacuum pressure, and the filters were immediately preserved in liquid nitrogen before being processed.The filters were soaked in 90% acetone for 24 h and refrigerated in darkness at 4 • C in the laboratory.Then, Chlorophyll-a concentration (Chl-a) was determined for each sample using Turner-Design 10 fluorometer which was calibrated using Chl-a standards (Sigma, Sigma-Aldrich, St. Louis, MO, USA) [28].

Measurement of Light Absorption of Seawater
For the estimation of absorption coefficient of CDOM (a g (λ), units: m −1 ), seawater was filtered through polycarbonate filters of 0.22 µm pore size under low vacuum pressure.Absorbance spectra of CDOM were measured using a dual-beam Shimadzu spectrophotometer (UV-1800, Shimadzu, Japan) over the spectral range of 280 to 900 nm with an interval of 2 nm.Baseline data were obtained by filling Milli-Q water both in the sample and reference cells, and a baseline correction was applied by subtracting the offset from each sample spectrum [29].
For the estimation of absorption coefficient of particles (a p (λ), units: m −1 ) and non-algal particles (a NAP (λ), units: m −1 ), an adequate amount of seawater (depending on the particle load) was filtered through glass fiber filters of 0.7 µm pore size (Whatman GF/F) under low vacuum pressure.The absorption spectrum of the particles retained on the filter was measured for wavelengths between 280 and 900 nm with 2 nm increments using a dual-beam Shimadzu spectrophotometer (UV-1800) according to the filter-pad technique [30].Non-algal spectral absorption was measured from the sample filter after bleaching of pigmented particles [31].The background signals were subtracted, and correction of path-length described in Tassan et al. (2000) [32] was applied.The absorption coefficient for phytoplankton pigments (a ph (λ), units: m −1 ) was then derived by subtracting a NAP (λ) from a p (λ).

Remote Sensing Reflectance
Inherent optical properties (IOPs) such as the absorption coefficient (a) and scattering coefficient (b), together with the directional structure of the ambient light field, determine the apparent optical properties (AOPs).As one of the important AOPs, the surface spectral reflectance derived from satellite sensors can be used to reflect the variation of absorption coefficient.To investigate the diurnal variability of water absorption in the larger area, the diurnal variation of remote sensing reflectance (R rs ) derived from GOCI (Geostationary Ocean Color Imager) images was analyzed.
The acquisition dates of the GOCI images with cloud-free or low cloud coverage were 22 September 2011, 16 October 2011, 25 October 2011 and 13 November 2011.The corresponding flow discharges of these days were 2670, 1150, 836 and 876 m 3 /s, respectively.According to the data of the National Center for Environmental Prediction (NCEP), strong winds began on 23 October, reached the highest speed of approximately 14.5 m/s at 6 p.m., and then gradually declined to 7.5 m/s on 25 October.The wind speed on 13 November was 9.81 m/s.Therefore, GOCI images (25 October and 13 November) can be used to map the spatio-temporal distributions of R rs after strong winds, which could reveal the effect of strong winds on the diurnal variation of R rs .R rs (680) was selected as an indicator of SPM in the analysis below, due to its significant sensitivity to SPM variation in the turbid waters [22].

Data Analysis
The absorption by CDOM can be generally expressed by the following exponential equation [33]: where a g (λ) and a g (λ 0 ) represent the absorption of CDOM at any wavelength λ and a reference wavelength λ 0 , respectively.In general, the choice of a reference absorption wavelength is arbitrary, and several values of λ 0 can be found in the literature [34].However, 440 nm is generally chosen to be a reference wavelength because it corresponds, approximately, to the peak absorption by phytoplankton and represents the CDOM concentration, since CDOM is often measured optically at that wavelength [29,34].S g is the spectral slope coefficient obtained using a non-linear fit within a certain wavelength range, which means how rapidly the absorption decreases with increasing wavelength.In this paper, we chose the wavelength range from 380 to 690 nm.It must to be noted that the choice of representative wavelength range is different in most previous studies, which make the inter-comparison between various studies difficult [29,[35][36][37].
Non-algal particles and CDOM often have similar spectral behavior, and have usually been modeled using the same exponential function as follows [38]: where λ 0 is the reference wavelength, a NAP (λ 0 ) is the measured absorption value at the reference wavelength, and exponent S NAP represents the spectral slope of the absorption spectrum.In this study, the non-algal particle absorption spectrum fitted by the least square technique ranged from 380 to 700 nm, and 440 nm was chosen as the reference wavelength.
The descriptive statistics for the diurnal variations of water absorption was applied, that included the mean, maximum, and minimum values with standard deviations (SDs).The relationship of the parameters was indicated by means of the correlation coefficient (R).Based on the assumptions about independence, normal distribution and homogeneity of variances, one-way ANOVA (analysis of variance) was used to determine the difference or correlation, which defined statistically significant at p < 0.05.The ANOVA F test compares the variability between stations to the variability within stations.The number of degrees of freedom (DF) was the number of values in the calculation of a statistic.

SPM and Chl-a Concentrations
The hourly variation of SPM and Chl-a concentration for three fixed stations is shown in Figure 2. It is noted that patterns of variations of SPM and Chl-a concentration at the P1, P2 and P3 stations were different.For Station P1, the SPM concentration was found to vary significantly with fluctuations within 26 h.According to the statistics shown in Table 1, a maximum value of 758.69 mg/L was observed at 9 a.m. on day one and a much lower minimum value of 11.92 mg/L was observed at 4 a.m. the next day with the difference between them of 1 order of magnitude.At Station P2, with low SPM concentration less than 30 mg/L, the SPM diurnal variation of SPM was much smaller compared to Station P1, and the SD value of SPM was only 5.68 mg/L.At Station P3, the SPM concentration was always high during observation and decreased continuously from 10 a.m. on day one.However, the variations of Chl-a were not large for three stations with the measurement range from 0.24 to 1.34 µg/L, and the variation patterns of Chl-a and SPM were not alike.
Figure 2 also shows the variation of sea surface salinity (SSS) during observation at Stations P1, P2 and P3.The SSS of the P1 station was less than 27 psu with the lowest salinity of 4.92 psu, which was affected by the water from the Yellow River with low salinity and low temperature.Moreover, just as the variation of SPM, the SSS varied significantly with obvious fluctuation within 26 h, but the trend of change of SSS was opposite to that of SPM except for the time from 16 to 21 with relatively high SSS which could be caused by strong vertical mixing of the low and high salinity water.The dramatic variation of SSS at Station P1 indicates the influence of the Yellow River freshwater by river plume as time goes by.The SSS at Station P2 varied between 25.25 psu and 29.65 psu, which was higher compared to the Station P1.It is indicated that the influence of the Yellow River decreased for Station P2.For Station P3, the variation of SSS was small with the maximum and minimum values of Figure 3 shows the wind speed at YRE coming from NCEP in November and December 2011 and in-situ wind speed observed at Stations P1, P2 and P3 during the observation period.The wind speed at Station P3 was higher than that at Stations P1 and P2, the maximum wind speed was up to 10.4 m/s during the observation period.In addition, strong winds occurred on 7 and 10 December 2011 before observation, with the highest wind speed of approximately 14.2 m/s and 12.0 m/s, respectively.
In summary, we can conclude that Station P1 represents a typical estuarine environment which is affected by the Yellow River discharge, Station P2 is likely to have intermediate properties of both the marine and freshwater end members which are affected weakly by the Yellow River discharge, and Station P3 represents the marine environment which is affected by re-suspended sediment mainly induced by the earlier strong wind.

Relationship between CDOM Absorption and Salinity or Temperature
The sources of CDOM could be deduced with absorption coefficients and their correlation with environmental factors [40].The relationship of a g (440) vs. SSS and SST (Sea Surface Temperature) for Stations P1, P2 and P3 is illustrated in Figure 5 influence linked to the Yellow River discharge, which indicates that the Yellow River discharge is a significant factor controlling the variation of CDOM absorption, especially for Station P1 with the SSS range of 4.9 to 26.9 psu.The terrestrial input is the primary source of natural CDOM to the ocean in the YRE area, where the Yellow River runoff mixes with seawater.However, the weak correlation between a g (440) and SSS at P3 station suggests that it is not the Yellow River discharge but the bottom material re-suspension caused by strong winds that plays an important role in the variability of CDOM optical properties, just like causing the variation of SPM.

Spectral Slope of CDOM Absorption
In this paper, the CDOM absorption in the range of 380 to 690 nm was used to calculate S g by exponential fitting with the correlation coefficient of more than 0.9, as the reference wavelength was set at 440 nm.The statistics of S g are provided in Table 3.The S g values at Station P1 and P2 have clear differences with the S g values at Station P3.ANOVA tests indicated significant differences in S g (DF = 46, F = 25.16,p < 0.001) between Stations P1 and P3, whereas no significant variability was observed in S g (DF = 47, F = 2.73, p = 0.11) between Stations P1 and P2.Although the diurnal variation of S g for each fixed station was found to be not as obvious as expected, it is apparent that the S g value of CDOM at Station P3 showed a distinguishable pattern by contrast to Station P1 and P2.We can see that S g which was affected by the re-suspended sediment was more variable and exhibited lower values than those derived from the CDOM affected by the Yellow River runoff.The variation of S g was not greatly different between Stations P1 and P2 as the absorption of CDOM was.This showed that the variation of S g could result from the change of composition, but not the concentration of CDOM [40].The spectral slopes were not significantly different at Stations P1 and P2, likely indicating that similar CDOM sources were influenced by the Yellow River discharge in spite of the difference in the CDOM concentration, by contrast there was a bigger difference of S g at Station P3 where the influence of wind was strong, which suggests a high degree of variability in the composition of CDOM in these waters.
Therefore, even in the YRE, it is difficult to choose one S g to represent the S g of the whole study area because there are large differences in S g which could be due to different sources of CDOM.This result demonstrates the potential for using S g to characterize the composition of CDOM [41,42].Some researchers found an inverse linear relationship between the logarithm of a g (440) and S g [41].In addition, according to [15], the absorption of CDOM at 440 nm and S g can be parameterized using a power law, i.e., S g = Y 0 a g (440 where Y 0 and Y 1 are undetermined coefficients.Figure 6 shows both kinds of relationship between S g and a g (440) mentioned above.It can be seen that there is a strong negative relationship between a g (440) and S g at Station P3 with the correlation coefficient greater than 0.92, and a weak negative relationship at Station P2 with the correlation coefficient higher than 0.24.However, the data from Station P1 do not exhibit the same behavior.
In Figure 6a,b, we can clearly identify two different CDOM sources which are from the terrestrial input and re-suspended sediments corresponding to Stations P1 and P3.These findings suggest that CDOM at Station P1 has different optical properties from those at Station P3, and this relationship between the spectral slope coefficient and absorption coefficient has potential to differentiate between terrestrially derived CDOM from marine source derived CDOM [34,43].

Total Particle Absorption
Figure 7 shows the spectral absorption coefficients for total particles (a p (λ)) and non-algal particles (a NAP (λ)) at Stations P1, P2 and P3.It can be seen that a p (λ) at Stations P1 and P3 were much higher than that at Station P2, which was accompanied by a greater diurnal variation.It is noted that the spectral absorption coefficient for total particles was nearly similar to that for non-algal particles, implying that the non-algal particles absorption could be the main contributor of a p (λ), which is different from that of open ocean water [44][45][46] and some eutrophic water [27], where the contribution of absorption by phytoplankton to the total particulate absorption is dominant.Figure 8 shows the variation of non-algal particles to total particle absorption ratio a NAP (λ)/a p (λ) at different wavelengths as a function of the SPM concentration for Stations P1, P2 and P3.It shows that the contributions of a NAP (λ) to a p (λ) were more than 56.4% and 60.7% at all main wavelengths (these wavelengths corresponded to channels of the GOCI instrument) for Stations P1 and P3, respectively, while such contribution was relatively low at Station P2 with more than 56.3% for 412 nm, 443 nm, 490 nm and 555 nm, and more than 38.8% for 660 nm and 680 nm which was caused by the peak in phytoplankton absorption around 676 nm.It is noted that the a NAP (λ)/a p (λ) values are approximately the same as were observed in the Yangtze Estuary (especially inside the river's mouth) [47], but significantly different from other areas.They were on average larger than those measured in the Mediterranean, Atlantic and subequatorial Pacific waters [46].This suggests that, in the YRE waters, not only the amplitude of a NAP but also its relative contribution to total particle absorption was significantly different from that in other areas.

Non-Algal Particle Absorption
Figure 9 shows hourly variation of non-algal particle absorption coefficient at 440 nm; large fluctuations of a NAP (440) were observed at Station P1, relatively lower values of a NAP (440) were observed at Station P2, and a gradual decline of a NAP (440) was observed at Station P3.The statistics of non-algal particle absorption coefficients at 440 nm at different stations are shown in Table 4.   Figure 10 shows a relationship between a NAP (440) and SPM for Stations P1, P2 and P3.There was a positive relationship between a NAP (440) and SPM for Stations P1, P2 and P3 with a correlation coefficient R 2 greater than 0.75, and this relationship can be described using a linear equation, as shown in Figure 10.It has been reported that organic particles may have a higher a NAP (443):SPM ratio than inorganic particles [38].We can note that the slope at station 2 (0.045 m 2 g −1 ) was higher than that at Stations P1 and P3 (0.031 m 2 g −1 and 0.032 m 2 g −1 ) which also supported this interpretation.Figure 10 shows that the variation of a NAP (440) was actually in accordance with the variation of SPM with high values at Stations P1 and P3, low values at Station P2, which indicated that SPM was composed mainly of non-algal particles in the YRE.

Total Water Absorption
Defining the total absorption in a water column as the sum of absorption by pure water (a w ), phytoplankton (a ph ), non-algal particles (a NAP ), and CDOM (a g ), the total absorption coefficient (a t ) at wavelength λ, can be expressed as [38,[48][49][50]: The relative contribution of a given component was calculated as the ratio of the absorption coefficient of that component to the total absorption coefficient of non-water constituents.Figure 11 illustrates the relative contribution of absorption by non-algal particles, phytoplankton pigments and CDOM derived from our data as triangular diagrams at six different wavelengths.As can be seen, the total absorption of non-water constituents in all our samples was almost dominated by non-algal particles at all wavelengths, except for only a few samples at Station P2 at longer wavelengths (660 nm and 680 nm) due to the absorption of phytoplankton.At all main wavelengths, the average contributions of non-algal particles (a NAP ) to the total non-water absorption (a ph + a NAP + a g ) were 75.3 ± 11.5%, 63.2 ± 11.6% and 76.5 ± 10.3% for Stations P1, P2 and P3, respectively.Only a small percentage of light in the estuary was absorbed by phytoplankton and CDOM, the spectrally averaged percentages for each material were: phytoplankton, 17.7 ± 9.7% (P1), 19.5 ± 8.3% (P2) and 18.0 ± 7.3% (P3), CDOM, 7.0 ± 6.2% (P1), 17.3 ± 9.0% (P2) and 5.5 ± 5.0% (P3).It is noted that the contributions of CDOM to the total absorption were less than 29.0% for P1, 46.7% for P2 and 21.9% P3.This result is different from the previous reports which indicated CDOM can contribute more than 50% to the light absorption budget at short wavelengths (e.g., 440 nm) at the surface of the ocean [14,38,46].It is also noted that the variations of the contributions of different components to the total absorption were related to the wavelength; with decreasing wavelength, the contributions of non-algal particles and CDOM increased.On the contrary, the contribution of phytoplankton was becoming the important absorber at the red end of the spectrum (660 nm and 680 nm), especially for Station P2.
Disregarding the absorption due to pure water, we found that the total non-water absorption was dominated by the particle absorption at almost all wavelengths, especially at Stations P1 and P3.Therefore, these results suggested that the water in the YRE was a typical Case 2 water, where non-algal particles were the dominant component affecting optical properties, and that the diurnal variability of absorption by water was principally due to non-algal particle matter rather than CDOM and Chl-a [11,42].

Diurnal Variation of Remote Sensing Reflectance
Figure 12 shows the diurnal variation of R rs at different GOCI bands at Stations P1, P2 and P3 in four different days.It shows that the magnitude of diurnal variation of R rs at Station P1 was larger than that at Stations P2 and P3 for all four days.According to the statistics shown in Table 6, the SDs in different days at Station P1 were the largest compared to those at Stations P2 and P3, and the values of SD were consistent with the flow discharges of Yellow River, with the maximum and minimum SDs occurring on 22 September and 13 November.The SDs at Station P2 were larger than the SDs at Station P3 due to the different distances to the Yellow River mouth.In addition, it was found that the SDs at Stations P1 and P2 on 25 October and 13 November were larger than they were on 22 September and 16 October, which could be related with the re-suspension of SPM caused by strong winds.As shown in Figure 13, using the GOCI 500-m data, we obtained hourly maps of R rs (680) in the YRE and its adjacent areas (including Bohai Bay and Laizhou Bay) on four different days.By the observations of hourly GOCI, the R rs (680) diurnal variations can be seen in the study area on these days.The SD images show that different regions exhibited different magnitudes of diurnal variations.The major regions of R rs (680) diurnal variation were generally located in near-coastal waters and in the estuary of the Yellow River.The diurnal variations of R rs (680) in the Bohai Bay and Laizhou Bay were higher than those in the center of Bohai Sea.With the effect of previous strong winds, the SDs of R rs (680) on 25 October and 13 November increased over a broad region, especially in the Bohai Bay and Laizhou Bay.Therefore, we can conclude that the flow discharge of Yellow River and strong winds have a significant impact on the diurnal variation of R rs in the YRE and its adjacent areas, which is consistent with the results of the light absorption study based on in-situ experimental investigation.

CDOM Absorption
Based on the observations at the P1 and P2 stations, it can be concluded that the diurnal variation of CDOM absorption decreased from the river estuary toward the offshore area.Meanwhile, observations at Station P3 show that the diurnal variation of CDOM absorption has a significant increase resulting from the re-suspended sediment.Considering the effect of river runoff and strong wind, the values of a g (440) found for Yellow River Estuary in this study were basically in the range observed for Bohai Sea by Zhou et al. [15], who reported that the range of a g (440) in the Bohai Sea was from 0.0958 to 0.7795 m −1 , and average value of a g (440) was 0.2311 m −1 .
CDOM usually has a much higher concentration in estuaries than in other marine systems [51], and most of it is introduced into the estuary by the way of river discharge.Generally, the component of CDOM is influenced by terrigenous sediment for the estuarine and coastal waters, and by marine sediment for the offshore waters.The high diurnal variation of CODM absorption at Station P1 was caused by terrigenous sediment that came with the Yellow River discharge, while this influence was reduced far away from the Yellow River plume at Station P2, where the variation of absorption of CODM was influenced by terrigenous sediment together with marine sediment.
In fact, the CDOM absorption shows a linear inverse correlation with salinity in many estuarine and coastal waters [5,[52][53][54][55][56][57], and based on this finding salinity fields are tracked through remote sensing-based observation of CDOM [10,57,58].Previous observations discovered that the negative relationship between SSS and a g (440) would begin to loosen as the SSS reached a certain extent, indicating a decreased influence of riverine input of CDOM [5,37,59].The present study demonstrates that this type of slope discontinuity breaks down at Station P3 with salinity more than 29.6 psu.
Previous studies in the Bohai Sea showed that the average value of S g derived for the wavelength range of 380-550 nm was 0.015 nm −1 and its standard deviation was 0.0039 nm −1 [15], which was similar to the average value of S g at all the three stations, but with a greater SD value due to the large study area.It also can be found that the S g values at Stations P1 and P2 except for Station P3 were nearly similar to the range of (0.013-0.018 nm −1 ) reported by Blough and Del Vecchio (2002) [53] for coastal waters influenced by river input and were close to the commonly assumed value of (0.015 nm −1 ) for CDOM in coastal waters [38].Zhou et al. (2015) [15] reported that weak negative correlations between S g and a g (440) were observed in the Bohai Sea and Northern Yellow Sea based on 895 stations surveyed during China's offshore marine optical investigation, but there was no clear relationship between S g and a g (440) in most of the Chinese offshore waters [15].In addition, this inverse relationship between the CDOM absorption at a certain wavelength and S g has been reported before in the literature [34,41,43,[60][61][62][63], whereas such a relationship was not found by other research reports [15,38,41,56].It is indicated that this negative relationship between S g and a g (440) can be associated with the source of CDOM, which was discussed by Xing et al. (2008) [64].

Particles Absorption
Although a large variability in the a NAP spectrum was observed in this study, the spectral slope S NAP varied within a narrow range which meant that the overall spectral shape was rather conservative for three stations.The low S NAP approaching 0.011 nm −1 was found for mineral particles in the North Sea [38] and Irish Sea [65], meanwhile the high S NAP approaching 0.013 nm −1 was observed in the Baltic Sea which is known for its high organic matter content [38].The range of S NAP values derived at all three stations was relatively narrow, with the mean S NAP less than 0.0113 nm −1 close to that found for mineral particles.Moreover, no obvious difference was observed for S NAP when comparing Station P1 and Station P3, and its average (0.0113 nm −1 and 0.0103 nm −1 ) was similar to the average found for mineral particles by Bowers et al. (1996) [65] and Ferreira et al. (2014) [56].In addition, the generally narrow range of variation in SNAP restricts this parameter used for purposes of monitoring [56], but allows a reliable modeling of a NAP due to the lack of dependence on the variables.

Total Water Absorption
The high variation of total absorption coefficient in a water column was caused by the variable contributions of the components other than pure water [11,38,45,66].The variability of absorption coefficient in open-ocean waters (Case 1 waters) has been well studied, and it has been reported that open ocean samples were dominated by phytoplankton absorption and all components except pure water are often assumed to covary with Chl-a concentration [11,[44][45][46]48,67].In contrast, the variability of absorption coefficient in coastal, estuarine and inland waters (Case 2 waters) were relatively poorly documented, and the components were considered to vary widely and independently (IOCCG, 2000).Moreover, the non-algal particles and CDOM make significant contributions to the total absorption in coastal and estuarine waters in addition to phytoplankton [38,56].

Factors Affecting Light Absorption Diurnal Variation
The YRE, well-known for its high SPM concentration in the world, is one of the most significant eco-regions that can be easily affected by its complicated hydrodynamic environment.In general, hydrological factors, which can cause a complicated water environment in and outside the estuary, mainly include river discharge, tidal currents, and wave forcing induced by wind [22,[68][69][70][71][72][73].

Strong Winds
As already analyzed for Station P3, under the strong wind influence, the water absorption experienced large variations within 26 h, especially with a rapid increase of particles and CDOM coming from re-suspended sediments, which cause the bio-optical properties much more complex.However, it has been reported that such a sharp increase, induced by re-suspended sediments at the sea surface, was episodic, had a short duration, and returned to normal quickly after the passage of a strong wind, which was also evidenced by satellite-retrieved data [22,72,74] and in-situ buoy data [74], respectively.However, the wind speed varies strongly with season, associated closely with monsoon activity in the study area.East Asian Monsoon has been considered to be a significant driving force for spatial and temporal sediment distribution in coastal waters [73,[75][76][77][78].The northerly winds in the winter are stronger than the dominant southerly winds in the summer, as is characteristic for the Bohai Sea.The events of maximum wind speed greater than 11 m/s from the northwest and northeast often occurred in the winter season [79].The submarine sediments are frequently re-suspended and transported by storm waves which are closely associated with the prevailing winds in the winter season [73,[80][81][82][83][84][85][86].Therefore, considering the shallow nature, climate and adequate terrestrial sediment supply from major rivers in the study area, winds in the winter season must play a major role in the diurnal variation of absorption by waters, resulting in significant seasonal differences of diurnal variations in water absorption.

Yellow River Discharge
The diurnal variation of water absorption in the YRE was affected by terrigenous matter that came with the Yellow River discharge, although this influence declined with the increase of distance from Yellow River plume.The degree of influence of terrigenous matter could be related with the river flow.In general, the water and sediment discharges of the Yellow River into the sea have a significant seasonal variability.Figure 14 shows the daily average discharge and sediment concentration of the Yellow River which were observed at the Lijin hydrological gauging station (obtained from the website http://www.mwr.gov.cn/) in 2011.In natural conditions, the flow discharge was larger in the wet season with high sediment concentration than in the dry season.However, apart from natural causes, the water and sediment discharges of Yellow River have been regulated since 2000 by the operation of the Xiaolangdi Reservoir, the largest reservoir in the mainstream, through the Project of "Artificial Regulation of the Yellow River Water and Sediment" [87].Thus, the flow discharge distribution had two high peaks from late June to early July and from late September to early October, respectively, which was due to reservoir's flood relief in the Yellow River upstream.Hence, the diurnal variation of water absorption in different seasons will definitely vary with the change of flow discharge and sediment concentration.Furthermore, considering the seasonal variation of monsoon, the strong wind could expand the influence of river discharge on the diurnal variation of water absorption at temporal and spatial scales.

Tidal Cycle
Aside from the strong wind-induced re-suspended sediment which caused the diurnal variation of water absorption, the tidal cycle may enhance vertical mixing and bottom sediment re-suspension, and thus can systematically and sustainably increase the variation of water absorption and change an area optically [56,72,74,[88][89][90].Cheng et al. (2016) [72] reported that macro-tidal currents play a significant role in the diurnal variation of suspended particle matter in the Yalu River Estuary of China, and there were usually two peaks of suspended particle matter in a tidal cycle corresponding to the maximum flood and ebb current.In addition, in Hangzhou Bay, He et al. (2013) [74] revealed that tide was the main driver of the diurnal spatial and temporal variability using in-situ data and GOCI-derived suspended particle matter concentration.Ferreira et al. (2014) [56] reported that strong tidal variations played a role in the distribution of suspended material and the relative contribution of organic and inorganic components in Santos Bay.
It is important to note that the various regions and tidal phases have different magnitudes of diurnal variations, which, to a great extent, could be dependent on water depth and tidal range.In general, the tidal cycle variations of re-suspended sediment often occur in the coastal region due to shallow water, while the tidal effect on the diurnal variation of re-suspended sediment is insignificant and even negligible in the open ocean [89].It has been demonstrated that the re-suspended sediment in the water column is correlated with the integration of tidal currents as well as the water depth in some studies [88,89,91,92].In addition, the magnitude of tide current was strongly correlated with the temporal and spatial variation of tidal range.A large magnitude of tidal current can induce a strong bottom shear stress and enhance stronger re-suspension of sediment, which could cause a high variation of water constituents in the coastal water [74].In contrast, a small tidal current with weak bottom shear stress and sediment re-suspension had limited influence on the variation of water constituents.With regard to the YRE area, the tide regime is complex and dominated by irregular semi-diurnal tides with an average tidal range of 0.6-0.8m, and that increases both southward and northward reaching 1.5-2.0m at Bohai Bay and Laizhou Bay [85,86,93].Therefore, we can conclude that tidal cycle is one of the hydrodynamic factors that can drive the diurnal variation of water absorption in the coastal region with shallow water of the YRE.However, this effect of tidal cycle on the diurnal variation of water absorption could be weak in the study area by comparison with macro-tidal estuary and coasts.In addition to tidal cycle, the current circulation was considered as the key hydrodynamic factor controlling the sediment distribution and transport in the region [85], so it may be a reasonable reason for elevated SPM and the factor affecting the light absorption diurnal variation to consider.

Implications of Diurnal Variation for Remote Sensing Products
The complex optical properties of Case 2 waters present distinct challenges to bio-optical modeling and remote sensing retrieval algorithms.In general, local empirical inversion algorithms were developed for the Case 2 water to solve its applicability problem caused by the temporal and spatial variation of optical properties.Our results indicate that bio-optical parameters could have significant diurnal variability caused by the hydrological and climatic factors.These short-term changes of optical properties challenge the applicability of local empirical inversion algorithms.To further improve the performance of ocean color inversion algorithms, the effect of temporal IOP variability needs to be fully considered [6, 94,95].
The validation of remote sensing products relies on a high-quality match-up dataset between satellite and in-situ data.It is difficult to perform measurements that are ideally synchronous with satellite overpass at specific sites for the majority of in-situ-satellite match-ups, so time differences are unavoidable in most cases.It was reported that a time window (e.g., 1-3 h) is usually accepted as a compromise [17,[96][97][98].According to our result, the errors induced by a 3 h time difference between in-situ measurements and satellite overpass could be great in certain areas or under certain climate conditions.Therefore, significant uncertainties introduced by time window due to the temporal variability and the effect of the relaxed match-up criterion on the assessment result should be considered in the validation activities, especially for the regions or conditions with significant diurnal variations.
In addition, our results indicated that short-term processes, such as the diurnal variation of IOPs caused by flow discharge of the River and strong winds, often occurred in the estuarine and coastal waters.Since the traditional polar-orbiting ocean color satellites may be insufficient to detect/monitor short-term processes for areas of significant temporal variability, the geostationary ocean color satellite with high-spatial and high-spectral resolution measurements would play a significant role in helping monitor short term and regional oceanic phenomena.

Conclusions
The diurnal variability of ocean optical absorption properties in the YRE is quantified, using in-situ measurements.The diurnal variation of absorption by water constituents, including CDOM and particles, can be affected by river discharge and wave forcing induced by strong winds.Terrestrial sources through the Yellow River discharge can cause a high diurnal variation of water absorption by the movement of river plume in the YRE.However, such an influence declines far away from the Yellow River plume.The diurnal variability of water absorption during a strong wind is found to be strengthened with a rapid increase of CDOM and particles coming from the re-suspended sediment induced by wave forcing.
We observe a general trend of decreasing a g (440) with increasing SSS in the YRE.However, such a negative relationship between CDOM and salinity begins to weaken as the SSS reaches a certain extent, which indicates a decreasing influence of riverine input of CDOM.It is demonstrated that S g could reflect the formation and component of CDOM, with weak relationship with the CDOM concentration.The diurnal variation of S g at one fixed station is found to be as not obvious as expected.However, considering the different sources of CDOM, S g which is affected by re-suspended sediment, is more variable and exhibits much lower values than those derived from CDOM affected by the Yellow River runoff.Consistent with previous studies, there exist an unstable negative relationship between S g and a g (440), which could be associated with the source of CDOM.
Different from the situation of open ocean waters and some eutrophic waters, the variability in particle absorption s controlled by non-algal particle absorption in the YRE, and the ratio of a NAP (λ)/a p (λ) at most wavelengths, except for 660 nm and 680 nm, is more than 56.4%.There was a positive relationship between a NAP (440) and SPM in the study area, which indicates that SPM was composed mainly of non-algal particles in the YRE.The diurnal variation of S NAP at all three stations is found to vary within a narrow range, although large variability in the a NAP spectrum is observed.
Disregarding the absorption of pure water, the total non-water absorption is dominated by the non-algal particle absorption at almost all wavelengths in the YRE, which suggests that the diurnal variability of water absorption is principally due to non-algal particle matter rather than CDOM and Chl-a.
By the observations of hourly GOCI, it is seen that the major regions of R rs (680) diurnal variation are generally located in near-coastal waters and in the estuary of the Yellow River.The flow discharge of the Yellow River and strong winds have a significant impact on the diurnal variation of R rs in the YRE and its adjacent areas, which is consistent with the results of the light absorption study based on in-situ experimental investigation.
Finally, we discuss the seasonal differences of diurnal variations in water absorption caused by strong wind and river discharge in the study area.Considering the characteristics of monsoon activities, frequent winds in the winter season must play a major role in the diurnal variation of water absorption.The diurnal variation of water absorption in different seasons varies with the change of flow discharge and sediment concentration in the area close to the river plume.In addition, the tidal cycle can increase the diurnal variation of water absorption by enhancing vertical mixing and bottom sediment re-suspension in the shallow coastal water.
Three continuous observations were conducted at 3 fixed stations in the YRE on 4 December 2011 (Station P1), 25 November 2011 (Station P2), and 11 December 2011 (Station P3), and the locations are shown in Figure 1.The average water depths at Stations P1, P2 and P3 are about 9.9 m, 16.5 m and 16.9 m, respectively.During the experiments, water samples were collected at the sea surface every hour from 8 a.m. to 9 a.m. the next day.Then, water samples were filtered and preserved immediately, and taken to the laboratory on shore to measure spectral absorption and biogeochemical properties.Salinity and temperature were measured with a CTD (Seabird, Seacat) during sampling.

Figure 2 .
Figure 2. Hourly variations at Stations P1 (red squares, N = 26), P2 (green circles, N = 26) and P3 (blue triangles, N = 26) of: SPM, SSS, and SST (a); and Chl-a (b).The variation of SSS is consistent with that of SST which can be used as an indicator of the Yellow River discharge.The change tendency of SPM for Station P1 is opposite to that of SSS except for the time from 16 to 21, which is not obvious for Station P2 and P3.

Figure 3 .
Figure 3.Time series of wind speed coming from NCEP around the Yellow River Estuary in: (a) November 2011; and (b) December 2011.(c) Time series of in-situ wind speed during the observation period at Stations P1 (open squares, N = 27), P2 (open circles, N = 27) and P3 (solid triangles, N = 27).The red lines indicate the dates when measurements were taken.Wind speed at Station P3 was higher than that at Stations P1 and P2, with the maximum wind speed up to 10.4 m/s.Strong winds occurred on 7 and 10 December before observation, with the highest wind speeds of approximately 14.2 m/s and 12.0 m/s, respectively.

Figure 4a shows the
Figure 4a shows the hourly spectral changes of CDOM absorption at Stations P1, P2 and P3.The diurnal variation of CDOM absorption was observed to have significant differences at Stations P1, P2 and P3.The hourly variations of absorption by CDOM at Station P1 and P3 were greater than that at Station P2. Figure 4b shows hourly variations of the CDOM absorption coefficient at 440 nm, large fluctuations of a g (440) were observed at the P1 and P3 stations.The statistics shown in Table 2 indicate that the difference of diurnal variation of a g (440) between Station P1 and Station P2 was significant (DF = 47, F = 16.29,p < 0.001) probably due to the different influence of terrigenous matter input from the Yellow River.The difference of diurnal variation of a g (440) between Station P3 and Station P2 was observed (DF = 48, F = 12.33, p < 0.001), which was attributed to the large variability of CDOM at Station P3 caused by re-suspended sediment induced by strong winds.It has been shown that strong winds can mix the shallow water column and release the CDOM trapped in sediments [39].

Figure 4 .
Figure 4. Spectral curve of CDOM absorption (a); and the CDOM absorption coefficient at 440 nm (b) at Stations P1 (red lines, N = 23), P2 (green lines, N = 25) and P3 (blue lines N = 24).The results show significant differences of diurnal variation of CDOM absorption at different stations.

Figure 5 .
Figure 5. Relationships of: a g (440) vs. SSS (a); and a g (440) vs. SST (b), at Stations P1 (red squares, N = 21), P2 (green circles, N = 24) and P3 (blue triangles, N = 24).The solid line represents the linear regression with SSS and SST for a g (440).A general trend of decreasing CDOM with increasing SSS and SST was observed for Station P1 and Station P2, but not observed for Station P3.

Figure 6 .
Figure 6.Power law relationship between a g (440) and S g (a), and inverse linear relationship between the logarithm of a g (440) and S g (b), at Stations P1 (red squares, N = 23), P2 (green circles, N = 25), and P3 (blue triangles, N = 24).The solid lines represent the fitted lines using the equations shown in the figure.The coefficients of determination R 2 are shown in the figure.There is a strong negative relationship between a g (440) and S g at Station P3, and a weak negative relationship at Station P2, but no observed relationship at Station P1.

Figure 7 .
Figure 7. Spectral curves of: total particle absorption (a); and non-algal particle absorption (b), at Stations P1 (red lines, N = 16), P2 (green lines, N = 22) and P3 (blue lines, N = 18).Note that the magnitude of the spectral absorption coefficient for total particles was nearly similar to that for non-algal particles.

Figure 8 .
Figure 8. Variations of the non-algal to total particle absorption ratio (a NAP /a p ) at different wavelengths as a function of the SPM concentration at Stations P1 (squares, N = 16), P2 (circles, N = 22) and P3 (triangles, N=18).Different colors represent different wavelengths.The contributions of a NAP (λ) to a p (λ) were more than 56.4% and 60.7% at all main wavelengths for Stations P1 and P3.

Figure 10 .
Figure 10.Relationship between a NAP (440) and SPM at Stations P1 (red squares, N=16), P2 (green circles, N = 22) and P3 (blue triangles, N = 18).The solid line represents the linear regression with SPM for a NAP (440), and the coefficients of determination R 2 are shown in the figure.Note that there was a positive linear relationship between a NAP (440) and SPM for Stations P1, P2 and P3 with a correlation coefficient R 2 greater than 0.75.

Figure 12 .
Figure 12.Remote sensing reflectance (R rs (λ)) at different GOCI bands at Stations P1, P2 and P3 derived from GOCI images on 22 September (N = 8), 16 October (N = 8), 25 October (N = 7) and 13 November 2011 (N = 6).Different colors represent different stations.Note that the magnitude of diurnal variation of R rs at Station P1 was larger than they were at Station P2 and P3 for all four days.

Figure 13 .
Figure 13.SD (standard deviation) maps of GOCI retrieved R rs (680) in the Yellow River Estuary and its adjacent areas on 22 September, 16 October, 25 October and 13 November 2011.Note that the major regions of R rs (680) diurnal variation were generally located in near-coastal waters and in the estuary of the Yellow River.

Figure 14 .
Figure 14.Time series of the Yellow River flow discharge and sediment concentration at Lijin hydrological gauging station in 2011.The flow discharge distribution had two high peaks from late June to early July and from late September to early October, respectively, which was mainly due to reservoir's flood relief in the Yellow River upstream.

Table 1 .
Statistics of SPM and Chl-a concentration measured at Stations P1, P2 and P3."Max.","Min.", and "SD" stand for maximum, minimum and standard deviation, respectively, and "N" is the number of samples.

Table 2 .
CDOM absorption coefficients at 440 nm in the Yellow River Estuary at different stations."Max.","Min.", and "SD" stand for maximum, minimum, and standard deviation, respectively, and "N" is the number of samples.

Table 5 .
Non-algal particle absorption spectral slopes S NAP at Stations P1, P2 and P3."Max.","Min.", and "SD" stand for maximum, minimum, and standard deviation, respectively, and "N" is the number of samples.