Sea surface temperature for climate from the along-track scanning radiometers

This thesis describes the construction of a sea surface temperature (SST) dataset from Along-Track Scanning Radiometer (ATSR) observations suitable for climate applications. The algorithms presented here are now used at ESA for reprocessing of historical ATSR data and will be the basis of the retrieval used on the forthcoming SLSTR instrument on ESA’s Sentinel-3 satellite. In order to ensure independence of ATSR SSTs from in situ measurements, the retrieval uses physics-based methods through the use of radiative transfer (RT) simulations. The RT simulations are based on the Reference ForwardModel line-by-line model linked to a new sea surface emissivity model which accounts for surface temperature, wind speed, viewing angle and salinity, and to a discrete ordinates scattering (DISORT) model to account for aerosol. An atmospheric profile dataset, based on full resolution ERA-40 numerical weather prediction (NWP) data, is defined and used as input to the RTmodel. Five atmospheric trace gases (N2O, CH4, HNO3, and CFC-11 and CFC-12) are identified as having temporal and geographical variability which have a significant (∼0.1K) impact on RT simulations. Several additional trace gases neglected in previous studies are included using fixed profiles contributing ∼0.04K to RT simulations. Comparison against ATSR-2 and AATSR observations indicates that RT model biases are reduced from 0.2–0.5K for previous studies to ∼0.1K. A new coefficient-based SST retrieval scheme is developed from the RT simulations. Coefficients are banded by total column water vapour (TCWV) from NWP analyses reducing simulated regional biases to <0.1K compared to ∼0.2K for global coefficients. An improved treatment of the instrument viewing geometry decreases simulated view-angle related biases from ∼0.1K to <0.005K for the day-time dual-view retrieval. To eliminate inter-algorithmbiases due to remaining RT model biases and uncertainty in the characterisation of the ATSR instruments the offset coefficient for each TCWV band is adjusted to match the results from a reference channel combination. As infrared radiometers are sensitive to the skin SST while in situ buoys measure SST at some depth below the surface an adjustment for the skin effect and diurnal stratification is included. The samemodel allows adjustment for the differing time of observation between ATSR-2 and AATSR to prevent the diurnal cycle being aliased into the final record. The RT simulations are harmonised between sensors using a double-difference technique eliminating discontinuities in the final SST record. Comparison against in situ drifting and tropical moored buoys shows the new SST dataset is of high quality. Systematic differences between ATSR retrieved SST and in situ drifters show zonal, regional, TCWV, and wind speed biases are less than 0.1K except for themost extreme cases (TCWV <5 kgm−2). The precision of ATSR retrieved SSTs is ∼0.15 K, lower than the precision ofmeasurement of the global ensemble of in situ drifting buoys. From 1995 onwards the ARC SSTs are stable with instability of less than 5mK year−1 to 95% confidence (demonstrated for tropical regions)

Adjusting for desert-dust-related biases in a climate data record of sea surface temperature

Atmospheric desert-dust aerosol, primarily from north Africa, causes negative biases in remotely sensed climate data records of sea surface temperature (SST). Here, large-scale bias adjustments are deduced and applied to the v2 climate data record of SST from the European Space Agency Climate Change Initiative (CCI). Unlike SST from infrared sensors, SST measured in situ is not prone to desert-dust bias. An in-situ-based SST analysis is combined with column dust mass from the Modern-Era Retrospective analysis for Research and Applications, Version 2 to deduce a monthly, large-scale adjustment to CCI analysis SSTs. Having reduced the dust-related biases, a further correction for some periods of anomalous satellite calibration is also derived. The corrections will increase the usability of the v2 CCI SST record for oceanographic and climate applications, such as understanding the role of Arabian Sea SSTs in the Indian monsoon. The corrections will also pave the way for a v3 climate data record with improved error characteristics with respect to atmospheric dust aerosol

Datasets related to in-land water for limnology and remote sensing applications: distance-to-land, distance-to-water, water-body identifier and lake-centre co-ordinates

Datasets containing information to locate and identify water bodies have been generated from data locating static-water-bodies with resolution of about 300 m (1/360 deg) recently released by the Land Cover Climate Change Initiative (LC CCI) of the European Space Agency. The LC CCI water-bodies dataset has been obtained from multi-temporal metrics based on time series of the backscattered intensity recorded by ASAR on Envisat between 2005 and 2010. The new derived datasets provide coherently: distance to land, distance to water, water-body identifiers and lake-centre locations. The water-body identifier dataset locates the water bodies assigning the identifiers of the Global Lakes and Wetlands Database (GLWD), and lake centres are defined for in-land waters for which GLWD IDs were determined. The new datasets therefore link recent lake/reservoir/wetlands extent to the GLWD, together with a set of coordinates which locates unambiguously the water bodies in the database. Information on distance-to-land for each water cell and the distance-to-water for each land cell has many potential applications in remote sensing, where the applicability of geophysical retrieval algorithms may be affected by the presence of water or land within a satellite field of view (image pixel). During the generation and validation of the datasets some limitations of the GLWD database and of the LC CCI water-bodies mask have been found. Some examples of the inaccuracies/limitations are presented and discussed. Temporal change in water-body extent is common. Future versions of the LC CCI dataset are planned to represent temporal variation, and this will permit these derived datasets to be updated

