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).

The Energy Exascale Earth System Model Atmosphere Model version 1, the atmospheric component of the Department of Energy's Energy Exascale Earth System Model is described. The model began as a fork of the well‐known Community Atmosphere Model, but it has evolved in new ways, and coding, performance, resolution, physical processes (primarily cloud and aerosols formulations), testing and development procedures now differ significantly. Vertical resolution was increased (from 30 to 72 layers), and the model top extended to 60 km (~0.1 hPa). A simple ozone photochemistry predicts stratospheric ozone, and the model now supports increased and more realistic variability in the upper troposphere and stratosphere. An optional improved treatment of light‐absorbing particle deposition to snowpack and ice is available, and stronger connections with Earth system biogeochemistry can be used for some science problems. Satellite and ground‐based cloud and aerosol simulators were implemented to facilitate evaluation of clouds, aerosols, and aerosol‐cloud interactions. Higher horizontal and vertical resolution, increased complexity, and more predicted and transported variables have increased the model computational cost and changed the simulations considerably. These changes required development of alternate strategies for tuning and evaluation as it was not feasible to “brute force” tune the high‐resolution configurations, so short‐term hindcasts, perturbed parameter ensemble simulations, and regionally refined simulations provided guidance on tuning and parameterization sensitivity to higher resolution. A brief overview of the model and model climate is provided. Model fidelity has generally improved compared to its predecessors and the CMIP5 generation of climate models.

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.

Eastern China has experienced severe and persistent winter haze episodes in recent years due to intensification of aerosol pollution. In addition to anthropogenic emissions, the winter aerosol pollution over eastern China is associated with unusual meteorological conditions, including weaker wind speeds. Here we show, based on model simulations, that during years with decreased wind speed, large decreases in dust emissions (29%) moderate the wintertime land–sea surface air temperature difference and further decrease winds by −0.06 (±0.05) m s⁻¹ averaged over eastern China. The dust-induced lower winds enhance stagnation of air and account for about 13% of increasing aerosol concentrations over eastern China. Although recent increases in anthropogenic emissions are the main factor causing haze over eastern China, we conclude that natural emissions also exert a significant influence on the increases in wintertime aerosol concentrations, with important implications that need to be taken into account by air quality studies.

To mitigate climate warming, many countries have committed to achieve carbon neutrality in the mid-21^st century. Here, we assess the global impacts of changing greenhouse gases (GHGs), aerosols, and tropospheric ozone (O₃) following a carbon neutrality pathway on climate and extreme weather events individually using the Community Earth System Model version 1 (CESM1). The results suggest that the future aerosol reductions significantly contribute to climate warming and increase the frequency and intensity of extreme weathers toward carbon neutrality and aerosol impacts far outweigh those of GHGs and tropospheric O₃. It reverses the knowledge that the changing GHGs dominate the future climate changes as predicted in the middle of the road pathway. Therefore, substantial reductions in GHGs and tropospheric O₃ are necessary to reach the 1.5 °C warming target and mitigate the harmful effects of concomitant aerosol reductions on climate and extreme weather events under carbon neutrality in the future.

As a major alkaline gas in the atmosphere, NH₃ significantly impacts atmospheric chemistry, ecological environment, and biodiversity. Gridded NH₃ emission inventories can significantly affect the accuracy of model concentrations and play a crucial role in the refinement of mitigation strategies. However, several uncertainties are still associated with existing NH₃ emission inventories in China. Therefore, in this study, we focused on improving fertilizer-application-related NH₃ emission inventories. We comprehensively evaluated the dates and times of fertilizer application to the major crops that are cultivated in China, improved the spatial allocation methods for NH₃ emissions from croplands with different rice types, and established a gridded NH₃ emission inventory for mainland China with a resolution of 5 min × 5 min in 2016. The results showed that the atmospheric NH₃ emissions in mainland China amounted to 12.11 Tg, with livestock waste (44.8 %) and fertilizer application (38.6 %) being the two main NH₃ emission sources in China. Obvious spatial variability in NH₃ emissions was also identified, and high emissions were predominantly concentrated in North China. Further, NH₃ emissions tended to be high in summer and low in winter, and the ratio for the July–January period was 3.08. Furthermore, maize and rice fertilization in summer was primarily responsible for the increase in NH₃ emissions in China, and the evaluation of the spatial and temporal accuracy of the NH₃ emission inventory established in this study using the WRF-Chem and ground-station- and satellite-based observations showed that it was more accurate than other inventories.