The ESA Climate Change Initiative: Satellite Data Records for Essential Climate Variables

Observations of Earth from space have been made for over 40 years and have contributed to advances in many aspects of climate science. However, attempts to exploit this wealth of data are often hampered by a lack of homogeneity and continuity and by insufficient understanding of the products and their uncertainties. There is, therefore, a need to reassess and reprocess satellite datasets to maximize their usefulness for climate science. The European Space Agency has responded to this need by establishing the Climate Change Initiative (CCI). The CCI will create new climate data records for (currently) 13 essential climate variables (ECVs) and make these open and easily accessible to all. Each ECV project works closely with users to produce time series from the available satellite observations relevant to users' needs. A climate modeling users' group provides a climate system perspective and a forum to bring the data and modeling communities together. This paper presents the CCI program. It outlines its benefit and presents approaches and challenges for each ECV project, covering clouds, aerosols, ozone, greenhouse gases, sea surface temperature, ocean color, sea level, sea ice, land cover, fire, glaciers, soil moisture, and ice sheets. It also discusses how the CCI approach may contribute to defining and shaping future developments in Earth observation for climate science.

Sun glint estimation in marine satellite images: a comparison of results from calculation and radiative transfer modeling

The intensity and location of Sun glint in two Medium Resolution Imaging Spectrometer (MERIS) images was modeled using a radiative transfer model that includes elevation features as well as the slope of the sea surface. The results are compared to estimates made using glint flagging and correction approaches used within standard atmospheric correction processing code. The model estimate gives a glint pattern with a similar width but lower peak level than any current method, or than that estimated by a radiative transfer model with surfaces that include slope but not height. The MERIS third reprocessing recently adopted a new slope statistics model for Sun glint correction; the results show that this model is an outlier with respect to both the elevation model and other slope statistics models and we recommend that its adoption should be reviewed.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed. populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (similar to 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Coccolithophore surface distributions in the North Atlantic and their modulation of the air-sea flux of CO2 from 10 years of satellite Earth observation data

Coccolithophores are the primary oceanic phytoplankton responsible for the production of calcium carbonate (CaCO3). These climatically important plankton play a key role in the oceanic carbon cycle as a major contributor of carbon to the open ocean carbonate pump (similar to 50 %) and their calcification can affect the atmosphere-to-ocean (air-sea) uptake of carbon dioxide (CO2) through increasing the seawater partial pressure of CO2 (pCO(2)). Here we document variations in the areal extent of surface blooms of the globally important coccolithophore, Emiliania huxleyi, in the North Atlantic over a 10-year period (1998-2007), using Earth observation data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). We calculate the annual mean sea surface areal coverage of E. huxleyi in the North Atlantic to be 474 000 +/- 104 000 km(2), which results in a net CaCO3 carbon (CaCO3-C) production of 0.14-1.71 Tg CaCO3-C per year. However, this surface coverage (and, thus, net production) can fluctuate inter-annually by -54/+81% about the mean value and is strongly correlated with the El Nino/Southern Oscillation (ENSO) climate oscillation index (r = 0.75, p < 0.02). Our analysis evaluates the spatial extent over which the E. huxleyi blooms in the North Atlantic can increase the pCO(2) and, thus, decrease the localised air-sea flux of atmospheric CO2. In regions where the blooms are prevalent, the average reduction in the monthly air-sea CO2 flux can reach 55%. The maximum reduction of the monthly air-sea CO2 flux in the time series is 155 %. This work suggests that the high variability, frequency and distribution of these calcifying plankton and their impact on pCO(2) should be considered if we are to fully understand the variability of the North Atlantic air-to-sea flux of CO2. We estimate that these blooms can reduce the annual N. Atlantic net sink atmospheric CO2 by between 3-28 %.

An intercomparison of bio-optical techniques for detecting dominant phytoplankton size class from satellite remote sensing

Satellite remote sensing of ocean colour is the only method currently available for synoptically measuring wide-area properties of ocean ecosystems, such as phytoplankton chlorophyll biomass. Recently, a variety of bio-optical and ecological methods have been established that use satellite data to identify and differentiate between either phytoplankton functional types (PFTs) or phytoplankton size classes (PSCs). In this study, several of these techniques were evaluated against in situ observations to determine their ability to detect dominant phytoplankton size classes (micro-, nano- and picoplankton). The techniques are applied to a 10-year ocean-colour data series from the SeaWiFS satellite sensor and compared with in situ data (6504 samples) from a variety of locations in the global ocean. Results show that spectral-response, ecological and abundance-based approaches can all perform with similar accuracy. Detection of microplankton and picoplankton were generally better than detection of nanoplankton. Abundance-based approaches were shown to provide better spatial retrieval of PSCs. Individual model performance varied according to PSC, input satellite data sources and in situ validation data types. Uncertainty in the comparison procedure and data sources was considered. Improved availability of in situ observations would aid ongoing research in this field.

