Extended optimal estimation techniques for sea surface temperature from the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI)

Sea surface temperature (SST) can be estimated from day and night observations of the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI) by optimal estimation (OE). We show that exploiting the 8.7 μm channel, in addition to the “traditional” wavelengths of 10.8 and 12.0 μm, improves OE SST retrieval statistics in validation. However, the main benefit is an improvement in the sensitivity of the SST estimate to variability in true SST. In a fair, single-pixel comparison, the 3-channel OE gives better results than the SST estimation technique presently operational within the Ocean and Sea Ice Satellite Application Facility. This operational technique is to use SST retrieval coefficients, followed by a bias-correction step informed by radiative transfer simulation. However, the operational technique has an additional “atmospheric correction smoothing”, which improves its noise performance, and hitherto had no analogue within the OE framework. Here, we propose an analogue to atmospheric correction smoothing, based on the expectation that atmospheric total column water vapour has a longer spatial correlation length scale than SST features. The approach extends the observations input to the OE to include the averaged brightness temperatures (BTs) of nearby clear-sky pixels, in addition to the BTs of the pixel for which SST is being retrieved. The retrieved quantities are then the single-pixel SST and the clear-sky total column water vapour averaged over the vicinity of the pixel. This reduces the noise in the retrieved SST significantly. The robust standard deviation of the new OE SST compared to matched drifting buoys becomes 0.39 K for all data. The smoothed OE gives SST sensitivity of 98% on average. This means that diurnal temperature variability and ocean frontal gradients are more faithfully estimated, and that the influence of the prior SST used is minimal (2%). This benefit is not available using traditional atmospheric correction smoothing.

Temperature trends at the Mauna Loa observatory, Hawaii

Observations at the Mauna Loa Observatory, Hawaii, established the systematic increase of anthropogenic CO2 in the atmosphere. For the same reasons that this site provides excellent globally averaged CO2 data, it may provide temperature data with global significance. Here, we examine hourly temperature records, averaged annually for 1977–2006, to determine linear trends as a function of time of day. For night-time data (22:00 to 06:00 LST (local standard time)) there is a near-uniform warming of 0.040 °C yr−1. During the day, the linear trend shows a slight cooling of −0.014 °C yr−1 at 12:00 LST (noon). Overall, at Mauna Loa Observatory, there is a mean warming trend of 0.021 °C yr−1. The dominance of night-time warming results in a relatively large annual decrease in the diurnal temperature range (DTR) of −0.050 °C yr−1 over the period 1977–2006. These trends are consistent with the observed increases in the concentrations of CO2 and its role as a greenhouse gas (demonstrated here by first-order radiative forcing calculations), and indicate the possible relevance of the Mauna Loa temperature measurements to global warming.

Understanding the rapid summer warming and changes in temperature extremes since the mid-1990s over Western Europe

Analysis of observations indicates that there was a rapid increase in summer (June-August, JJA) mean surface air temperature (SAT) since the mid-1990s over Western Europe. Accompanying this rapid warming are significant increases in summer mean daily maximum temperature, daily minimum temperature, annual hottest day temperature and warmest night temperature, and an increase in frequency of summer days and tropical nights, while the change in the diurnal temperature range (DTR) is small. This study focuses on understanding causes of the rapid summer warming and associated temperature extreme changes. A set of experiments using the atmospheric component of the state-of-the-art HadGEM3 global climate model have been carried out to quantify relative roles of changes in sea surface temperature (SST)/sea ice extent (SIE), anthropogenic greenhouse gases (GHGs), and anthropogenic aerosols (AAer). Results indicate that the model forced by changes in all forcings reproduces many of the observed changes since the mid-1990s over Western Europe. Changes in SST/SIE explain 62.2% ± 13.0% of the area averaged seasonal mean warming signal over Western Europe, with the remaining 37.8% ± 13.6% of the warming explained by the direct impact of changes in GHGs and AAer. Results further indicate that the direct impact of the reduction of AAer precursor emissions over Europe, mainly through aerosol-radiation interaction with additional contributions from aerosol-cloud interaction and coupled atmosphere-land surface feedbacks, is a key factor for increases in annual hottest day temperature and in frequency of summer days. It explains 45.5% ± 17.6% and 40.9% ± 18.4% of area averaged signals for these temperature extremes. The direct impact of the reduction of AAer precursor emissions over Europe acts to increase DTR locally, but the change in DTR is countered by the direct impact of GHGs forcing. In the next few decades, greenhouse gas concentrations will continue to rise and AAer precursor emissions over Europe and North America will continue to decline. Our results suggest that the changes in summer seasonal mean SAT and temperature extremes over Western Europe since the mid-1990s are most likely to be sustained or amplified in the near term, unless other factors intervene.

How many facets are needed to represent the surface energy balance of an urban area?

