Host model uncertainties in aerosol radiative forcing estimates: results from the AeroCom Prescribed intercomparison study

Simulated multi-model “diversity” in aerosol direct radiative forcing estimates is often perceived as a measure of aerosol uncertainty. However, current models used for aerosol radiative forcing calculations vary considerably in model components relevant for forcing calculations and the associated “host-model uncertainties” are generally convoluted with the actual aerosol uncertainty. In this AeroCom Prescribed intercomparison study we systematically isolate and quantify host model uncertainties on aerosol forcing experiments through prescription of identical aerosol radiative properties in twelve participating models. Even with prescribed aerosol radiative properties, simulated clear-sky and all-sky aerosol radiative forcings show significant diversity. For a purely scattering case with globally constant optical depth of 0.2, the global-mean all-sky top-of-atmosphere radiative forcing is −4.47Wm−2 and the inter-model standard deviation is 0.55Wm−2, corresponding to a relative standard deviation of 12 %. For a case with partially absorbing aerosol with an aerosol optical depth of 0.2 and single scattering albedo of 0.8, the forcing changes to 1.04Wm−2, and the standard deviation increases to 1.01W−2, corresponding to a significant relative standard deviation of 97 %. However, the top-of-atmosphere forcing variability owing to absorption (subtracting the scattering case from the case with scattering and absorption) is low, with absolute (relative) standard deviations of 0.45Wm−2 (8 %) clear-sky and 0.62Wm−2 (11 %) all-sky. Scaling the forcing standard deviation for a purely scattering case to match the sulfate radiative forcing in the Aero- Com Direct Effect experiment demonstrates that host model uncertainties could explain about 36% of the overall sulfate forcing diversity of 0.11Wm−2 in the AeroCom Direct Radiative Effect experiment.

A review of measurement-based assessments of the aerosol direct radiative effect and forcing

Aerosols affect the Earth's energy budget directly by scattering and absorbing radiation and indirectly by acting as cloud condensation nuclei and, thereby, affecting cloud properties. However, large uncertainties exist in current estimates of aerosol forcing because of incomplete knowledge concerning the distribution and the physical and chemical properties of aerosols as well as aerosol-cloud interactions. In recent years, a great deal of effort has gone into improving measurements and datasets. It is thus feasible to shift the estimates of aerosol forcing from largely model-based to increasingly measurement-based. Our goal is to assess current observational capabilities and identify uncertainties in the aerosol direct forcing through comparisons of different methods with independent sources of uncertainties. Here we assess the aerosol optical depth (τ), direct radiative effect (DRE) by natural and anthropogenic aerosols, and direct climate forcing (DCF) by anthropogenic aerosols, focusing on satellite and ground-based measurements supplemented by global chemical transport model (CTM) simulations. The multi-spectral MODIS measures global distributions of aerosol optical depth (τ) on a daily scale, with a high accuracy of ±0.03±0.05τ over ocean. The annual average τ is about 0.14 over global ocean, of which about 21%±7% is contributed by human activities, as estimated by MODIS fine-mode fraction. The multi-angle MISR derives an annual average AOD of 0.23 over global land with an uncertainty of ~20% or ±0.05. These high-accuracy aerosol products and broadband flux measurements from CERES make it feasible to obtain observational constraints for the aerosol direct effect, especially over global the ocean. A number of measurement-based approaches estimate the clear-sky DRE (on solar radiation) at the top-of-atmosphere (TOA) to be about -5.5±0.2 Wm-2 (median ± standard error from various methods) over the global ocean. Accounting for thin cirrus contamination of the satellite derived aerosol field will reduce the TOA DRE to -5.0 Wm-2. Because of a lack of measurements of aerosol absorption and difficulty in characterizing land surface reflection, estimates of DRE over land and at the ocean surface are currently realized through a combination of satellite retrievals, surface measurements, and model simulations, and are less constrained. Over the oceans the surface DRE is estimated to be -8.8±0.7 Wm-2. Over land, an integration of satellite retrievals and model simulations derives a DRE of -4.9±0.7 Wm-2 and -11.8±1.9 Wm-2 at the TOA and surface, respectively. CTM simulations derive a wide range of DRE estimates that on average are smaller than the measurement-based DRE by about 30-40%, even after accounting for thin cirrus and cloud contamination. A number of issues remain. Current estimates of the aerosol direct effect over land are poorly constrained. Uncertainties of DRE estimates are also larger on regional scales than on a global scale and large discrepancies exist between different approaches. The characterization of aerosol absorption and vertical distribution remains challenging. The aerosol direct effect in the thermal infrared range and in cloudy conditions remains relatively unexplored and quite uncertain, because of a lack of global systematic aerosol vertical profile measurements. A coordinated research strategy needs to be developed for integration and assimilation of satellite measurements into models to constrain model simulations. Enhanced measurement capabilities in the next few years and high-level scientific cooperation will further advance our knowledge.

