Examination of Relationships between Clear-Sky Longwave Radiation and Aspects of the Atmospheric Hydrological Cycle in Climate Models, Reanalyses, and Observations

Relationships between clear-sky longwave radiation and aspects of the atmospheric hydrological cycle are quantified in models, reanalyses, and observations over the period 1980-2000. The robust sensitivity of clear-sky surface net longwave radiation (SNLc) to column-integrated water vapor (CWV) of 1-1.5 Wm(-2) mm(-1) combined with the positive relationship between CWV and surface temperature (T-s) explains substantial increases in clear-sky longwave radiative cooling of the atmosphere (Q(LWc)) to the surface over the period. Clear-sky outgoing longwave radiation (OLRc) is highly sensitive to changes in aerosol and greenhouse gas concentrations in addition to temperature and humidity. Over tropical ocean regions of mean descent, Q(LWc) increases with T-s at similar to 3.5-5.5 W m(-2) K-1 for reanalyses, estimates derived from satellite data, and models without volcanic forcing included. Increased Q(LWc) with warming across the tropical oceans helps to explain model ensemble mean increases in precipitation of 0.1-0.15 mm day(-1) K-1, which are primarily determined by ascent regions where precipitation increases at the rate expected from the Clausius-Clapeyron equation. The implications for future projections in the atmospheric hydrological cycle are discussed

Estimating surface CO2 fluxes from space-borne dry air mole fraction observations using an Ensemble Kalman filter

We have developed an ensemble Kalman Filter (EnKF) to estimate 8-day regional surface fluxes of CO2 from space-borne CO2 dry-air mole fraction observations (XCO2) and evaluate the approach using a series of synthetic experiments, in preparation for data from the NASA Orbiting Carbon Observatory (OCO). The 32-day duty cycle of OCO alternates every 16 days between nadir and glint measurements of backscattered solar radiation at short-wave infrared wavelengths. The EnKF uses an ensemble of states to represent the error covariances to estimate 8-day CO2 surface fluxes over 144 geographical regions. We use a 12×8-day lag window, recognising that XCO2 measurements include surface flux information from prior time windows. The observation operator that relates surface CO2 fluxes to atmospheric distributions of XCO2 includes: a) the GEOS-Chem transport model that relates surface fluxes to global 3-D distributions of CO2 concentrations, which are sampled at the time and location of OCO measurements that are cloud-free and have aerosol optical depths <0.3; and b) scene-dependent averaging kernels that relate the CO2 profiles to XCO2, accounting for differences between nadir and glint measurements, and the associated scene-dependent observation errors. We show that OCO XCO2 measurements significantly reduce the uncertainties of surface CO2 flux estimates. Glint measurements are generally better at constraining ocean CO2 flux estimates. Nadir XCO2 measurements over the terrestrial tropics are sparse throughout the year because of either clouds or smoke. Glint measurements provide the most effective constraint for estimating tropical terrestrial CO2 fluxes by accurately sampling fresh continental outflow over neighbouring oceans. We also present results from sensitivity experiments that investigate how flux estimates change with 1) bias and unbiased errors, 2) alternative duty cycles, 3) measurement density and correlations, 4) the spatial resolution of estimated flux estimates, and 5) reducing the length of the lag window and the size of the ensemble. At the revision stage of this manuscript, the OCO instrument failed to reach its orbit after it was launched on 24 February 2009. The EnKF formulation presented here is also applicable to GOSAT measurements of CO2 and CH4.

Simulation of the Earth's radiation budget by the European Centre for Medium-Range Weather Forecasts 40-year reanalysis (ERA40)

The radiation budget simulated by the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) is evaluated for the period 1979–2001 using independent satellite data and additional model data. This provides information on the quality of the radiation products and indirect evaluation of other aspects of the climate produced by ERA40. The climatology of clear-sky outgoing longwave radiation (OLR) is well captured by ERA40. Underestimations of about 10 W m−2 in clear-sky OLR over tropical convective regions by ERA40 compared to satellite data are substantially reduced when the satellite sampling is taken into account. The climatology of column-integrated water vapor is well simulated by ERA40 compared to satellite data over the ocean, indicating that the simulation of downward clear-sky longwave fluxes at the surface is likely to be good. Clear-sky absorbed solar radiation (ASR) and clear-sky OLR are overestimated by ERA40 over north Africa and high-latitude land regions. The observed interannual changes in low-latitude means are not well reproduced. Using ERA40 to analyze trends and climate feedbacks globally is therefore not recommended. The all-sky radiation budget is poorly simulated by ERA40. OLR is overestimated by around 10 W m−2 over much of the globe. ASR is underestimated by around 30 W m−2 over tropical ocean regions. Away from marine stratocumulus regions, where cloud fraction is underestimated by ERA40, the poor radiation simulation by ERA40 appears to be related to inaccurate radiative properties of cloud rather than inaccurate cloud distributions.

