A reprocessing for climate of sea surface temperature from the Along-Track Scanning Radiometers: initial validation, accounting for skin and diurnal variability

An initial validation of the Along Track Scanning Radiometer (ATSR) Reprocessing for Climate (ARC) retrievals of sea surface temperature (SST) is presented. ATSR-2 and Advanced ATSR (AATSR) SST estimates are compared to drifting buoy and moored buoy observations over the period 1995 to 2008. The primary ATSR estimates are of skin SST, whereas buoys measure SST below the surface. Adjustment is therefore made for the skin effect, for diurnal stratification and for differences in buoy–satellite observation time. With such adjustments, satellite-in situ differences are consistent between day and night within ~ 0.01 K. Satellite-in situ differences are correlated with differences in observation time, because of the diurnal warming and cooling of the ocean. The data are used to verify the average behaviour of physical and empirical models of the warming/cooling rates. Systematic differences between adjusted AATSR and in-situ SSTs against latitude, total column water vapour (TCWV), and wind speed are less than 0.1 K, for all except the most extreme cases (TCWV < 5 kg m–2, TCWV > 60 kg m–2). For all types of retrieval except the nadir-only two-channel (N2), regional biases are less than 0.1 K for 80% of the ocean. Global comparison against drifting buoys shows night time dual-view two-channel (D2) SSTs are warm by 0.06 ± 0.23 K and dual-view three-channel (D3) SSTs are warm by 0.06 ± 0.21 K (day-time D2: 0.07 ± 0.23 K). Nadir-only results are N2: 0.03 ± 0.33 K and N3: 0.03 ± 0.19 K showing the improved inter-algorithm consistency to ~ 0.02 K. This represents a marked improvement from the existing operational retrieval algorithms for which inter-algorithm inconsistency is > 0.5 K. Comparison against tropical moored buoys, which are more accurate than drifting buoys, gives lower error estimates (N3: 0.02 ± 0.13 K, D2: 0.03 ± 0.18 K). Comparable results are obtained for ATSR-2, except that the ATSR-2 SSTs are around 0.1 K warm compared to AATSR

A reprocessing for climate of sea surface temperature from the Along-Track Scanning Radiometers: a new retrieval scheme

We present a new coefficient-based retrieval scheme for estimation of sea surface temperature (SST) from the Along Track Scanning Radiometer (ATSR) instruments. The new coefficients are banded by total column water vapour (TCWV), obtained from numerical weather prediction analyses. TCWV banding reduces simulated regional retrieval biases to < 0.1 K compared to biases ~ 0.2 K for global coefficients. Further, detailed treatment of the instrumental viewing geometry reduces simulated view-angle related biases from ~ 0.1 K down to < 0.005 K for dual-view retrievals using channels at 11 and 12 μm. A novel analysis of trade-offs related to the assumed noise level when defining coefficients is undertaken, and we conclude that adding a small nominal level of noise (0.01 K) is optimal for our purposes. When applied to ATSR observations, some inter-algorithm biases appear as TCWV-related differences in SSTs estimated from different channel combinations. The final step in coefficient determination is to adjust the offset coefficient in each TCWV band to match results from a reference algorithm. This reference uses the dual-view observations of 3.7 and 11 μm. The adjustment is independent of in situ measurements, preserving independence of the retrievals. The choice of reference is partly motivated by uncertainty in the calibration of the 12 μm of Advanced ATSR. Lastly, we model the sensitivities of the new retrievals to changes to TCWV and changes in true SST, confirming that dual-view SSTs are most appropriate for climatological applications

A reprocessing for climate of sea surface temperature from the Along-Track Scanning Radiometers: basis in radiative transfer

We present new radiative transfer simulations to support determination of sea surface temperature (SST) from Along Track Scanning Radiometer (ATSR) imagery. The simulations are to be used within the ATSR Reprocessing for Climate project. The simulations are based on the “Reference Forward Model” line-by-line model linked with a sea surface emissivity model that accounts for wind speed and temperature, and with a discrete ordinates scattering model (DISORT). Input to the forward model is a revised atmospheric profile dataset, based on full resolution ERA-40, with a wider range of high-latitude profiles to address known retrieval biases in those regions. Analysis of the radiative impacts of atmospheric trace gases shows that geographical and temporal variation of N2O, CH4, HNO3, and CFC-11 and CFC-12 have effects of order 0.05, 0.2, 0.1 K on the 3.7, 11, 12 μm channels respectively. In addition several trace gases, neglected in previous studies, are included using fixed profiles contributing ~ 0.04 K to top-of-atmosphere BTs. Comparison against observations for ATSR2 and AATSR indicates that forward model biases have been reduced from 0.2 to 0.5 K for previous simulations to ~ 0.1 K.

