Ocean colour is recognised as an Essential Climate Variable (ECV) by the Global Climate Observing System (GCOS); and spectrally-resolved water-leaving radiances (or remote-sensing reflectances) in the visible domain, and chlorophyll-a concentration are identified as required ECV products. Time series of the products at the global scale and at high spatial resolution, derived from ocean-colour data, are key to studying the dynamics of phytoplankton at seasonal and inter-annual scales; their role in marine biogeochemistry; the global carbon cycle; the modulation of how phytoplankton distribute solar-induced heat in the upper layers of the ocean; and the response of the marine ecosystem to climate variability and change. However, generating a long time series of these products from ocean-colour data is not a trivial task: algorithms that are best suited for climate studies have to be selected from a number that are available for atmospheric correction of the satellite signal and for retrieval of chlorophyll-a concentration; since satellites have a finite life span, data from multiple sensors have to be merged to create a single time series, and any uncorrected inter-sensor biases could introduce artefacts in the series, e.g., different sensors monitor radiances at different wavebands such that producing a consistent time series of reflectances is not straightforward. Another requirement is that the products have to be validated against in situ observations. Furthermore, the uncertainties in the products have to be quantified, ideally on a pixel-by-pixel basis, to facilitate applications and interpretations that are consistent with the quality of the data. This paper outlines an approach that was adopted for generating an ocean-colour time series for climate studies, using data from the MERIS (MEdium spectral Resolution Imaging Spectrometer) sensor of the European Space Agency; the SeaWiFS (Sea-viewing Wide-Field-of-view Sensor) and MODIS-Aqua (Moderate-resolution Imaging Spectroradiometer-Aqua) sensors from the National Aeronautics and Space Administration (USA); and VIIRS (Visible and Infrared Imaging Radiometer Suite) from the National Oceanic and Atmospheric Administration (USA). The time series now covers the period from late 1997 to end of 2018. To ensure that the products meet, as well as possible, the requirements of the user community, marine-ecosystem modellers, and remote-sensing scientists were consulted at the outset on their immediate and longer-term requirements as well as on their expectations of ocean-colour data for use in climate research. Taking the user requirements into account, a series of objective criteria were established, against which available algorithms for processing ocean-colour data were evaluated and ranked. The algorithms that performed best with respect to the climate user requirements were selected to process data from the satellite sensors. Remote-sensing reflectance data from MODIS-Aqua, MERIS, and VIIRS were band-shifted to match the wavebands of SeaWiFS. Overlapping data were used to correct for mean biases between sensors at every pixel. The remote-sensing reflectance data derived from the sensors were merged, and the selected in-water algorithm was applied to the merged data to generate maps of chlorophyll concentration, inherent optical properties at SeaWiFS wavelengths, and the diffuse attenuation coefficient at 490 nm. The merged products were validated against in situ observations. The uncertainties established on the basis of comparisons with in situ data were combined with an optical classification of the remote-sensing reflectance data using a fuzzy-logic approach, and were used to generate uncertainties (root mean square difference and bias) for each product at each pixel.

