Ocean color remote sensing: choosing the correct depth weighting function

Combined Earth observations reveal the sequence of conditions leading to a large algal bloom in Lake Geneva

Freshwater algae exhibit complex dynamics, particularly in meso-oligotrophic lakes with sudden and dramatic increases in algal biomass following long periods of low background concentration. While the fundamental prerequisites for algal blooms, namely light and nutrient availability, are well-known, their specific causation involves an intricate chain of conditions. Here we examine a recent massive Uroglena bloom in Lake Geneva (Switzerland/France). We show that a certain sequence of meteorological conditions triggered this specific algal bloom event: heavy rainfall promoting excessive organic matter and nutrients loading, followed by wind-induced coastal upwelling, and a prolonged period of warm, calm weather. The combination of satellite remote sensing, in-situ measurements, ad-hoc biogeochemical analyses, and three-dimensional modeling proved invaluable in unraveling the complex dynamics of algal blooms highlighting the substantial role of littoral-pelagic connectivities in large low-nutrient lakes. These findings underscore the advantages of state-of-the-art multidisciplinary approaches for an improved understanding of dynamic systems as a whole.

Transmission across the water-air interface: resolving the impact of multiple interactions at the sea surface

A series of Monte Carlo and HydroLight radiative transfer simulations are used to demonstrate that the traditional form of the Fresnel transmission across the water-air interface is accurate. This contradicts assertions to the contrary in a recent paper [Opt. Express25, 27086 (2017)10.1364/OE.25.027086] that suggested that the impact of multiple surface interactions had previously been ignored and that the transmission factor was dependent upon the turbidity of the water.

Evaluation of forward reflectance models and empirical algorithms for chlorophyll concentration of stratified waters

For waters with stratified chlorophyll concentration (Chl), numerical simulations were carried out to gain insight into the forward models of subsurface reflectance and empirical algorithms for Chl from the ocean color. It is found that the Gordon and Clark (1980) forward model for reflectance using an equivalent homogeneous water with a weighted average Chl (⟨Chl⟩) as the input works well, but depending on the contribution of gelbstoff, the difference in reflectance between stratified and the equivalent homogeneous water can be more than 10%. Further, the attenuation of upward light is better approximated as ∼1.5times that of the diffuse attenuation coefficient of downwelling irradiance. On the other hand, although the forward model for reflectance developed in Zaneveld et al. [Opt. Express13, 9052 (2005)] using equivalent homogeneous water with a weighted average of the backscattering to absorption ratio as the input also works well, this model cannot be used to obtain equivalent ⟨Chl⟩ for reflectance. Further, for empirical Chl algorithms designed for "Case 1" waters, it has been discovered that, for surface Chl in a range of ∼0.06-22.0mg/m³, the predictability of surface Chl is basically the same as that of ⟨Chl⟩ from the blue-green band ratio or the band difference of reflectance. Because ⟨Chl⟩ is wavelength and weighting-formula dependent, and it is required to have profiles of both Chl and the optical properties, these results emphasize that for empirical Chl algorithms, it is easier, less ambiguous, and certainly more straightforward and simple to use surface Chl for algorithm development and then its evaluation, rather than to use ⟨Chl⟩, regardless of whether or not the water is stratified.

Evaluation of Weighting Average Functions as a Simplification of the Radiative Transfer Simulation in Vertically Inhomogeneous Eutrophic Waters

Current water color remote sensing algorithms typically do not consider the vertical variations of phytoplankton. Ecolight with a radiative transfer program was used to simulate the underwater light field of vertical inhomogeneous waters based on the optical properties of a eutrophic lake (i.e., Lake Chaohu, China). Results showed that the vertical distribution of chlorophyll-a (Chla(z)) can considerably affect spectrum shape and magnitude of apparent optical properties (AOPs), including subsurface remote sensing reflectance in water (rrs(λ, z)) and the diffuse attenuation coefficient (Kx(λ, z)). The vertical variations of Chla(z) changed the spectrum shapes of rrs(λ, z) at the green and red wavelengths with a maximum value at approximately 590 nm, and changed the Kx(λ, z) from blue to red wavelength range with no obvious spectral variation. The difference between rrs(λ, z) at depth z m and its asymptotic value (Δrrs(λ, z)) could reach to ~78% in highly stratified waters. Diffuse attenuation coefficient of downwelling plane irradiance (Kd(λ, z)) had larger vertical variations, especially near water surface, in highly stratified waters. Three weighting average functions performed well in less stratified waters, and the weighting average function proposed by Zaneveld et al., (2005) performed best in highly stratified waters. The total contribution of the first three layers to rrs(λ, 0−) was approximately 90%, but the contribution of each layer in the water column to rrs(λ, 0−) varied with wavelength, vertical distribution of Chla(z) profiles, concentration of suspended particulate inorganic matter (SPIM), and colored dissolved organic matter (CDOM). A simple stratified remote sensing reflectance model considering the vertical distribution of phytoplankton was built based on the contribution of each layer to rrs(λ, 0−).