Estimation of sea surface temperature from the Spinning Enhanced Visible and Infra Red Imager, improved using numerical weather prediction

Most of the operational Sea Surface Temperature (SST) products derived from satellite infrared radiometry use multi-spectral algorithms. They show, in general, reasonable performances with root mean square (RMS) residuals around 0.5 K when validated against buoy measurements, but have limitations, particularly a component of the retrieval error that relates to such algorithms' limited ability to cope with the full variability of atmospheric absorption and emission. We propose to use forecast atmospheric profiles and a radiative transfer model to simulate the algorithmic errors of multi-spectral algorithms. In the practical case of SST derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation (MSG), we demonstrate that simulated algorithmic errors do explain a significant component of the actual errors observed for the non linear (NL) split window algorithm in operational use at the Centre de Météorologie Spatiale (CMS). The simulated errors, used as correction terms, reduce significantly the regional biases of the NL algorithm as well as the standard deviation of the differences with drifting buoy measurements. The availability of atmospheric profiles associated with observed satellite-buoy differences allows us to analyze the origins of the main algorithmic errors observed in the SEVIRI field of view: a negative bias in the inter-tropical zone, and a mid-latitude positive bias. We demonstrate how these errors are explained by the sensitivity of observed brightness temperatures to the vertical distribution of water vapour, propagated through the SST retrieval algorithm.

The impact of diurnal variability in sea surface temperature on the atlantic air-sea CO2 flux

The effect of diurnal variations in sea surface temperature (SST) on the air-sea flux of CO2 over the central Atlantic ocean and Mediterranean Sea (60 S–60 N, 60 W–45 E) is evaluated for 2005–2006. We use high spatial resolution hourly satellite ocean skin temperature data to determine the diurnal warming (ΔSST). The CO2 flux is then computed using three different temperature fields – a foundation temperature (Tf, measured at a depth where there is no diurnal variation), Tf, plus the hourly ΔSST and Tf, plus the monthly average of the ΔSSTs. This is done in conjunction with a physically-based parameterisation for the gas transfer velocity (NOAA-COARE). The differences between the fluxes evaluated for these three different temperature fields quantify the effects of both diurnal warming and diurnal covariations. We find that including diurnal warming increases the CO2 flux out of this region of the Atlantic for 2005–2006 from 9.6 Tg C a−1 to 30.4 Tg C a−1 (hourly ΔSST) and 31.2 Tg C a−1 (monthly average of ΔSST measurements). Diurnal warming in this region, therefore, has a large impact on the annual net CO2 flux but diurnal covariations are negligible. However, in this region of the Atlantic the uptake and outgassing of CO2 is approximately balanced over the annual cycle, so although we find diurnal warming has a very large effect here, the Atlantic as a whole is a very strong carbon sink (e.g. −920 Tg C a−1 Takahashi et al., 2002) making this is a small contribution to the Atlantic carbon budget.

Sea-surface temperature records of Termination 1 in the Gulf of California: Challenges for seasonal and interannual analogues of tropical Pacific climate change

Centennial-scale records of sea-surface temperature and opal composition spanning the Last Glacial Maximum and Termination 1 (circa 25–6 ka) are presented here from Guaymas Basin in the Gulf of California. Through the application of two organic geochemistry proxies, the U37K′ index and the TEX86H index, we present evidence for rapid, stepped changes in temperatures during deglaciation. These occur in both temperature proxies at 13 ka (∼3°C increase in 270 years), 10.0 ka (∼2°C decrease over ∼250 years) and at 8.2 ka (3°C increase in <200 years). An additional rapid warming step is also observed in TEX86H at 11.5 ka. In comparing the two temperature proxies and opal content, we consider the potential for upwelling intensity to be recorded and link this millennial-scale variability to shifting Intertropical Convergence Zone position and variations in the strength of the Subtropical High. The onset of the deglacial warming from 17 to 18 ka is comparable to a “southern hemisphere” signal, although the opal record mimics the ice-rafting events of the north Atlantic (Heinrich events). Neither the modern seasonal cycle nor El Niño/Southern Oscillation patterns provide valid analogues for the trends we observe in comparison with other regional records. Fully coupled climate model simulations confirm this result, and in combination we question whether the seasonal or interannual climate variations of the modern climate are valid analogues for the glacial and deglacial tropical Pacific.

