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.

Mighty small: Observing and modeling individual microbes becomes big science

Progress in microbiology has always been driven by technological advances, ever since Antonie van Leeuwenhoek discovered bacteria by making an improved compound microscope. However, until very recently we have not been able to identify microbes and record their mostly invisible activities, such as nutrient consumption or toxin production on the level of the single cell, not even in the laboratory. This is now changing with the rapid rise of exciting new technologies for single-cell microbiology (1, 2), which enable microbiologists to do what plant and animal ecologists have been doing for a long time: observe who does what, when, where, and next to whom. Single cells taken from the environment can be identified and even their genomes sequenced. Ex situ, their size, elemental, and biochemical composition, as well as other characteristics can be measured with high-throughput and cells sorted accordingly. Even better, individual microbes can be observed in situ with a range of novel microscopic and spectroscopic methods, enabling localization, identification, or functional characterization of cells in a natural sample, combined with detecting uptake of labeled compounds. Alternatively, they can be placed into fabricated microfluidic environments, where they can be positioned, exposed to stimuli, monitored, and their interactions controlled “in microfluido.” By introducing genetically engineered reporter cells into a fabricated landscape or a microcosm taken from nature, their reproductive success or activity can be followed, or their sensing of their local environment recorded.

Climate change impacts on sea–air fluxes of CO2 in three Arctic seas: a sensitivity study using Earth observation

We applied coincident Earth observation data collected during 2008 and 2009 from multiple sensors (RA2, AATSR and MERIS, mounted on the European Space Agency satellite Envisat) to characterise environmental conditions and integrated sea-air fluxes of CO2 in three Arctic seas (Greenland, Barents, Kara). We assessed net CO2 sink sensitivity due to changes in temperature, salinity and sea ice duration arising from future climate scenarios. During the study period the Greenland and Barents seas were net sinks for atmospheric CO2, with integrated sea-air fluxes of -36 +/- 14 and -11 +/- 5 Tg C yr(-1), respectively, and the Kara Sea was a weak net CO2 source with an integrated sea-air flux of +2.2 +/- 1.4 TgC yr(-1). The combined integrated CO2 sea-air flux from all three was -45 +/- 18 TgC yr(-1). In a sensitivity analysis we varied temperature, salinity and sea ice duration. Variations in temperature and salinity led to modification of the transfer velocity, solubility and partial pressure of CO2 taking into account the resultant variations in alkalinity and dissolved organic carbon (DOC). Our results showed that warming had a strong positive effect on the annual integrated sea-air flux of CO2 (i.e. reducing the sink), freshening had a strong negative effect and reduced sea ice duration had a small but measurable positive effect. In the climate change scenario examined, the effects of warming in just over a decade of climate change up to 2020 outweighed the combined effects of freshening and reduced sea ice duration. Collectively these effects gave an integrated sea-air flux change of +4.0 TgC in the Greenland Sea, +6.0 Tg C in the Barents Sea and +1.7 Tg C in the Kara Sea, reducing the Greenland and Barents sinks by 11% and 53 %, respectively, and increasing the weak Kara Sea source by 81 %. Overall, the regional integrated flux changed by +11.7 Tg C, which is a 26% reduction in the regional sink. In terms of CO2 sink strength, we conclude that the Barents Sea is the most susceptible of the three regions to the climate changes examined. Our results imply that the region will cease to be a net CO2 sink in the 2050s.

Above-water reflectance for the evaluation of adjacency effects in Earth observation data: initial results and methods comparison for near-coastal waters in the Western Channel, UK

Un-supervised hyperspectral remote-sensing reflectance data (<15 km from the shore) were collected from a moving research vessel. Two different processing methods were compared. The results were similar to concurrent Aqua-MODIS and Suomi-NPP-VIIRS satellite data.

Space-borne study of seasonal, multi-year, and decadal phytoplankton dynamics in the Bay of Biscay

Seasonal and inter-annual variations in phytoplankton community abundance in the Bay of Biscay are studied. Preliminarily processed by the National Aeronautics and Space Administration (NASA) to yield normalized water-leaving radiance and the top-of-the-atmosphere solar radiance, Sea-viewing Wide Field-of-View Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Coastal Zone Color Scanner (CZCS) data are further supplied to our dedicated retrieval algorithms to infer the sought for parameters. By applying the National Oceanic and Atmospheric Administration's (NOAA's) Advanced Very High Resolution Radiometer (AVHRR) data, the surface reflection coefficient in the only band in the visible spectrum is derived and employed for analysis. Decadal bridged time series of variations of diatom-dominated phytoplankton and green dinoflagellate Lepidodinium chlorophorum within the shelf zone and the coccolithophore Emiliania huxleyi in the pelagic area of the Bay are documented and analysed in terms of impacts of some biogeochemical and geophysical forcing factors.

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.

Surface Ocean CO2 Atlas (SOCAT) gridded data products

As a response to public demand for a well-documented, quality controlled, publically available, global surface ocean carbon dioxide (CO2) data set, the international marine carbon science community developed the Surface Ocean CO2 Atlas (SOCAT). The first SOCAT product is a collection of 6.3 million quality controlled surface CO2 data from the global oceans and coastal seas, spanning four decades (1968–2007). The SOCAT gridded data presented here is the second data product to come from the SOCAT project. Recognizing that some groups may have trouble working with millions of measurements, the SOCAT gridded product was generated to provide a robust, regularly spaced CO2 fugacity (fCO2) product with minimal spatial and temporal interpolation, which should be easier to work with for many applications. Gridded SOCAT is rich with information that has not been fully explored yet (e.g., regional differences in the seasonal cycles), but also contains biases and limitations that the user needs to recognize and address (e.g., local influences on values in some coastal regions).

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 %.

