This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). 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 with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, 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 Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. 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.

A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011–2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of −24 % and −35 % for particles with dry diameters >50 and >120 nm, as well as −36 % and −34 % for CCN at supersaturations of 0.2 % and 1.0 %, respectively. However, they seem to behave differently for particles activating at very low supersaturations (<0.1 %) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N₃ (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2 % (CCN_0.2) compared to that for N₃, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120 nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40 % during winter and 20 % in summer. In contrast to the large spread in simulated aerosol particle and CCN number concentrations, the CDNC derived from simulated CCN spectra is less diverse and in better agreement with CDNC estimates consistently derived from the observations (average NMB −13 % and −22 % for updraft velocities 0.3 and 0.6 m s⁻¹, respectively). In addition, simulated CDNC is in slightly better agreement with observationally derived values at lower than at higher updraft velocities (index of agreement 0.64 vs. 0.65). The reduced spread of CDNC compared to that of CCN is attributed to the sublinear response of CDNC to aerosol particle number variations and the negative correlation between the sensitivities of CDNC to aerosol particle number concentration ($\partial N_{\mathrm{d}}/\partial N_{\mathrm{a}}$ ) and to updraft velocity ($\partial N_{\mathrm{d}}/\partial w$ ). Overall, we find that while CCN is controlled by both aerosol particle number and composition, CDNC is sensitive to CCN at low and moderate CCN concentrations and to the updraft velocity when CCN levels are high. Discrepancies are found in sensitivities $\partial N_{\mathrm{d}}/\partial N_{\mathrm{a}}$ and $\partial N_{\mathrm{d}}/\partial w$ ; models may be predisposed to be too “aerosol sensitive” or “aerosol insensitive” in aerosol–cloud–climate interaction studies, even if they may capture average droplet numbers well. This is a subtle but profound finding that only the sensitivities can clearly reveal and may explain inter-model biases on the aerosol indirect effect.

A recent development in the representation of aerosols in climate models is the realization that some components of organic aerosol (OA), emitted from biomass and biofuel burning, can have a significant contribution to shortwave radiation absorption in the atmosphere. The absorbing fraction of OA is referred to as brown carbon (BrC). This study introduces one of the first implementations of BrC into the Community Atmosphere Model version 5 (CAM5), using a parameterization for BrC absorptivity described in Saleh et al. (2014). Nine-year experiments are run (2003–2011) with prescribed emissions and sea surface temperatures to analyze the effect of BrC in the atmosphere. Model validation is conducted via model comparison to single-scatter albedo and aerosol optical depth from the Aerosol Robotic Network (AERONET). This comparison reveals a model underestimation of single scattering albedo (SSA) in biomass burning regions for both default and BrC model runs, while a comparison between AERONET and the model absorption Ångström exponent shows a marked improvement with BrC implementation. Global annual average radiative effects are calculated due to aerosol–radiation interaction (REari; 0.13±0.01 W m⁻²) and aerosol–cloud interaction (REaci; 0.01±0.04 W m⁻²). REari is similar to other studies' estimations of BrC direct radiative effect, while REaci indicates a global reduction in low clouds due to the BrC semi-direct effect. The mechanisms for these physical changes are investigated and found to correspond with changes in global circulation patterns. Comparisons of BrC implementation approaches find that this implementation predicts a lower BrC REari in the Arctic regions than previous studies with CAM5. Implementation of BrC bleaching effect shows a significant reduction in REari (0.06±0.008 W m⁻²). Also, variations in OA density can lead to differences in REari and REaci, indicating the importance of specifying this property when estimating the BrC radiative effects and when comparing similar studies.

Herein, we present a description of the Mechanism of Intermediate complexity for Modelling Iron (MIMI v1.0). This iron processing module was developed for use within Earth system models and has been updated within a modal aerosol framework from the original implementation in a bulk aerosol model. MIMI simulates the emission and atmospheric processing of two main sources of iron in aerosol prior to deposition: mineral dust and combustion processes. Atmospheric dissolution of insoluble to soluble iron is parameterized by an acidic interstitial aerosol reaction and a separate in-cloud aerosol reaction scheme based on observations of enhanced aerosol iron solubility in the presence of oxalate. Updates include a more comprehensive treatment of combustion iron emissions, improvements to the iron dissolution scheme, and an improved physical dust mobilization scheme. An extensive dataset consisting predominantly of cruise-based observations was compiled to compare to the model. The annual mean modelled concentration of surface-level total iron compared well with observations but less so in the soluble fraction (iron solubility) for which observations are much more variable in space and time. Comparing model and observational data is sensitive to the definition of the average as well as the temporal and spatial range over which it is calculated. Through statistical analysis and examples, we show that a median or log-normal distribution is preferred when comparing with soluble iron observations. The iron solubility calculated at each model time step versus that calculated based on a ratio of the monthly mean values, which is routinely presented in aerosol studies and used in ocean biogeochemistry models, is on average globally one-third (34 %) higher. We redefined ocean deposition regions based on dominant iron emission sources and found that the daily variability in soluble iron simulated by MIMI was larger than that of previous model simulations. MIMI simulated a general increase in soluble iron deposition to Southern Hemisphere oceans by a factor of 2 to 4 compared with the previous version, which has implications for our understanding of the ocean biogeochemistry of these predominantly iron-limited ocean regions.

Dust particles, serving as ice-nucleating particles (INPs), may impact the Arctic surface energy budget and regional climate by modulating the mixed-phase cloud properties and lifetime. In addition to long-range transport from low-latitude deserts, dust particles in the Arctic can originate from local sources. However, the importance of high-latitude dust (HLD) as a source of Arctic INPs (compared to low-latitude dust, LLD) and its effects on Arctic mixed-phase clouds are overlooked. In this study, we evaluate the contribution to Arctic dust loading and INP population from HLD and six LLD source regions by implementing a source-tagging technique for dust aerosols in version 1 of the US Department of Energy's Energy Exascale Earth System Model (E3SMv1). Our results show that HLD is responsible for 30.7 % of the total dust burden in the Arctic, whereas LLD from Asia and North Africa contributes 44.2 % and 24.2 %, respectively. Due to its limited vertical transport as a result of stable boundary layers, HLD contributes more in the lower troposphere, especially in boreal summer and autumn when the HLD emissions are stronger. LLD from North Africa and East Asia dominates the dust loading in the upper troposphere with peak contributions in boreal spring and winter. The modeled INP concentrations show better agreement with both ground and aircraft INP measurements in the Arctic when including HLD INPs. The HLD INPs are found to induce a net cooling effect (−0.24 W m⁻² above 60^∘ N) on the Arctic surface downwelling radiative flux by changing the cloud phase of the Arctic mixed-phase clouds. The magnitude of this cooling is larger than that induced by North African and East Asian dust (0.08 and −0.06 W m⁻², respectively), mainly due to different seasonalities of HLD and LLD. Uncertainties of this study are discussed, which highlights the importance of further constraining the HLD emissions.

Dust aerosol is important in modulating the climate system at local and global scales, yet its spatiotemporal distributions simulated by global climate models (GCMs) are highly uncertain. In this study, we evaluate the spatiotemporal variations of dust extinction profiles and dust optical depth (DOD) simulated by the Community Earth System Model version 1 (CESM1) and version 2 (CESM2), the Energy Exascale Earth System Model version 1 (E3SMv1), and the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) against satellite retrievals from Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), Moderate Resolution Imaging Spectroradiometer (MODIS), and Multi-angle Imaging SpectroRadiometer (MISR). We find that CESM1, CESM2, and E3SMv1 underestimate dust transport to remote regions. E3SMv1 performs better than CESM1 and CESM2 in simulating dust transport and the northern hemispheric DOD due to its higher mass fraction of fine dust. CESM2 performs the worst in the Northern Hemisphere due to its lower dust emission than in the other two models but has a better dust simulation over the Southern Ocean due to the overestimation of dust emission in the Southern Hemisphere. DOD from MERRA-2 agrees well with CALIOP DOD in remote regions due to its higher mass fraction of fine dust and the assimilation of aerosol optical depth. The large disagreements in the dust extinction profiles and DOD among CALIOP, MODIS, and MISR retrievals make the model evaluation of dust spatial distributions challenging. Our study indicates the importance of representing dust emission, dry/wet deposition, and size distribution in GCMs in correctly simulating dust spatiotemporal distributions.