We investigate the question of how many facets are needed to represent the energy balance of an urban area by developing simplified 3-, 2- and 1-facet versions of a 4-facet energy balance model of two-dimensional streets and buildings. The 3-facet model simplifies the 4-facet model by averaging over the canyon orientation, which results in similar net shortwave and longwave balances for both wall facets, but maintains the asymmetry in the heat fluxes within the street canyon. For the 2-facet model, on the assumption that the wall and road temperatures are equal, the road and wall facets can be combined mathematically into a single street-canyon facet with effective values of the heat transfer coefficient, albedo, emissivity and thermodynamic properties, without further approximation. The 1-facet model requires the additional assumption that the roof temperature is also equal to the road and wall temperatures. Idealised simulations show that the geometry and material properties of the walls and road lead to a large heat capacity of the combined street canyon, whereas the roof behaves like a flat surface with low heat capacity. This means that the magnitude of the diurnal temperature variation of the street-canyon facets are broadly similar and much smaller than the diurnal temperature variation of the roof facets. Consequently, the approximation that the street-canyon facets have similar temperatures is sound, and the road and walls can be combined into a single facet. The roof behaves very differently and a separate roof facet is required. Consequently, the 2-facet model performs similarly to the 4-facet model, while the 1-facet model does not. The models are compared with previously published observations collected in Mexico City. Although the 3- and 2-facet models perform better than the 1-facet model, the present models are unable to represent the phase of the sensible heat flux. This result is consistent with previous model comparisons, and we argue that this feature of the data cannot be produced by a single column model. We conclude that a 2-facet model is necessary, and for numerical weather prediction sufficient, to model an urban surface, and that this conclusion is robust and therefore applicable to more general geometries.

Evidence of trends in daily climate extremes over southern and west Africa

[ 1] There has been a paucity of information on trends in daily climate and climate extremes, especially from developing countries. We report the results of the analysis of daily temperature ( maximum and minimum) and precipitation data from 14 south and west African countries over the period 1961 - 2000. Data were subject to quality control and processing into indices of climate extremes for release to the global community. Temperature extremes show patterns consistent with warming over most of the regions analyzed, with a large proportion of stations showing statistically significant trends for all temperature indices. Over 1961 to 2000, the regionally averaged occurrence of extreme cold ( fifth percentile) days and nights has decreased by - 3.7 and - 6.0 days/decade, respectively. Over the same period, the occurrence of extreme hot (95th percentile) days and nights has increased by 8.2 and 8.6 days/decade, respectively. The average duration of warm ( cold) has increased ( decreased) by 2.4 (0.5) days/decade and warm spells. Overall, it appears that the hot tails of the distributions of daily maximum temperature have changed more than the cold tails; for minimum temperatures, hot tails show greater changes in the NW of the region, while cold tails have changed more in the SE and east. The diurnal temperature range (DTR) does not exhibit a consistent trend across the region, with many neighboring stations showing opposite trends. However, the DTR shows consistent increases in a zone across Namibia, Botswana, Zambia, and Mozambique, coinciding with more rapid increases in maximum temperature than minimum temperature extremes. Most precipitation indices do not exhibit consistent or statistically significant trends across the region. Regionally averaged total precipitation has decreased but is not statistically significant. At the same time, there has been a statistically significant increase in regionally averaged daily rainfall intensity and dry spell duration. While the majority of stations also show increasing trends for these two indices, only a few of these are statistically significant. There are increasing trends in regionally averaged rainfall on extreme precipitation days and in maximum annual 5-day and 1-day rainfall, but only trends for the latter are statistically significant.

Causal relationships vs. emergent patterns in the global controls of fire frequency

Global controls on month-by-month fractional burnt area (2000–2005) were investigated by fitting a generalised linear model (GLM) to Global Fire Emissions Database (GFED) data, with 11 predictor variables representing vegetation, climate, land use and potential ignition sources. Burnt area is shown to increase with annual net primary production (NPP), number of dry days, maximum temperature, grazing-land area, grass/shrub cover and diurnal temperature range, and to decrease with soil moisture, cropland area and population density. Lightning showed an apparent (weak) negative influence, but this disappeared when pure seasonal-cycle effects were taken into account. The model predicts observed geographic and seasonal patterns, as well as the emergent relationships seen when burnt area is plotted against each variable separately. Unimodal relationships with mean annual temperature and precipitation, population density and gross domestic product (GDP) are reproduced too, and are thus shown to be secondary consequences of correlations between different controls (e.g. high NPP with high precipitation; low NPP with low population density and GDP). These findings have major implications for the design of global fire models, as several assumptions in current models – most notably, the widely assumed dependence of fire frequency on ignition rates – are evidently incorrect.

The diurnal course of soil moisture as measured by various dielectric sensors: Effects of soil temperature and the implications for evaporation estimates