A multimodel assessment of the influence of regional anthropogenic emission reductions on aerosol direct radiative forcing and the role of intercontinental transport

In this study, we assess changes of aerosol optical depth (AOD) and direct radiative forcing (DRF) in response to the reduction of anthropogenic emissions in four major pollution regions in the Northern Hemisphere by using results from nine global models in the framework of the Hemispheric Transport of Air Pollution (HTAP). DRF at top of atmosphere (TOA) and surface is estimated based on AOD results from the HTAP models and AOD-normalized DRF (NDRF) from a chemical transport model. The multimodel results show that, on average, a 20% reduction of anthropogenic emissions in North America, Europe, East Asia, and South Asia lowers the global mean AOD (all-sky TOA DRF) by 9.2% (9.0%), 3.5% (3.0%), and 9.4% (10.0%) for sulfate, particulate organic matter (POM), and black carbon (BC), respectively. Global annual average TOA all-sky forcing efficiency relative to particle or gaseous precursor emissions from the four regions (expressed as multimodel mean ± one standard deviation) is  ±3.5 ±0.8,  ±4.0 ±1.7, and 29.5 ±18.1mWm ±2 per Tg for sulfate (relative to SO2), POM, and BC, respectively. The impacts of the regional emission reductions on AOD and DRF extend well beyond the source regions because of intercontinental transport (ICT). On an annual basis, ICT accounts for 11 ±5% to 31 ±9% of AOD and DRF in a receptor region at continental or subcontinental scale, with domestic emissions accounting for the remainder, depending on regions and species. For sulfate AOD, the largest ICT contribution of 31 ±9% occurs in South Asia, which is dominated by the emissions from Europe. For BC AOD, the largest ICT contribution of 28 ±18% occurs in North America, which is dominated by the emissions from East Asia. The large spreads among models highlight the need to improve aerosol processes in models, and evaluate and constrain models with observations.

Radiative forcing of climate by hydrochlorofluorocarbons and hydrofluorocarbons

We measure infrared absorption spectra of 18 hydrochlorofluorocarbons and hydrofluorocarbons, seven of which do not yet appear in the literature. The spectra are used in a narrowband model of the terrestrial infrared radiation to calculate radiative forcing and global warming potentials. We investigate the sensitivity of the radiative forcing to the absorption spectrum temperature dependence, halocarbon vertical profile, stratospheric adjustment, cloudiness, spectral overlap, and latitude, and we make some recommendations for the reporting of radiative forcings that would help to resolve discrepancies between assessments. We investigate simple methods of estimating instantaneous radiative forcing directly from a molecule's absorption spectrum and we present a new method that agrees to within 0.3% with our narrowband model results.

Radiative forcing and global warming potentials of 11 halogenated compounds

Radiative forcing values have been calculated for 11 halogenated compounds which are in current use or which have been suggested as possible replacements for the chlorofluorocarbons. Absorption cross-sections measured over a range of atmospheric temperature and pressure conditions as part of a multi-laboratory programme have been used together with a narrow band radiative transfer model. We provide a “best estimate” radiative forcing taking into account the likely vertical profile of the gas in each case. The Global Warming Potential over a variety of time horizons has also been calculated where the lifetime is available. We present the first such information for 1,2-dichloroethane. For chloroform our radiative forcing is 5 times higher than the value used in previous assessments, possibly because these ignored the effect of absorption outside the 800–1200 cm−1 “window”. For several of the other compounds considered here, our forcing is between 10 and 30% lower than previous assessments. The perfluorocarbons have been found to have large global warming potentials, many times that of CFC-11, due to both strong absorption and long lifetimes. The importance of absorption features at wavenumbers below 800 cm−1 and the effect of temperature variations in absorption cross-section on the radiative forcing are also investigated.

Updated radiative forcing estimates of 65 halocarbons and nonmethane hydrocarbons