The effect of dispersed Petrobaltic oil droplet size on photosynthetically active radiation in marine environment

Oil pollution in seawater, primarily visible on sea surface, becomes dispersed as an effect of wave mixing as well as chemical dispersant treatment, and forms spherical oil droplets. In this study, we examined the influence of oil droplet size of highly dispersed Petrobaltic crude on the underwater visible light flux and the inherent optical properties (IOPs) of seawater, including absorption, scattering, backscattering and attenuation coefficients. On the basis of measured data and Mie theory, we calculated the IOPs of dispersed Petrobaltic crude oil in constant concentration, but different log-normal size distributions. We also performed a radiative transfer analysis, in order to evaluate the influence on the downwelling irradiance Ed, remote sensing reflectance Rrs and diffuse reflectance R, using in situ data from the Baltic Sea. We found that during dispersion, there occurs a boundary size distribution characterized by a peak diameter d0  = 0.3 μm causing a maximum E d increase of 40% within 0.5-m depth, and the maximum Ed decrease of 100% at depths below 5 m. Moreover, we showed that the impact of size distribution on the "blue to green" ratios of Rrs and R varies from 24% increase to 27% decrease at the same crude oil concentration.

Application of remote sensing for the optimization of in-situ sampling for monitoring of phytoplankton abundance in a large lake

Directives and legislations worldwide aim at representatively and continuously monitoring the ecological status of surface waters. In many countries, chlorophyll-a concentrations (CHL) are used as an indicator of phytoplankton abundance and the trophic level of lakes or reservoirs. In-situ measurements of water quality parameters, however, are time-consuming, costly and of unknown but naturally limited spatial representativeness. In addition, the variety of the involved lab and field measurement methods and instruments complicates comparability and reproducibility. Taking Lake Geneva as an example, 1234 satellite images from the MERIS sensor on the Envisat satellite from 2002 to 2012 are used to quantify the spatial and temporal variations of CHL concentrations. Based on histograms of spring, summer and autumn CHL estimates, the spatial representativeness of two existing in-situ measurement locations is analysed. Appropriate sampling frequencies to capture CHL peaks are examined by means of statistical resampling. The approaches proposed allow determining optimal in-situ sampling locations and frequencies. Their generic nature allows for adaptation to other lakes, especially to establish new survey programmes where no previous records are available.

Modeling the effects of near-surface plumes of suspended particulate matter on remote-sensing reflectance of coastal waters

A radiative transfer model was applied to examine the effects of vertically stratified inherent optical properties of the water column associated with near-surface plumes of suspended particulate matter on spectral remote-sensing reflectance, R(rs)(λ), of coastal marine environments. The simulations for nonuniform ocean consisting of two layers with different concentrations of suspended particulate matter (SPM) are compared with simulations for a reference homogeneous ocean whose SPM is identical to the surface SPM of the two-layer cases. The near-surface plumes of particles are shown to exert significant influence on R(rs)(λ). The sensitivity of R(rs)(λ) to vertical profile of SPM is dependent on the optical beam attenuation coefficient within the top layer, c(1)(λ), thickness of the top layer, z(1), and the ratio of SPM in the underlying layer to that in the top layer, SPM(2)/SPM(1), as well as the wavelength of light, λ. We defined a dimensionless spectral parameter, P(λ)=c(1)(λ)×z(1)×(SPM(2)/SPM(1)), to quantify and examine the effects of these characteristics of the two-layer profile of SPM on the magnitude and spectral shape of R(rs)(λ). In general, the difference of R(rs)(λ) between the two-layer and uniform ocean decreases to zero with an increase in P(λ). For the interpretation of ocean color measurements of water column influenced by near-surface plumes of particles, another dimensionless parameter P'(λ) was introduced, which is a product of terms representing homogenous ocean and a change caused by the two-layer structure of SPM. Based on the analysis of this parameter, we found that for the two-layer ocean there is a good relationship between R(rs)(λ) in the red and near-infrared spectral regions and the parameters describing the SPM(z) profile, i.e., SPM(1), SPM(2), and z(1).

