Demonstration of a Modular Prototype End-to-End Simulator for Aquatic Remote Sensing Applications.

This study introduces a prototype end-to-end Simulator software tool for simulating two-dimensional satellite multispectral imagery for a variety of satellite instrument models in aquatic environments. Using case studies, the impact of variable sensor configurations on the performance of value-added products for challenging applications, such as coral reefs and cyanobacterial algal blooms, is assessed. This demonstrates how decisions regarding satellite sensor design, driven by cost constraints, directly influence the quality of value-added remote sensing products. Furthermore, the Simulator is used to identify situations where retrieval algorithms require further parameterization before application to unsimulated satellite data, where error sources cannot always be identified or isolated. The application of the Simulator can verify whether a given instrument design meets the performance requirements of end-users before build and launch, critically allowing for the justification of the cost and specifications for planned and future sensors. It is hoped that the Simulator will enable engineers and scientists to understand important design trade-offs in phase 0/A studies easily, quickly, reliably, and accurately in future Earth observation satellites and systems.

Diel changes of the optical backscattering coefficient of oceanic particulate matter determined from diel changes in apparent optical properties: a case study in the Mediterranean Sea (BOUSSOLE site).

Using in situ measurements of radiometric quantities and of the optical backscattering coefficient of particulate matter (bbp) at an oceanic site, we show that diel cycles of bbp are large enough to generate measurable diel variability of the ocean reflectance. This means that biogeochemical quantities such as net phytoplankton primary production, which are derivable from the diel bbp signal, can be potentially derived also from the diel variability of ocean color radiometry (OCR). This is a promising avenue for basin-scale quantification of such quantities because OCR is now performed from geostationary platforms that enable quantification of diel changes in the ocean reflectance over large ocean expanses. To assess the feasibility of this inversion, we applied three numerical inversion algorithms to derive bbp from measured reflectance data. The uncertainty in deriving bbp transfers to the retrieval of its diel cycle, making the performance of the inversion better in the green part of the spectrum (555 nm), with correlation coefficients >0.75 and a variability of 40% between the observed and derived bbp diel changes. While the results are encouraging, they also emphasize the inherent limitation of current inversion algorithms in deriving diel changes of bbp, which essentially stems from the empirical parameterizations that such algorithms include.

Uncertainty assessment of unattended above-water radiometric data collection from research vessels with the Dynamic Above-water Radiance (L) and Irradiance (E) Collector (DALEC).

We used above- and below-water radiometry measurements collected during a research voyage in the eastern Indian Ocean to assess uncertainties in deriving the remote sensing reflectance, Rrs, from unattended above-water radiometric data collection with the In-Situ Marine Optics Pty. Ltd. (IMO) Dynamic Above-water Radiance (L) and Irradiance (E) Collector (DALEC). To achieve this, the Rrs values derived from using the latest version of this hyperspectral radiometer were compared to values obtained from two in-water profiling radiometer systems of rather general use in the ocean optics research community, i.e., the Biospherical Instruments Inc. Compact Optical Profiling System (C-OPS) and the Seabird HyperPro II. Our results show that unattended, carefully quality-controlled, DALEC measurements provide Rrs for wavelengths < 600 nm that match those derived from the in-water systems with no bias and a dispersion of about 8%, provided that the appropriate technique is used to quantify the contribution of sky light reflection to the measured signal. The dispersion is larger (25-50%) for red bands, which is expected for clear oligotrophic waters as encountered during the voyage, where ∼2 10-5 < Rrs < ∼2 10-4 sr-1. For comparison, the two in-water systems provided Rrs in agreement within 4% for wavelengths < 600 nm.

High latitude Southern Ocean phytoplankton have distinctive bio-optical properties.

