Changes in Atmospheric Constituents and in Radiative Forcing

Global tropospheric ozone distributions, budgets, and radiative forcings from an ensemble of 26 state-of-the-art atmospheric chemistry models have been intercompared and synthesized as part of a wider study into both the air quality and climate roles of ozone. Results from three 2030 emissions scenarios, broadly representing optimistic, likely, and pessimistic options, are compared to a base year 2000 simulation. This base case realistically represents the current global distribution of tropospheric ozone. A further set of simulations considers the influence of climate change over the same time period by forcing the central emissions scenario with a surface warming of around 0.7K. The use of a large multimodel ensemble allows us to identify key areas of uncertainty and improves the robustness of the results. Ensemble mean changes in tropospheric ozone burden between 2000 and 2030 for the 3 scenarios range from a 5% decrease, through a 6% increase, to a 15% increase. The intermodel uncertainty (±1 standard deviation) associated with these values is about ±25%. Model outliers have no significant influence on the ensemble mean results. Combining ozone and methane changes, the three scenarios produce radiative forcings of -50, 180, and 300 mW m-2, compared to a CO 2 forcing over the same time period of 800-1100 mW m-2. These values indicate the importance of air pollution emissions in short- to medium-term climate forcing and the potential for stringent/lax control measures to improve/worsen future climate forcing. The model sensitivity of ozone to imposed climate change varies between models but modulates zonal mean mixing ratios by ±5 ppbv via a variety of feedback mechanisms, in particular those involving water vapor and stratosphere-troposphere exchange. This level of climate change also reduces the methane lifetime by around 4%. The ensemble mean year 2000 tropospheric ozone budget indicates chemical production, chemical destruction, dry deposition and stratospheric input fluxes of 5100, 4650, 1000 and 550 Tg(O 3 ) yr-1, respectively. These values are significantly different to the mean budget documented by the Intergovernmental Panel on Climate Change (IPCC) Third Assessment Report (TAR). The mean ozone burden (340 Tg(O 3 )) is 10% larger than the IPCC TAR estimate, while the mean ozone lifetime (22 days) is 10% shorter. Results from individual models show a correlation between ozone burden and lifetime, and each model's ozone burden and lifetime respond in similar ways across the emissions scenarios. The response to climate change is much less consistent. Models show more variability in the tropics compared to midlatitudes. Some of the most uncertain areas of the models include treatments of deep tropical convection, including lightning NO x production; isoprene emissions from vegetation and isoprene's degradation chemistry; stratosphere-troposphere exchange; biomass burning; and water vapor concentrations. Copyright 2006 by the American Geophysical Union.

/pdf/multimodel-ensemble-simulations-of-present-day-and-near-tbq5w9fgf4.pdf

Multimodel ensemble simulations of present-day and near-future tropospheric ozone

Simulation results of global aerosol models have been assembled in the framework of the AeroCom intercomparison exercise. In this paper, we analyze the life cycles of dust, sea salt, sulfate, black carbon and particulate organic matter as simulated by sixteen global aerosol models. The differences among the results (model diversities) for sources and sinks, burdens, particle sizes, water uptakes, and spatial dispersals have been established. These diversities have large consequences for the calculated radiative forcing and the aerosol concentrations at the surface. Processes and parameters are identified which deserve further research. The AeroCom all-models-average emissions are dominated by the mass of sea salt (SS), followed by dust (DU), sulfate (SO 4 ), particulate organic matter (POM), and finally black carbon (BC). Interactive parameterizations of the emissions and contrasting particles sizes of SS and DU lead generally to higher diversities of these species, and for total aerosol. The lower diversity of the emissions of the fine aerosols, BC, POM, and SO 4 , is due to the use of similar emission inventories, and does therefore not necessarily indicate a better understanding of their sources. The diversity of SO 4 -sources is mainly caused by the disagreement on depositional loss of precursor gases and on chemical production. The diversities of the emissions are passed on to the burdens, but the latter are also strongly affected by the model-specific treatments of transport and aerosol processes. The burdens of dry masses decrease from largest to smallest: DU, SS, SO 4 , POM, and BC. The all-models-average residence time is shortest for SS with about half a day, followed by SO 4 and DU with four days, and POM and BC with six and seven days, respectively. The wet deposition rate is controlled by the solubility and increases from DU, BC, POM to SO 4 and SS. It is the dominant sink for SO 4 , BC, and POM, and contributes about one third to the total removal of SS and DU species. For SS and DU we find high diversities for the removal rate coefficients and deposition pathways. Models do neither agree on the split between wet and dry deposition, nor on that between sedimentation and other dry deposition processes. We diagnose an extremely high diversity for the uptake of ambient water vapor that influences the particle size and thus the sink rate coefficients. Furthermore, we find little agreement among the model results for the partitioning of wet removal into scavenging by convective and stratiform rain. Large differences exist for aerosol dispersal both in the vertical and in the horizontal direction. In some models, a minimum of total aerosol concentration is simulated at the surface. Aerosol dispersal is most pronounced for SO 4 and BC and lowest for SS. Diversities are higher for meridional than for vertical dispersal, they are similar for the individual species and highest for SS and DU. For these two components we do not find a correlation between vertical and meridional aerosol dispersal. In addition the degree of dispersals of SS and DU is not related to their residence times. SO 4 , BC, and POM, however, show increased meridional dispersal in models with larger vertical dispersal, and dispersal is larger for longer simulated residence times.