The direct radiative forcing of 65 chlorofluorocarbons, hydrochlorofluorocarbons, hydrofluorocarbons, hydrofluoroethers, halons, iodoalkanes, chloroalkanes, bromoalkanes, perfluorocarbons and nonmethane hydrocarbons has been evaluated using a consistent set of infrared absorption cross sections. For the radiative transfer models, both line-by-line and random band model approaches were employed for each gas. The line-by-line model was first validated against measurements taken by the Airborne Research Interferometer Evaluation System (ARIES) of the U.K. Meteorological Office; the computed spectrally integrated radiance of agreed to within 2% with experimental measurements. Three model atmospheres, derived from a three-dimensional climatology, were used in the radiative forcing calculations to more accurately represent hemispheric differences in water vapor, ozone concentrations, and cloud cover. Instantaneous, clear-sky radiative forcing values calculated by the line-by-line and band models were in close agreement. The band model values were subsequently modified to ensure exact agreement with the line-by-line model values. Calibrated band model radiative forcing values, for atmospheric profiles with clouds and using stratospheric adjustment, are reported and compared with previous literature values. Fourteen of the 65 molecules have forcings that differ by more than 15% from those in the World Meteorological Organization [1999] compilation. Eleven of the molecules have not been reported previously. The 65-molecule data set reported here is the most comprehensive and consistent database yet available to evaluate the relative impact of halocarbons and hydrocarbons on climate change.

Sensitivity of direct radiative forcing by mineral dust to particle characteristics

Airborne dust affects the Earth's energy balance — an impact that is measured in terms of the implied change in net radiation (or radiative forcing, in W m-2) at the top of the atmosphere. There remains considerable uncertainty in the magnitude and sign of direct forcing by airborne dust under current climate. Much of this uncertainty stems from simplified assumptions about mineral dust-particle size, composition and shape, which are applied in remote sensing retrievals of dust characteristics and dust-cycle models. Improved estimates of direct radiative forcing by dust will require improved characterization of the spatial variability in particle characteristics to provide reliable information dust optical properties. This includes constraints on: (1) particle-size distribution, including discrimination of particle subpopulations and quantification of the amount of dust in the sub-10 µm to <0.1 µm mass fraction; (2) particle composition, specifically the abundance of iron oxides, and whether particles consist of single or multi-mineral grains; (3) particle shape, including degree of sphericity and surface roughness, as a function of size and mineralogy; and (4) the degree to which dust particles are aggregated together. The use of techniques that measure the size, composition and shape of individual particles will provide a better basis for optical modelling.

Radiative forcing effect of ice-age dust

During glacial periods, dust deposition rates and inferred atmospheric concentrations were globally much higher than present. According to recent model results, the large enhancement of atmospheric dust content at the last glacial maximum (LGM) can be explained only if increases in the potential dust source areas are taken into account. Such increases are to be expected, due to effects of low precipitation and low atmospheric (CO2) on plant growth. Here the modelled three-dimensional dust fields from Mahowald et al. and modelled seasonally varying surface-albedo fields derived in a parallel manner, are used to quantify the mean radiative forcing due to modern (non-anthropogenic) and LGM dust. The effect of mineralogical provenance on the radiative properties of the dust is taken into account, as is the range of optical properties associated with uncertainties about the mixing state of the dust particles. The high-latitude (poleward of 45°) mean change in forcing (LGM minus modern) is estimated to be small (–0.9 to +0.2 W m–2), especially when compared to nearly –20 W m–2 due to reflection from the extended ice sheets. Although the net effect of dust over ice sheets is a positive forcing (warming), much of the simulated high-latitude dust was not over the ice sheets, but over unglaciated regions close to the expanded dust source region in central Asia. In the tropics the change in forcing is estimated to be overall negative, and of similarly large magnitude (–2.2 to –3.2 W m–2) to the radiative cooling effect of low atmospheric (CO2). Thus, the largest long-term climatic effect of the LGM dust is likely to have been a cooling of the tropics. Low tropical sea-surface temperatures, low atmospheric (CO2) and high atmospheric dust loading may be mutually reinforcing due to multiple positive feedbacks, including the negative radiative forcing effect of dust.

Changes in global net radiative imbalance 1985-2012

Combining satellite data, atmospheric reanalyses and climate model simulations, variability in the net downward radiative flux imbalance at the top of Earth's atmosphere (N) is reconstructed and linked to recent climate change. Over the 1985-1999 period mean N (0.34 ± 0.67 Wm–2) is lower than for the 2000-2012 period (0.62 ± 0.43 Wm–2, uncertainties at 90% confidence level) despite the slower rate of surface temperature rise since 2000. While the precise magnitude of N remains uncertain, the reconstruction captures interannual variability which is dominated by the eruption of Mt. Pinatubo in 1991 and the El Niño Southern Oscillation. Monthly deseasonalized interannual variability in N generated by an ensemble of 9 climate model simulations using prescribed sea surface temperature and radiative forcings and from the satellite-based reconstruction is significantly correlated (r ∼ 0.6) over the 1985-2012 period.