Examining reliability of seasonal to decadal sea surface temperature forecasts: the role of ensemble dispersion

Useful probabilistic climate forecasts on decadal timescales should be reliable (i.e. forecast probabilities match the observed relative frequencies) but this is seldom examined. This paper assesses a necessary condition for reliability, that the ratio of ensemble spread to forecast error being close to one, for seasonal to decadal sea surface temperature retrospective forecasts from the Met Office Decadal Prediction System (DePreSys). Factors which may affect reliability are diagnosed by comparing this spread-error ratio for an initial condition ensemble and two perturbed physics ensembles for initialized and uninitialized predictions. At lead times less than 2 years, the initialized ensembles tend to be under-dispersed, and hence produce overconfident and hence unreliable forecasts. For longer lead times, all three ensembles are predominantly over-dispersed. Such over-dispersion is primarily related to excessive inter-annual variability in the climate model. These findings highlight the need to carefully evaluate simulated variability in seasonal and decadal prediction systems.Useful probabilistic climate forecasts on decadal timescales should be reliable (i.e. forecast probabilities match the observed relative frequencies) but this is seldom examined. This paper assesses a necessary condition for reliability, that the ratio of ensemble spread to forecast error being close to one, for seasonal to decadal sea surface temperature retrospective forecasts from the Met Office Decadal Prediction System (DePreSys). Factors which may affect reliability are diagnosed by comparing this spread-error ratio for an initial condition ensemble and two perturbed physics ensembles for initialized and uninitialized predictions. At lead times less than 2 years, the initialized ensembles tend to be under-dispersed, and hence produce overconfident and hence unreliable forecasts. For longer lead times, all three ensembles are predominantly over-dispersed. Such over-dispersion is primarily related to excessive inter-annual variability in the climate model. These findings highlight the need to carefully evaluate simulated variability in seasonal and decadal prediction systems.

Examining reliability of seasonal to decadal sea surface temperature forecasts: the role of ensemble dispersion

Useful probabilistic climate forecasts on decadal timescales should be reliable (i.e. forecast probabilities match the observed relative frequencies) but this is seldom examined. This paper assesses a necessary condition for reliability, that the ratio of ensemble spread to forecast error being close to one, for seasonal to decadal sea surface temperature retrospective forecasts from the Met Office Decadal Prediction System (DePreSys). Factors which may affect reliability are diagnosed by comparing this spread-error ratio for an initial condition ensemble and two perturbed physics ensembles for initialized and uninitialized predictions. At lead times less than 2 years, the initialized ensembles tend to be under-dispersed, and hence produce overconfident and hence unreliable forecasts. For longer lead times, all three ensembles are predominantly over-dispersed. Such over-dispersion is primarily related to excessive inter-annual variability in the climate model. These findings highlight the need to carefully evaluate simulated variability in seasonal and decadal prediction systems.Useful probabilistic climate forecasts on decadal timescales should be reliable (i.e. forecast probabilities match the observed relative frequencies) but this is seldom examined. This paper assesses a necessary condition for reliability, that the ratio of ensemble spread to forecast error being close to one, for seasonal to decadal sea surface temperature retrospective forecasts from the Met Office Decadal Prediction System (DePreSys). Factors which may affect reliability are diagnosed by comparing this spread-error ratio for an initial condition ensemble and two perturbed physics ensembles for initialized and uninitialized predictions. At lead times less than 2 years, the initialized ensembles tend to be under-dispersed, and hence produce overconfident and hence unreliable forecasts. For longer lead times, all three ensembles are predominantly over-dispersed. Such over-dispersion is primarily related to excessive inter-annual variability in the climate model. These findings highlight the need to carefully evaluate simulated variability in seasonal and decadal prediction systems.

Sea surface temperature datasets for climate applications from Phase 1 of the European Space Agency Climate Change Initiative (SST CCI)