Can desert dust explain the outgoing longwave radiation anomaly over the Sahara during July 2003?

Measurements of the top‐of‐the‐atmosphere outgoing longwave radiation (OLR) for July 2003 from Meteosat‐7 are used to assess the performance of the numerical weather prediction version of the Met Office Unified Model. A significant difference is found over desert regions of northern Africa where the model emits too much OLR by up to 35 Wm−2 in the monthly mean. By cloud‐screening the data we find an error of up to 50 Wm−2 associated with cloud‐free areas, which suggests an error in the model surface temperature, surface emissivity, or atmospheric transmission. By building up a physical model of the radiative properties of mineral dust based on in situ, and surface‐based and satellite remote sensing observations we show that the most plausible explanation for the discrepancy in OLR is due to the neglect of mineral dust in the model. The calculations suggest that mineral dust can exert a longwave radiative forcing by as much as 50 Wm−2 in the monthly mean for 1200 UTC in cloud‐free regions, which accounts for the discrepancy between the model and the Meteosat‐7 observations. This suggests that inclusion of the radiative effects of mineral dust will lead to a significant improvement in the radiation balance of numerical weather prediction models with subsequent improvements in performance.

Exploitation of Geostationary Earth Radiation Budget data using simulations from a numerical weather prediction model: Methodology and data validation

We describe a new methodology for comparing satellite radiation budget data with a numerical weather prediction (NWP) model. This is applied to data from the Geostationary Earth Radiation Budget (GERB) instrument on Meteosat-8. The methodology brings together, in near-real time, GERB broadband shortwave and longwave fluxes with simulations based on analyses produced by the Met Office global NWP model. Results for the period May 2003 to February 2005 illustrate the progressive improvements in the data products as various initial problems were resolved. In most areas the comparisons reveal systematic errors in the model's representation of surface properties and clouds, which are discussed elsewhere. However, for clear-sky regions over the oceans the model simulations are believed to be sufficiently accurate to allow the quality of the GERB fluxes themselves to be assessed and any changes in time of the performance of the instrument to be identified. Using model and radiosonde profiles of temperature and humidity as input to a single-column version of the model's radiation code, we conduct sensitivity experiments which provide estimates of the expected model errors over the ocean of about ±5–10 W m−2 in clear-sky outgoing longwave radiation (OLR) and ±0.01 in clear-sky albedo. For the more recent data the differences between the observed and modeled OLR and albedo are well within these error estimates. The close agreement between the observed and modeled values, particularly for the most recent period, illustrates the value of the methodology. It also contributes to the validation of the GERB products and increases confidence in the quality of the data, prior to their release.

The Geostationary Earth Radiation Budget Project

This paper reports on a new satellite sensor, the Geostationary Earth Radiation Budget (GERB) experiment. GERB is designed to make the first measurements of the Earth's radiation budget from geostationary orbit. Measurements at high absolute accuracy of the reflected sunlight from the Earth, and the thermal radiation emitted by the Earth are made every 15 min, with a spatial resolution at the subsatellite point of 44.6 km (north–south) by 39.3 km (east–west). With knowledge of the incoming solar constant, this gives the primary forcing and response components of the top-of-atmosphere radiation. The first GERB instrument is an instrument of opportunity on Meteosat-8, a new spin-stabilized spacecraft platform also carrying the Spinning Enhanced Visible and Infrared (SEVIRI) sensor, which is currently positioned over the equator at 3.5°W. This overview of the project includes a description of the instrument design and its preflight and in-flight calibration. An evaluation of the instrument performance after its first year in orbit, including comparisons with data from the Clouds and the Earth's Radiant Energy System (CERES) satellite sensors and with output from numerical models, are also presented. After a brief summary of the data processing system and data products, some of the scientific studies that are being undertaken using these early data are described. This marks the beginning of a decade or more of observations from GERB, as subsequent models will fly on each of the four Meteosat Second Generation satellites.

Vertical current flow through extensive layer clouds

The global atmospheric electrical circuit sustains a vertical current density between the ionosphere and the Earth's surface, the existence of which is well-established from measurements made in fair-weather conditions. In overcast, but non-thunderstorm, non-precipitating conditions, the current travels through the cloud present, despite cloud layers having low electrical conductivity. For extensive layer clouds, this leads to space charge at the upper and lower cloud boundaries. Using a combination of atmospheric electricity and solar radiation measurements at three UK sites, vertical current measurements have been categorised into clear, broken, and overcast cloud conditions. This approach shows that the vertical “fair weather” current is maintained despite the presence of cloud. In fully overcast conditions with thick cloud, the vertical current is reduced compared to thin cloud overcast conditions, associated with the cloud's resistance contributions. Contribution of cloud to the columnar resistance depends both on cloud thickness, and the cloud's height.