Soil moisture content, theta, of a bare and vegetated UK gravelly sandy loam soil (in situ and repacked in small lysimeters) was measured using various dielectric instruments (single-sensor ThetaProbes, multi-sensor Profile Probes, and Aquaflex Sensors), at depths ranging between 0.03 and I m, during the summers of 2001 (in situ soil) and 2002 (mini-lysimeters). Half-hourly values of evaporation, E, were calculated from diurnal changes in total soil profile water content, using the soil water balance equation. For the bare soil field, Profile Probes and ML2x ThetaProbes indicated a diurnal course of theta that did not concur with typical soil physical observations: surface layer soil moisture content increased from early morning until about midday, after which theta declined, generally until the early evening. The unexpected course of theta was positively correlated to soil temperature, T-s, also at deeper depths. Aquaflex and ML1 ThetaProbe (older models) outputs, however, reflected common observations: 0 increased slightly during the night (capillary rise) and decreased from the morning until late afternoon (as a result of evaporation). For the vegetated plot, the spurious diurnal theta fluctuations were less obvious, because canopy shading resulted in lower amplitudes of T-s. The unrealistic theta profiles measured for the bare and vegetated field sites caused diurnal estimates of E to attain downward daytime and upward night-time values. In the mini-lysimeters, at medium to high moisture contents, theta values measured by (ML2x) ThetaProbes followed a relatively realistic course, and predictions of E from diurnal changes in vertically integrated theta generally compared well with lysimeter estimates of E. However, time courses of theta and E became comparable to those observed for the field plots when the soil in the lysimeters reached relatively low values of theta. Attempts to correct measured theta for fluctuations in T, revealed that no generally applicable formula could be derived. (c) 2005 Elsevier B.V. All rights reserved.

Diurnal cycle of land surface temperature in a desert encroachment zone as observed from satellites

Climate variability in the African Soudano-Sahel savanna zone has attracted much attention because of the persistence of anomalously low rainfall. Past efforts to monitor the climate of this region have focused on rainfall and vegetation conditions, while land surface temperature (LST) has received less attention. Remote sensing of LST is feasible and possible at global scale. Most remotely sensed estimates of LST are based on the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) that are limited in their ability to capture the full diurnal cycle. Although more frequent observations are available from past geostationary satellites, their spatial resolution is coarser than that of polar orbiting satellites. In this study, the improved capabilities of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the METEOSAT Second Generation (MSG) instrument are used to remotely sense the LST in the African Soudano-Sahel savanna zone at a resolution of 3 km and 15 minutes. In support of the Radiative Atmospheric Divergence using the ARM Mobile Facility (AMF), GERB and AMMA Stations (RADAGAST) project, African Monsoon Multidisciplinary Analyses (AMMA) project and the Department of Energy's Atmospheric Radiation Measurement (ARM) program, the ARM Mobile Facility was deployed during 2006 in this climatically sensitive region, thereby providing a unique opportunity to evaluate remotely sensed algorithms for deriving LST.

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

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.

Observations of the diurnal cycle of outgoing longwave radiation from the Geostationary Earth Radiation Budget instrument

The Geostationary Earth Radiation Budget instrument on Meteosat-8, located over Africa, provides unprecedented temporal sampling (~17 minutes) of the broadband emitted thermal and reflected solar radiances. We analyse the diurnal cycle of the outgoing longwave radiation (OLR) fluxes derived from the thermal radiances for July 2006. Principal component (PC) analysis separates the signals of the surface temperature response to solar heating and of the development of convective clouds. The first two PCs explain most of the OLR variations: PC1 (surface heating) explains 82.3% of the total variance and PC2 (cloud development) explains 12.8% of the variance. Convection is initiated preferentially over mountainous regions and the cloud then advects downstream in the ambient flow. Diurnal variations are much weaker over the oceans, but a coherent signal over the Gulf of Guinea suggests that the cloudiness is modulated by the diurnally varying contrast between the Gulf and the adjacent land mass.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

The diurnal cycle in the tropics

A global archive of high-resolution (3-hourly, 0.58 latitude–longitude grid) window (11–12 mm) brightness temperature (Tb) data from multiple satellites is being developed by the European Union Cloud Archive User Service (CLAUS) project. It has been used to construct a climatology of the diurnal cycle in convection, cloudiness, and surface temperature for all regions of the Tropics. An example of the application of the climatology to the evaluation of the climate version of the U.K. Met. Office Unified Model (UM), version HadAM3, is presented. The characteristics of the diurnal cycle described by the CLAUS data agree with previous observational studies, demonstrating the universality of the characteristics of the diurnal cycle for land versus ocean, clear sky versus convective regimes. It is shown that oceanic deep convection tends to reach its maximum in the early morning. Continental convection generally peaks in the evening, although there are interesting regional variations, indicative of the effects of complex land–sea and mountain–valley breezes, as well as the life cycle of mesoscale convective systems. A striking result from the analysis of the CLAUS data has been the extent to which the strong diurnal signal over land is spread out over the adjacent oceans, probably through gravity waves of varying depths. These coherent signals can be seen for several hundred kilometers and in some instances, such as over the Bay of Bengal, can lead to substantial diurnal variations in convection and precipitation. The example of the use of the CLAUS data in the evaluation of the Met. Office UM has demonstrated that the model has considerable difficulty in capturing the observed phase of the diurnal cycle in convection, which suggests some fundamental difficulties in the model’s physical parameterizations. Analysis of the diurnal cycle represents a powerful tool for identifying and correcting model deficiencies.