The physical properties of the atmosphere in the new Hadley Centre Global Environmental Model (HadGEM1). Part II: Aspects of variability and regional climate

The performance of the atmospheric component of the new Hadley Centre Global Environmental Model (HadGEM1) is assessed in terms of its ability to represent a selection of key aspects of variability in the Tropics and extratropics. These include midlatitude storm tracks and blocking activity, synoptic variability over Europe, and the North Atlantic Oscillation together with tropical convection, the Madden-Julian oscillation, and the Asian summer monsoon. Comparisons with the previous model, the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3), demonstrate that there has been a considerable increase in the transient eddy kinetic energy (EKE), bringing HadGEM1 into closer agreement with current reanalyses. This increase in EKE results from the increased horizontal resolution and, in combination with the improved physical parameterizations, leads to improvements in the representation of Northern Hemisphere storm tracks and blocking. The simulation of synoptic weather regimes over Europe is also greatly improved compared to HadCM3, again due to both increased resolution and other model developments. The variability of convection in the equatorial region is generally stronger and closer to observations than in HadCM3. There is, however, still limited convective variance coincident with several of the observed equatorial wave modes. Simulation of the Madden-Julian oscillation is improved in HadGEM1: both the activity and interannual variability are increased and the eastward propagation, although slower than observed, is much better simulated. While some aspects of the climatology of the Asian summer monsoon are improved in HadGEM1, the upper-level winds are too weak and the simulation of precipitation deteriorates. The dominant modes of monsoon interannual variability are similar in the two models, although in HadCM3 this is linked to SST forcing, while in HadGEM1 internal variability dominates. Overall, analysis of the phenomena considered here indicates that HadGEM1 performs well and, in many important respects, improves upon HadCM3. Together with the improved representation of the mean climate, this improvement in the simulation of atmospheric variability suggests that HadGEM1 provides a sound basis for future studies of climate and climate change.

It is premature to include non-CO2 effects of aviation in emission trading schemes

