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.

As global temperatures increase, sea ice loss will increasingly enable commercial shipping traffic to cross the Arctic Ocean, where the ships' gas and particulate emissions may have strong regional effects. Here we investigate impacts of shipping emissions on Arctic climate using a fully coupled Earth system model (CESM 1.2.2) and a suite of newly developed projections of 21st-century trans-Arctic shipping emissions. We find that trans-Arctic shipping will reduce Arctic warming by nearly 1 °C by 2099, due to sulfate-driven liquid water cloud formation. Cloud fraction and liquid water path exhibit significant positive trends, cooling the lower atmosphere and surface. Positive feedbacks from sea ice growth-induced albedo increases and decreased downwelling longwave radiation due to reduced water vapor content amplify the cooling relative to the shipping-free Arctic. Our findings thus point to the complexity in Arctic climate responses to increased shipping traffic, justifying further study and policy considerations as trade routes open.

Desert dust simulations generated by the National Center for Atmospheric Research's Community Climate System Model for the current climate are shown to be consistent with present day satellite and deposition data. The response of the dust cycle to last glacial maximum, preindustrial, modern, and doubled-carbon dioxide climates is analyzed. Only natural (non-land use related) dust sources are included in this simulation. Similar to some previous studies, dust production mainly responds to changes in the source areas from vegetation changes, not from winds or soil moisture changes alone. This model simulates a +92%, +33%, and −60% change in dust loading for the last glacial maximum, preindustrial, and doubled-carbon dioxide climate, respectively, when impacts of carbon dioxide fertilization on vegetation are included in the model. Terrestrial sediment records from the last glacial maximum compiled here indicate a large underestimate of deposition in continental regions, probably due to the lack of simulation of glaciogenic dust sources. In order to include the glaciogenic dust sources as a first approximation, we designate the location of these sources, and infer the size of the sources using an inversion method that best matches the available data. The inclusion of these inferred glaciogenic dust sources increases our dust flux in the last glacial maximum from 2.1 to 3.3 times current deposition.

Atmospheric rivers (ARs), intrusions of warm and moist air, can effectively drive weather extremes over the Arctic and trigger subsequent impact on sea ice and climate. What controls the observed multi-decadal Arctic AR trends remains unclear. Here, using multiple sources of observations and model experiments, we find that, contrary to the uniform positive trend in climate simulations, the observed Arctic AR frequency increases by twice as much over the Atlantic sector compared to the Pacific sector in 1981-2021. This discrepancy can be reconciled by the observed positive-to-negative phase shift of Interdecadal Pacific Oscillation (IPO) and the negative-to-positive phase shift of Atlantic Multidecadal Oscillation (AMO), which increase and reduce Arctic ARs over the Atlantic and Pacific sectors, respectively. Removing the influence of the IPO and AMO can reduce the projection uncertainties in near-future Arctic AR trends by about 24%, which is important for constraining projection of Arctic warming and the timing of an ice-free Arctic.

This study examines multiyear climate variability associated with sea salt aerosols and their contribution to the variability of shortwave cloud forcing (SWCF) using a 150 year simulation for preindustrial conditions of the Community Earth System Model version 1.0. The results suggest that changes in sea salt and related cloud and radiative properties on interannual timescales are dominated by the El Niño-Southern Oscillation cycle. Sea salt variability on longer (interdecadal) timescales is associated with low-frequency variability in the Pacific Ocean similar to the Interdecadal Pacific Oscillation but does not show a statistically significant spectral peak. A multivariate regression suggests that sea salt aerosol variability may contribute to SWCF variability in the tropical Pacific, explaining up to 20-30% of the variance in that region. Elsewhere, there is only a small sea salt aerosol influence on SWCF through modifying cloud droplet number and liquid water path that contributes to the change of cloud effective radius and cloud optical depth (and hence cloud albedo), producing a multiyear aerosol-cloud-wind interaction.

We use three 150 year preindustrial simulations of the Community Earth System Model to quantify the impacts of El Niño-Southern Oscillation (ENSO) events on shortwave and longwave cloud radiative effects (CRESW and CRELW). Compared to recent observations from the Clouds and the Earth’s Radiant Energy System data set, the model simulation successfully reproduces larger variations of CRESW and CRELW over the tropics. The ENSO cycle is found to dominate interannual variations of cloud radiative effects. Simulated cooling (warming) effects from CRESW (CRELW) are strongest over the tropical western and central Pacific Ocean during warm ENSO events, with the largest difference between 20 and 60Wm-2, with weaker effects of 10-40Wm-2 over Indonesian regions and the subtropical Pacific Ocean. Sensitivity tests show that variations of cloud radiative effects are mainly driven by ENSO-related changes in cloud fraction. The variations in midlevel and high cloud fractions each account for approximately 20-50% of the interannual variations of CRESW over the tropics and almost all of the variations of CRELW between 60°S and 60°N. The variation of low cloud fraction contributes to most of the variations of CRESW over the midlatitude oceans. Variations in natural aerosol concentrations explained 10-30% of the variations of both CRESW and CRELW over the tropical Pacific, Indonesian regions, and the tropical Indian Ocean. Changes in natural aerosol emissions and concentrations enhance 3-5% and 1-3% of the variations of cloud radiative effects averaged over the tropics.

The effect of an increase in atmospheric aerosol concentrations on the distribution and radiative properties of Earth's clouds is the most uncertain component of the overall global radiative forcing from preindustrial time. General circulation models (GCMs) are the tool for predicting future climate, but the treatment of aerosols, clouds, and aerosol-cloud radiative effects carries large uncertainties that directly affect GCM predictions, such as climate sensitivity. Predictions are hampered by the large range of scales of interaction between various components that need to be captured. Observation systems (remote sensing, in situ) are increasingly being used to constrain predictions, but significant challenges exist, to some extent because of the large range of scales and the fact that the various measuring systems tend to address different scales. Fine-scale models represent clouds, aerosols, and aerosol-cloud interactions with high fidelity but do not include interactions with the larger scale and are therefore limited from a climatic point of view. We suggest strategies for improving estimates of aerosol-cloud relationships in climate models, for new remote sensing and in situ measurements, and for quantifying and reducing model uncertainty.

Numerical weather, climate, or Earth system models involve the coupling of components. At a broad level, these components can be classified as the resolved fluid dynamics, unresolved fluid dynamical aspects (i.e., those represented by physical parameterizations such as subgrid-scale mixing), and nonfluid dynamical aspects such as radiation and microphysical processes. Typically, each component is developed, at least initially, independently.Once development ismature, the components are coupled to deliver a model of the required complexity. The implementation of the coupling can have a significant impact on the model.As the error associated with each component decreases, the errors introduced by the coupling will eventually dominate. Hence, any improvement in one of the components is unlikely to improve the performance of the overall system. The challenges associated with combining the components to create a coherentmodel are here termed physics-dynamics coupling. The issue goes beyond the coupling between the parameterizations and the resolved fluid dynamics. This paper highlights recent progress and some of the current challenges. It focuses on three objectives: to illustrate the phenomenology of the coupling problemwith references to examples in the literature, to show howthe problem can be analyzed, and to create awareness of the issue across the disciplines and specializations. The topics addressed are different ways of advancing full models in time, approaches to understanding the role of the coupling and evaluation of approaches, coupling ocean and atmosphere models, thermodynamic compatibility between model components, and emerging issues such as those that arise as model resolutions increase and/ormodels use variable resolutions.