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.

Inferring inherent optical properties and water constituent profiles from apparent optical properties.

The BP09 experiment conducted by the Centre for Maritime Research and Experimentation in the Ligurian Sea in March 2009 provided paired vertical profiles of nadir-viewing radiances L(u)(z) and downward irradiances E(d)(z) and inherent optical properties (IOPs, absorption, scattering and backscattering coefficients). An inversion algorithm was implemented to retrieve IOPs from apparent optical properties (AOPs, radiance reflectance R(L), irradiance reflectance R(E) and diffuse attenuation coefficient K(d)) derived from the radiometric measurements. Then another inversion algorithm was developed to infer vertical profiles of water constituent concentrations, including chlorophyll-a concentration, non-algal particle concentration, and colored dissolved organic matter from the retrieved IOPs based on a bio-optical model. The algorithm was tested on a synthetic dataset and found to give reliable results with an accuracy better than 1%. When the algorithm was applied to the BP09 dataset it was found that good retrievals of IOPs could be obtained for sufficiently deep waters, i.e. for L(u)(z) and E(d)(z) measurements conducted to depths of 50 m or more. This requirement needs to be satisfied in order to obtain a good estimation of the backscattering coefficient. For such radiometric measurements a correlation of 0.88, 0.96 and 0.93 was found between retrieved and measured absorption, scattering and backscattering coefficients, respectively. A comparison between water constituent values derived from the measured IOPs and in-situ measured values, yielded a correlation of 0.80, 0.78, and 0.73 for chlorophyll-a concentration, non-algal particle concentration, and absorption coefficient of colored dissolved organic matter at 443 nm, respectively. This comparison indicates that adjustments to the bio-optical model are needed in order to obtain a better match between inferred and measured water constituent values in the Ligurian Sea using the methodology developed in this paper.

Detailed validation of the bidirectional effect in various Case I and Case II waters.

Simulated bidirectional reflectance distribution functions (BRDF) were compared with measurements made just beneath the water's surface. In Case I water, the set of simulations that varied the particle scattering phase function depending on chlorophyll concentration agreed more closely with the data than other models. In Case II water, however, the simulations using fixed phase functions agreed well with the data and were nearly indistinguishable from each other, on average. The results suggest that BRDF corrections in Case II water are feasible using single, average, particle scattering phase functions, but that the existing approach using variable particle scattering phase functions is still warranted in Case I water.