Quantifying the effects of perturbing the physics of an interactive sulfur scheme using an ensemble of GCMs on the climateprediction.net platform

Aerosols from anthropogenic and natural sources have been recognized as having an important impact on the climate system. However, the small size of aerosol particles (ranging from 0.01 to more than 10 μm in diameter) and their influence on solar and terrestrial radiation makes them difficult to represent within the coarse resolution of general circulation models (GCMs) such that small-scale processes, for example, sulfate formation and conversion, need parameterizing. It is the parameterization of emissions, conversion, and deposition and the radiative effects of aerosol particles that causes uncertainty in their representation within GCMs. The aim of this study was to perturb aspects of a sulfur cycle scheme used within a GCM to represent the climatological impacts of sulfate aerosol derived from natural and anthropogenic sulfur sources. It was found that perturbing volcanic SO2 emissions and the scavenging rate of SO2 by precipitation had the largest influence on the sulfate burden. When these parameters were perturbed the sulfate burden ranged from 0.73 to 1.17 TgS for 2050 sulfur emissions (A2 Special Report on Emissions Scenarios (SRES)), comparable with the range in sulfate burden across all the Intergovernmental Panel on Climate Change SRESs. Thus, the results here suggest that the range in sulfate burden due to model uncertainty is comparable with scenario uncertainty. Despite the large range in sulfate burden there was little influence on the climate sensitivity, which had a range of less than 0.5 K across the ensemble. We hypothesize that this small effect was partly associated with high sulfate loadings in the control phase of the experiment.

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.

Jupiter's and Saturn's convectively driven banded jets in the laboratory

The banded patterns of cloud and wind are among the most striking features of the atmospheres of Jupiter and Saturn, but their dynamical origin remains poorly understood. Most approaches towards understanding zonation so far (also in the terrestrial oceans) have used highly idealized models to show that it might originate from dynamical anisotropy in a shallow turbulent fluid layer due to the planetary β-effect. Here we report the results of laboratory experiments, conducted on a 14-m diameter turntable, which quantitatively confirm that multiple zonal jets may indeed be generated and maintained by this mechanism in the presence of deep convection and a topographic β-effect. At the very small values of Ekman number (≤2 × 10−5) and large local Reynolds numbers (≥2000, based on jet scales) achieved, the kinetic energy spectra suggest the presence of both energy-cascading and enstrophy-cascading inertial ranges in addition to the zonation near twice the Rhines wave number.

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.

Comparison of balloon-carried atmospheric motion sensors with Doppler lidar turbulence measurements

Magnetic sensors have been added to a standard weather balloon radiosonde package to detect motion in turbulent air. These measure the terrestrial magnetic field and return data over the standard uhf radio telemetry. Variability in the magnetic sensor data is caused by motion of the instrument package. A series of radiosonde ascents carrying these sensors has been made near a Doppler lidar measuring atmospheric properties. Lidar-retrieved quantities include vertical velocity (w) profile and its standard deviation (w). w determined over 1 h is compared with the radiosonde motion variability at the same heights. Vertical motion in the radiosonde is found to be robustly increased when w>0.75 m s−1 and is linearly proportional to w. ©2009 American Institute of Physics

Evidence for global circuit current flow through water droplet layers

Abstract Foggy air and clear air have appreciably different electrical conductivities. The conductivity gradient at horizontal droplet boundaries causes droplet charging, as a result of vertical current flow in the global atmospheric electrical circuit. The charging is poorly known, as both the current flow through atmospheric water droplet layers and the air conductivity are poorly characterised experimentally. Surface measurements during three days of continuous fog using new instrument techniques show that a shallow (of order 100 m deep) fog layer still permits the vertical conduction current to pass. Further, the conductivity in the fog is estimated to be approximately 20% lower than in clear air. Assuming a fog transition thickness of one metre, this implies a vertical conductivity gradient of order 10 fS m−2 at the boundary. The actual vertical conductivity gradient at a cloud boundary would probably be greater, due to the presence of larger droplets in clouds compared to fog, and cleaner, more conductive clear air aloft.

An Integrated Assessment of changes in the thermohaline circulation