Monitoring changes in phytoplankton abundance and composition in the Northwest Atlantic: a comparison of results obtained by Continuous Plankton Recorder sampling and colour satellite imagery

Phytoplankton abundance in the NW Atlantic was measured by continuous plankton recorder (CPR) sampling along tracks between Iceland and the western Scotian Shelf from 1998 to 2006, when sea-surface chlorophyll (SSChl) measurements were also being made by ocean colour satellite imagery using the SeaWiFS sensor. Seasonal and inter-annual changes in phytoplankton abundance were examined using data collected by both techniques, averaged over each of four shelf regions and four deep ocean regions. CPR sampling had gaps (missing months) in all regions and in the four deep ocean regions satellite observations were too sparse between November and February to be of use. Average seasonal cycles of SSChl were similar to those of total diatom abundance in seven regions, to those of the phytoplankton colour index in six regions, but were not similar to those of total dinoflagellate abundance anywhere. Large inter-annual changes in spring bloom dynamics were captured by both samplers in shelf regions. Changes in annual (or 8 months) averages of SSChl did not generally follow those of the CPR indices within regions and multi-year averages of SSChl, and the three CPR indices were generally higher in shelf than in deep ocean regions. Remote sensing and CPR sampling provide complementary ways of monitoring phytoplankton in the ocean: the former has superior temporal and spatial coverage and temporal resolution, and the latter provides better taxonomic information.

From silk to satellite: half a century of ocean colour anomalies in the Northeast Atlantic

Changes in phytoplankton dynamics influence marine biogeochemical cycles, climate processes, and food webs, with substantial social and economic consequences. Large-scale estimation of phytoplankton biomass was possible via ocean colour measurements from two remote sensing satellites – the Coastal Zone Color Scanner (CZCS, 1979-1986) and the Sea-viewing Wide Field-of-view Sensor (SeaWiFS, 1998-2010). Due to the large gap between the two satellite eras and differences in sensor characteristics, comparison of the absolute values retrieved from the two instruments remains challenging. Using a unique in situ ocean colour dataset that spans more than half a century, the two satellite-derived chlorophyll-a (Chl-a) eras are linked to assess concurrent changes in phytoplankton variability and bloom timing over the Northeast Atlantic Ocean and North Sea. Results from this unique re-analysis reflect a clear increasing pattern of Chl-a, a merging of the two seasonal phytoplankton blooms producing a longer growing season and higher seasonal biomass, since the mid-1980s. The broader climate plays a key role in Chl-a variability as the ocean colour anomalies parallel the oscillations of the Northern Hemisphere Temperature (NHT) since 1948.

Macroscale factors affecting diatom abundance: a synergistic use of Continuous Plankton Recorder and satellite remote sensing data

Diatoms exist in almost every aquatic regime; they are responsible for 20% of global carbon fixation and 25% of global primary production, and are regarded as a key food for copepods, which are subsequently consumed by larger predators such as fish and marine mammals. A decreasing abundance and a vulnerability to climatic change in the North Atlantic Ocean have been reported in the literature. In the present work, a data matrix composed of concurrent satellite remote sensing and Continuous Plankton Recorder (CPR) in situ measurements was collated for the same spatial and temporal coverage in the Northeast Atlantic. Artificial neural networks (ANNs) were applied to recognize and learn the complex non-monotonic and non-linear relationships between diatom abundance and spatiotemporal environmental factors. Because of their ability to mimic non-linear systems, ANNs proved far more effective in modelling the diatom distribution in the marine ecosystem. The results of this study reveal that diatoms have a regular seasonal cycle, with their abundance most strongly influenced by sea surface temperature (SST) and light intensity. The models indicate that extreme positive SSTs decrease diatom abundances regardless of other climatic conditions. These results provide information on the ecology of diatoms that may advance our understanding of the potential response of diatoms to climatic change.

Coccolithophore surface distributions in the North Atlantic and their modulation of the air-sea flux of CO2 from 10 years of satellite Earth observation data

Coccolithophores are the primary oceanic phytoplankton responsible for the production of calcium carbonate (CaCO3). These climatically important plankton play a key role in the oceanic carbon cycle as a major contributor of carbon to the open ocean 5 carbonate pump (�50%) and their formation can affect the atmosphere-to-ocean (airsea) uptake of carbon dioxide (CO2) through increasing the seawater partial pressure of CO2 (pCO2). Here we document variations in the areal extent of surface blooms of the globally important coccolithophore, Emiliania huxleyi, in the North Atlantic over a 10-year period (1998–2007), using Earth observation data from the Sea-viewing Wide 10 Field of view Sensor (SeaWiFS).We calculate the annual mean surface areal coverage of E. huxleyi in the North Atlantic to be 474 000±119 000km2 yr−1, which results in a net CaCO3 production of 0.62±0.15 Tg CaCO3 carbon per year. However, this surface coverage and net production can fluctuate by −54/+81% about these mean values and are strongly correlated with the El Ni˜no/Southern Oscillation (ENSO) climate os15 cillation index (r =0.75, p<0.02). Our analysis evaluates the spatial extent over which the E. huxleyi blooms in the North Atlantic can increase the pCO2 and thus decrease the localised sink of atmospheric CO2. In regions where the blooms are prevalent, the average reduction in the monthly CO2 sink can reach 12 %. The maximum reduction of the monthly CO2 sink in the time series is 32 %. This work suggests that the high 20 variability, frequency and distribution of these calcifying plankton and their impact on pCO2 should be considered within modelling studies of the North Atlantic if we are to fully understand the variability of its air-to-sea CO2 flux.

