Interpreting Differential Temperature Trends at the Surface and in the Lower Troposphere

Estimated global-scale temperature trends at Earth's surface (as recorded by thermometers) and in the lower troposphere (as monitored by satellites) diverge by up to 0.14°C per decade over the period 1979 to 1998. Accounting for differences in the spatial coverage of satellite and surface measurements reduces this differential, but still leaves a statistically significant residual of roughly 0.1°C per decade. Natural internal climate variability alone, as simulated in three state-of-the-art coupled atmosphere-ocean models, cannot completely explain this residual trend difference. A model forced by a combination of anthropogenic factors and volcanic aerosols yields surface-troposphere temperature trend differences closest to those observed.

Why is the global warming proceeding much slower than expected?

Upper air observations from radiosondes and microwave satellite instruments does not indicate any global warming during the last 19 years, contrary to surface measurements, where a warming trend is supposedly being found. This result is somewhat difficult to reconcile, since climate model experiments do indicate a reverse trend, namely, that upper tropospheric air should warm faster than the surface. To contribute toward an understanding of this difficulty, we have here undertaken some specific experiments to study the effect on climate due to the decrease in stratospheric ozone and the Mount Pinatubo eruption in 1991. The associated forcing was added to the forcing from greenhouse gases, sulfate aerosols (direct and indirect effect), and tropospheric ozone, which was investigated in a separate series of experiments. Furthermore, we have undertaken an ensemble study in order to explore the natural variability of an advanced climate model exposed to such a forcing over 19 years. The result shows that the reduction of stratospheric ozone cools not only the lower stratosphere but also the troposphere, in particular, the upper and middle part. In the upper troposphere the cooling from stratospheric ozone leads to a significant reduction of greenhouse warming. The modeled stratospheric aerosols from Mount Pinatubo generate a climate response (stratospheric warming and tropospheric cooling) in good agreement with microwave satellite measurements. Finally, analysis of a series of experiments with both stratospheric ozone and the Mount Pinatubo effect shows considerable variability in climate response, suggesting that an evolution having no warming in the period is as likely as another evolution showing modest warming. However, the observed trend of no warming in the midtroposphere and clear warming at the surface is not found in the model simulations.

Transient climate change simulations with a coupled atmosphere–ocean GCM including the tropospheric sulfur cycle

The time-dependent climate response to changing concentrations of greenhouse gases and sulfate aerosols is studied using a coupled general circulation model of the atmosphere and the ocean (ECHAM4/OPYC3). The concentrations of the well-mixed greenhouse gases like CO2, CH4, N2O, and CFCs are prescribed for the past (1860–1990) and projected into the future according to International Panel on Climate Change (IPCC) scenario IS92a. In addition, the space–time distribution of tropospheric ozone is prescribed, and the tropospheric sulfur cycle is calculated within the coupled model using sulfur emissions of the past and projected into the future (IS92a). The radiative impact of the aerosols is considered via both the direct and the indirect (i.e., through cloud albedo) effect. It is shown that the simulated trend in sulfate deposition since the end of the last century is broadly consistent with ice core measurements, and the calculated radiative forcings from preindustrial to present time are within the uncertainty range estimated by IPCC. Three climate perturbation experiments are performed, applying different forcing mechanisms, and the results are compared with those obtained from a 300-yr unforced control experiment. As in previous experiments, the climate response is similar, but weaker, if aerosol effects are included in addition to greenhouse gases. One notable difference to previous experiments is that the strength of the Indian summer monsoon is not fundamentally affected by the inclusion of aerosol effects. Although the monsoon is damped compared to a greenhouse gas only experiment, it is still more vigorous than in the control experiment. This different behavior, compared to previous studies, is the result of the different land–sea distribution of aerosol forcing. Somewhat unexpected, the intensity of the global hydrological cycle becomes weaker in a warmer climate if both direct and indirect aerosol effects are included in addition to the greenhouse gases. This can be related to anomalous net radiative cooling of the earth’s surface through aerosols, which is balanced by reduced turbulent transfer of both sensible and latent heat from the surface to the atmosphere.

Detection of anthropogenic climate change using a fingerprint method

A fingerprint method for detecting anthropogenic climate change is applied to new simulations with a coupled ocean-atmosphere general circulation model (CGCM) forced by increasing concentrations of greenhouse gases and aerosols covering the years 1880 to 2050. In addition to the anthropogenic climate change signal, the space-time structure of the natural climate variability for near-surface temperatures is estimated from instrumental data over the last 134 years and two 1000 year simulations with CGCMs. The estimates are compared with paleoclimate data over 570 years. The space-time information on both the signal and the noise is used to maximize the signal-to-noise ratio of a detection variable obtained by applying an optimal filter (fingerprint) to the observed data. The inclusion of aerosols slows the predicted future warming. The probability that the observed increase in near-surface temperatures in recent decades is of natural origin is estimated to be less than 5%. However, this number is dependent on the estimated natural variability level, which is still subject to some uncertainty.

