This work describes the implementation and evaluation of the Slab Ocean Model com16 ponent of the Energy Exascale Earth System Model version 2 (E3SMv2-SOM) and its application to understanding the climate sensitivity to ocean heat transports (OHTs) and CO2 forcing. E3SMv2-SOM reproduces the baseline climate and Equilibrium Climate Sensitivity (ECS) of the fully coupled E3SMv2 experiments reasonably well, with a pattern correlation close to 1 and a global mean bias of less than 1% of the fully coupled surface temperature and precipitation. Sea ice extent and volume are also well reproduced in the SOM. Consistent with general model behavior, the ECS estimated from the SOM (4.5K) exceeds the effective climate sensitivity obtained from extrapolation to equilibrium in the fully coupled model (4.0K). The E3SMv2 baseline climate also shows a large sensitivity to OHT strengths, with a global surface temperature difference of about 4.0◦ C between high-/low-OHT experiments with prescribed forcings derived from fully coupled experiments with realistic/weak ocean circulation strengths. Similar to their forc ng pattern, the surface temperature response occurs mainly over the subpolar regions in both hemispheres. However, the Southern Ocean shows more surface temperature sensitivity to high/low-OHT forcing due to a positive/negative shortwave cloud radiative effect caused by decreases/increases in mid-latitude marine low-level clouds. This large temperature sensitivity also causes an overcompensation between the prescribed OHTs and atmosphere heat transports. The SOM’s ECS estimate is also sensitive to the prescribed OHT and the associated baseline climate it is initialized from; the high-OHT ECS is 0.5K lower than the low-OHT ECS.

Diagnosing the role of internal variability over recent decades is critically important for both model validation and projections of future warming. Recent research suggests that for 1980-2022 internal variability manifested as Global Cooling and Arctic Warming (i-GCAW), leading to enhanced Arctic Amplification (AA) and suppressed global warming over this period. Here we show that the observationally derived i-GCAW is rare in CMIP6 large ensembles, but simulations that do produce similar i-GCAW exhibit a unique and robust internally driven global surface air temperature (SAT) trend pattern. This unique pattern of SAT change features enhanced warming in Barents and Kara Sea and cooling in the tropical Eastern Pacific and Southern Ocean. Given that these features are imprinted in the observed record over recent decades, this work suggests that internal variability makes a crucial contribution to the discrepancy between model-simulated forced SAT trend pattern and observations.

Process modeling of aerosol-cloud interaction is essential to bridging gaps between observational analysis and climate modeling of aerosol effects in the Earth system and eventually reducing climate projection uncertainties. In this study, we examine aerosol-cloud interaction in summertime precipitating shallow cumuli observed during the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE). Aerosols and precipitating shallow cumuli were extensively observed with in-situ and remote-sensing instruments during two research flight cases on 02 June and 07 June, respectively, during the ACTIVATE summer 2021 deployment phase. We perform observational analysis and large-eddy simulation (LES) of aerosol effect on precipitating cumulus in these two cases. Given the measured aerosol size distributions and meteorological conditions, LES is able to reproduce the observed cloud properties by aircraft such as liquid water content (LWC), cloud droplet number concentration (Nc) and effective radius reff. However, it produces smaller liquid water path (LWP) and larger Nc compared to the satellite retrievals. Both 02 and 07 June cases are over warm waters of the Gulf Stream and have a cloud top height over 3 km, but the 07 June case is more polluted and has larger LWC. We find that the aerosol-induced LWP adjustment is dominated by precipitation and is anticorrelated with cloud-top entrainment for both cases. A negative cloud fraction adjustment due to an increase of aerosol number concentration is also shown in the simulations.

A document by Yiling Huo. Click on the document to view its contents.

This work documents version two of the Department of Energy’s Energy Exascale Earth System Model (E3SM). E3SM version 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid-latitudes and 30 km at the equator and poles. The model performance is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima (DECK) simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate is generally realistic, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Programme (WCRP) assessment. However, E3SMv2 significantly underestimates the global mean surface temperature in the second half of the historical record. An analysis of single-forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol-related forcing.