A next-generation in-water profiler designed to measure the apparent optical properties (AOPs) of seawater was developed and validated across a wide dynamic range of in-water properties. The new free-falling instrument, the Compact-Optical Profiling System (C-OPS), was based on sensors built with a cluster of 19 state-of-the-art microradiometers spanning 320–780 nm and a novel kite-shaped backplane. The new backplane includes tunable ballast, a hydrobaric buoyancy chamber, plus pitch and roll adjustments, to provide unprecedented stability and vertical resolution in near-surface waters. A unique data set was collected as part of the development activity plus the first major field campaign that used the new instrument, the Malina expedition to the Beaufort Sea in the vicinity of the Mackenzie River outflow. The data were of sufficient resolution and quality to show that errors – more correctly, uncertainties – in the execution of data sampling protocols were measurable at the 1% and 1 cm level with C-OPS. A theoretical sensitivity analysis as a function of three water types established by the peak in the remote sensing reflectance spectrum, R_rs(λ), revealed which water types and which parts of the spectrum were the most sensitive to data acquisition uncertainties. Shallow riverine waters were the most sensitive water type, and the ultraviolet and near-infrared spectral end members, which are critical to next-generation satellite missions, were the most sensitive parts of the spectrum. The sensitivity analysis also showed how the use of data products based on band ratios significantly mitigated the influence of data acquisition uncertainties. The unprecedented vertical resolution provided high-quality data products, which supported an alternative classification capability based on the spectral diffuse attenuation coefficient, K_d(λ). The K_d(320) and K_d(780) data showed how complex coastal systems can be distinguished two-dimensionally and how near-ice water masses are different from the neighboring open ocean. Finally, an algorithm for predicting the spectral absorption due to colored dissolved organic matter (CDOM), denoted a_CDOM(λ), was developed using the K_d(320) / K_d(780) ratio, which was based on a linear relationship with respect to a_CDOM(440). The robustness of the approach was established by expanding the use of the algorithm to include a geographically different coastal environment, the Southern Mid-Atlantic Bight, with no significant change in accuracy (approximately 98% of the variance explained). Alternative spectral end members reminiscent of next-generation (340 and 710 nm) as well as legacy satellite missions (412 and 670 nm) were also used to accurately derive a_CDOM(440) from K_d(λ) ratios.

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end-members, i.e., a short- and long-wavelength pair. The global perspective is established by sampling “extremely” clear water plus a generalized extent in turbidity and optical properties that each span 3 decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the western coast of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses (i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent) plus 15 subcategories of more complex water masses that were not necessarily evolving conservatively. The dependence on optical complexity was confirmed with an objective methodology based on a cluster analysis technique. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four algorithms based on different validation quality end-members were validated with accuracies of 1.2 %–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing nonconservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1 mm (with pressure transducer precision of 0.03–0.08 mm) and a boundary constraint to mitigate wave-focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an in-water constituent in keeping with the accuracy obtained in the laboratory.

Presently, operational ocean color satellite sensors are designed with a legacy perspective for sampling the open ocean primarily in the visible domain, while high spatial resolution sensors such as Sentinel-2, Sentinel-3, and Landsat8 are increasingly used for observations of coastal and inland water quality. Next-generation satellites such as the NASA Plankton, Aerosol, Cloud and ocean Ecosystem (PACE) and Surface Biology and Geology (SBG) sensors are anticipated to increase spatial and/or spectral resolution. An important consideration is determining the minimum signal-to-noise ratio (SNR) needed to retrieve typical biogeochemical products, such as biomass, in aquatic systems, and whether legacy sensors can be used for algorithm development. Here, we evaluate SNR and remote-sensing reflectance (R_rs) uncertainty for representative bright and dim targets in coastal California, USA. The majority of existing sensors fail to meet proposed criteria. Despite these limitations, uncertainties in retrieved biomass as chlorophyll or normalized difference vegetation index (NDVI) remain well below a proposed threshold of 17.5%, suggesting that existing sensors can be used in coastal systems. Existing commercially available in-water and airborne instrument suites can exceed all proposed thresholds for SNR and R_rs uncertainty, providing a path forward for collection of calibration and validation data for future satellite missions.

The launch of the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) and the Surface Biology and Geology (SBG) satellite sensors will provide increased spectral resolution compared to existing platforms. These new sensors will require robust calibration and validation datasets, but existing field-based instrumentation is limited in its availability and potential for geographic coverage, particularly for coastal and inland waters, where optical complexity is substantially greater than in the open ocean. The minimum signal-to-noise ratio (SNR) is an important metric for assessing the reliability of derived biogeochemical products and their subsequent use as proxies, such as for biomass, in aquatic systems. The SNR can provide insight into whether legacy sensors can be used for algorithm development as well as calibration and validation activities for next-generation platforms. We extend our previous evaluation of SNR and associated uncertainties for representative coastal and inland targets to include the imaging sensors PRISM and AVIRIS-NG, the airborne-deployed C-AIR radiometers, and the shipboard HydroRad and HyperSAS radiometers, which were not included in the original analysis. Nearly all the assessed hyperspectral sensors fail to meet proposed criteria for SNR or uncertainty in remote sensing reflectance (R_rs) for some part of the spectrum, with the most common failures (>20% uncertainty) below 400 nm, but all the sensors were below the proposed 17.5% uncertainty for derived chlorophyll-a. Instrument suites for both in-water and airborne platforms that are capable of exceeding all the proposed thresholds for SNR and R_rs uncertainty are commercially available. Thus, there is a straightforward path to obtaining calibration and validation data for current and next-generation sensors, but the availability of suitable high spectral resolution sensors is limited.