Sea surface temperature (SST) datasets have been generated from satellite observations for the period 1991–2010, intended for use in climate science applications. Attributes of the datasets specifically relevant to climate applications are: first, independence from in situ observations; second, effort to ensure homogeneity and stability through the time-series; third, context-specific uncertainty estimates attached to each SST value; and, fourth, provision of estimates of both skin SST (the fundamental measure- ment, relevant to air-sea fluxes) and SST at standard depth and local time (partly model mediated, enabling comparison with his- torical in situ datasets). These attributes in part reflect requirements solicited from climate data users prior to and during the project. Datasets consisting of SSTs on satellite swaths are derived from the Along-Track Scanning Radiometers (ATSRs) and Advanced Very High Resolution Radiometers (AVHRRs). These are then used as sole SST inputs to a daily, spatially complete, analysis SST product, with a latitude-longitude resolution of 0.05°C and good discrimination of ocean surface thermal features. A product user guide is available, linking to reports describing the datasets’ algorithmic basis, validation results, format, uncer- tainty information and experimental use in trial climate applications. Future versions of the datasets will span at least 1982–2015, better addressing the need in many climate applications for stable records of global SST that are at least 30 years in length.

Spatial Patterns in Seasonal and Interannual Variability of Chlorophyll and Sea Surface Temperature in the California Current

Physical forcing and biological response within the California Current System (CCS) are highly variable over a wide range of scales. Satellite remote sensing offers the only feasible means of quantifying this variability over the full extent of the CCS. Using six years (1997-2003) of daily SST and chlorophyll imagery, we map the spatial dependence of dominant temporal variability at resolutions sufficient to identify recurrent mesoscale circulation and local pattern associated with coastal topography. Here we describe mean seasonal cycles and interannual variation; intraseasonal variability is left to a companion paper ( K. R. Legaard and A. C. Thomas, manuscript in preparation, 2006). Coastal upwelling dictates seasonality along north-central California, where weak cycles of SST fluctuate between spring minima and late summer maxima and chlorophyll peaks in early summer. Off northern California, chlorophyll maxima are bounded offshore by the seasonally recurrent upwelling jet. Seasonal cycles differ across higher latitudes and in the midlatitude Southern California Bight, where upwelling winds are less vigorous and/or persistent. Seasonality along south-central Baja is strongly affected by processes other than upwelling, despite year-round upwelling-favorable winds. Interannual variation is generally dominated by El Nino and La Nina conditions. Interannual SST variance is greatest along south-central Baja, although interannual variability constitutes a greater fraction of total variance inshore along southern Oregon and much of California. Patterns of interannual chlorophyll variance are consistent with dominant forcing through the widespread depression and elevation of the nutricline during El Nino and La Nina, respectively. Interannual variability constitutes a greater fraction of total chlorophyll variance offshore.

Spatial Patterns of Intraseasonal Variability of Chlorophyll and Sea Surface Temperature in the California Current

Six years of daily satellite data are used to quantify and map intraseasonal variability of chlorophyll and sea surface temperature (SST) in the California Current. We define intraseasonal variability as temporal variation remaining after removal of interannual variability and stationary seasonal cycles. Semivariograms are used to quantify the temporal structure of residual time series. Empirical orthogonal function (EOF) analyses of semivariograms calculated across the region isolate dominant scales and corresponding spatial patterns of intraseasonal variability. The mode 1 EOFs for both chlorophyll and SST semivariograms indicate a dominant timescale of similar to 60 days. Spatial amplitudes and patterns of intraseasonal variance derived from mode 1 suggest dominant forcing of intraseasonal variability through distortion of large scale chlorophyll and SST gradients by mesoscale circulation. Intraseasonal SST variance is greatest off southern Baja and along southern Oregon and northern California. Chlorophyll variance is greatest over the shelf and slope, with elevated values closely confined to the Baja shelf and extending farthest from shore off California and the Pacific Northwest. Intraseasonal contributions to total SST variability are strongest near upwelling centers off southern Oregon and northern California, where seasonal contributions are weak. Intraseasonal variability accounts for the majority of total chlorophyll variance in most inshore areas save for southern Baja, where seasonal cycles dominate. Contributions of higher EOF modes to semivariogram structure indicate the degree to which intraseasonal variability is shifted to shorter timescales in certain areas. Comparisons of satellite-derived SST semivariograms to those calculated from co-located and concurrent buoy SST time series show similar features.

Sea surface temperature and sea ice variability in the subpolar North Atlantic from explosive volcanism of the late thirteenth century