Inhalation exposure to particulate matter in rooms with underfloor air distribution

This work investigated the personal exposure to indoor particulate matters using the intake fraction metric and provided a possible way to trace the particle inhaled from an indoor particle source. A turbulence model validated by the particle measurements in a room with underfloor air distribution (UFAD) system was used to predict the indoor particle concentrations. Inhalation intake fraction of indoor particles was defined and evaluated in two rooms equipped with the UFAD, i.e., the experimental room and a small office. According to the exposure characteristics and a typical respiratory rate, the intake fraction was determined in two rooms with a continuous and episodic (human cough) source of particles, respectively. The findings showed that the well-mixing assumption of indoor air failed to give an accurate estimation of inhalation exposure and the average concentration at return outlet or within the overall room could not relate well the intake fraction to the amount of particle emitted from an indoor source.

The influence of anthropogenic aerosol on multi-decadal variations of historical global climate

Analysis of single forcing runs from CMIP5 (the fifth Coupled Model Intercomparison Project) simulations shows that the mid-twentieth century temperature hiatus, and the coincident decrease in precipitation, is likely to have been influenced strongly by anthropogenic aerosol forcing. Models that include a representation of the indirect effect of aerosol better reproduce inter-decadal variability in historical global-mean near-surface temperatures, particularly the cooling in the 1950s and 1960s, compared to models with representation of the aerosol direct effect only. Models with the indirect effect also show a more pronounced decrease in precipitation during this period, which is in better agreement with observations, and greater inter-decadal variability in the inter-hemispheric temperature difference. This study demonstrates the importance of representing aerosols, and their indirect effects, in general circulation models, and suggests that inter-model diversity in aerosol burden and representation of aerosol–cloud interaction can produce substantial variation in simulations of climate variability on multi decadal timescales.

Ozone climatology using interactive chemistry: results from the Canadian Middle Atmosphere Model

The climatology of ozone produced by the Canadian Middle Atmosphere Model (CMAM) is presented. This three-dimensional global model incorporates the radiative feedbacks of ozone and water vapor calculated on-line with a photochemical module. This module includes a comprehensive gas-phase reaction set and a limited set of heterogeneous reactions to account for processes occurring on background sulphate aerosols. While transport is global, photochemistry is solved from about 400 hPa to the top of the model at ∼95 km. This approach provides a complete and comprehensive representation of transport, emission, and photochemistry of various constituents from the surface to the mesopause region. A comparison of model results with observations indicates that the ozone distribution and variability are in agreement with observations throughout most of the model domain. Column ozone annual variation is represented to within 5–10% of the observations except in the Southern Hemisphere for springtime high latitudes. The vertical ozone distribution is generally well represented by the model up to the mesopause region. Nevertheless, in the upper stratosphere, the model generally underestimates the amount of ozone as well as the latitudinal tilting of ozone isopleths at high latitude. Ozone variability is analyzed and compared with measurements. The comparison shows that the phase and amplitude of the seasonal variation as well as shorter timescale variations are well represented by the model at various latitudes and heights. Finally, the impact of incorporating ozone radiative feedback on the model climatology is isolated. It is found that the incorporation of ozone radiative feedback results in a cooling of ∼8 K in the summer stratopause region, which corrects a warm bias that results when climatological ozone is used.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

High-latitude climate sensitivity to ice-sheet forcing over the last 120kyr

Interpretation of ice-core records is currently limited by paucity of modelling at adequate temporal and spatial resolutions. Several key questions relate to mechanisms of polar amplification and inter-hemispheric coupling on glacial/interglacial timescales. Here, we present the first results from a large set of global ocean–atmosphere climate model ‘snap-shot’ simulations covering the last 120 000 years using the Hadley Centre climate model (HadCM3) at up to 1 kyr temporal resolution. Two sets of simulations were performed in order to examine the roles of orbit and greenhouse gases versus ice-sheet forcing of orbital-scale climate change. A series of idealised Heinrich events were also simulated, but no changes to aerosols or vegetation were prescribed. This paper focuses on high latitudes and inter-hemispheric linkages. The simulations reproduce polar temperature trends well compared to ice-core reconstructions, although the magnitude is underestimated. Polar amplification varies with obliquity, but this variability is dampened by including variations in land ice coverage, while the overall amplification factor increases. The relatively constant amplification of Antarctic temperatures (with ice-sheet forcing included) suggests it is possible to use Antarctic temperature reconstructions to estimate global changes (which are roughly half the magnitude). Atlantic Ocean overturning circulation varies considerably only with the introduction of Northern Hemisphere ice sheets, but only weakens in the North Atlantic in the deep glacial, when ocean–sea-ice feedbacks result in the movement of the region of deep convection to lower latitudes and with the introduction of freshwater to the surface North Atlantic in order to simulate Heinrich events.

