Measured and modeled radiometric quantities in coastal waters: toward a closure

Sea-surface reflectance factor: replicability of computed values

The sea-surface reflectance factor ρ required for the determination of the water- leaving radiance from above-water radiometric measurements is derived from radiative transfer simulations relying on models of the sky-radiance distribution and sea-surface statistics. This work primarily investigates the impact on ρ of various sky-radiance and sea-surface models. A specific replicability analysis, restricted to the 550 nm wavelength, has been performed with the Monte Carlo code for Ocean Color Simulations (so-called MOX) accounting for the measurement geometry recommended in protocols for the validation of satellite ocean color data and commonly applied for operational measurements. Results indicate that the variability of ρ increases with wind speed and reaches the largest values for sun elevations close to the zenith or approaching the horizon. In particular, a variability up to about 2% is observed for wind speeds below 4 ms^-1 and sun zenith angles larger than 20°. Finally, the benchmark of the ρ values from this study with those formally determined with the Hydrolight code and widely utilized by the ocean color community, exhibits systematic differences. The source of these differences is discussed and the implications for field measurements are addressed.

Adjacency radiance around a small island: implications for system vicarious calibrations

The adjacency radiance field surrounding a small island (i.e., the Lampedusa Island in the Central Mediterranean Sea) was theoretically analyzed to address implications on a hypothetical nearby system vicarious calibration (SVC) infrastructure for satellite ocean color sensors. Simulations, performed in the visible and near-infrared regions for the Ocean Land Color Instrument (OLCI) operated onboard Sentinel-3 satellites, show different patterns of adjacency effects (AE) around the island. In the direction of the reflected sunbeam (i.e., in the north-western region), AE mainly originate by missing glint contributions from the sea surface masked by the island. These AE are mainly negative, decrease with wavelength, and strongly depend on sea surface anisotropy (i.e., sea state) and illumination conditions; this hinders the capability to provide a general unique description of their features. In the remaining marine regions, AE are positive and do not exceed the radiometric sensitivity of OLCI data beyond approximately 14 km from the coast. At shorter distances, uncertainties in satellite radiance due to AE would hence not allow fulfilling requirements for SVC.

Characterization of the Light Field and Apparent Optical Properties in the Ocean Euphotic Layer Based on Hyperspectral Measurements of Irradiance Quartet

Although the light fields and apparent optical properties (AOPs) within the ocean euphotic layer have been studied for many decades through extensive measurements and theoretical modeling, there is virtually a lack of simultaneous high spectral resolution measurements of plane and scalar downwelling and upwelling irradiances (the so-called irradiance quartet). We describe a unique dataset of hyperspectral irradiance quartet, which was acquired under a broad range of environmental conditions within the water column from the near-surface depths to about 80 m in the Gulf of California. This dataset enabled the characterization of a comprehensive suite of AOPs for realistic non-uniform vertical distributions of seawater inherent optical properties (IOPs) and chlorophyll-a concentration (Chl) in the common presence of inelastic radiative processes within the water column, in particular Raman scattering by water molecules and chlorophyll-a fluorescence. In the blue and green spectral regions, the vertical patterns of AOPs are driven primarily by IOPs of seawater with weak or no discernible effects of inelastic processes. In the red, the light field and AOPs are strongly affected or totally dominated by inelastic processes of Raman scattering by water molecules, and additionally by chlorophyll-a fluorescence within the fluorescence emission band. The strongest effects occur in the chlorophyll-a fluorescence band within the chlorophyll-a maximum layer, where the average cosines of the light field approach the values of uniform light field, irradiance reflectance is exceptionally high approaching 1, and the diffuse attenuation coefficients for various irradiances are exceptionally low, including the negative values for the attenuation of upwelling plane and scalar irradiances. We established the empirical relationships describing the vertical patterns of some AOPs in the red spectral region as well as the relationships between some AOPs which can be useful in common experimental situations when only the downwelling plane irradiance measurements are available. We also demonstrated the applicability of irradiance quartet data in conjunction with Gershun’s equation for estimating the absorption coefficient of seawater in the blue-green spectral region, in which the effects of inelastic processes are weak or negligible.

On the minimization of adjacency effects in SeaWiFS primary data products from coastal areas

