Black carbon aerosol plays a unique and important role in Earth's climate system. Black carbon is a type of carbonaceous material with a unique combination of physical properties. This assessment provides an evaluation of black‐carbon climate forcing that is comprehensive in its inclusion of all known and relevant processes and that is quantitative in providing best estimates and uncertainties of the main forcing terms: direct solar absorption; influence on liquid, mixed phase, and ice clouds; and deposition on snow and ice. These effects are calculated with climate models, but when possible, they are evaluated with both microphysical measurements and field observations. Predominant sources are combustion related, namely, fossil fuels for transportation, solid fuels for industrial and residential uses, and open burning of biomass. Total global emissions of black carbon using bottom‐up inventory methods are 7500 Gg yr⁻¹in the year 2000 with an uncertainty range of 2000 to 29000. However, global atmospheric absorption attributable to black carbon is too low in many models and should be increased by a factor of almost 3. After this scaling, the best estimate for the industrial‐era (1750 to 2005) direct radiative forcing of atmospheric black carbon is +0.71 W m⁻²with 90% uncertainty bounds of (+0.08, +1.27) W m⁻². Total direct forcing by all black carbon sources, without subtracting the preindustrial background, is estimated as +0.88 (+0.17, +1.48) W m⁻². Direct radiative forcing alone does not capture important rapid adjustment mechanisms. A framework is described and used for quantifying climate forcings, including rapid adjustments. The best estimate of industrial‐era climate forcing of black carbon through all forcing mechanisms, including clouds and cryosphere forcing, is +1.1 W m⁻²with 90% uncertainty bounds of +0.17 to +2.1 W m⁻². Thus, there is a very high probability that black carbon emissions, independent of co‐emitted species, have a positive forcing and warm the climate. We estimate that black carbon, with a total climate forcing of +1.1 W m⁻², is the second most important human emission in terms of its climate forcing in the present‐day atmosphere; only carbon dioxide is estimated to have a greater forcing. Sources that emit black carbon also emit other short‐lived species that may either cool or warm climate. Climate forcings from co‐emitted species are estimated and used in the framework described herein. When the principal effects of short‐lived co‐emissions, including cooling agents such as sulfur dioxide, are included in net forcing, energy‐related sources (fossil fuel and biofuel) have an industrial‐era climate forcing of +0.22 (−0.50 to +1.08) W m⁻²during the first year after emission. For a few of these sources, such as diesel engines and possibly residential biofuels, warming is strong enough that eliminating all short‐lived emissions from these sources would reduce net climate forcing (i.e., produce cooling). When open burning emissions, which emit high levels of organic matter, are included in the total, the best estimate of net industrial‐era climate forcing by all short‐lived species from black‐carbon‐rich sources becomes slightly negative (−0.06 W m⁻²with 90% uncertainty bounds of −1.45 to +1.29 W m⁻²). The uncertainties in net climate forcing from black‐carbon‐rich sources are substantial, largely due to lack of knowledge about cloud interactions with both black carbon and co‐emitted organic carbon. In prioritizing potential black‐carbon mitigation actions, non‐science factors, such as technical feasibility, costs, policy design, and implementation feasibility play important roles. The major sources of black carbon are presently in different stages with regard to the feasibility for near‐term mitigation. This assessment, by evaluating the large number and complexity of the associated physical and radiative processes in black‐carbon climate forcing, sets a baseline from which to improve future climate forcing estimates.

We apply our Snow, Ice, and Aerosol Radiative (SNICAR) model, coupled to a general circulation model with prognostic carbon aerosol transport, to improve understanding of climate forcing and response from black carbon (BC) in snow. Building on two previous studies, we account for interannually varying biomass burning BC emissions, snow aging, and aerosol scavenging by snow meltwater. We assess uncertainty in forcing estimates from these factors, as well as BC optical properties and snow cover fraction. BC emissions are the largest source of uncertainty, followed by snow aging. The rate of snow aging determines snowpack effective radius (r_e), which directly controls snow reflectance and the magnitude of albedo change caused by BC. For a reasonable r_e range, reflectance reduction from BC varies threefold. Inefficient meltwater scavenging keeps hydrophobic impurities near the surface during melt and enhances forcing. Applying biomass burning BC emission inventories for a strong (1998) and weak (2001) boreal fire year, we estimate global annual mean BC/snow surface radiative forcing from all sources (fossil fuel, biofuel, and biomass burning) of +0.054 (0.007–0.13) and +0.049 (0.007–0.12) W m⁻², respectively. Snow forcing from only fossil fuel + biofuel sources is +0.043 W m⁻² (forcing from only fossil fuels is +0.033 W m⁻²), suggesting that the anthropogenic contribution to total forcing is at least 80%. The 1998 global land and sea‐ice snowpack absorbed 0.60 and 0.23 W m⁻², respectively, because of direct BC/snow forcing. The forcing is maximum coincidentally with snowmelt onset, triggering strong snow‐albedo feedback in local springtime. Consequently, the “efficacy” of BC/snow forcing is more than three times greater than forcing by CO₂. The 1998 and 2001 land snowmelt rates north of 50°N are 28% and 19% greater in the month preceding maximum melt of control simulations without BC in snow. With climate feedbacks, global annual mean 2‐meter air temperature warms 0.15 and 0.10°C, when BC is included in snow, whereas annual arctic warming is 1.61 and 0.50°C. Stronger high‐latitude climate response in 1998 than 2001 is at least partially caused by boreal fires, which account for nearly all of the 35% biomass burning contribution to 1998 arctic forcing. Efficacy was anomalously large in this experiment, however, and more research is required to elucidate the role of boreal fires, which we suggest have maximum arctic BC/snow forcing potential during April–June. Model BC concentrations in snow agree reasonably well (r = 0.78) with a set of 23 observations from various locations, spanning nearly 4 orders of magnitude. We predict concentrations in excess of 1000 ng g⁻¹ for snow in northeast China, enough to lower snow albedo by more than 0.13. The greatest instantaneous forcing is over the Tibetan Plateau, exceeding 20 W m⁻² in some places during spring. These results indicate that snow darkening is an important component of carbon aerosol climate forcing.