Harmonization of space-borne infra-red sensors measuring sea surface temperature

Sea surface temperature (SST) is observed by a constellation of sensors, and SST retrievals are commonly combined into gridded SST analyses and climate data records (CDRs). Differential biases between SSTs from different sensors cause errors in such products, including feature artefacts. We introduce a new method for reducing differential biases across the SST constellation, by reconciling the brightness temperature (BT) calibration and SST retrieval parameters between sensors. We use the Advanced Along-Track Scanning Radiometer (AATSR) and the Sea and Land Surface Temperature Radiometer (SLSTR) as reference sensors, and the Advanced Very High Resolution Radiometer (AVHRR) of the MetOp-A mission to bridge the gap between these references. Observations across a range of AVHRR zenith angles are matched with dual-view three-channel skin SST retrievals from the AATSR and SLSTR. These skin SSTs act as the harmonization reference for AVHRR retrievals by optimal estimation (OE). Parameters for the harmonized AVHRR OE are iteratively determined, including BT bias corrections and observation error covariance matrices as functions of water-vapor path. The OE SSTs obtained from AVHRR are shown to be closely consistent with the reference sensor SSTs. Independent validation against drifting buoy SSTs shows that the AVHRR OE retrieval is stable across the reference-sensor gap. We discuss that this method is suitable to improve consistency across the whole constellation of SST sensors. The approach will help stabilize and reduce errors in future SST CDRs, as well as having application to other domains of remote sensing

Stability assessment of the (A)ATSR sea surface temperature climate dataset from the European Space Agency Climate Change Initiative

Sea surface temperature is a key component of the climate record, with multiple independent records giving confidence in observed changes. As part of the European Space Agencies (ESA) Climate Change Initiative (CCI) the satellite archives have been reprocessed with the aim of creating a new dataset that is independent of the in situ observations, and stable with no artificial drift (<0.1 K decade−1 globally) or step changes. We present a method to assess the satellite sea surface temperature (SST) record for step changes using the Penalized Maximal t Test (PMT) applied to aggregate time series. We demonstrated the application of the method using data from version EXP1.8 of the ESA SST CCI dataset averaged on a 7 km grid and in situ observations from moored buoys, drifting buoys and Argo floats. The CCI dataset was shown to be stable after ~1994, with minimal divergence (~0.01 K decade−1) between the CCI data and in situ observations. Two steps were identified due to the failure of a gyroscope on the ERS-2 satellite, and subsequent correction mechanisms applied. These had minimal impact on the stability due to having equal magnitudes but opposite signs. The statistical power and false alarm rate of the method were assessed

Independent uncertainty estimates for coefficient based sea surface temperature retrieval from the Along-Track Scanning Radiometer instruments

We establish a methodology for calculating uncertainties in sea surface temperature estimates from coefficient based satellite retrievals. The uncertainty estimates are derived independently of in-situ data. This enables validation of both the retrieved SSTs and their uncertainty estimate using in-situ data records. The total uncertainty budget is comprised of a number of components, arising from uncorrelated (eg. noise), locally systematic (eg. atmospheric), large scale systematic and sampling effects (for gridded products). The importance of distinguishing these components arises in propagating uncertainty across spatio-temporal scales. We apply the method to SST data retrieved from the Advanced Along Track Scanning Radiometer (AATSR) and validate the results for two different SST retrieval algorithms, both at a per pixel level and for gridded data. We find good agreement between our estimated uncertainties and validation data. This approach to calculating uncertainties in SST retrievals has a wider application to data from other instruments and retrieval of other geophysical variables

The role of Advanced Microwave Scanning Radiometer 2 channels within an optimal estimation scheme for sea surface temperature

We present an analysis of information content for sea surface temperature (SST) retrieval from the Advanced Microwave Scanning Radiometer 2 (AMSR2). We find that SST uncertainty of ∼0.37 K can be achieved within an optimal estimation framework in the presence of wind, water vapour and cloud liquid water effects, given appropriate assumptions for instrumental uncertainty and prior knowledge, and using all channels. We test all possible combinations of AMSR2 channels and demonstrate the importance of including cloud liquid water in the retrieval vector. The channel combinations, with the minimum number of channels, that carry most SST information content are calculated, since in practice calibration error drives a trade-off between retrieved SST uncertainty and the number of channels used. The most informative set of five channels is 6.9 V, 6.9 H, 7.3 V, 10.7 V and 36.5 H and these are suitable for optimal estimation retrievals. We discuss the relevance of microwave SSTs and issues related to them compared to SSTs derived from infra-red observations