Studying the biogeochemistry of the Southern Ocean using remote sensing relies on accurate interpretation of ocean colour through bio-optical and biogeochemical relationships between quantities and properties of interest. During the Antarctic Circumnavigation Expedition of the 2016/2017 Austral Summer, we collected a spatially comprehensive dataset of phytoplankton pigment concentrations, particulate absorption and particle size distribution and compared simple bio-optical and particle property relationships as a function of chlorophyll a. Similar to previous studies we find that the chlorophyll-specific phytoplankton absorption coefficient is significantly lower than in other oceans at comparable chlorophyll concentrations. This appears to be driven in part by lower concentrations of accessory pigments per unit chlorophyll a as well as increased pigment packaging due to relatively larger sized phytoplankton at low chlorophyll a than is typically observed in other oceans. We find that the contribution of microphytoplankton (>20 µm size) to chlorophyll a estimates of phytoplankton biomass is significantly higher than expected for the given chlorophyll a concentration, especially in higher latitudes south of the Southern Antarctic Circumpolar Current Front. Phytoplankton pigments are more packaged in larger cells, which resulted in a flattening of phytoplankton spectra as measured in these samples when compared to other ocean regions with similar chlorophyll a concentration. Additionally, we find that at high latitude locations in the Southern Ocean, pheopigment concentrations can exceed mono-vinyl chlorophyll a concentrations. Finally, we observed very different relationships between particle volume and chlorophyll a concentrations in high and low latitude Southern Ocean waters, driven by differences in phytoplankton community composition and acclimation to environmental conditions and varying contribution of non-algal particles to the particulate matter. Our data confirm that, as previously suggested, the relationships between bio-optical properties and chlorophyll a in the Southern Ocean are different to other oceans. In addition, distinct bio-optical properties were evident between high and low latitude regions of the Southern Ocean basin. Here we provide a region-specific set of power law functions describing the phytoplankton absorption spectrum as a function of chlorophyll a.

Diurnal variations of the optical properties of phytoplankton in a laboratory experiment and their implication for using inherent optical properties to measure biomass.

Diurnal variations of phytoplankton size distributions, chlorophyll, carbon and nitrogen content, in vivo fluorescence and associated optical absorption and scattering properties were observed in the laboratory to help understand in situ and spatial observations. We grew triplicate semi-continuous cultures of T. pseudonana, D. tertiolecta, P. tricornutum and E. huxleyi under a sinusoidal light regime. We observed diurnal variations in the particulate absorption (ap), scattering (bp), attenuation (cp), and backscattering coefficients (bbp), which correlate with carbon and Chl concentrations. Relative variations from sunrise of bbp are slightly lower than those of cp, suggesting that bbp diurnal increases observed in nature are partly caused by phytoplankton. Non-concurrent changes of carbon and Chl-specific backscattering and scattering coefficients and optical cross-sections however indicates that using backscattering to infer scattering or biomass must be done with care.

Evaluation of sea-surface photosynthetically available radiation algorithms under various sky conditions and solar elevations.

In this study, we report on the performance of satellite-based photosynthetically available radiation (PAR) algorithms used in published oceanic primary production models. The performance of these algorithms was evaluated using buoy observations under clear and cloudy skies, and for the particular case of low sun angles typically encountered at high latitudes or at moderate latitudes in winter. The PAR models consisted of (i) the standard one from the NASA-Ocean Biology Processing Group (OBPG), (ii) the Gregg and Carder (GC) semi-analytical clear-sky model, and (iii) look-up-tables based on the Santa Barbara DISORT atmospheric radiative transfer (SBDART) model. Various combinations of atmospheric inputs, empirical cloud corrections, and semi-analytical irradiance models yielded a total of 13 (11 + 2 developed in this study) different PAR products, which were compared with in situ measurements collected at high frequency (15 min) at a buoy site in the Mediterranean Sea (the "BOUée pour l'acquiSition d'une Série Optique à Long termE," or, "BOUSSOLE" site). An objective ranking method applied to the algorithm results indicated that seven PAR products out of 13 were well in agreement with the in situ measurements. Specifically, the OBPG method showed the best overall performance with a root mean square difference (RMSD) (bias) of 19.7% (6.6%) and 10% (6.3%) followed by the look-up-table method with a RMSD (bias) of 25.5% (6.8%) and 9.6% (2.6%) at daily and monthly scales, respectively. Among the four methods based on clear-sky PAR empirically corrected for cloud cover, the Dobson and Smith method consistently underestimated daily PAR while the Budyko formulation overestimated daily PAR. Empirically cloud-corrected methods using cloud fraction (CF) performed better under quasi-clear skies (CF<0.3) with an RMSD (bias) of 9.7%-14.8% (3.6%-11.3%) than under partially clear to cloudy skies (0.30.7), however, all methods showed larger RMSD differences (biases) ranging between 32% and 80.6% (-54.5%-8.7%). Finally, three methods tested for low sun elevations revealed systematic overestimation, and one method showed a systematic underestimation of daily PAR, with relative RMSDs as large as 50% under all sky conditions. Under partially clear to overcast conditions all the methods underestimated PAR. Model uncertainties predominantly depend on which cloud products were used.

Warmer, deeper, and greener mixed layers in the North Atlantic subpolar gyre over the last 50 years.