Comparison of chlorophyll a concentration detected by remote sensors and other chlorophyll indices in inhomogeneous turbid waters

A new analytical approach for retrieval of the vertically weighted chlorophyll a concentration (<linear span>Chl<linear span>(rs)) detected by remote sensors is presented. Model calculations were carried out for the turbid waters of Lake Kinneret, Israel, and showed that <linear span>Chl<linear span>(rs) may be replaced by the average chlorophyll a concentration (<linear span>Chl<linear span>(p)) within the upper "penetration layer" 0-Z(p). The study also showed a high correlation between <linear span>Chl<linear span>(rs) and Chl concentration averaged in the other depth layers, namely, the 0-1 m layer, the euphotic layer (0-Z(e)), and the production layer (0-Z(pr)). Our findings are closely related to models developed for the world ocean, with the exception of periods when the dinoflagellate Peridinium gatunense blooms in the lake. We showed the effect of the pattern of vertical Chl distributions within the penetration layer on the difference between <linear span>Chl<linear span>(rs) and other Chl indices was conspicuous when the Chl maximum was in the uppermost 0- m layer of the water column. We assume that the presented approaches are instrumental for further development of optimal, locally adapted algorithms for remote sensing of Chl in any type of natural waters.

Light scattering by coccoliths detached from Emiliania huxleyi

We used in situ radiance/irradiance profiles to retrieve profiles of the spectral backscattering coefficient for all particles in an E. huxleyi coccolithophore bloom off the coast of Plymouth, UK. At high detached coccolith concentrations the spectra of backscattering all showed a minimum near approximately 550 to 600 nm. Using flow cytometry estimates of the detached coccolith concentration, and assuming all of the backscattering (over and above the backscattering by the water itself) was due to detached coccoliths, we determined the upper limit of the backscattering cross section (sigma(b)) of individual coccoliths to be 0.123+/-0.039 microm(2)/coccolith at 500 nm. Physical models of detached coccoliths were then developed and the discrete dipole approximation was used to compute their average backscattering cross section in random orientation. The result was 0.092 microm(2) at 500 nm, with the computed sigma(b) displaying a spectral shape similar to the measurements, but with less apparent increase in backscattering toward the red. When sigma(b) is computed on a per mole of calcite, rather than a per coccolith basis, it agreed reasonably well with that determined for acid-labile backscattering at 632 nm averaged over several species of cultured calcifying algae. Intact coccolithophore cells were taken into account by arguing that coccoliths attached to coccolithophore cells (forming a "coccosphere") backscatter in a manner similar to free coccoliths in random orientation. Estimating the number of coccoliths per coccosphere and using the observed number of coccolithophore cells resulted is an apparent backscattering cross section at 500 nm of 0.114+/-0.013 microm(2)/coccolith, in satisfactory agreement with the measured backscattering.

Small-scale effects of underwater bubble clouds on ocean reflectance: 3-D modeling results

We examined the effect of individual bubble clouds on remote-sensing reflectance of the ocean with a 3-D Monte Carlo model of radiative transfer. The concentrations and size distribution of bubbles were defined based on acoustical measurements of bubbles in the surface ocean. The light scattering properties of bubbles for various void fractions were calculated using Mie scattering theory. We show how the spatial pattern, magnitude, and spectral behavior of remote-sensing reflectance produced by modeled bubble clouds change due to variations in their geometric and optical properties as well as the background optical properties of the ambient water. We also determined that for realistic sizes of bubble clouds, a plane-parallel horizontally homogeneous geometry (1-D radiative transfer model) is inadequate for modeling water-leaving radiance above the cloud.