The minimization of adjacency effects (AE) in SeaWiFS primary products at the Aqua Alta Oceanographic Tower (AAOT) was investigated using sample images concurrent with in situ measurements. The validation exercise was performed with the NASA SeaDAS processing scheme ingesting original SeaWiFS data and alternatively SeaWiFS top-of-atmosphere data corrected for AE, and additionally including and excluding the default turbid water (TW) correction algorithm. Results show overestimates of the TW contributions partially compensating for AE. The analysis also suggests that intra-annual biases observed in SeaWiFS radiometric products at the AAOT may result from a misinterpretation of the NIR atmospheric signal as water contribution in data acquired in winter, and from uncompensated AE in data acquired in summer.

On the detectability of adjacency effects in ocean color remote sensing of mid-latitude coastal environments by SeaWiFS, MODIS-A, MERIS, OLCI, OLI and MSI

The detectability of adjacency effects (AE) in ocean color remote sensing by SeaWiFS, MODIS-A, MERIS, OLCI, OLI and MSI is theoretically assessed for typical observation conditions up to 36 km offshore (20 km for MSI). The methodology detailed in Bulgarelli et al. (2014) is applied to expand previous investigations to the wide range of terrestrial land covers and water types usually encountered in mid-latitude coastal environments. Simulations fully account for multiple scattering within a stratified atmosphere bounded by a non-uniform reflecting surface, sea surface roughness, sun position and off-nadir sensor view. A harmonized comparison of AE is ensured by adjusting the radiometric sensitivity of each sensor to the same input radiance. Results show that average AE in data from MODIS-A, and from MERIS and OLCI in reduced spatial resolution, are still above the sensor noise level (NL) at 36 km offshore, except for AE caused by green vegetation at the red wavelengths. Conversely, in data from the less sensitive SeaWiFS, OLI and MSI sensors, as well as from MERIS and OLCI in full spatial resolution, sole AE caused by highly reflecting land covers (such as snow, dry vegetation, white sand and concrete) are above the sensor NL throughout the transect, while AE originated from green vegetation and bare soil at visible wavelengths may become lower than NL at close distance from the coast. Such a distance increases with the radiometric resolution of the sensor. It is finally observed that AE are slightly sensitive to the water type only at the blue wavelengths. Notably, for an atmospheric correction scheme inferring the aerosol properties from NIR data, perturbations induced by AE at NIR and visible wavelengths might compensate each other. As a consequence, biases induced by AE on radiometric products (e.g., the water-leaving radiance) are not strictly correlated to the intensity of the reflectance of the nearby land.

Re-examining the effect of particle phase functions on the remote-sensing reflectance

Even though it is well known that both the magnitude and detailed angular shape of scattering (phase function, PF), particularly in the backward angles, affect the color of the ocean, the current remote-sensing reflectance (R_rs) models typically account for the effect of its magnitude only through the backscattering coefficient (b_b). Using 116 volume scattering function (VSF) measurements previously collected in three coastal waters around the U.S. and in the water of the North Atlantic Ocean, we re-examined the effect of particle PF on R_rs in four scenarios. In each scenario, the magnitude of particle backscattering (i.e., b_bp) is known, but the knowledge on the angular shape of particle backscattering is assumed to increase from knowing nothing about the shape of particle PFs to partially knowing the particle backscattering ratio (B_p), the exact backscattering shape as defined by β˜_p(γ≥90°) (particle VSF normalized by the particle total scattering coefficient), and the exact backscattering shape as defined by the χ_p factor (particle VSF normalized by the particle backscattering coefficient). At sun zenith angle=30°, the nadir-viewed R_rs would vary up to 65%, 35%, 20%, and 10%, respectively, as the constraints on the shape of particle backscattering become increasingly stringent from scenarios 1 to 4. In all four scenarios, the R_rs variations increase with both viewing and sun angles and are most prominent in the direction opposite the sun. Our results show a greater impact of the measured particle PFs on R_rs than previously found, mainly because our VSF data show a much greater variability in B_p, β˜_p(γ≥90°), and χ_p than previously known. Among the uncertainties in R_rs due to the particle PFs, about 97% can be explained by χ_p, 90% by β˜_p(γ≥90°), and 27% by B_p. The results indicate that the uncertainty in ocean color remote sensing can be significantly constrained by accounting for χ_p of the VSFs.

Adjacency effects in satellite radiometric products from coastal waters: a theoretical analysis for the northern Adriatic Sea