This paper discusses the risks of a shutdown of the thermohaline circulation (THC) for the climate system, for ecosystems in and around the North Atlantic as well as for fisheries and agriculture by way of an Integrated Assessment. The climate model simulations are based on greenhouse gas scenarios for the 21st century and beyond. A shutdown of the THC, complete by 2150, is triggered if increased freshwater input from inland ice melt or enhanced runoff is assumed. The shutdown retards the greenhouse gas-induced atmospheric warming trend in the Northern Hemisphere, but does not lead to a persistent net cooling. Due to the simulated THC shutdown the sea level at the North Atlantic shores rises by up to 80 cm by 2150, in addition to the global sea level rise. This could potentially be a serious impact that requires expensive coastal protection measures. A reduction of marine net primary productivity is associated with the impacts of warming rather than a THC shutdown. Regional shifts in the currents in the Nordic Seas could strongly deteriorate survival chances for cod larvae and juveniles. This could lead to cod fisheries becoming unprofitable by the end of the 21st century. While regional socioeconomic impacts might be large, damages would be probably small in relation to the respective gross national products. Terrestrial ecosystem productivity is affected much more by the fertilization from the increasing CO2 concentration than by a THC shutdown. In addition, the level of warming in the 22nd to 24th century favours crop production in northern Europe a lot, no matter whether the THC shuts down or not. CO2 emissions corridors aimed at limiting the risk of a THC breakdown to 10% or less are narrow, requiring departure from business-as-usual in the next few decades. The uncertainty about THC risks is still high. This is seen in model analyses as well as in the experts’ views that were elicited. The overview of results presented here is the outcome of the Integrated Assessment project INTEGRATION.

Modeling the Stream Water Nitrate Dynamics in 60,000-km(2) European Catchment, the Garonne, Southwest France

The spatial and temporal dynamics in the stream water NO3-N concentrations in a major European river-system, the Garonne (62,700 km(2)), are described and related to variations in climate, land management, and effluent point-sources using multivariate statistics. Building on this, the Hydrologiska Byrans Vattenbalansavdelning (HBV) rainfall-runoff model and the Integrated Catchment Model of Nitrogen (INCA-N) are applied to simulate the observed flow and N dynamics. This is done to help us to understand which factors and processes control the flow and N dynamics in different climate zones and to assess the relative inputs from diffuse and point sources across the catchment. This is the first application of the linked HBV and INCA-N models to a major European river system commensurate with the largest basins to be managed tinder the Water Framework Directive. The simulations suggest that in the lowlands, seasonal patterns in the stream water NO3-N concentrations emerge and are dominated by diffuse agricultural inputs, with an estimated 75% of the river load in the lowlands derived from arable farming. The results confirm earlier European catchment studies. Namely, current semi-distrubuted catchment-scale dynamic models, which integrate variations in land cover, climate, and a simple representation of the terrestrial and in-stream N cycle, are able to simulate seasonal NO3-N patterns at large spatial (> 300 km(2)) and temporal (>= monthly) scales using available national datasets.

Over-parameterised, uncertain 'mathematical marionettes' - How can we best use catchment water quality models? An example of an 80-year catchment-scale nutrient balance

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the River Lambourn, a Chalk river-system in southern England. The model's abilities to simulate the long-term trend and seasonal patterns in observed stream water nitrate concentrations from 1920 to 2003 were tested. This is the first time a semi-distributed, daily time-step model has been applied to simulate such a long time period and then used to calculate detailed catchment nutrient budgets which span the conversion of pasture to arable during the late 1930s and 1940s. Thus, this work goes beyond source apportionment and looks to demonstrate how such simulations can be used to assess the state of the catchment and develop an understanding of system behaviour. The mass-balance results from 1921, 1922, 1991, 2001 and 2002 are presented and those for 1991 are compared to other modelled and literature values of loads associated with nitrogen soil processes and export. The variations highlighted the problem of comparing modelled fluxes with point measurements but proved useful for identifying the most poorly understood inputs and processes thereby providing an assessment of input data and model structural uncertainty. The modelled terrestrial and instream mass-balances also highlight the importance of the hydrological conditions in pollutant transport. Between 1922 and 2002, increased inputs of nitrogen from fertiliser, livestock and deposition have altered the nitrogen balance with a shift from possible reduction in soil fertility but little environmental impact in 1922, to a situation of nitrogen accumulation in the soil, groundwater and instream biota in 2002. In 1922 and 2002 it was estimated that approximately 2 and 18 kg N ha(-1) yr(-1) respectively were exported from the land to the stream. The utility of the approach and further considerations for the best use of models are discussed. (C) 2008 Elsevier B.V. All rights reserved.

Cropping systems, carbon sequestration and erosion in Brazil, a review