The recent G8 Gleneagles climate statement signed on 8 July 2005 specifically mentions a determination to lessen the impact of aviation on climate [Gleneagles, 2005. The Gleneagles communique: climate change, energy and sustainable development. http://www.fco.gov.uk/Files/kfile/PostG8_Gleneagles_Communique.pdf]. In January 2005 the European Union Emission Trading Scheme (ETS) commenced operation as the largest multi-country, multi-sector ETS in the world, albeit currently limited only to CO2 emissions. At present the scheme makes no provision for aircraft emissions. However, the UK Government would like to see aircraft included in the ETS and plans to use its Presidencies of both the EU and G8 in 2005 to implement these schemes within the EU and perhaps internationally. Non-CO2 effects have been included in some policy-orientated studies of the impact of aviation but we argue that the inclusion of such effects in any such ETS scheme is premature; we specifically argue that use of the Radiative Forcing Index for comparing emissions from different sources is inappropriate and that there is currently no metric for such a purpose that is likely to enable their inclusion in the near future. (c) 2005 Elsevier Ltd. All rights reserved.

An estimate of the global impact of multiple scattering by clouds on outgoing long-wave radiation

Although the potential importance of scattering of long-wave radiation by clouds has been recognised, most studies have concentrated on the impact of high clouds and few estimates of the global impact of scattering have been presented. This study shows that scattering in low clouds has a significant impact on outgoing long-wave radiation (OLR) in regions of marine stratocumulus (-3.5 W m(-2) for overcast conditions) where the column water vapour is relatively low. This corresponds to an enhancement of the greenhouse effect of such clouds by 10%. The near-global impact of scattering on OLR is estimated to be -3.0 W m(-2), with low clouds contributing -0.9 W m(-2), mid-level cloud -0.7 W m(-2) and high clouds -1.4 W m(-2). Although this effect appears small compared to the global mean OLR of 240 W m(-2), it indicates that neglect of scattering will lead to an error in cloud long-wave forcing of about 10% and an error in net cloud forcing of about 20%.

Mesoscale simulations of organized convection: importance of convective equilibrium

The validity of convective parametrization breaks down at the resolution of mesoscale models, and the success of parametrized versus explicit treatments of convection is likely to depend on the large-scale environment. In this paper we examine the hypothesis that a key feature determining the sensitivity to the environment is whether the forcing of convection is sufficiently homogeneous and slowly varying that the convection can be considered to be in equilibrium. Two case studies of mesoscale convective systems over the UK, one where equilibrium conditions are expected and one where equilibrium is unlikely, are simulated using a mesoscale forecasting model. The time evolution of area-average convective available potential energy and the time evolution and magnitude of the timescale of convective adjustment are consistent with the hypothesis of equilibrium for case 1 and non-equilibrium for case 2. For each case, three experiments are performed with different partitionings between parametrized and explicit convection: fully parametrized convection, fully explicit convection and a simulation with significant amounts of both. In the equilibrium case, bulk properties of the convection such as area-integrated rain rates are insensitive to the treatment of convection. However, the detailed structure of the precipitation field changes; the simulation with parametrized convection behaves well and produces a smooth field that follows the forcing region, and the simulation with explicit convection has a small number of localized intense regions of precipitation that track with the mid-levelflow. For the non-equilibrium case, bulk properties of the convection such as area-integrated rain rates are sensitive to the treatment of convection. The simulation with explicit convection behaves similarly to the equilibrium case with a few localized precipitation regions. In contrast, the cumulus parametrization fails dramatically and develops intense propagating bows of precipitation that were not observed. The simulations with both parametrized and explicit convection follow the pattern seen in the other experiments, with a transition over the duration of the run from parametrized to explicit precipitation. The impact of convection on the large-scaleflow, as measured by upper-level wind and potential-vorticity perturbations, is very sensitive to the partitioning of convection for both cases. © Royal Meteorological Society, 2006. Contributions by P. A. Clark and M. E. B. Gray are Crown Copyright.

Boundary layer structure and decoupling from synoptic scale flow during NAMBLEX

This paper presents an overview of the meteorology and planetary boundary layer structure observed during the NAMBLEX field campaign to aid interpretation of the chemical and aerosol measurements. The campaign has been separated into five periods corresponding to the prevailing synoptic condition. Comparisons between meteorological measurements ( UHF wind profiler, Doppler sodar, sonic aneometers mounted on a tower at varying heights and a standard anemometer) and the ECMWF analysis at 10 m and 1100 m identified days when the internal boundary layer was decoupled from the synoptic flow aloft. Generally the agreement was remarkably good apart from during period one and on a few days during period four when the diurnal swing in wind direction implies a sea/land breeze circulation near the surface. During these periods the origin of air sampled at Mace Head would not be accurately represented by back trajectories following the winds resolved in ECMWF analyses. The wind profiler observations give a detailed record of boundary layer structure including an indication of its depth, average wind speed and direction. Turbulence statistics have been used to assess the height to which the developing internal boundary layer, caused by the increased surface drag at the coast, reaches the sampling location under a wide range of marine conditions. Sampling conducted below 10 m will be impacted by emission sources at the shoreline in all wind directions and tidal conditions, whereas sampling above 15 m is unlikely to be affected in any of the wind directions and tidal heights sampled during the experiment.

The climate sensitivity and its components diagnosed from Earth Radiation Budget data

One of the major uncertainties in the ability to predict future climate change, and hence its impacts, is the lack of knowledge of the earth's climate sensitivity. Here, data are combined from the 1985-96 Earth Radiation Budget Experiment (ERBE) with surface temperature change information and estimates of radiative forcing to diagnose the climate sensitivity. Importantly, the estimate is completely independent of climate model results. A climate feedback parameter of 2.3 +/- 1.4 W m(-2) K-1 is found. This corresponds to a 1.0-4.1-K range for the equilibrium warming due to a doubling of carbon dioxide (assuming Gaussian errors in observable parameters, which is approximately equivalent to a uniform "prior" in feedback parameter). The uncertainty range is due to a combination of the short time period for the analysis as well as uncertainties in the surface temperature time series and radiative forcing time series, mostly the former. Radiative forcings may not all be fully accounted for; however, all argument is presented that the estimate of climate sensitivity is still likely to be representative of longer-term climate change. The methodology can be used to 1) retrieve shortwave and longwave components of climate feedback and 2) suggest clear-sky and cloud feedback terms. There is preliminary evidence of a neutral or even negative longwave feedback in the observations, suggesting that current climate models may not be representing some processes correctly if they give a net positive longwave feedback.

The meteorology of the Western Indian Ocean, and the influence of the east African highlands

This paper reviews the meteorology of the Western Indian Ocean and uses a state–of–the–art atmospheric general circulation model to investigate the influence of the East African Highlands on the climate of the Indian Ocean and its surrounding regions. The new 44–year re–analysis produced by the European Centre for Medium range Weather Forecasts (ECMWF) has been used to construct a new climatology of the Western Indian Ocean. A brief overview of the seasonal cycle of the Western Indian Ocean is presented which emphasizes the importance of the geography of the Indian Ocean basin for controlling the meteorology of the Western Indian Ocean. The principal modes of inter–annual variability are described, associated with El Niño and the Indian Ocean Dipole or Zonal Mode, and the basic characteristics of the subseasonal weather over the Western Indian Ocean are presented, including new statistics on cyclone tracks derived from the ECMWF re–analyses. Sensitivity experiments, in which the orographic effects of East Africa are removed, have shown that the East African Highlands, although not very high, play a significant role in the climate of Africa, India and Southeast Asia, and in the heat, salinity and momentum forcing of the Western Indian Ocean. The hydrological cycle over Africa is systematically enhanced in all seasons by the presence of the East African Highlands, and during the Asian summer monsoon there is a major redistribution of the rainfall across India and Southeast Asia. The implied impact of the East African Highlands on the ocean is substantial. The East African Highlands systematically freshen the tropical Indian Ocean, and act to focus the monsoon winds along the coast, leading to greater upwelling and cooler sea–surface temperatures.

On the initiation of surface waves by turbulent shear flow

An analytical model is developed for the initial stage of surface wave generation at an air-water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417-445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate F, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer. For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as proportional to 1/Gamma), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.