Monitoring of a quasi-stationary eddy in the Bay of Biscay by means of satellite, in situ and model results

The presence of a quasi-stationary anticyclonic eddy within the southeastern Bay of Biscay (centred around 44°30′N-4°W) has been reported on various occasions in the bibliography. The analysis made in this study for the period 2003–2010, by using in situ and remote sensing measurements and model results shows that this mesoscale coherent structure is present almost every year from the end of winter-beginning of spring, to the beginning of fall. During this period it remains in an area limited to the east by the Landes Plateau, to the west by Le Danois Bank and Torrelavega canyon and to the northwest by the Jovellanos seamount. All the observations and analysis made in this contribution, suggest that this structure is generated between Capbreton and Torrelavega canyons. Detailed monitoring from in situ and remote sensing data of an anticyclonic quasi-stationary eddy, in 2008, shows the origin of this structure from a warm water current located around 43°42′N-3°30′W in mid-January. This coherent structure is monitored until August around the same area, where it has a marked influence on the Sea Level Anomaly, Sea Surface Temperature and surface Chlorophyll-a concentration. An eddy tracking method, applied to the outputs of a numerical model, shows that the model is able to reproduce this type of eddy, with similar 2D characteristics and lifetimes to that suggested by the observations and previous works. This is the case, for instance, of the simulated MAY04 eddy, which was generated in May 2004 around Torrelavega canyon and remained quasi-stationary in the area for 4 months. The diameter of this eddy ranged from 40 to 60 km, its azimuthal velocity was less than 20 cm s−1, its vertical extension reached 3000–3500 m depth during April and May and it was observed to interact with other coherent structures.

Satellite-detected fluorescence reveals global physiology of ocean phytoplankton

Phytoplankton photosynthesis links global ocean biology and climate-driven fluctuations in the physical environment. These interactions are largely expressed through changes in phytoplankton physiology, but physiological status has proven extremely challenging to characterize globally. Phytoplankton fluorescence does provide a rich source of physiological information long exploited in laboratory and field studies, and is now observed from space. Here we evaluate the physiological underpinnings of global variations in satellite-based phytoplankton chlorophyll fluorescence. The three dominant factors influencing fluorescence distributions are chlorophyll concentration, pigment packaging effects on light absorption, and light-dependent energy-quenching processes. After accounting for these three factors, resultant global distributions of quenching-corrected fluorescence quantum yields reveal a striking consistency with anticipated patterns of iron availability. High fluorescence quantum yields are typically found in low iron waters, while low quantum yields dominate regions where other environmental factors are most limiting to phytoplankton growth. Specific properties of photosynthetic membranes are discussed that provide a mechanistic view linking iron stress to satellite-detected fluorescence. Our results present satellite-based fluorescence as a valuable tool for evaluating nutrient stress predictions in ocean ecosystem models and give the first synoptic observational evidence that iron plays an important role in seasonal phytoplankton dynamics of the Indian Ocean. Satellite fluorescence may also provide a path for monitoring climate-phytoplankton physiology interactions and improving descriptions of phytoplankton light use efficiencies in ocean productivity models.

Global assessment of ocean carbon export by combining satellite observations and food-web models

The export of organic carbon from the surface ocean by sinking particles is an important, yet highly uncertain, component of the global carbon cycle. Here we introduce a mechanistic assessment of the global ocean carbon export using satellite observations, including determinations of net primary production and the slope of the particle size spectrum, to drive a food-web model that estimates the production of sinking zooplankton feces and algal aggregates comprising the sinking particle flux at the base of the euphotic zone. The synthesis of observations and models reveals fundamentally different and ecologically consistent regional-scale patterns in export and export efficiency not found in previous global carbon export assessments. The model reproduces regional-scale particle export field observations and predicts a climatological mean global carbon export from the euphotic zone of ~6 Pg C yr−1. Global export estimates show small variation (typically < 10%) to factor of 2 changes in model parameter values. The model is also robust to the choices of the satellite data products used and enables interannual changes to be quantified. The present synthesis of observations and models provides a path for quantifying the ocean's biological pump.