Soils represent a large carbon pool, approximately 1500 Gt, which is equivalent to almost three times the quantity stored in terrestrial biomass and twice the amount stored in the atmosphere. Any modification of land use or land management can induce variations in soil carbon stocks, even in agricultural systems that are perceived to be in a steady state. Tillage practices often induce soil aerobic conditions that are favourable to microbial activity and may lead to a degradation of soil structure. As a result, mineralisation of soil organic matter increases in the long term. The adoption of no-tillage systems and the maintenance of a permanent vegetation cover using Direct seeding Mulch-based Cropping system or DMC, may increase carbon levels in the topsoil. In Brazil, no-tillage practices (mainly DMC), were introduced approximately 30 years ago in the south in the Parana state, primarily as a means of reducing erosion. Subsequently, research has begun to study the management of the crop waste products and their effects on soil fertility, either in terms of phosphorus management, as a means of controlling soil acidity, or determining how manures can be applied in a more localised manner. The spread of no-till in Brazil has involved a large amount of extension work. The area under no-tillage is still increasing in the centre and north of the country and currently occupies ca. 20 million hectares, covering a diversity of environmental conditions, cropping systems and management practices. Most studies of Brazilian soils give rates of carbon storage in the top 40 cm of the soil of 0.4 to 1.7 t C ha(-1) per year, with the highest rates in the Cerrado region. However, caution must be taken when analysing DMC systems in terms of carbon sequestration. Comparisons should include changes in trace gas fluxes and should not be limited to a consideration of carbon storage in the soil alone if the full implications for global warming are to be assessed.

Pyrohydrolysis-IRMS determination of silicate chlorine stable isotope compositions. Application to oceanic crust and meteorite samples

This contribution describes the optimization of chlorine extraction from silicate samples by pyrohydrolysis prior to the precise determination of Cl stable-isotope compositions (637 Cl) by gas source, dual inlet Isotope Ratio Mass Spectrometry (IRMS) on CH(3)Clg. The complete method was checked on three international reference materials for Cl-content and two laboratory glass standards. Whole procedure blanks are lower than 0. 5 mu mol, corresponding to less than 10 wt.% of most of the sample chloride analysed. In the absence of international chlorine isotope rock, we report here Cl extracted compared to accepted Cl contents and reproducibilities on Cl and delta Cl-37 measurements for the standard rocks. After extraction, the Cl contents of the three international references compared within error with the accepted values (mean yield = 94 +/-10%) with reproducibilities better than 12% (10). The laboratory glass standards - andesite SO100DS92 and phonolite S9(2) - were used specifically to test the effect of chloride amount on the measurements. They gave Cl extraction yields of 100 +/-6% (1 sigma-; n = 15) and 105 +/- 8% (1 sigma-; n = 7), respectively, with delta Cl-37 values of -0.51 0.14%o and -0.39 0.17%o (1g). In summary, for silicate samples with Cl contents between 39 and 9042 ppm, the Pyrohydrolysis/HPLC method leads to overall CI extraction yields of 100 8%, reproducibilities on Cl contents of 7% and on delta Cl-37 measurements of 0.12%o (all 1 sigma). The method was further applied to ten silicate rocks of various mineralogy and chemistry (meteorite, fresh MORB glasses, altered basalts and setpentinized peridotites) chosen for their large range of Cl contents (70-2156 ppm) and their geological significance. delta Cl-37 values range between -2.33 and -0.50%o. These strictly negative values contrast with the large range and mainly positive values previously reported for comparable silicate samples and shown here to be affected by analytical problems. Thus we propose a preliminary, revised terrestrial CI cycle, mainly dominated by negative and zero delta Cl-37 values. (C) 2007 Elsevier B.V. All rights reserved.

Magnitude of climate variability during middle Pliocene warmth: a palaeoclimate modelling study

The nature and magnitude of climatic variability during the period of middle Pliocene warmth (ca 3.29–2.97 Ma) is poorly understood. We present a suite of palaeoclimate modelling experiments incorporating an advanced atmospheric general circulation model (GCM), coupled to a Q-flux ocean model for 3.29, 3.12 and 2.97 Ma BP. Astronomical solutions for the periods in question were derived from the Berger and Loutre BL2 astronomical solution. Boundary conditions, excluding sea surface temperatures (SSTs) which were predicted by the slab-ocean model, were provided from the USGS PRISM2 2°×2° digital data set. The model results indicate that little annual variation (0.5°C) in SSTs, relative to a ‘control’ experiment, occurred during the middle Pliocene in response to the altered orbital configurations. Annual surface air temperatures also displayed little variation. Seasonally, surface air temperatures displayed a trend of cooler temperatures during December, January and February, and warmer temperatures during June, July and August. This pattern is consistent with altered seasonality resulting from the prescribed orbital configurations. Precipitation changes follow the seasonal trend observed for surface air temperature. Compared to present-day, surface wind strength and wind stress over the North Atlantic, North Pacific and Southern Ocean remained greater in each of the Pliocene experiments. This suggests that wind-driven gyral circulation may have been consistently greater during the middle Pliocene. The trend of climatic variability predicted by the GCM for the middle Pliocene accords with geological data. However, it is unclear if the model correctly simulates the magnitude of the variation. This uncertainty is derived from, (a) the relative insensitivity of the GCM to perturbation in the imposed boundary conditions, (b) a lack of detailed time series data concerning changes to terrestrial ice cover and greenhouse gas concentrations for the middle Pliocene and (c) difficulties in representing the effects of ‘climatic history’ in snap-shot GCM experiments.