Concepts for a geostationary-like polar missions

An evidence-led scientific case for development of a space-based polar remote sensing platform at geostationary-like (GEO-like) altitudes is developed through methods including a data user survey. Whilst a GEO platform provides a near static perspective, multiple platforms are required to provide circumferential coverage. Systems for achieving GEO-like polar observation likewise require multiple platforms however the perspective is non-stationery. A key choice is between designs that provide complete polar view from a single platform at any given instant, and designs where this is obtained by compositing partial views from multiple sensors. Users foresee an increased challenge in extracting geophysical information from composite images and consider the use of non-composited images advantageous. Users also find the placement of apogee over the pole to be preferable to the alternative scenarios. Thus, a clear majority of data users find the “Taranis” orbit concept to be better than a critical inclination orbit, due to the improved perspective offered. The geophysical products that would benefit from a GEO-like polar platform are mainly estimated from radiances in the visible/near infrared and thermal parts of the electromagnetic spectrum, which is consistent with currently proven technologies from GEO. Based on the survey results, needs analysis, and current technology proven from GEO, scientific and observation requirements are developed along with two instrument concepts with eight and four channels, based on Flexible Combined Imager heritage. It is found that an operational system could, mostly likely, be deployed from an Ariane 5 ES to a 16-hour orbit, while a proof-of-concept system could be deployed from a Soyuz launch to the same orbit

Bayesian cloud detection over land for climate data records

Cloud detection is a necessary step in the generation of land surface temperature (LST) climate data records (CDRs) and affects data quality and uncertainty. We present here a sensor- independent Bayesian cloud detection algorithm and show that it is suitable for use in the production of LST CDRs. We evaluate the performance of the cloud detection with reference to two man- ually masked datasets for the Advanced Along-Track Scanning Radiometer (AATSR) and find a 7.9% increase in the hit rate and 4.9% decrease in the false alarm rate when compared to the opera- tional cloud mask. We then apply the algorithm to four instruments aboard polar-orbiting satellites, which together can produce a global, 25-year LST CDR: the second Along-Track Scanning Radiometer (ATSR-2), AATSR, the Moderate Resolution Spectroradiometer (MODIS Terra) and the Sea and Land Surface Temperature Radiometer (SLSTR-A). The Bayesian cloud detection hit rate is assessed with respect to in situ ceilometer measurements for periods of overlap between sensors. The consistency of the hit rate is assessed between sensors, with mean differences in the cloud hit rate of 4.5% for ATSR-2 vs. AATSR, 4.9% for AATSR vs. MODIS, and 2.5% for MODIS vs. SLSTR-A. This is important because consistent cloud detection performance is needed for the observational stability of a CDR. The application of a sensor-independent cloud detection scheme in the production of CDRs is thus shown to be a viable approach to achieving LST observational stability over time

Sea surface temperature in global analyses: gains from the copernicus imaging microwave radiometer

Sea surface temperatures (SSTs) derived from passive microwave (PMW) observations benefit global ocean and SST analyses because of their near-all-weather availability. Present PMW SSTs have a real aperture-limited spatial resolution in excess of 50 km, limiting the spatial fidelity with which SST features, reflecting ocean dynamics, can be captured. This contrasts with the target resolution of global analyses of 5 to 10 km. The Copernicus Imaging Microwave Radiometer (CIMR) is a mission concept under consideration as a high-priority candidate mission for the expansion of the Copernicus space programme. This instrument would be capable of real aperture resolution < 15 km with low total uncertainties in the range 0.4–0.8 K for channels between 1.4 and 36.5 GHz, and a dual-view arrangement that further reduces noise. This paper provides a comparative study of SST uncertainty and feature resolution with and without the availability of CIMR in the future SST-observing satellite constellation based on a detailed simulation of CIMR plus infrared observations and the processing of global SST analyses with 0.05◦ final grid resolution. Simulations of CIMR data including structured errors were added to an observing system consisting of the Sea and Land Surface Temperature Radiometer (SLSTR) on Sentinel-3A and the Advanced Very High Resolution Radiometer (AVHRR) on MetOp-A. This resulted in a large improvement in the global root-mean-square error (RMSE) for SST from 0.37 K to 0.21 K for January and 0.40 K to 0.25 K for July. There was a particularly noticeable improvement in the performance of the analysis, as measured by the reduction in RMSE, for dynamical and persistently cloudy areas. Of these, the Aghulas Current showed an improvement of 43% in January and 48% in July, the Gulf Stream showed 70% and 44% improvements, the Southern Ocean showed 57% and 74% improvements, and the Maritime Continent showed 50% and 40% improvements, respectively