The Upper Valanginian (Early Cretaceous) positive carbon-isotope event recorded in terrestrial plants

Our understanding of the ancient ocean-atmosphere system has focused on oceanic proxies. However, the study of terrestrial proxies is equally necessary to constrain our understanding of ancient climates and linkages between the terrestrial and oceanic carbon reservoirs. We have analyzed carbon-isotope ratios from fossil plant material through the Valanginian and Lower Hauterivian from a shallow-marine, ammonite-constrained succession in the Crimean Peninsula of the southern Ukraine in order to determine if the Upper Valanginian positive carbon-isotope excursion is expressed in the atmosphere. delta(13)C(plant) values fluctuate around -23% to -22% for the Valanginian-Hauterivian, except during the Upper Valanginian where delta(13)C(plant) values record a positive excursion to similar to-18%. Based upon ammonite biostratigraphy from Crimea, and in conjunction with a composite Tethyan marine delta(13)C(carb) curve, several conclusions can be drawn: (1) the delta(13)C(plant) record indicates that the atmospheric carbon reservoir was affected; (2) the defined ammonite correlations between Europe and Crimea are synchronous; and (3) a change in photosynthetic carbon-isotope fractionation, caused by a decrease in atmospheric PCO2, occurred during the Upper Valanginian Positive delta(13)C excursion. Our new data, combined with other paleoenvironmental and paleoclimatic information, indicate that the Upper Valanginian was a cool period (icehouse) and highlights that the Cretaceous period was interrupted by periods of cooling and was not an equable climate as previously thought. (C) 2005 Elsevier B.V. All rights reserved.

Terrestrial implications for the maritime geoarchaeological resource: A view from the Lower Palaeolithic

Stone tools and faunal remains have been recovered from the English Channel and the North Sea through trawling, dredging for aggregates, channel clearance, and coring. These finds highlight the potential for a maritime Lower Palaeolithic archaeological resource. It is proposed here that any Lower Palaeolithic artefacts, faunal remains, and sediments deposited in the maritime zone during dry, low-stand phases were once (and may still be) contextually similar to their counterparts in the terrestrial Lower Palaeolithic records of north-western Europe. Given these similarities, can interpretive models and analytical frameworks developed for terrestrial archaeology be profitably applied to an assessment of the potential value of any maritime resource? The terrestrial geoarchaeological resource for the Lower Palaeolithic is dominated by artefacts and ecofacts that have been fluvially reworked. The spatio-temporal resolution of these data varies from entire river valleys and marine isotope stages to river channel gravel bar surfaces and decadal timescales, thus supporting a variety of questions and approaches. However, the structure of the terrestrial resource also highlights two fundamental limitations in current maritime knowledge that can restrict the application of terrestrial approaches to any potential maritime resource: (i) how have the repetitive transgressions and regressions of the Middle and Late Pleistocene modified the terrace landforms and sediments associated with the river systems of the English Channel and southern North Sea basins?; and (ii) do the surviving submerged terrace landforms and fluvial sedimentary deposits support robust geochronological models, as is the case with the classical terrestrial terrace sequences? This paper highlights potential approaches to these questions, and concludes that the fluvial palaeogeography, Pleistocene fossils, and potential Lower Palaeolithic artefacts of the maritime geoarchaeological resource can be profitably investigated in future as derived, low-resolution data sets, facilitating questions of colonisation, occupation, demography, and material culture.

Laboratory simulation of terrestrial meteorite weathering using the Bensour (LL6) ordinary chondrite