Estimating the direct aerosol radiative perturbation: impact of ocean surface representation and aerosol non-sphericity

Atmospheric aerosols are now actively studied, in particular because of their radiative and climate impacts. Estimations of the direct aerosol radiative perturbation, caused by extinction of incident solar radiation, usually rely on radiative transfer codes and involve simplifying hypotheses. This paper addresses two approximations which are widely used for the sake of simplicity and limiting the computational cost of the calculations. Firstly, it is shown that using a Lambertian albedo instead of the more rigorous bidirectional reflectance distribution function (BRDF) to model the ocean surface radiative properties leads to large relative errors in the instantaneous aerosol radiative perturbation. When averaging over the day, these errors cancel out to acceptable levels of less than 3% (except in the northern hemisphere winter). The other scope of this study is to address aerosol non-sphericity effects. Comparing an experimental phase function with an equivalent Mie-calculated phase function, we found acceptable relative errors if the aerosol radiative perturbation calculated for a given optical thickness is daily averaged. However, retrieval of the optical thickness of non-spherical aerosols assuming spherical particles can lead to significant errors. This is due to significant differences between the spherical and non-spherical phase functions. Discrepancies in aerosol radiative perturbation between the spherical and non-spherical cases are sometimes reduced and sometimes enhanced if the aerosol optical thickness for the spherical case is adjusted to fit the simulated radiance of the non-spherical case.

Global estimate of aerosol direct radiative forcing from satellite measurements

Atmospheric aerosols cause scattering and absorption of incoming solar radiation. Additional anthropogenic aerosols released into the atmosphere thus exert a direct radiative forcing on the climate system1. The degree of present-day aerosol forcing is estimated from global models that incorporate a representation of the aerosol cycles1–3. Although the models are compared and validated against observations, these estimates remain uncertain. Previous satellite measurements of the direct effect of aerosols contained limited information about aerosol type, and were confined to oceans only4,5. Here we use state-of-the-art satellitebased measurements of aerosols6–8 and surface wind speed9 to estimate the clear-sky direct radiative forcing for 2002, incorporating measurements over land and ocean. We use a Monte Carlo approach to account for uncertainties in aerosol measurements and in the algorithm used. Probability density functions obtained for the direct radiative forcing at the top of the atmosphere give a clear-sky, global, annual average of 21.9Wm22 with standard deviation, 60.3Wm22. These results suggest that present-day direct radiative forcing is stronger than present model estimates, implying future atmospheric warming greater than is presently predicted, as aerosol emissions continue to decline10.

Estimates of aerosol radiative forcing from the MACC re-analysis

The European Centre for Medium-range Weather Forecast (ECMWF) provides an aerosol re-analysis starting from year 2003 for the Monitoring Atmospheric Composition and Climate (MACC) project. The re-analysis assimilates total aerosol optical depth retrieved by the Moderate Resolution Imaging Spectroradiometer (MODIS) to correct for model departures from observed aerosols. The reanalysis therefore combines satellite retrievals with the full spatial coverage of a numerical model. Re-analysed products are used here to estimate the shortwave direct and first indirect radiative forcing of anthropogenic aerosols over the period 2003–2010, using methods previously applied to satellite retrievals of aerosols and clouds. The best estimate of globally-averaged, all-sky direct radiative forcing is −0.7±0.3Wm−2. The standard deviation is obtained by a Monte-Carlo analysis of uncertainties, which accounts for uncertainties in the aerosol anthropogenic fraction, aerosol absorption, and cloudy-sky effects. Further accounting for differences between the present-day natural and pre-industrial aerosols provides a direct radiative forcing estimate of −0.4±0.3Wm−2. The best estimate of globally-averaged, all-sky first indirect radiative forcing is −0.6±0.4Wm−2. Its standard deviation accounts for uncertainties in the aerosol anthropogenic fraction, and in cloud albedo and cloud droplet number concentration susceptibilities to aerosol changes. The distribution of first indirect radiative forcing is asymmetric and is bounded by −0.1 and −2.0Wm−2. In order to decrease uncertainty ranges, better observational constraints on aerosol absorption and sensitivity of cloud droplet number concentrations to aerosol changes are required.

Impact of the modal aerosol scheme GLOMAP-mode on aerosol forcing in the Hadley Centre Global Environmental Model