/pdf/analysis-and-quantification-of-the-diversities-of-aerosol-nam6kgi1wz.pdf

Analysis and quantification of the diversities of aerosol life cycles within AeroCom

Monthly and 3-hourly precipitation data from twentieth-century climate simulations by the newest generation of 18 coupled climate system models are analyzed and compared with available observations. The characteristics examined include the mean spatial patterns, intraseasonal-to-interannual and ENSO-related variability, convective versus stratiform precipitation ratio, precipitation frequency and intensity for different precipitation categories, and diurnal cycle. Although most models reproduce the observed broad patterns of precipitation amount and year-to-year variability, models without flux corrections still show an unrealistic double-ITCZ pattern over the tropical Pacific, whereas the flux-corrected models, especially the Meteorological Research Institute (MRI) Coupled Global Climate Model (CGCM; version 2.3.2a), produce realistic rainfall patterns at low latitudes. As in previous generations of coupled models, the rainfall double ITCZs are related to westward expansion of the cold tongue of...

/pdf/precipitation-characteristics-in-eighteen-coupled-climate-3cb39n92yh.pdf

Precipitation Characteristics in Eighteen Coupled Climate Models

A new two-moment stratiform cloud microphysics scheme in a general circulation model is described. Prognostic variables include cloud droplet and cloud ice mass mixing ratios and number concentrations. The scheme treats several microphysical processes, including hydrometeor collection, condensation/ evaporation, freezing, melting, and sedimentation. The activation of droplets on aerosol is physically based and coupled to a subgrid vertical velocity. Unique aspects of the scheme, relative to existing two-moment schemes developed for general circulation models, are the diagnostic treatment of rain and snow number concentration and mixing ratio and the explicit treatment of subgrid cloud water variability for calculation of the microphysical process rates. Numerical aspects of the scheme are described in detail using idealized one-dimensional offline tests of the microphysics. Sensitivity of the scheme to time step, vertical resolution, and numerical method for diagnostic precipitation is investigated over a range of conditions. It is found that, in general, two substeps are required for numerical stability and reasonably small time truncation errors using a time step of 20 min; however, substepping is only required for the precipitation microphysical processes rather than the entire scheme. A new numerical approach for the diagnostic rain and snow produces reasonable results compared to a benchmark simulation, especially at low vertical resolution. Part II of this study details results of the scheme in single-column and global simulations, including comparison with observations.

/pdf/a-new-two-moment-bulk-stratiform-cloud-microphysics-scheme-2qgapbfkah.pdf

A New Two-Moment Bulk Stratiform Cloud Microphysics Scheme in the Community Atmosphere Model, Version 3 (CAM3). Part I: Description and Numerical Tests

An efficient new prognostic cloud water parameterization designed for use in global climate models is described. The scheme allows for life cycle effects in stratiform clouds and permits cloud optical properties to be determined interactively. The parameterization contains representations of all important microphysical processes, including autoconversion, accretion, Bergeron–Findeisen diffusional growth, and cloud/rain water evaporation. Small-scale dynamical processes, including detrainment of convective condensate, cloud-top entrainment instability, and stability-dependent cloud physical thickness variations, are also taken into account. Cloud optical thickness is calculated from the predicted liquid/ice water path and a variable droplet effective radius estimated by assuming constant droplet number concentration. Microphysical and radiative properties are assumed to be different for liquid and ice clouds, and for liquid clouds over land and ocean. The parameterization is validated in several s...

Kenneth K. W. Lo

Papers

A Prognostic Cloud Water Parameterization for Global Climate Models