Snow and ice albedo reduction due to deposition of absorbing particles (i.e., snow darkening effect (SDE)) warms the Earth system and is largely attributed to black carbon (BC) and dust. Absorbing organic aerosol (BrC) also contributes to SDE but has received less attention due to uncertainty and challenges in model representation. This work incorporates the SDE of absorbing organic aerosol (BrC) from biomass burning and biofuel sources into the Snow Ice and Aerosol Radiative (SNICAR) model within a variant of the Community Earth System Model (CESM). Additionally, 12 different emission regions of BrC and BC from biomass burning and biofuel sources are tagged to quantify the relative contribution to global and regional SDE. BrC global SDE (0.021–0.056 Wm-2) is larger than other model estimates, corresponding to 37%–98% of the SDE from BC. When compared to observations, BrC simulations have a range in median bias (-2.5%–+21%), with better agreement in the simulations that include BrC photochemical bleaching. The largest relative contributions to global BrC SDE are traced to Northern Asia (23%–31%), Southeast Asia (16%–21%), and South Africa (13%–17%). Transport from Southeast Asia contributes nearly half of the regional BrC SDE in Antarctica (0.084–0.3 Wm-2), which is the largest regional input to global BrC SDE. Lower latitude BrC SDE is correlated with snowmelt, in-snow BrC concentrations, and snow cover fraction, while polar BrC SDE is correlated with surface insolation and snowmelt. This indicates the importance of in-snow processes and snow feedbacks on modeled BrC SDE.

14 15 Key Points: 16 (1) E3SMv1 captures spatial and temporal variability in the observed dust aerosol optical 17 depth, but underestimates long-range transport. 18 (2) The net direct radiative effect of dust simulated by E3SMv1 is-0.42 Wm-2 with a smaller 19 longwave warming than other recent studies. 20 (3) In addition to emission, dry removal of dust are highly sensitive to the increase of 21 horizontal or vertical model resolution. 22 23 24

Nitrate aerosol plays an important role in affecting regional air quality as well as Earth’s climate. However, it is not well represented or even neglected in many global climate models. In this study, we couple the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) module with the four-mode version of the Modal Aerosol Module (MAM4) in DOE’s Energy Exascale Earth System Model version 2 (E3SMv2) to treat nitrate aerosol and its radiative effects. We find that nitrate aerosol simulated by E3SMv2-MAM4-MOSAIC is sensitive to the treatment of gaseous HNO3 transfer to/from interstitial particles related to accommodation coefficients of HNO3 (αHNO3) on dust and non-dust particles. We compare three different treatments of HNO3 transfer: 1) a treatment (MTC_SLOW) that uses a low αHNO3 in the mass transfer coefficient (MTC) calculation; 2) a dust-weighted MTC treatment (MTC_WGT) that uses a high αHNO3 on non-dust particles; and 3) a dust-weighted MTC treatment that also splits coarse mode aerosols into the coarse dust and sea salt sub-modes in MOSAIC (MTC_SPLC). MTC_WGT and MTC_SPLC increase the global annual mean (2005-2014) nitrate burden from 0.096 (MTC_SLOW) to 0.237 and 0.185 Tg N, respectively, mostly in the coarse mode. They also produce stronger nitrate direct radiative forcing (–0.048 and –0.051 W m–2, respectively) and indirect forcing (–0.33 and –0.35 W m–2, respectively) than MTC_SLOW (–0.021 and –0.24 W m–2). All three treatments overestimate nitrate surface concentrations compared with ground-based observations. MTC_WGT and MTC_SPLC improve the vertical profiles of nitrate concentrations against aircraft measurements below 400 hPa.

The warm Gulf Stream sea surface temperatures (SSTs) strongly impact the evolution of winter clouds behind atmospheric cold fronts. Such cloud evolution remains challenging to model. The Gulf Stream is too wide within the ERA5 and MERRA2 reanalyses, affecting the turbulent surface fluxes. Known problems within the ERA5 boundary layer (too-dry and too-cool with too strong westerlies), ascertained primarily from ACTIVATE 2020 campaign aircraft dropsondes and secondarily from older buoy measurements, reinforce surface flux biases. In contrast, MERRA2 winter surface winds and air-sea temperature/humidity differences are slightly too weak, producing surface fluxes that are too low. Reanalyses boundary layer heights in the strongly-forced winter cold-air-outbreak regime are realistic, whereas late-summer quiescent stable boundary layers are too shallow. Nevertheless, the reanalysis biases are small, and reanalyses adequately support their use for initializing higher-resolution cloud process modeling studies of cold-air outbreaks.