Shifts in global climate resonate in plankton dynamics, biogeochemical cycles, and marine food webs. We studied these linkages in the North Atlantic subpolar gyre (NASG), which hosts extensive phytoplankton blooms. We show that phytoplankton abundance increased since the 1960s in parallel to a deepening of the mixed layer and a strengthening of winds and heat losses from the ocean, as driven by the low frequency of the North Atlantic Oscillation (NAO). In parallel to these bottom-up processes, the top-down control of phytoplankton by copepods decreased over the same time period in the western NASG, following sea surface temperature changes typical of the Atlantic Multi-decadal Oscillation (AMO). While previous studies have hypothesized that climate-driven warming would facilitate seasonal stratification of surface waters and long-term phytoplankton increase in subpolar regions, here we show that deeper mixed layers in the NASG can be warmer and host a higher phytoplankton biomass. These results emphasize that different modes of climate variability regulate bottom-up (NAO control) and top-down (AMO control) forcing on phytoplankton at decadal timescales. As a consequence, different relationships between phytoplankton, zooplankton, and their physical environment appear subject to the disparate temporal scale of the observations (seasonal, interannual, or decadal). The prediction of phytoplankton response to climate change should be built upon what is learnt from observations at the longest timescales.

Multivariate approach for the retrieval of phytoplankton size structure from measured light absorption spectra in the Mediterranean Sea (BOUSSOLE site).

Models based on the multivariate partial least squares (PLS) regression technique are developed for the retrieval of phytoplankton size structure from measured light absorption spectra (BOUSSOLE site, northwestern Mediterranean Sea). PLS-models trained with data from the Mediterranean Sea showed good accuracy in retrieving, over the nine-year BOUSSOLE time series, the concentrations of total chlorophyll a [Tchl a], of the sum of seven diagnostic pigments and of pigments associated with micro, nano, and picophytoplankton size classes separately. PLS-models trained using either total particle or phytoplankton absorption spectra performed similarly, and both reproduced seasonal variations of biomass and size classes derived by high performance liquid chromatography. Satisfactory retrievals were also obtained using PLS-models trained with a data set including various locations of the world's oceans, with however a lower accuracy. These results open the way to an application of this method to absorption spectra derived from hyperspectral and field satellite radiance measurements.

Generalized ocean color inversion model for retrieving marine inherent optical properties.

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

Relationship of left ventricular outflow tract velocity time integral to treatment strategy in submassive and massive pulmonary embolism.

Pulmonary embolism is associated with high rates of mortality and morbidity. It is important to understand direct comparisons of current interventions to differentiate favorable outcomes and complications. The objective of this study was to compare ultrasound-accelerated thrombolysis versus systemic thrombolysis versus anticoagulation alone and their effect on left ventricular outflow tract velocity time integral. This was a retrospective cohort study of subjects ≥18 years of age with a diagnosis of submassive or massive pulmonary embolism. The primary outcome was the percent change in left ventricular outflow tract velocity time integral between pre- and post-treatment echocardiograms. Ultrasound-accelerated thrombolysis compared to anticoagulation had a greater improvement in left ventricular outflow tract velocity time integral, measured by percent change. No significant change was noted between the ultrasound-accelerated thrombolysis and systemic thrombolysis nor systemic thrombolysis and anticoagulation groups. Pulmonary artery systolic pressure only showed a significant reduction in the ultrasound-accelerated thrombolysis versus anticoagulation group. The percent change of right ventricular to left ventricular ratios was improved when systemic thrombolysis was compared to both ultrasound-accelerated thrombolysis and anticoagulation. In this retrospective study of submassive or massive pulmonary embolisms, left ventricular outflow tract velocity time integral demonstrated greater improvement in patients treated with ultrasound-accelerated thrombolysis as compared to anticoagulation alone, a finding not seen with systemic thrombolysis. While this improvement in left ventricular outflow tract velocity time integral parallels the trend seen in mortality outcomes across the three groups, it only correlates with changes seen in pulmonary artery systolic pressure, not in other markers of echocardiographic right ventricular dysfunction (tricuspid annular plane systolic excursion and right ventricular to left ventricular ratios). Changes in left ventricular outflow tract velocity time integral, rather than echocardiographic markers of right ventricular dysfunction, may be considered a more useful prognostic marker of both dysfunction and improvement after reperfusion therapy.

An assessment of phytoplankton primary productivity in the Arctic Ocean from satellite ocean color/in situ chlorophyll-a based models.

