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

Worldwide alteration of lake mixing regimes in response to climate change

Lakes hold much of Earth’s accessible liquid freshwater, support biodiversity and provide key ecosystem services to people around the world. However, they are vulnerable to climate change, for example through shorter durations of ice cover, or through rising lake surface temperatures. Here we use a one-dimensional numerical lake model to assess climate change impacts on mixing regimes in 635 lakes worldwide. We run the lake model with input data from four state-of-the-art model projections of twenty-first-century climate under two emissions scenarios. Under the scenario with higher emissions (Representative Concentration Pathway 6.0), many lakes are projected to have reduced ice cover; about one-quarter of seasonally ice-covered lakes are projected to be permanently ice-free by 2080–2100. Surface waters are projected to warm, with a median warming across lakes of about 2.5 °C, and the most extreme warming about 5.5 °C. Our simulations suggest that around 100 of the stud- ied lakes are projected to undergo changes in their mixing regimes. About one-quarter of these 100 lakes are currently clas- sified as monomictic—undergoing one mixing event in most years— and will become permanently stratified systems. About one-sixth of these are currently dimictic—mixing twice per year—and will become monomictic. We conclude that many lakes will mix less frequently in response to climate change

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

Error correlations in High-Resolution Infrared Radiation Sounder (HIRS) Radiances

The High-resolution Infrared Radiation Sounder (HIRS) has been flown on 17 polar-orbiting satellites between the late 1970s and the present day. HIRS applications require accurate characterisation of uncertainties and inter-channel error correlations, which has so far been lacking. Here, we calculate error correlation matrices by accumulating count deviations for sequential sets of calibration measurements, and then correlating deviations between channels (for a fixed view) or views (for a fixed channel). The inter-channel error covariance is usually assumed to be diagonal, but we show that large error correlations, both positive and negative, exist between channels and between views close in time. We show that correlated error exists for all HIRS and that the degree of correlation varies markedly on both short and long timescales. Error correlations in excess of 0.5 are not unusual. Correlations between calibration observations taken sequentially in time arise from periodic error affecting both calibration and Earth counts. A Fourier spectral analysis shows that, for some HIRS instruments, this instrumental effect dominates at some or all spatial frequencies. These findings are significant for application of HIRS data in various applications, and related information will be made available as part of an upcoming Fundamental Climate Data Record covering all HIRS channels and satellites

Applying principles of metrology to historical Earth observations from satellites

Approaches from metrology can assist Earth Observation (EO) practitioners to develop quantitative characterisation of uncertainty in EO data. This is necessary for the credibility of statements based on Earth observations in relation to topics of public concern, particularly climate and environmental change. This paper presents the application of metrological uncertainty analysis to historical Earth observations from satellites, and is intended to aid mutual understanding of metrology and EO. The nature of satellite observations is summarised for different EO data processing levels, and key metrological nomenclature and principles for uncertainty characterisation are reviewed. We then address metrological approaches to developing estimates of uncertainty that are traceable from the satellite sensor, through levels of data processing, to products describing the evolution of the geophysical state of the Earth. EO radiances have errors with complex error correlation structures that are significant when performing common higher-level transformations of EO imagery. Principles of measurement-function-centred uncertainty analysis are described that apply sequentially to each EO data processing level. Practical tools for organising and traceably documenting uncertainty analysis are presented. We illustrate these principles and tools with examples including some specific sources of error seen in EO satellite data as well as with an example of the estimation of sea surface temperature from satellite infra-red imagery. This includes a simulation-based estimate for the error distribution of clear-sky infra-red brightness temperature (BT) in which calibration uncertainty and digitisation are found to dominate. The propagation of these errors to sea surface temperature is then presented, illustrating the relevance of the approach to derivation of EO-based climate datasets. We conclude with a discussion arguing that there is broad scope and need for improvement in EO practice as a measurement science. EO practitioners and metrologists willing to extend and adapt their disciplinary knowledge to meet this need can make valuable contributions to EO

Determining lake surface water temperatures worldwide using a tuned one-dimensional lake model (FLake, v1)

A tuning method for FLake, a one-dimensional (1-D) freshwater lake model, is applied for the individual tuning of 244 globally distributed large lakes using observed lake surface water temperatures (LSWTs) derived from along-track scanning radiometers (ATSRs). The model, which was tuned using only three lake properties (lake depth, snow and ice albedo and light extinction coefficient), substantially improves the measured mean differences in various features of the LSWT annual cycle, including the LSWTs of saline and high altitude lakes, when compared to the observed LSWTs. Lakes whose lake-mean LSWT persists below 1 °C for part of the annual cycle are considered to be seasonally ice-covered. For trial seasonally ice-covered lakes (21 lakes), the daily mean and standard deviation (2σ) of absolute differences between the modelled and observed LSWTs are reduced from 3.07 °C ± 2.25 °C to 0.84 °C ± 0.51 °C by tuning the model. For all other trial lakes (14 non-ice-covered lakes), the improvement is from 3.55 °C ± 3.20 °C to 0.96 °C ± 0.63 °C. The post tuning results for the 35 trial lakes (21 seasonally ice-covered lakes and 14 non-ice-covered lakes) are highly representative of the post-tuning results of the 244 lakes. For the 21 seasonally ice-covered lakes, the modelled response of the summer LSWTs to changes in snow and ice albedo is found to be statistically related to lake depth and latitude, which together explain 0.50 (R2adj, p = 0.001) of the inter-lake variance in summer LSWTs. Lake depth alone explains 0.35 (p = 0.003) of the variance. Lake characteristic information (snow and ice albedo and light extinction coefficient) is not available for many lakes. The approach taken to tune the model, bypasses the need to acquire detailed lake characteristic values. Furthermore, the tuned values for lake depth, snow and ice albedo and light extinction coefficient for the 244 lakes provide some guidance on improving FLake LSWT modelling

High performance software framework for the calculation of satellite-to-satellite data matchups (MMS version 1.2)

We present a Multisensor Matchup System (MMS) that allows systematic detection of satellite based sensor-to- sensor matchups and the extraction of local subsets of satellite data around matchup locations. The software system implements a generic matchup-detection approach and is currently being used for validation and sensor harmonisation purposes. An overview of the flexible and highly configurable software architecture and the target processing environments is given. We discuss improvements implemented with respect to heritage systems, and present some performance comparisons. A detailed description of the intersection algorithm is given which allows a fast matchup detection in geometry and time

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

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

Cloud clearing techniques over land for land surface temperature retrieval from the Advanced Along Track Scanning Radiometer

We present five new cloud detection algorithms over land based on dynamic threshold or Bayesian techniques, applicable to the Advanced Along Track Scanning Radiometer (AATSR) instrument and compare these with the standard threshold based SADIST cloud detection scheme. We use a manually classified dataset as a reference to assess algorithm performance and quantify the impact of each cloud detection scheme on land surface temperature (LST) retrieval. The use of probabilistic Bayesian cloud detection methods improves algorithm true skill scores by 8-9 % over SADIST (maximum score of 77.93 % compared to 69.27 %). We present an assessment of the impact of imperfect cloud masking, in relation to the reference cloud mask, on the retrieved AATSR LST imposing a 2 K tolerance over a 3x3 pixel domain. We find an increase of 5-7 % in the observations falling within this tolerance when using Bayesian methods (maximum of 92.02 % compared to 85.69 %). We also demonstrate that the use of dynamic thresholds in the tests employed by SADIST can significantly improve performance, applicable to cloud-test data to provided by the Sea and Land Surface Temperature Radiometer (SLSTR) due to be launched on the Sentinel 3 mission (estimated 2014)