Laboratory dissolution experiments using the LL6 ordinary chondrite Bensour demonstrate that meteoritic minerals readily react with distilled water at low temperatures, liberating ions into solution and forming reaction products. Three experiments were performed, all for 68 days and at atmospheric fO(2) but using a range of water/rock ratios and different ternperatures. Experiments I and 2 were batch experiments and undertaken at room temperature, whereas in experiment 3, condensed boiling water was dripped onto meteorite subsamples within a Soxhlet extractor. Solutions from experiment 1 were chemically analyzed at the end of the experiment, whereas aliquots were extracted from experiments 2 and 3 for analysis at regular intervals. In all three experiments, a very significant proportion of the Na, Cl, and K within the Bensour subsamples entered solution, demonstrating that chlorapatite and feldspar were especially susceptible to dissolution. Concentrations of Mg, Al, Si, Ca, and Fe in solution were strongly affected by the precipitation of reaction products and Mg and Ca may also have been removed by sorption. Calculations predict saturation of experimental solutions with respect to Al hydroxides, Fe oxides, and Fe (oxy)hydroxides, which would have frequently been accompanied by hydrous aluminosilicates. Some reaction products were identified and include silica, a Mg-rich silicate, Fe oxides, and Fe (oxy)hydroxides. The implications of these results are that even very short periods of subaerial exposure of ordinary chondrites will lead to dissolution of primary minerals and crystallization of weathering products that are likely to include aluminosilicates and silicates, Mg-Ca carbonates, and sulfates in addition to the ubiquitous Fe oxides and (oxy)hydroxides.

IntCal04 terrestrial radiocarbon age calibration, 0-26 cal kyr BP

A new calibration curve for the conversion of radiocarbon ages to calibrated (cal) ages has been constructed and internationally ratified to replace ImCal98, which extended from 0-24 cal kyr BP (Before Present, 0 cal BP = AD 1950). The new calibration data set for terrestrial samples extends from 0-26 cal kyr BP, but with much higher resolution beyond 11.4 cal kyr BP than ImCal98. Dendrochronologically-dated tree-ring samples cover the period from 0-12.4 cal kyr BP. Beyond the end of the tree rings, data from marine records (corals and foraminifera) are converted to the atmospheric equivalent with a site-specific marine reservoir correction to provide terrestrial calibration from 12.4-26.0 cal kyr BP. A substantial enhancement relative to ImCal98 is the introduction of a coherent statistical approach based on a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The tree-ring data sets, sources of uncertainty, and regional offsets are discussed here. The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine 04) are discussed in brief, but details are presented in Hughen et al. (this issue a). We do not make a recommendation for calibration beyond 26 cal kyr BP at this time; however, potential calibration data sets are compared in another paper (van der Plicht et al., this issue).

Estimating uncertainty in terrestrial critical loads and their exceedances at four sites in the UK

Critical loads are the basis for policies controlling emissions of acidic substances in Europe and elsewhere. They are assessed by several elaborate and ingenious models, each of which requires many parameters, and have to be applied on a spatially-distributed basis. Often the values of the input parameters are poorly known, calling into question the validity of the calculated critical loads. This paper attempts to quantify the uncertainty in the critical loads due to this "parameter uncertainty", using examples from the UK. Models used for calculating critical loads for deposition of acidity and nitrogen in forest and heathland ecosystems were tested at four contrasting sites. Uncertainty was assessed by Monte Carlo methods. Each input parameter or variable was assigned a value, range and distribution in an objective a fashion as possible. Each model was run 5000 times at each site using parameters sampled from these input distributions. Output distributions of various critical load parameters were calculated. The results were surprising. Confidence limits of the calculated critical loads were typically considerably narrower than those of most of the input parameters. This may be due to a "compensation of errors" mechanism. The range of possible critical load values at a given site is however rather wide, and the tails of the distributions are typically long. The deposition reductions required for a high level of confidence that the critical load is not exceeded are thus likely to be large. The implication for pollutant regulation is that requiring a high probability of non-exceedance is likely to carry high costs. The relative contribution of the input variables to critical load uncertainty varied from site to site: any input variable could be important, and thus it was not possible to identify variables as likely targets for research into narrowing uncertainties. Sites where a number of good measurements of input parameters were available had lower uncertainties, so use of in situ measurement could be a valuable way of reducing critical load uncertainty at particularly valuable or disputed sites. From a restricted number of samples, uncertainties in heathland critical loads appear comparable to those of coniferous forest, and nutrient nitrogen critical loads to those of acidity. It was important to include correlations between input variables in the Monte Carlo analysis, but choice of statistical distribution type was of lesser importance. Overall, the analysis provided objective support for the continued use of critical loads in policy development. (c) 2007 Elsevier B.V. All rights reserved.

Estimating surface net solar radiation by use of Landsat-5 TM and digital elevation models

The ground surface net solar radiation is the energy that drives physical and chemical processes at the ground surface. In this paper, multi-spectral data from the Landsat-5 TM, topographic data from a gridded digital elevation model, field measurements, and the atmosphere model LOWTRAN 7 are used to estimate surface net solar radiation over the FIFE site. Firstly an improved method is presented and used for calculating total surface incoming radiation. Then, surface albedo is integrated from surface reflectance factors derived from remotely sensed data from Landsat-5 TM. Finally, surface net solar radiation is calculated by subtracting surface upwelling radiation from the total surface incoming radiation.