The Hadley Centre Global Environmental Model (HadGEM) includes two aerosol schemes: the Coupled Large-scale Aerosol Simulator for Studies in Climate (CLASSIC), and the new Global Model of Aerosol Processes (GLOMAP-mode). GLOMAP-mode is a modal aerosol microphysics scheme that simulates not only aerosol mass but also aerosol number, represents internally-mixed particles, and includes aerosol microphysical processes such as nucleation. In this study, both schemes provide hindcast simulations of natural and anthropogenic aerosol species for the period 2000–2006. HadGEM simulations of the aerosol optical depth using GLOMAP-mode compare better than CLASSIC against a data-assimilated aerosol re-analysis and aerosol ground-based observations. Because of differences in wet deposition rates, GLOMAP-mode sulphate aerosol residence time is two days longer than CLASSIC sulphate aerosols, whereas black carbon residence time is much shorter. As a result, CLASSIC underestimates aerosol optical depths in continental regions of the Northern Hemisphere and likely overestimates absorption in remote regions. Aerosol direct and first indirect radiative forcings are computed from simulations of aerosols with emissions for the year 1850 and 2000. In 1850, GLOMAP-mode predicts lower aerosol optical depths and higher cloud droplet number concentrations than CLASSIC. Consequently, simulated clouds are much less susceptible to natural and anthropogenic aerosol changes when the microphysical scheme is used. In particular, the response of cloud condensation nuclei to an increase in dimethyl sulphide emissions becomes a factor of four smaller. The combined effect of different 1850 baselines, residence times, and abilities to affect cloud droplet number, leads to substantial differences in the aerosol forcings simulated by the two schemes. GLOMAP-mode finds a presentday direct aerosol forcing of −0.49Wm−2 on a global average, 72% stronger than the corresponding forcing from CLASSIC. This difference is compensated by changes in first indirect aerosol forcing: the forcing of −1.17Wm−2 obtained with GLOMAP-mode is 20% weaker than with CLASSIC. Results suggest that mass-based schemes such as CLASSIC lack the necessary sophistication to provide realistic input to aerosol-cloud interaction schemes. Furthermore, the importance of the 1850 baseline highlights how model skill in predicting present-day aerosol does not guarantee reliable forcing estimates. Those findings suggest that the more complex representation of aerosol processes in microphysical schemes improves the fidelity of simulated aerosol forcings.

Observations of the eruption of the Sarychev volcano and simulations using the HadGEM2 climate model

In June 2009 the Sarychev volcano located in the Kuril Islands to the northeast of Japan erupted explosively, injecting ash and an estimated 1.2 ± 0.2 Tg of sulfur dioxide into the upper troposphere and lower stratosphere, making it arguably one of the 10 largest stratospheric injections in the last 50 years. During the period immediately after the eruption, we show that the sulfur dioxide (SO2) cloud was clearly detected by retrievals developed for the Infrared Atmospheric Sounding Interferometer (IASI) satellite instrument and that the resultant stratospheric sulfate aerosol was detected by the Optical Spectrograph and Infrared Imaging System (OSIRIS) limb sounder and CALIPSO lidar. Additional surface‐based instrumentation allows assessment of the impact of the eruption on the stratospheric aerosol optical depth. We use a nudged version of the HadGEM2 climate model to investigate how well this state‐of‐the‐science climate model can replicate the distributions of SO2 and sulfate aerosol. The model simulations and OSIRIS measurements suggest that in the Northern Hemisphere the stratospheric aerosol optical depth was enhanced by around a factor of 3 (0.01 at 550 nm), with resultant impacts upon the radiation budget. The simulations indicate that, in the Northern Hemisphere for July 2009, the magnitude of the mean radiative impact from the volcanic aerosols is more than 60% of the direct radiative forcing of all anthropogenic aerosols put together. While the cooling induced by the eruption will likely not be detectable in the observational record, the combination of modeling and measurements would provide an ideal framework for simulating future larger volcanic eruptions.

Highly contrasting effects of different climate forcing agents on terrestrial ecosystem services

Many atmospheric constituents besides carbon dioxide (CO2) contribute to global warming, and it is common to compare their influence on climate in terms of radiative forcing, which measures their impact on the planetary energy budget. A number of recent studies have shown that many radiatively active constituents also have important impacts on the physiological functioning of ecosystems, and thus the ‘ecosystem services’ that humankind relies upon. CO2 increases have most probably increased river runoff and had generally positive impacts on plant growth where nutrients are non-limiting, whereas increases in near-surface ozone (O3) are very detrimental to plant productivity. Atmospheric aerosols increase the fraction of surface diffuse light, which is beneficial for plant growth. To illustrate these differences, we present the impact on net primary productivity and runoff of higher CO2, higher near-surface O3, and lower sulphate aerosols, and for equivalent changes in radiative forcing.We compare this with the impact of climate change alone, arising, for example, from a physiologically inactive gas such as methane (CH4). For equivalent levels of change in radiative forcing, we show that the combined climate and physiological impacts of these individual agents vary markedly and in some cases actually differ in sign. This study highlights the need to develop more informative metrics of the impact of changing atmospheric constituents that go beyond simple radiative forcing.