In this study, we use IP and alkenone biomarker proxies to document the subdecadal variations of sea ice and sea surface temperature in the subpolar North Atlantic induced by the decadally paced explosive tropical volcanic eruptions of the second half of the thirteenth century. The short- and long-term evolutions of both variables were investigated by cross analysis with a simulation of the IPSL-CM5A LR model. Our results show short-term ocean cooling and sea ice expansion in response to each volcanic eruption. They also highlight that the long response time of the ocean leads to cumulative surface cooling and subsurface heat buildup due to sea ice capping. As volcanic forcing relaxes, the surface ocean rapidly warms, likely amplified by subsurface heat, and remains almost ice free for several decades

A four-year record of UK'37- and TEX86-derived sea surface temperature estimates from sinking particles in the filamentous upwelling region off Cape Blanc, Mauritania

Lipid biomarker records from sinking particles collected by sediment traps are excellent tools to study the seasonality of biomarker production as well as processes of particle formation and settling, ultimately leading to the preservation of the biomarkers in sediments. Here we present records of the biomarker indices UK'37 based on alkenones and TEX86 based on isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs), both used for the reconstruction of sea surface temperatures (SST). These records were obtained from sinking particles collected using a sediment trap moored in the filamentous upwelling zone off Cape Blanc, Mauritania, at approximately 1300 water depth during a four-year time interval between 2003 and 2007. Mass and lipid fluxes are highest during peak upwelling periods between October and June. The alkenone and GDGT records both display pronounced seasonal variability. Sinking velocities calculated from the time lag between measured SST maxima and minima and corresponding index maxima and minima in the trap samples are higher for particles containing alkenones (14-59 m/d) than for GDGTs (9-17 m/d). It is suggested that GDGTs are predominantly exported from shallow waters by incorporation in opal-rich particles. SST estimates based on the UK'37 index faithfully record observed fluctuations in SST during the study period. Temperature estimates based on TEX86 show smaller seasonal amplitudes, which can be explained with either predominant production of GDGTs during the warm season, or a contribution of GDGTs exported from deep waters carrying GDGTs in a distribution that translates to a high TEX86 signal.

(Figure 2) Organic carbon content, stable isotope ratios, TEX86 index and sea surface temperature reconstruction for the early Aptian of ODP Hole 198-1207B

The Cretaceous has long been recognized as a time when greenhouse conditions were fueled by elevated atmospheric CO2 and accompanied by perturbations of the global carbon cycle described as oceanic anoxic events (OAEs). Yet, the magnitude and frequency of temperature change during this interval of warm and equable climate are poorly constrained. Here we present a high-resolution record of sea-surface temperatures (SSTs) reconstructed using the TEX86 paleothermometer for a sequence of early Aptian organic-rich sediments deposited during the first Cretaceous OAE (OAE1a) at Shatsky Rise in the tropical Pacific. SSTs range from ~30 to ~36 °C and include two prominent cooling episodes of ~4 °C. The cooler temperatures reflect significant temperature instability in the tropics likely triggered by changes in carbon cycling induced by enhanced burial of organic matter. SST instability recorded during the early Aptian in the Pacific is comparable to that reported for the late Albian-early Cenomanian in the Atlantic, suggesting that such climate perturbations may have recurred during the Cretaceous with concomitant consequences for biota and the marine environment.

Magnesium/Calcium ratios and sea surface temperature estimation for sediments of the Timor Sea

We use a multiproxy approach to monitor changes in the vertical profile of the Indonesian Throughflow as well as monsoonal wind and precipitation patterns in the Timor Sea on glacial-interglacial, precessional, and suborbital timescales. We focus on an interval of extreme climate change and sea level variation: marine isotope (MIS) 6 to MIS 5e. Paleoproductivity fluctuations in the Timor Sea follow a precessional beat related to the intensity of the Australian (NW) monsoon. Paired Mg/Ca and d18O measurements of surface- and thermocline-dwelling planktonic foraminifers (G. ruber and P. obliquiloculata) indicate an increase of >4°C in both surface and thermocline water temperatures during Termination II. Tropical sea surface temperature changed synchronously with ice volume (benthic d18O) during deglaciation, implying a direct coupling of high- and low-latitude climate via atmospheric and/or upper ocean circulation. Substantial cooling and freshening of thermocline waters occurred toward the end of Termination II and during MIS 5e, indicating a change in the vertical profile of the Indonesian Throughflow from surface- to thermocline-dominated flow.