We describe a model for predicting the size‐resolved distribution of atmospheric dust for climate and chemistry‐related studies. The dust distribution from 1990 to 1999 is simulated with our mineral aerosol entrainment and deposition module embedded in a chemical transport model. Mobilization processes include entrainment thresholds for saltation, moisture inhibition, drag partitioning, and saltation feedback. For mobilization we assume that soil texture is globally uniform and is replete with saltators. Soil erodibility is prescribed by a new physically based geomorphic index that is proportional to the runoff area upstream of each source region. Dry deposition processes include sedimentation and turbulent mix‐out. Nucleation scavenging and size‐resolved washout in both stratiform and convective cloud types are represented. Simulations of the 1990s broadly agree with station observations and satellite‐inferred dust distributions. Without invoking anthropogenic mechanisms the model captures the seasonal migration of the transatlantic African dust plume, and it captures the spring maximum in Asian dust outflow and concentration over the Pacific. We estimate the 1990s global annual mean and variability of D < 10 μm dust to be the following: emissions, 1490 ± 160 Tg yr⁻¹; burden, 17 ± 2 Tg; and optical depth at 0.63 μm, 0.030 ± 0.004. This emission, burden, and optical depth are significantly lower than some recent estimates. The model underestimates transport and deposition of East Asian and Australian dust to some regions of the Pacific Ocean. An underestimate of long‐range transport of particles larger than 3 μm contributes to this bias. Our experiments support the hypothesis that dust emission “hot spots” exist in regions where alluvial sediments have accumulated and may be disturbed.

This study presents the results of a broad intercomparison of a total of 15 global aerosol models within the AeroCom project. Each model is compared to observations related to desert dust aerosols, their direct radiative effect, and their impact on the biogeochemical cycle, i.e., aerosol optical depth (AOD) and dust deposition. Additional comparisons to Angström exponent (AE), coarse mode AOD and dust surface concentrations are included to extend the assessment of model performance and to identify common biases present in models. These data comprise a benchmark dataset that is proposed for model inspection and future dust model development. There are large differences among the global models that simulate the dust cycle and its impact on climate. In general, models simulate the climatology of vertically integrated parameters (AOD and AE) within a factor of two whereas the total deposition and surface concentration are reproduced within a factor of 10. In addition, smaller mean normalized bias and root mean square errors are obtained for the climatology of AOD and AE than for total deposition and surface concentration. Characteristics of the datasets used and their uncertainties may influence these differences. Large uncertainties still exist with respect to the deposition fluxes in the southern oceans. Further measurements and model studies are necessary to assess the general model performance to reproduce dust deposition in ocean regions sensible to iron contributions. Models overestimate the wet deposition in regions dominated by dry deposition. They generally simulate more realistic surface concentration at stations downwind of the main sources than at remote ones. Most models simulate the gradient in AOD and AE between the different dusty regions. However the seasonality and magnitude of both variables is better simulated at African stations than Middle East ones. The models simulate the offshore transport of West Africa throughout the year but they overestimate the AOD and they transport too fine particles. The models also reproduce the dust transport across the Atlantic in the summer in terms of both AOD and AE but not so well in winter-spring nor the southward displacement of the dust cloud that is responsible of the dust transport into South America. Based on the dependency of AOD on aerosol burden and size distribution we use model bias with respect to AOD and AE to infer the bias of the dust emissions in Africa and the Middle East. According to this analysis we suggest that a range of possible emissions for North Africa is 400 to 2200 Tg yr⁻¹ and in the Middle East 26 to 526 Tg yr⁻¹.

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission‐relevant water cycle questions. Its components include atmosphere and land (110‐km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO₂ forcing simulations. The model performs well overall with biases typical of other CMIP‐class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP‐class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two‐layer energy balance model, we attribute this divergence to the model's strong aerosol‐related effective radiative forcing (ERF_ari+aci = −1.65 W/m²) and high equilibrium climate sensitivity (ECS = 5.3 K).

We report measurements and analysis of a boreal forest fire, integrating the effects of greenhouse gases, aerosols, black carbon deposition on snow and sea ice, and postfire changes in surface albedo. The net effect of all agents was to increase radiative forcing during the first year (34 ± 31 Watts per square meter of burned area), but to decrease radiative forcing when averaged over an 80-year fire cycle (–2.3 ± 2.2 Watts per square meter) because multidecadal increases in surface albedo had a larger impact than fire-emitted greenhouse gases. This result implies that future increases in boreal fire may not accelerate climate warming.