Infrared absorption spectra, radiative efficiencies, and global warming potentials of newly-detected halogenated compounds: CFC-113a, CFC-112 and HCFC-133a

CFC-113a (CF3CCl3), CFC-112 (CFCl2CFCl2) and HCFC-133a (CF3CH2Cl) are three newly detected molecules in the atmosphere that are almost certainly emitted as a result of human activity. It is important to characterise the possible contribution of these gases to radiative forcing of climate change and also to provide information on the CO2-equivalence of their emissions. We report new laboratory measurements of absorption cross-sections of these three compounds at a resolution of 0.01 cm−1 for two temperatures 250 K and 295 K in the spectral range of 600–1730 cm−1. These spectra are then used to calculate the radiative efficiencies and global warming potentials (GWP). The radiative efficiencies are found to be between 0.15 and 0.3 W∙m−2∙ppbv−1. The GWP for a 100 year time horizon, relative to carbon dioxide, ranges from 340 for the relatively short-lived HCFC-133a to 3840 for the longer-lived CFC-112. At current (2012) concentrations, these gases make a trivial contribution to total radiative forcing; however, the concentrations of CFC-113a and HCFC-133a are continuing to increase. The 2012 CO2-equivalent emissions, using the GWP (100), are estimated to be about 4% of the current global CO2-equivalent emissions of HFC-134a

Global radiative and climate effect of the water vapour continuum at visible and near-infrared wavelengths

Recent laboratory measurements show that absorption by the water vapour continuum in near-infrared windows may be about an order of magnitude higher than assumed in many radiation codes. The radiative impact of the continuum at visible and near-infrared wavelengths is examined for the present day and for a possible future warmer climate (with a global-mean total column water increase of 33%). The calculations use a continuum model frequently used in climate models (‘CKD’) and a continuum model where absorption is enhanced at wavelengths greater than 1 µm based on recent measurements (‘CAVIAR’). The continuum predominantly changes the partitioning between solar radiation absorbed by the surface and the atmosphere; changes in top-of-atmosphere net irradiances are smaller. The global-mean clear-sky atmospheric absorption is enhanced by 1.5 W m−2 (about 2%) and 2.8 W m−2 (about 3.5%) for CKD and CAVIAR respectively, relative to a hypothetical no-continuum case, with all-sky enhancements about 80% of these values. The continuum is, in relative terms, more important for radiation budget changes between the present day and a possible future climate. Relative to the no-continuum case, the increase in global-mean clear-sky absorption is 8% higher using CKD and almost 20% higher using CAVIAR; all-sky enhancements are about half these values. The effect of the continuum is estimated for the solar component of the water vapour feedback, the reduction in downward surface irradiance and precipitation change in a warmer world. For CKD and CAVIAR respectively, and relative to the no-continuum case, the solar component of the water vapour feedback is enhanced by about 4 and 9%, the change in clear-sky downward surface irradiance is 7 and 18% more negative, and the global-mean precipitation response decreases by 1 and 4%. There is a continued need for improved continuum measurements, especially at atmospheric temperatures and at wavelengths below 2 µm.

The dependence of radiative forcing and feedback on evolving patterns of surface temperature change in climate models

Experiments with CO2 instantaneously quadrupled and then held constant are used to show that the relationship between the global-mean net heat input to the climate system and the global-mean surface-air-temperature change is nonlinear in Coupled Model Intercomparison Project phase 5 (CMIP5) Atmosphere-Ocean General Circulation Models (AOGCMs). The nonlinearity is shown to arise from a change in strength of climate feedbacks driven by an evolving pattern of surface warming. In 23 out of the 27 AOGCMs examined the climate feedback parameter becomes significantly (95% confidence) less negative – i.e. the effective climate sensitivity increases – as time passes. Cloud feedback parameters show the largest changes. In the AOGCM-mean approximately 60% of the change in feedback parameter comes from the topics (30N-30S). An important region involved is the tropical Pacific where the surface warming intensifies in the east after a few decades. The dependence of climate feedbacks on an evolving pattern of surface warming is confirmed using the HadGEM2 and HadCM3 atmosphere GCMs (AGCMs). With monthly evolving sea-surface-temperatures and sea-ice prescribed from its AOGCM counterpart each AGCM reproduces the time-varying feedbacks, but when a fixed pattern of warming is prescribed the radiative response is linear with global temperature change or nearly so. We also demonstrate that the regression and fixed-SST methods for evaluating effective radiative forcing are in principle different, because rapid SST adjustment when CO2 is changed can produce a pattern of surface temperature change with zero global mean but non-zero change in net radiation at the top of the atmosphere (~ -0.5 Wm-2 in HadCM3).