Indian Ocean-monsoon coupled interactions and impending monsoon droughts

Monsoon droughts over the Indian subcontinent emanate from failures in the seasonal (June-September) monsoon rains. While prolonged dry-spells ("monsoon-breaks'') pervade on sub-seasonal/intra-seasonal time-scales, the underlying causes for these long-lasting anomalies remain elusive. Based on analyses of a suite of observed data sets, we report an ocean-atmosphere dynamical coupling on intra-seasonal time-scales, in the tropical Indian Ocean, which is pivotal in forcing extended monsoon-breaks and causing droughts over the subcontinent. This coupling involves a feedback between the monsoonal flow and thermocline depth in the Equatorial Eastern Indian Ocean (EEIO), in which an anomaly of the summer monsoon circulation induces downwelling and maintains a higher-than-normal heat-content. The near-equatorial anomalies induce strong and sustained suppression of monsoon rainfall over the subcontinent. It is concluded that the intra-seasonal evolution of the ocean-monsoon coupled system is a vital key to unlocking the dynamics of monsoon droughts.

Krakatoa lives: The effect of volcanic eruptions on ocean heat content and thermal expansion

A suite of climate model experiments indicates that 20th Century increases in ocean heat content and sea-level ( via thermal expansion) were substantially reduced by the 1883 eruption of Krakatoa. The volcanically-induced cooling of the ocean surface is subducted into deeper ocean layers, where it persists for decades. Temporary reductions in ocean heat content associated with the comparable eruptions of El Chichon ( 1982) and Pinatubo ( 1991) were much shorter lived because they occurred relative to a non-stationary background of large, anthropogenically-forced ocean warming. Our results suggest that inclusion of the effects of Krakatoa ( and perhaps even earlier eruptions) is important for reliable simulation of 20th century ocean heat uptake and thermal expansion. Inter-model differences in the oceanic thermal response to Krakatoa are large and arise from differences in external forcing, model physics, and experimental design. Systematic experimentation is required to quantify the relative importance of these factors. The next generation of historical forcing experiments may require more careful treatment of pre-industrial volcanic aerosol loadings.

Mechanisms of ocean heat uptake in a coupled climate model and the implications for tracer based predictions of ocean heat uptake

The distribution of tracers in the ocean is often taken as an indication of the ventilation pathways for oceanic water masses. It has been suggested that under anthropogenic forcing heat will be taken up into the interior of the ocean along isopycnal ventilation pathways. This notion is investigated by examining distributions of potential temperature and a passive anomaly temperature tracer in a coupled climate experiment where CO2 is increased at a rate of 2% per year. We show that interior temperature changes cannot be explained solely by passive tracer transport along isopycnals. Heat uptake is strongly affected by changes in circulation and has a substantial diapycnal component.

Changes in the global sulfate burden due to perturbations in global CO2 concentrations

A large ensemble of general circulation model (GCM) integrations coupled to a fully interactive sulfur cycle scheme were run on the climateprediction.net platform to investigate the uncertainty in the climate response to sulfate aerosol and carbon dioxide (CO2) forcing. The sulfate burden within the model (and the atmosphere) depends on the balance between formation processes and deposition (wet and dry). The wet removal processes for sulfate aerosol are much faster than dry removal and so any changes in atmospheric circulation, cloud cover, and precipitation will feed back on the sulfate burden. When CO2 is doubled in the Hadley Centre Slab Ocean Model (HadSM3), global mean precipitation increased by 5%; however, the global mean sulfate burden increased by 10%. Despite the global mean increase in precipitation, there were large areas of the model showing decreases in precipitation (and cloud cover) in the Northern Hemisphere during June–August, which reduced wet deposition and allowed the sulfate burden to increase. Further experiments were also undertaken with and without doubling CO2 while including a future anthropogenic sulfur emissions scenario. Doubling CO2 further enhanced the increases in sulfate burden associated with increased anthropogenic sulfur emissions as observed in the doubled CO2-only experiment. The implications are that the climate response to doubling CO2 can influence the amount of sulfate within the atmosphere and, despite increases in global mean precipitation, may act to increase it.

The significance of volcanic eruption strength and frequency for climate

A simple physical model of the atmospheric effects of large explosive volcanic eruptions is developed. Using only one input parameter - the initial amount of sulphur dioxide injected into the stratosphere - the global-average stratospheric optical-depth perturbation and surface temperature response are modelled. The simplicity of this model avoids issues of incomplete data (applicable to more comprehensive models), making it a powerful and useful tool for atmospheric diagnostics of this climate forcing mechanism. It may also provide a computationally inexpensive and accurate way of introducing volcanic activity into larger climate models. The modelled surface temperature response for an initial sulphur-dioxide injection, coupled with emission-history statistics, is used to demonstrate that the most climatically significant volcanic eruptions are those of sufficient explosivity to just reach into the stratosphere (and achieve longevity). This study also highlights the fact that this measure of significance is highly sensitive to the representation of the climatic response and the frequency data used, and that we are far from producing a definitive history of explosive volcanism for at least the past 1000 years. Given this high degree of uncertainty, these results suggest that eruptions that release around and above 0.1 Mt SO2 into the stratosphere have the maximum climatic impact.

Factors influencing anthropogenic carbon dioxide uptake in the North Atlantic in models of the ocean carbon cycle

The uptake and storage of anthropogenic carbon in the North Atlantic is investigated using different configurations of ocean general circulation/carbon cycle models. We investigate how different representations of the ocean physics in the models, which represent the range of models currently in use, affect the evolution of CO2 uptake in the North Atlantic. The buffer effect of the ocean carbon system would be expected to reduce ocean CO2 uptake as the ocean absorbs increasing amounts of CO2. We find that the strength of the buffer effect is very dependent on the model ocean state, as it affects both the magnitude and timing of the changes in uptake. The timescale over which uptake of CO2 in the North Atlantic drops to below preindustrial levels is particularly sensitive to the ocean state which sets the degree of buffering; it is less sensitive to the choice of atmospheric CO2 forcing scenario. Neglecting physical climate change effects, North Atlantic CO2 uptake drops below preindustrial levels between 50 and 300 years after stabilisation of atmospheric CO2 in different model configurations. Storage of anthropogenic carbon in the North Atlantic varies much less among the different model configurations, as differences in ocean transport of dissolved inorganic carbon and uptake of CO2 compensate each other. This supports the idea that measured inventories of anthropogenic carbon in the real ocean cannot be used to constrain the surface uptake. Including physical climate change effects reduces anthropogenic CO2 uptake and storage in the North Atlantic further, due to the combined effects of surface warming, increased freshwater input, and a slowdown of the meridional overturning circulation. The timescale over which North Atlantic CO2 uptake drops to below preindustrial levels is reduced by about one-third, leading to an estimate of this timescale for the real world of about 50 years after the stabilisation of atmospheric CO2. In the climate change experiment, a shallowing of the mixed layer depths in the North Atlantic results in a significant reduction in primary production, reducing the potential role for biology in drawing down anthropogenic CO2.

Direct radiative effect of aerosols emitted by transport from road, shipping and aviation

Aerosols and their precursors are emitted abundantly by transport activities. Transportation constitutes one of the fastest growing activities and its growth is predicted to increase significantly in the future. Previous studies have estimated the aerosol direct radiative forcing from one transport sub-sector, but only one study to our knowledge estimated the range of radiative forcing from the main aerosol components (sulphate, black carbon (BC) and organic carbon) for the whole transportation sector. In this study, we compare results from two different chemical transport models and three radiation codes under different hypothesis of mixing: internal and external mixing using emission inventories for the year 2000. The main results from this study consist of a positive direct radiative forcing for aerosols emitted by road traffic of +20±11 mW m−2 for an externally mixed aerosol, and of +32±13 mW m−2 when BC is internally mixed. These direct radiative forcings are much higher than the previously published estimate of +3±11 mW m−2. For transport activities from shipping, the net direct aerosol radiative forcing is negative. This forcing is dominated by the contribution of the sulphate. For both an external and an internal mixture, the radiative forcing from shipping is estimated at −26±4 mW m−2. These estimates are in very good agreement with the range of a previously published one (from −46 to −13 mW m−2) but with a much narrower range. By contrast, the direct aerosol forcing from aviation is estimated to be small, and in the range −0.9 to +0.3 mW m−2.