The properties of ocean color remote sensors utilized in this paper.
Analysis of in-situ measurements during a spring cruise survey in the Pearl River Estuary (PRE) reveals that, controlled by the two-layer gravitational circulation, chromophoric dissolved organic matter (CDOM) absorption shows a clear horizontal distribution pattern at both water surface and bottom, with higher CDOM absorption and lower spectral slope in the northwestern estuary, and a reversed pattern in the southeastern estuary and near the Hong Kong waters. The surface CDOM has higher absorption and lower spectral slope than the bottom. Horizontal transport is suggested to be the dominant hydrodynamic mechanism affecting CDOM distribution pattern in the PRE. With a regional algorithm tailored for the PRE CDOM absorption retrieval, a time series of CDOM absorption and spectral slope in the PRE and the Hong Kong waters in spring from 2012 to 2018 is produced based on satellite images obtained by four sensors with different spatial and spectral resolutions: the Visible Infrared Imaging Radiometer Suite (VIIRS), the Ocean and Land Colour Instrument (OLCI), the Hyperspectral Imager for the Coastal Ocean (HICO), and the Operational Land Imager (OLI). A correlation is revealed between the multi-temporal CDOM absorption and the monthly averaged river discharge, indicating the capability of CDOM ocean color products in identifying hydrodynamic processes in estuaries and coastal waters.
- chromophoric dissolved organic matter
- ocean color
- time series
- Pearl River Estuary
The chromophoric (or colored) dissolved organic matter (CDOM) is the light-absorbing part of the dissolved organic matter (DOM) in natural waters, which may cause a yellow or brown water color at high concentrations . Chemically, CDOM refers to an ensemble of structures including carboxylic acids and carboxyl-rich alicyclic molecules, substituted phenols, ketones, aldehydes, quinones, carbohydrates, saturated and unsaturated hydrocarbons, and nitrogenous material .
The major sources of CDOM in the marine environment are terrestrial-derived, oceanic-produced, or sediment-introduced . The terrestrial CDOM is originally synthesized by land and water plants, subsequently processed and modified in limnic systems and eventually exported to coastal waters. The oceanic CDOM is originally fixed by marine plants or phytoplankton, produced by heterotrophic and autotrophic organisms, or formed by photooxidation of colorless DOM [3, 4, 5]. The sediment CDOM is often observed in coastal waters and shelf seas where sediment resuspension, hypoxia events or hydrothermal events occur [6, 7, 8].
The processes involved in the removal of CDOM in nature waters include photodegradation and microbial activities. The photochemical reaction triggered by CDOM absorption of high-energy (low wavelength) light can have great impacts on biogeochemical processes and water ecology [9, 10, 11]. Heterotrophic microbes either incorporate or respire organic matter and modify CDOM into labile DOM, which are rapidly degraded in the process .
Estuaries and coastal waters are very productive systems, where high loading of terrestrial CDOM and high local production are mixed, processed and exported to shelf seas . Therefore, a better understanding of CDOM variation in estuaries and coastal waters can help estimating the oceanic carbon budget and evaluating the anthropogenic impacts on marine environments and global climate change. Furthermore, in estuarine and coastal waters, CDOM absorption usually co-varies with salinity. Variation of CDOM absorption can be used as a tracer of water mass mixing in near shore waters . Investigating and explaining CDOM variations are therefore crucial for understanding various processes in the aquatic environment.
The Pearl River Estuary (PRE) is located on the southern coast of Guangdong, China (22–22.75°N, 113.5–114°N) (Figure 1), adjacent to the Northern South China Sea (NSCS). Its major part, the Lingdingyang Estuary, has a trumpet-like shape, with four gates (Humen, Jiaomen, Hongqili, and Hengmen) on the west side of the upper estuary, discharging freshwater of 6.83 × 1010 m3 into the NSCS every year . The eastern side, however, has two deep channels, along which the coastal sea water can flood into . The various water masses with complex properties make it difficult to interpret the CDOM variations in the estuary.
Furthermore, as China’s second largest river in terms of water discharge, the Pearl River has an annual mean discharge of 10,000 m3 s−1 [15, 16, 17]. A significant seasonal variation exists in the Pearl River discharge. Therefore, the optical properties of the estuarine water can be strongly influenced by river discharge, especially during the wet season , which leads to distinctive characteristics when interpreting and quantifying CDOM variations in the PRE in different seasons.
Studies on CDOM absorption in the PRE last for decades. Previous studies of CDOM optical properties have covered the entire estuary in different seasons of the year, but detailed picture of CDOM variation in the estuary remains poorly understood, the dynamics of CDOM optical properties under control of multiple hydrodynamic and biogeochemical processes are not yet to be revealed satisfactorily, and the ocean color algorithms for CDOM optical property retrieval with high reliability are still needed to be developed, as well as the CDOM ocean color products with high spatial and temporal resolutions. Therefore, a comprehensive investigation of CDOM optical properties synthesizing advantages of in-situ observations and ocean color interpretations is still in necessity.
The aim of this research is to enhance the understanding of the spatial and temporal variations of CDOM optical properties in the PRE through analyzing in-situ measurements and interpreting satellite ocean color observations. The horizontal and vertical variations of CDOM absorption in the PRE are depicted based on dense and detailed in-situ observations and the dominant driven forces affecting the variations are discriminated. Spatial and temporal variations of CDOM optical properties in the PRE are analyzed based on the CDOM products derived from multi-source ocean color data with complemented spatial and temporal resolutions.
2. Materials and methods
2.1 In-situ observations
Generally, the absorption coefficient (ag(λ), m−1) is a direct measure of CDOM absorptivity. The spectrum of ag usually spans from ultraviolet band to red band. The spectral shape is unstructured and typically decreases with increasing wavelength in an exponential fashion (Eq. (1)) [18, 19, 20, 21, 22, 23].
where λ0 is the reference wavelength, and Sg is the spectral slope, describing the decreasing rate of the spectrum.
The absorptivity of CDOM is first measured on a spectrophotometer as absorbance (A(λ)), a unit-less ratio of spectral radiant power transmitted through the sample across a path length (L) . A(λ) is then converted to the (Napierian) absorption coefficient (ag(λ)) according to Eq. (2):
Typically, measurements are performed using a long-path length quartz cuvette (e.g., 0.1 m) or a liquid core waveguide (0.5–5 m) with submicron-filtered (0.2 or ~0.7 μm) seawater. In low-CDOM waters, measurement of the absorption spectrum of CDOM using a conventional spectrophotometer is challenging. Due to the exponential decline of absorption with increasing wavelength, absorbance values can fall below the detection threshold of the detectors used, typically 0.03–0.06 m−1 . Devices with longer effective path lengths, such as reflective tube absorption meters (up to 25 cm ), liquid waveguide cells (up to 200 cm [27, 28]), and integrating cavity meters (up to 25 m ) attempt to resolve this problem.
A cruise was implemented in the PRE previously from 2 to 12 May 2014. A total of 148 water samples were collected (Figure 2) from the water surface (with sampling depth of 0 and 1 m) and the bottom (with sampling depth 1 m above the bottom). CDOM absorbance (A(λ)) was measured in a 10-cm quartz cell using a Shimadzu UV-2550 spectrophotometer. The absorption coefficient (ag(λ)) is calculated at 1 nm interval from 250 to 700 nm according to Eq. (2). Spectral slope, Sg, over wavelength range of UV (250–400 nm), noted as Sg(250–400), is used as a representative of CDOM spectral slope when analyzing the spatial heterogeneity of CDOM .
In addition to the spectroscopic analysis of CDOM absorption, in-situ measurements of remote sensing reflectance in visible range (400–900 nm) are also conducted.
2.2 Satellite imagery
The in-situ measurements provide detailed observations on CDOM optical properties. However, the spatial and temporal resolutions of in-situ observations are always limited by high costs in the survey implementation. The satellite imagery can offer better temporal and spatial coverage with high efficiency and low cost. Therefore, space-born and airborne sensors have been used to monitor spectral properties of natural waters for over four decades, and ocean color remote sensing has become an important technology to study water environment.
Nevertheless, the applicability of satellite observation depends on the reliability of CDOM ocean color algorithms. The spectral properties of water bodies may change substantially with geographical locations and time. The color (spectrum) of water can be used to estimate concentrations of optically active constituents in water, such as the phytoplankton, suspended solids, and CDOM based on effective ocean color algorithms.
A regional algorithm for CDOM absorption estimation from satellite ocean color imagery is developed for the spring PRE (Eqs. (3)–(5)), which derives CDOM UV absorption and spectral slope from the visible remote sensing reflectance [Rrs()], which is established based on the in-situ measured above water remote sensing reflectance (Figure 3). The algorithm is with an overall mean absolute percentage difference (MAPD) of ~30 and ~5% for the estimation of CDOM absorption and the spectral slope over 250–450 nm, respectively.
Using all available satellite image data obtained by four ocean color sensors with different spatial and spectral resolutions: the Visible Infrared Imaging Radiometer Suite (VIIRS) on board Suomi National Polar-orbiting Partnership (Suomi NPP), the Ocean and Land Colour Instrument (OLCI) on board Sentinel-3A (S3A), the Hyperspectral Imager for the Coastal Ocean (HICO) integrated in the International Space Station (ISS) Window Observational Research Facility (WORF), and the Operational Land Imager (OLI) on board Landsat 8 (LS8) (Table 1), a time series of CDOM absorption and spectral slope in the PRE and the Hong Kong waters in spring from 2012 to 2018 is produced. Relevant factors related to the temporal variation of CDOM absorption and spectral slope are analyzed.
|Data period||Jan 2010 to Now||Oct 2016 to Now||Sep 2009 to Sep 2014||Feb 2013 to Now|
|Number of bands||22||21||128||9|
|Spectral ranges (nm)||402–12,490||400–1020||353–1080||435–2294|
|Spatial resolution (m)||375/750/1000||300||100||15/30|
To match the season of the in-situ observations in this study, a set of satellite images acquired in spring (April and May) from 2012 to 2018 is selected for the application of the developed algorithm. For the OLI data with a 16-day revisit cycle, the time constraint is relaxed to March and June. The study area covers 22° and 23°N, 113.5° and 114.5°E. The images with cloud coverage below 20% are downloaded (see Tables 2 and 3 for detailed information of the available images).
|Data period||Jan 2010 to Now||Oct 2016 to Now||Sep 2009 to Sep 2014||Feb 2013 to Now|
|Months searched||Apr to May||Apr to May||Apr to May||Mar to Jun|
Eq. (3) is applied to derive ag(290) from the selected satellite images. The Rrs(596) is matched to the bands of the four sensors by the criteria listed in Table 4. When retrieving Sg(250–400), considering the available spectral range of the satellite imagery and the performance of atmospheric correction, the lower ends for Rrs gradient calculation (the λmin in Eq. (4)) are set as 445, 415.5, 400, and 443 nm for the VIIRS, HICO, OLCI and OLI data, respectively. The range for the Rrs maximum is limited below 700 nm for all sensors. Sg(250–400) is afterward derived from Rrs gradient by Eq. (5).
|Sensor||Spectral bands in visible range (nm)||Rrs(596)|
|VIIRS||412/445/488/555/672||[0.66B(555) + 0.34B(672)]/2|
|OLCI||400/412.5/442.5/490/510/560/620/665/673.75/681.25||[B(560) + B(620)]/2|
|HICO||404–696 (52 bands with interval of ~5.7 nm)||[B(593) + B(599)]/2|
Figures 4 and 5 present time-series products of ocean color-retrieved CDOM absorption (ag(290)) and spectral slopes (Sg(250–400)) from 2012 to 2018, which are derived from the satellite images listed in Table 3 using Eqs. (3)–(5). All the products can capture the general distribution pattern of CDOM in the PRE, showing the gradual decrease of absorption and the increase of spectral slope along the northwest to southeast direction.
In addition to the clear and consistent trend illustrated in Figures 4 and 5, the quantitative variations of CDOM absorption and spectral slope with time are also remarkable. The temporal fluctuation is especially evident in the upper and western PRE, where the water properties are greatly influenced by freshwater discharge, especially during the flooding season.
Source and removal are the two major aspects that control and balance the CDOM absorption budget in nature waters. In open oceans, phytoplankton production and upwelling of deep water can bring new CDOM into upper layer to elevate CDOM absorption in surface water, while the photobleaching and high stratification can strongly decrease the absorptivity of surface CDOM and increase the spectral slope. In estuaries and coastal waters, the most distinctive condition different from the open oceans is the terrestrial input, which is a significant source of organic matter in surface water, and leads to the highest variation of CDOM concentration in the region . Therefore, in estuaries and coastal waters, river flow is always an important factor influencing the temporal variation of CDOM absorption during the wet season.
A positive correlation between the monthly averaged river discharge and the monthly averaged ag(290) maxima from the CDOM ocean color products is shown in the PRE (Figure 6), suggesting the high CDOM absorption in the PRE is always associated with high river flow. This pattern is especially typical in spring, when the Pearl River enters the flood season, and the river discharge can dramatically increase after a heavy rainfall. As a contrast to the high river flow, the existing CDOM in surface water of the PRE can be rapidly photodegraded or microbial consumed under the low river discharge with limited inputs of new CDOM, resulting in low level of CDOM absorption and large value of spectral slopes. In this point of view, the time-series of CDOM absorption can be well correlated with river discharge, and therefore is a good indicator of estuarine hydrodynamic conditions.
For dissolved matter in water, the variation of CDOM is greatly influenced by hydrodynamic conditions, such as horizontal transport and vertical mixing. The PRE circulation is dominated by the gravitational exchange current with the surface freshwater flowing seaward in the surface, mainly on the west side, and the compensated seawater in the lower layer flowing landward along the east coast [31, 32].
As illustrated by in-situ observations and satellite interpretation, CDOM shows the strongest absorption and lowest spectral slope at the head of the estuary (Station A1, A2 and A3) (Figures 4 and 5). Water in this region comes from the discharge of the Humen and Jiaomen Gates (Figure 2), carrying abundant dissolved organic materials. These CDOM components with complex composition are then transported seaward along the west coast, causing the western estuary to have stronger CDOM absorption and lower spectral slope than the eastern estuary, where the intrusive seawater brings large amount of marine CDOM with lower absorption and higher spectral slope.
It has been reported that terrestrial CDOM is much more sensitive to photochemical reactions than marine CDOM . DOM originating from upper streams during high-flow events is abundant in aromatic compounds, and therefore is highly susceptible to solar irradiation, which is in agreement with laboratory experiments which reveal a high degree of riverine DOM photoreaction . Thus, photodegradation can be a possible reason for the removal of terrigenous CDOM in the surface water of the PRE, therefore the decrease of CDOM absorption. Furthermore, photobleaching can convert high molecular weight CDOM species to lower molecular weight species, result in a great decrease in ag(295) than in ag(275), and consequently cause the spectral slope increase [30, 35, 36].
With respect to the ocean color applications, the spatial distribution pattern of CDOM absorption in upper layer is successfully captured by the algorithm, and the time-series CDOM variation which collaboratively changes with river discharge, is also reflected from satellite multi-senor imagery. Furthermore, the variation of CDOM spectral slope (Sg) within the PRE is depicted with very high spatial resolution, by retrieving Sg from satellite imagery pixel by pixel using Eq. (5). This can lead to much more delicate CDOM absorption products (Figures 4 and 5), because previous algorithms to derive CDOM absorption spectra generally use a uniform Sg for the entire water area of interest to calculate ag , which is evidently unreasonable for estuarine and coastal waters, where the CDOM of different sources may have different levels of Sg.
This research analyzes the variation of the CDOM optical properties (absorption) in the PRE based on the in-situ measurements during a spring cruise survey and the ocean color interpretations. Based on a UV to visible scheme to retrieve CDOM absorption from remote sensing reflectance, the time-series product of CDOM absorption in the PRE waters is produced from satellite images with multiple spatial and spectral resolutions, and an evident correlation with the temporal variation of river discharge is captured in the time-series variation of CDOM absorption in the PRE in spring.
This study was supported by the National Natural Science Foundation of China (grant number 41376035), the General Research Fund of Hong Kong Research Grants Council (RGC) (grant numbers CUHK 14303818, 402912, and 403113), and the talent startup fund of Jiangxi Normal University.
Appendices and nomenclature
absorbance absorption coefficient of CDOM chromophoric dissolved organic matter dissolved organic carbon dissolved organic matter hyperspectral imager for the coastal ocean International Space Station mean absolute percent difference Ocean and Land Colour Instrument Operational Land Imager ultraviolet Visible Infrared Imaging Radiometer Suite
absorption coefficient of CDOM
chromophoric dissolved organic matter
dissolved organic carbon
dissolved organic matter
hyperspectral imager for the coastal ocean
International Space Station
mean absolute percent difference
Ocean and Land Colour Instrument
Operational Land Imager
Visible Infrared Imaging Radiometer Suite