Biases induced by land perturbations in satellite-derived water-leaving radiance are theoretically estimated for typical observation conditions in a coastal area of the northern Adriatic Sea hosting the Aqua Alta Oceanographic Tower (AAOT) validation site. Two different correction procedures are considered: not deriving (AC-1) or alternatively deriving (AC-2) the atmospheric properties from the remote sensing data. In both cases, biases due to adjacency effects largely increase by approaching the coast and with the satellite viewing angle. Conversely, the seasonal and spectral dependence of biases significantly differ between AC-1 and AC-2 schemes. For AC-1 schemes average biases are within ±5% throughout the transect at yellow-green wavelengths, but at the coast they can reach -21% and 34% at 412 and 670 nm, respectively, and exceed 100% at 865 nm. For AC-2 schemes, adjacency effects at those wavelengths from which atmospheric properties are inferred add significant perturbations. For the specific case of a correction scheme determining the atmospheric properties from the near-infrared region and by adopting a power-law spectral extrapolation of adjacency perturbations on the derived atmospheric radiance, average biases become all negative with values up to -60% and -74% at 412 and 670 nm at the coast, respectively. The seasonal trend of estimated biases at the AAOT is consistent with intra-annual variation of biases from match-ups between in situ and satellite products derived with SeaDAS from SeaWiFS and MODIS data. Nevertheless, estimated biases at blue wavelengths exceed systematic differences determined from match-up analysis. This may be explained by uncertainties and approximations in the simulation procedure, and by mechanisms of compensation introduced by the turbid water correction algorithm implemented in SeaDAS.

Absorption correction and phase function shape effects on the closure of apparent optical properties

We present a closure experiment between new inherent optical properties (IOPs: absorption a, scattering b, backscattering b_b) and apparent optical properties (AOPs: remote-sensing reflectance R_rs, irradiance reflectance R, and anisotropic factor at nadir Q_n) data of Ionian and Adriatic seawaters, from very clear to turbid waters, ranging across one order of magnitude in R_rs. The internal consistency of the IOP-AOP matchups was investigated though radiative transfer closure. Using the in situ IOPs, we predicted the AOPs with the commercial radiative transfer solver Hydrolight. Closure was limited by two unresolved issues, one regarding processing of in situ data and the other related to radiative transfer modeling. First, different correction methods of the absorption data measured by the Wetlabs ac-s produced high variations in simulated reflectances, reaching 40% for the highest reflectances in our dataset. Second, the lack of detailed volume scattering function measurements forces us to adopt analytical functions that are consistent with a given particle backscattering ratio. The analytical phase functions named Fournier-Forand and two-term Kopelevich presented reasonable angular shapes with respect to measurements at a few backward angles. Between these phase functions, induced changes were within 4% for R_rs, within 11% for R, and within 10% for Q_n. Additionally, closure of Q_n was generally not successful considering radiometric uncertainties. Simulated Q_n overestimated low values and underestimated high values, especially at 665 nm, where Hydrolight was unable to predict measured Q_n values greater than 6 sr. The physical nature of Q_n makes this mismatch almost independent of the measured IOPs, thus precluding Q_n tuning by varying the former. The non-closure of Q_n might be caused by an inaccurate phase function and, to a lesser extent, by the modeling of the incoming radiance. For the future, this remains the task of accurate absorption and phase function measurements, especially at red wavelengths.

Optical closure in marine waters from in situ inherent optical property measurements

Optical closure using radiative transfer simulations can be used to determine the consistency of in situ measurements of inherent optical properties (IOPs) and radiometry. Three scattering corrections are applied to in situ absorption and attenuation profile data for a range of coastal and oceanic waters, but are found to have only very limited impact on subsequent closure attempts for these stations. Best-fit regressions on log-transformed measured and modelled downwards irradiance, Ed, and upwards radiance, Lu, profiles have median slopes between 0.92 - 1.24, revealing a tendency to underestimate Ed and Lu with depth. This is only partly explained by non-inclusion of fluorescence emission from CDOM and chlorophyll in the simulations. There are several stations where multiple volume scattering function related data processing steps perform poorly which suggests the potential existence of unresolved features in the modelling of the angular distribution of scattered photons. General optical closure therefore remains problematic, even though there are many cases in the data set where the match between measured and modelled radiometric data is within 25% RMS%E. These results are significant for applications that rely on optical closure e.g. assimilating ocean colour data into coupled physical-ecosystem models.

