Seasonal variations of the physical and optical characteristics of Saharan Dust: results from the Dust outflow and Deposition to the Ocean (DODO) experiment

North African dust is important for climate through its direct radiative effect on solar and terrestrial radiation and its role in the biogeochemical system. The Dust Outflow and Deposition to the Ocean project (DODO) aimed to characterize the physical and optical properties of airborne North African dust in two seasons and to use these observations to constrain model simulations, with the ultimate aim of being able to quantify the deposition of iron to the North Atlantic Ocean. The in situ properties of dust from airborne campaigns measured during February and August 2006, based at Dakar, Senegal, are presented here. Average values of the single scattering albedo (0.99, 0.98), mass specific extinction (0.85 m^2 g^-1 , 1.14 m^2 g^-1 ), asymmetry parameter (0.68, 0.68), and refractive index (1.53--0.0005i,1.53--0.0014i) for the accumulation mode were found to differ by varying degrees between the dry and wet season, respectively. It is hypothesized that these differences are due to different source regions and transport processes which also differ between the DODO campaigns. Elemental ratios of Ca/Al were found to differ between the dry and wet season (1.1 and 0.5, respectively). Differences in vertical profiles are found between seasons and between land and ocean locations and reflect the different dynamics of the seasons. Using measurements of the coarse mode size distribution and illustrative Mie calculations, the optical properties are found to be very sensitive to the presence and amount of coarse mode of mineral dust, and the importance of accurate measurements of the coarse mode of dust is highlighted.

An overview of Earth's global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

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.