We investigated 32 net primary productivity (NPP) models by assessing skills to reproduce integrated NPP in the Arctic Ocean. The models were provided with two sources each of surface chlorophyll-a concentration (chlorophyll), photosynthetically available radiation (PAR), sea surface temperature (SST), and mixed-layer depth (MLD). The models were most sensitive to uncertainties in surface chlorophyll, generally performing better with in situ chlorophyll than with satellite-derived values. They were much less sensitive to uncertainties in PAR, SST, and MLD, possibly due to relatively narrow ranges of input data and/or relatively little difference between input data sources. Regardless of type or complexity, most of the models were not able to fully reproduce the variability of in situ NPP, whereas some of them exhibited almost no bias (i.e., reproduced the mean of in situ NPP). The models performed relatively well in low-productivity seasons as well as in sea ice-covered/deep-water regions. Depth-resolved models correlated more with in situ NPP than other model types, but had a greater tendency to overestimate mean NPP whereas absorption-based models exhibited the lowest bias associated with weaker correlation. The models performed better when a subsurface chlorophyll-a maximum (SCM) was absent. As a group, the models overestimated mean NPP, however this was partly offset by some models underestimating NPP when a SCM was present. Our study suggests that NPP models need to be carefully tuned for the Arctic Ocean because most of the models performing relatively well were those that used Arctic-relevant parameters.

Shedding light on the sea: André Morel's legacy to optical oceanography.

André Morel (1933-2012) was a prominent pioneer of modern optical oceanography, enabling significant advances in this field. Through his forward thinking and research over more than 40 years, he made key contributions that this field needed to grow and to reach its current status. This article first summarizes his career and then successively covers different aspects of optical oceanography where he made significant contributions, from fundamental work on optical properties of water and particles to global oceanographic applications using satellite ocean color observations. At the end, we share our views on André's legacy to our research field and scientific community.

Climate-driven basin-scale decadal oscillations of oceanic phytoplankton.

Phytoplankton--the microalgae that populate the upper lit layers of the ocean--fuel the oceanic food web and affect oceanic and atmospheric carbon dioxide levels through photosynthetic carbon fixation. Here, we show that multidecadal changes in global phytoplankton abundances are related to basin-scale oscillations of the physical ocean, specifically the Pacific Decadal Oscillation and the Atlantic Multidecadal Oscillation. This relationship is revealed in approximately 20 years of satellite observations of chlorophyll and sea surface temperature. Interaction between the main pycnocline and the upper ocean seasonal mixed layer is one mechanism behind this correlation. Our findings provide a context for the interpretation of contemporary changes in global phytoplankton and should improve predictions of their future evolution with climate change.

Sources and assumptions for the vicarious calibration of ocean color satellite observations.

Spaceborne ocean color sensors require vicarious calibration to sea-truth data to achieve accurate water-leaving radiance retrievals. The assumed requirements of an in situ data set necessary to achieve accurate vicarious calibration were set forth in a series of papers and reports developed nearly a decade ago, which were embodied in the development and site location of the Marine Optical BuoY (MOBY). Since that time, NASA has successfully used data collected by MOBY as the sole source of sea-truth data for vicarious calibration of the Sea-viewing Wide field-of-view Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer instruments. In this paper, we make use of the 10-year, global time series of SeaWiFS measurements to test the sensitivity of vicarious calibration to the assumptions inherent in the in situ requirements (e.g., very low chlorophyll waters, hyperspectral measurements). Our study utilized field measurements from a variety of sources with sufficient diversity in data collection methods and geophysical variability to challenge those in situ restrictions. We found that some requirements could be relaxed without compromising the ability to vicariously calibrate to the level required for accurate water-leaving radiance retrievals from satellite-based sensors.

Bidirectional reflectance of oceanic waters: accounting for Raman emission and varying particle scattering phase function.

The bidirectionality of the upward radiance field in oceanic case 1 waters has been reinvestigated by incorporation of revised parameterizations of inherent optical properties as a function of the chlorophyll concentration (Chl), considering Raman scattering and making the particle phase function shape (beta(rho)) continuously varying along with the Chl. Internal consistency is thus reached, as the decrease in backscattering probability (for increasing Chl) translates into a correlative change in beta(rho). The single particle phase function (previously used) precluded a realistic assessment of bidirectionality for waters with Chl > 1 mg m(-3). This limitation is now removed. For low Chl, Raman emissions significantly affect the radiance field. For moderate Chl (0.1-1 mg m(-3)), new and previous bidirectional parameters remain close. The ocean reflectance anisotropy has implications in ocean color remote-sensing problems, in derivation of coherent water-leaving radiances, in associated calibration-validation activities, and in the merging of data obtained under various geometrical configurations.