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.

Land cover is one of the essential climate variables of the ESA Climate Change Initiative (CCI). In this context, the Land Cover CCI (LC CCI) project aims at building global land cover maps suitable for climate modeling based on Earth observation by satellite sensors. The challenge is to generate a set of successive maps that are both accurate and consistent over time. To do so, operational methods for the automated classification of optical images are investigated. The proposed approach consists of a locally trained classification using an automated selection of training samples from existing, but outdated land cover information. Combinations of local extraction (based on spatial criteria) and self-cleaning of training samples (based on spectral criteria) are quantitatively assessed. Two large study areas, one in Eurasia and the other in South America, are considered. The proposed morphological cleaning of the training samples leads to higher accuracies than the statistical outlier removal in the spectral domain. An optimal neighborhood has been identified for the local sample extraction. The results are coherent for the two test areas, showing an improvement of the overall accuracy compared with the original reference datasets and a significant reduction of macroscopic errors. More importantly, the proposed method partly controls the reliability of existing land cover maps as sources of training samples for supervised classification.

Accurate maps of surface water extent are of paramount importance for water management, satellite data processing and climate modeling. Several maps of water bodies based on remote sensing data have been released during the last decade. Nonetheless, none has a truly (90 ${}^{\circ }$ N/90 ${}^{\circ }$ S) global coverage while being thoroughly validated. This paper describes a global, spatially-complete (void-free) and accurate mask of inland/ocean water for the 2000–2012 period, built in the framework of the European Space Agency (ESA) Climate Change Initiative (CCI). This map results from the synergistic combination of multiple individual SAR and optical water body and auxiliary datasets. A key aspect of this work is the original and rigorous stratified random sampling designed for the quality assessment of binary classifications where one class is marginally distributed. Input and consolidated products were assessed qualitatively and quantitatively against a reference validation database of 2110 samples spread throughout the globe. Using all samples, overall accuracy was always very high among all products, between $98\%$ and $100\%$ . The CCI global map of open water bodies provided the best water class representation (F-score of $89\%$ ) compared to its constitutive inputs. When focusing on the challenging areas for water bodies’ mapping, such as shorelines, lakes and river banks, all products yielded substantially lower accuracy figures with overall accuracies ranging between $74\%$ and $89\%$ . The inland water area of the CCI global map of open water bodies was estimated to be 3.17 million km ${}^2$ ± 0.24 million km ${}^2$ . The dataset is freely available through the ESA CCI Land Cover viewer.

The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400–1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies—especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.

The OLIVE (On Line Interactive Validation Exercise) platform is dedicated to the validation of global biophysical products such as LAI (Leaf Area Index) and FAPAR (Fraction of Absorbed Photosynthetically Active Radiation). It was developed under the framework of the CEOS (Committee on Earth Observation Satellites) Land Product Validation (LPV) sub-group. OLIVE has three main objectives: (i) to provide a consistent and centralized information on the definition of the biophysical variables, as well as a description of the main available products and their performances (ii) to provide transparency and traceability by an online validation procedure compliant with the CEOS LPV and QA4EO (Quality Assurance for Earth Observation) recommendations (iii) and finally, to provide a tool to benchmark new products, update product validation results and host new ground measurement sites for accuracy assessment. The functionalities and algorithms of OLIVE are described to provide full transparency of its procedures to the community. The validation process and typical results are illustrated for three FAPAR products: GEOV1 (VEGETATION sensor), MGVIo (MERIS sensor) and MODIS collection 5 FPAR. OLIVE is available on the European Space Agency CAL/VAL portal), including full documentation, validation exercise results, and product extracts.

Shallow and deep lakes receive and recycle organic and inorganic substances from within the confines of these lakes, their watershed and beyond. Hence, a large range in absorption and scattering and extreme differences in optical variability can be found between and within global lakes. This poses a challenge for atmospheric correction and bio-optical algorithms applied to optical remote sensing for water quality monitoring applications. To optimize these applications for the wide variety of lake optical conditions, we adapted a spectral classification scheme based on the concept of optical water types. The optical water types were defined through a cluster analysis of in situ hyperspectral remote sensing reflectance spectra collected by partners and advisors of the European Union 7th Framework Programme (FP7) Global Lakes Sentinel Services (GLaSS) project. The method has been integrated in the Envisat-BEAM software and the Sentinel Application Platform (SNAP) and generates maps of water types from image data. Two variations of water type classification are provided: one based on area-normalized spectral reflectance focusing on spectral shape (6CN, six-class normalized) and one that retains magnitude with no modification to the reflectance signal (6C). This resulted in a protocol, or processing scheme, that can also be applied or adapted for Sentinel-3 Ocean and Land Colour Imager (OLCI) datasets. We apply both treatments to MERIS imagery of a variety of European lakes to demonstrate its applicability. The studied target lakes cover a range of biophysical types, from shallow turbid to deep and clear, as well as eutrophic and dark absorbing waters, rich in colored dissolved organic matter (CDOM). In shallow, high-reflecting Dutch and Estonian lakes with high sediment load, 6C performed better, while in deep, low-reflecting clear Italian and Swedish lakes, 6CN performed better. The 6CN classification of in situ data is promising for very dark, high CDOM, absorbing lakes, but we show that our atmospheric correction of the imagery was insufficient to corroborate this. We anticipate that the application of the protocol to other lakes with unknown in-water characterization, but with comparable biophysical properties will suggest similar atmospheric correction (AC) and in-water retrieval algorithms for global lakes.