Marine04 marine radiocarbon age calibration, 0-26 cal kyr BP

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0-26 cal kyr BP (Before Present, 0 cal. BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0-10.5 call kyr BR Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific C-14 reservoir age information to provide a single global marine mixed-layer calibration from 10.5-26.0 cal kyr BR A substantial enhancement relative to IntCal98 is the introduction of 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 marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

National and sub-national assessments of soil organic carbon stocks and changes: The GEFSOC modelling system

Soil organic carbon (SOC) plays a vital role in ecosystem function, determining soil fertility, water holding capacity and susceptibility to land degradation. In addition, SOC is related to atmospheric CO, levels with soils having the potential for C release or sequestration, depending on land use, land management and climate. The United Nations Convention on Climate Change and its Kyoto Protocol, and other United Nations Conventions to Combat Desertification and on Biodiversity all recognize the importance of SOC and point to the need for quantification of SOC stocks and changes. An understanding of SOC stocks and changes at the national and regional scale is necessary to further our understanding of the global C cycle, to assess the responses of terrestrial ecosystems to climate change and to aid policy makers in making land use/management decisions. Several studies have considered SOC stocks at the plot scale, but these are site specific and of limited value in making inferences about larger areas. Some studies have used empirical methods to estimate SOC stocks and changes at the regional scale, but such studies are limited in their ability to project future changes, and most have been carried out using temperate data sets. The computational method outlined by the Intergovernmental Panel on Climate Change (IPCC) has been used to estimate SOC stock changes at the regional scale in several studies, including a recent study considering five contrasting eco regions. This 'one step' approach fails to account for the dynamic manner in which SOC changes are likely to occur following changes in land use and land management. A dynamic modelling approach allows estimates to be made in a manner that accounts for the underlying processes leading to SOC change. Ecosystem models, designed for site scale applications can be linked to spatial databases, giving spatially explicit results that allow geographic areas of change in SOC stocks to be identified. Some studies have used variations on this approach to estimate SOC stock changes at the sub-national and national scale for areas of the USA and Europe and at the watershed scale for areas of Mexico and Cuba. However, a need remained for a national and regional scale, spatially explicit system that is generically applicable and can be applied to as wide a range of soil types, climates and land uses as possible. The Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System was developed in response to this need. The GEFSOC system allows estimates of SOC stocks and changes to be made for diverse conditions, providing essential information for countries wishing to take part in an emerging C market, and bringing us closer to an understanding of the future role of soils in the global C cycle. (C) 2007 Elsevier B.V. All rights reserved.

Results from the Covered Catchment Experiment at Gårdsjön, Sweden, after ten years of clean precipitation treatment

The Covered Catchment Experiment at Gordsjon is a large scale forest ecosystem manipulation, where acid precipitation was intercepted by a 7000 m(2) plastic roof and replaced by 'clean precipitation' sprinkled below the roof for ten years between 1991 and 2001. The treatment resulted in a strong positive response of runoff quality. The runoff sulphate, inorganic aluminium and base cations decreased, while there was a strong increase in runoff ANC and a moderate increase in pH. The runoff continued to improve over the whole duration of the experiment. The achieved quality was, however, after ten years still considerably worse than estimated pre-industrial runoff at the site. Stable isotopes of sulphur were analysed to study the soil sulphur cycling. At the initial years of the experiment, the desorption of SO4 from the mineral soil appeared to control the runoff SO4 concentration. However, as the experiment proceeded, there was growing evidence that net mineralisation of soil organic sulphur in the humus layer was an additional source of SO4 in runoff. This might provide a challenge to current acidification models. The experiment convincingly demonstrated on a catchment scale, that reduction in acid deposition causes an immediate improvement of surface water quality even at heavily acidified sites. The improvement of the runoff appeared to be largely a result of cation exchange processes in the soil due to decreasing concentrations of the soil solution, while any potential change in soil base saturation seemed to be less important for the runoff chemistry over the short time period of one decade. These findings should be considered when interpreting and extrapolating regional trends in surface water chemistry to the terrestrial parts of ecosystems.