Coccolithophore species as indicators of surface oceanographic conditions in the vicinity of Azores islands

During summer 2008 and spring 2009, surface oceanographic surveys were carried out around three islands of the Azores archipelago (Terceira, Sao Miguel and Santa Maria) to assess the phytoplankton distribution and associated physico-chemical processes. The Azores archipelago is a major feature in the biogeochemical North Atlantic Subtropical Gyre (NAST) province although its influence on the productivity of the surrounding ocean is poorly known. Surface phytoplankton was studied by microscopy and HPLC (High Precision Liquid Chromatography). The mean values for biomass proxy Chlorophyll a (Chla) ranged from 0.04 to 0.55 mu g L-1 (Chla maximum = 0.86 mu g L-1) and coccolithophores were the most abundant group, followed by small flagellates, Cyanobacteria, diatoms and dinoflagellates being the least abundant group. The distribution of phytoplankton and coccolithophore species in particular presented seasonal differences and was consistent with the nearshore influence of warm subtropical waters from the south Azores current and colder subpolar waters from the north. The satellite-derived circulation patterns showed southward cold water intrusions off Terceira and northward warm water intrusions off Santa Maria. The warmer waters signal was confirmed by the subtropical coccolithophore assemblage, being Discosphaera tubifera a constant presence under these conditions. The regions of enhanced biomass, either resulting from northern cooler waters or from island induced processes, were characterized by the presence of Emiliania huxleyi. Diatoms and dinoflagellates indicated coastal and regional processes of nutrient enrichment and areas of physical stability, respectively.

Variation in ocean colour may help predict cod and haddock recruitment

Characteristics of the spring and fall phytoplankton blooms in spawning areas on the Scotian Shelf, Canada, were estimated from remote sensing data. These blooms, along with anomalies in the North Atlantic Oscillation, were used to explain variation in the recruitment of 4 populations of cod and haddock. We tested the effects of the timing of the bloom using the chlorophyll a (chl a) signal, the maximum amount of chl a, the timing of the diatom bloom, and the maximum relative dominance of diatoms on the recruitment (to Age 1) of cod and haddock on the Scotian Shelf. Models were run separately for the effects of the spring and fall blooms. Only 3 of 10 models tested (0-lag) explained significant (80 to 92%) variation in recruitment. However, the performance of these models was not consistent across populations or species, suggesting that generalities about how spring and fall phytoplankton blooms affect recruitment cannot yet be made. The differences among models suggest that fish larvae are probably adapted locally to food production and thus indirectly to the characteristics of the phytoplankton bloom, which in turn are influenced by regional (meso-scale) oceanographic conditions.

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.

Real-time sea surface air temperature retrieval over Southern Ocean using AMSR-E data

The AMSR-E satellite data and in-situ data were applied to retrieve sea surface air temperature (Ta) over the Southern Ocean. The in-situ data were obtained from the 24~(th) -26~(th) Chinese Antarctic Expeditions during 2008-2010. First, Ta was used to analyze the relativity with the bright temperature (Tb) from the twelve channels of AMSR-E, and no high relativity was found between Ta and Tb from any of the channels. The highest relativity was 0.38 (with 23.8 GHz). The dataset for the modeling was obtained by using in-situ data to match up with Tb, and two methods were applied to build the retrieval model. In multi-parameters regression method, the Tbs from 12 channels were used to the model and the region was divided into two parts according to the latitude of 50°S. The retrieval results were compared with the in-situ data. The Root Mean Square Error (RMS) and relativity of high latitude zone were 0.96℃and 0.93, respectively. And those of low latitude zone were 1.29 ℃ and 0.96, respectively. Artificial neural network (ANN) method was applied to retrieve Ta.The RMS and relativity were 1.26 ℃ and 0.98, respectively.

The response of phytoplankton to climate variability associated with the North Atlantic Oscillation

Inter-annual variability in the timing of phytoplankton spring bloom and phytoplankton community structure in the central North Atlantic Ocean was quantified using ocean color data and continuous plankton recorder (CPR) data. This variability was related to the North Atlantic Oscillation using correlation analysis and multivariate auto-regression models. The initiation of the spring bloom derived from CPR phytoplankton color index data is similar to that derived from satellite chlorophyll, and exhibits a nominal correlation with the sea surface temperature (SST) and the North Atlantic Oscillation (NAO). The extrapolated spring bloom timing suggested later initiation of blooms in the mid-1980s and earlier initiation of blooms in the 1990s. The climatological phytoplankton community structure in the central North Atlantic is dominated by diatoms, except for a shift in community composition favoring dinoflagellates in August. The ratio of diatoms to total phytoplankton abundance and the ratio of dinoflagellates to total phytoplankton abundance are both closely correlated with the NAO and SST. The extended time series of phytoplankton community structure between 1985 and 2009, deduced from the time series of SST and NAO over the same interval, showed a decadal shift away from diatoms towards dinoflagellates. The linkages between the NAO, and changes in stratification and phytoplankton processes occur over a larger scale than previously observed.

Joint distributions of waves and rain

The transfer of gases between the atmosphere and ocean is affected by a number of processes, of which wave action and rainfall are two of potential significance. Efforts have been made to quantify separately their contributions; however such assessments neglect the interaction of these phenomena. Here we look at the correlation statistics of waves and rain to note which regions display a strong association between rainfall and the local sea state. The conditional probability of rain varies from ~0.5% to ~15%, with most of the equatorial belt (which contains the ITCZ) showing a greater likelihood of rain at the lowest sea states. In contrast the occurrence of rain is independent of wave height in the Southern Ocean. The 1997/98 El Niño enhances the frequency of rain in some Pacific regions, with this change showing some association with wave conditions.