Simulation and analysis of adjacency effects in coastal waters: a case study

A methodology has been developed and applied to accurately quantify and analyze adjacency effects in satellite ocean color data for a set of realistic and representative observation conditions in the northern Adriatic Sea. The procedure properly accounts for sea surface reflectance anisotropy, off-nadir views, coastal morphology, and atmospheric multiple scattering. The study further includes a sensitivity analysis on commonly applied approximations. Results indicate that, within the accuracy limits defined by the radiometric resolution of ocean color sensors, adjacency effects in coastal waters might be significant at both visible and near-infrared wavelengths up to several kilometers off the coast. These results additionally highlight a significant dependence on the angle of observation, on the directional reflectance properties of the sea surface, and on the atmospheric multiple scattering.

Monte Carlo code for high spatial resolution ocean color simulations

A Monte Carlo code for ocean color simulations has been developed to model in-water radiometric fields of downward and upward irradiance (E(d) and E(u)), and upwelling radiance (L(u)) in a two-dimensional domain with a high spatial resolution. The efficiency of the code has been optimized by applying state-of-the-art computing solutions, while the accuracy of simulation results has been quantified through benchmark with the widely used Hydrolight code for various values of seawater inherent optical properties and different illumination conditions. Considering a seawater single scattering albedo of 0.9, as well as surface waves of 5 m width and 0.5 m height, the study has shown that the number of photons required to quantify uncertainties induced by wave focusing effects on E(d), E(u), and L(u) data products is of the order of 10(6), 10(9), and 10(10), respectively. On this basis, the effects of sea-surface geometries on radiometric quantities have been investigated for different surface gravity waves. Data products from simulated radiometric profiles have finally been analyzed as a function of the deployment speed and sampling frequency of current free-fall systems in view of providing recommendations to improve measurement protocols.

Measurements and modeling of the volume scattering function in the coastal northern Adriatic Sea

We performed measurements of the volume scattering function (VSF) between 0.5 degrees and 179 degrees with an angular resolution of 0.3 degrees in the northern Adriatic Sea onboard an oceanographic platform during three different seasons, using the multispectral volume scattering meter (MVSM) instrument. We observed important differences with respect to Petzold's commonly used functions, whereas the Fournier-Forand's analytical formulation provided a rather good description of the measured VSF. The comparison of the derived scattering, b(p)(lambda) and backscattering, b(bp)(lambda) coefficients for particles with the measurements performed with the classical AC-9 and Hydroscat-6 showed agreement to within 20%. The use of an empirical relationship for the derivation of b(b)(lambda) from beta(psi,lambda) at psi=140 degrees was validated for this coastal site although psi=118 degrees was confirmed to be the most appropriate angle. The low value of the factor used to convert beta(psi,lambda) into b(b)(lambda) within the Hydroscat-6 processing partially contributed to the underestimation of b(b)(lambda) with respect to the MVSM. Finally, use of the Kopelevich model together with a measurement of b(p)(lambda) at lambda=555 nm allowed us to reconstruct the VSF with average rms percent differences between 8 and 15%.

Effects of cosine error in irradiance measurements from field ocean color radiometers

The cosine error of in situ seven-channel radiometers designed to measure the in-air downward irradiance for ocean color applications was investigated in the 412-683 nm spectral range with a sample of three instruments. The interchannel variability of cosine errors showed values generally lower than +/-3% below 50 degrees incidence angle with extreme values of approximately 4-20% (absolute) at 50-80 degrees for the channels at 412 and 443 nm. The intrachannel variability, estimated from the standard deviation of the cosine errors of different sensors for each center wavelength, displayed values generally lower than 2% for incidence angles up to 50 degrees and occasionally increasing up to 6% at 80 degrees. Simulations of total downward irradiance measurements, accounting for average angular responses of the investigated radiometers, were made with an accurate radiative transfer code. The estimated errors showed a significant dependence on wavelength, sun zenith, and aerosol optical thickness. For a clear sky maritime atmosphere, these errors displayed values spectrally varying and generally within +/-3%, with extreme values of approximately 4-10% (absolute) at 40-80 degrees sun zenith for the channels at 412 and 443 nm. Schemes for minimizing the cosine errors have also been proposed and discussed.