Intercomparison of shortwave radiative transfer schemes in global aerosol modeling: results from the AeroCom Radiative Transfer Experiment

In this study we examine the performance of 31 global model radiative transfer schemes in cloud-free conditions with prescribed gaseous absorbers and no aerosols (Rayleigh atmosphere), with prescribed scattering-only aerosols, and with more absorbing aerosols. Results are compared to benchmark results from high-resolution, multi-angular line-by-line radiation models. For purely scattering aerosols, model bias relative to the line-by-line models in the top-of-the atmosphere aerosol radiative forcing ranges from roughly −10 to 20%, with over- and underestimates of radiative cooling at lower and higher solar zenith angle, respectively. Inter-model diversity (relative standard deviation) increases from ~10 to 15% as solar zenith angle decreases. Inter-model diversity in atmospheric and surface forcing decreases with increased aerosol absorption, indicating that the treatment of multiple-scattering is more variable than aerosol absorption in the models considered. Aerosol radiative forcing results from multi-stream models are generally in better agreement with the line-by-line results than the simpler two-stream schemes. Considering radiative fluxes, model performance is generally the same or slightly better than results from previous radiation scheme intercomparisons. However, the inter-model diversity in aerosol radiative forcing remains large, primarily as a result of the treatment of multiple-scattering. Results indicate that global models that estimate aerosol radiative forcing with two-stream radiation schemes may be subject to persistent biases introduced by these schemes, particularly for regional aerosol forcing.

Regional and monthly and clear-sky aerosol direct radiative effect (and forcing) derived from the GlobAEROSOL-AATSR satellite aerosol product

Using the GlobAEROSOL-AATSR dataset, estimates of the instantaneous, clear-sky, direct aerosol radiative effect and radiative forcing have been produced for the year 2006. Aerosol Robotic Network sun-photometer measurements have been used to characterise the random and systematic error in the GlobAEROSOL product for 22 regions covering the globe. Representative aerosol properties for each region were derived from the results of a wide range of literature sources and, along with the de-biased GlobAEROSOL AODs, were used to drive an offline version of the Met Office unified model radiation scheme. In addition to the mean AOD, best-estimate run of the radiation scheme, a range of additional calculations were done to propagate uncertainty estimates in the AOD, optical properties, surface albedo and errors due to the temporal and spatial averaging of the AOD fields. This analysis produced monthly, regional estimates of the clear-sky aerosol radiative effect and its uncertainty, which were combined to produce annual, global mean values of (−6.7±3.9)Wm−2 at the top of atmosphere (TOA) and (−12±6)Wm−2 at the surface. These results were then used to give estimates of regional, clear-sky aerosol direct radiative forcing, using modelled pre-industrial AOD fields for the year 1750 calculated for the AEROCOM PRE experiment. However, as it was not possible to quantify the uncertainty in the pre-industrial aerosol loading, these figures can only be taken as indicative and their uncertainties as lower bounds on the likely errors. Although the uncertainty on aerosol radiative effect presented here is considerably larger than most previous estimates, the explicit inclusion of the major sources of error in the calculations suggest that they are closer to the true constraint on this figure from similar methodologies, and point to the need for more, improved estimates of both global aerosol loading and aerosol optical properties.

Aerosol direct radiative effect of smoke over clouds over the southeast Atlantic Ocean from 2006 to 2009

The aerosol direct radiative effect (DRE) of African smoke was analyzed in cloud scenes over the southeast Atlantic Ocean, using Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) satellite observations and Hadley Centre Global Environmental Model version 2 (HadGEM2) climate model simulations. The observed mean DRE was about 30–35 W m−2 in August and September 2006–2009. In some years, short episodes of high-aerosol DRE can be observed, due to high-aerosol loadings, while in other years the loadings are lower but more prolonged. Climate models that use evenly distributed monthly averaged emission fields will not reproduce these high-aerosol loadings. Furthermore, the simulated monthly mean aerosol DRE in HadGEM2 is only about 6 W m−2 in August. The difference with SCIAMACHY mean observations can be partly explained by an underestimation of the aerosol absorption Ångström exponent in the ultraviolet. However, the subsequent increase of aerosol DRE simulation by about 20% is not enough to explain the observed discrepancy between simulations and observations.