An overview of the Community Earth System Model Version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long preindustrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low‐top” (40 km, with limited chemistry) and “high‐top” (130 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low‐latitude precipitation and shortwave cloud forcing biases; better representation of the Madden‐Julian Oscillation; better El Niño‐Southern Oscillation‐related teleconnections; and a global land carbon accumulation trend that agrees well with observationally based estimates. Most tropospheric and surface features of the low‐ and high‐top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1–5.3 °C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2's transient climate response of 1.9–2.0 °C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and they represent CESM2's contributions to the Coupled Model Intercomparison Project Phase 6.

Anthropogenic emissions and land use changes have modified atmospheric aerosol concentrations and size distributions over time. Understanding preindustrial conditions and changes in organic aerosol due to anthropogenic activities is important because these features (1) influence estimates of aerosol radiative forcing and (2) can confound estimates of the historical response of climate to increases in greenhouse gases. Secondary organic aerosol (SOA), formed in the atmosphere by oxidation of organic gases, represents a major fraction of global submicron‐sized atmospheric organic aerosol. Over the past decade, significant advances in understanding SOA properties and formation mechanisms have occurred through measurements, yet current climate models typically do not comprehensively include all important processes. This review summarizes some of the important developments during the past decade in understanding SOA formation. We highlight the importance of some processes that influence the growth of SOA particles to sizes relevant for clouds and radiative forcing, including formation of extremely low volatility organics in the gas phase, acid‐catalyzed multiphase chemistry of isoprene epoxydiols, particle‐phase oligomerization, and physical properties such as volatility and viscosity. Several SOA processes highlighted in this review are complex and interdependent and have nonlinear effects on the properties, formation, and evolution of SOA. Current global models neglect this complexity and nonlinearity and thus are less likely to accurately predict the climate forcing of SOA and project future climate sensitivity to greenhouse gases. Efforts are also needed to rank the most influential processes and nonlinear process‐related interactions, so that these processes can be accurately represented in atmospheric chemistry‐climate models.

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

Atmospheric carbonaceous aerosols play an important role in the climate system by influencing the Earth's radiation budgets and modifying the cloud properties. Despite the importance, their representations in large-scale atmospheric models are still crude, which can influence model simulated burden, lifetime, physical, chemical and optical properties, and the climate forcing of carbonaceous aerosols. In this study, we improve the current three-mode version of the Modal Aerosol Module (MAM3) in the Community Atmosphere Model version 5 (CAM5) by introducing an additional primary carbon mode to explicitly account for the microphysical ageing of primary carbonaceous aerosols in the atmosphere. Compared to MAM3, the four-mode version of MAM (MAM4) significantly increases the column burdens of primary particulate organic matter (POM) and black carbon (BC) by up to 40 % in many remote regions, where in-cloud scavenging plays an important role in determining the aerosol concentrations. Differences in the column burdens for other types of aerosol (e.g., sulfate, secondary organic aerosols, mineral dust, sea salt) are less than 1 %. Evaluating the MAM4 simulation against in situ surface and aircraft observations, we find that MAM4 significantly improves the simulation of seasonal variation of near-surface BC concentrations in the polar regions, by increasing the BC concentrations in all seasons and particularly in cold seasons. However, it exacerbates the overestimation of modeled BC concentrations in the upper troposphere in the Pacific regions. The comparisons suggest that, to address the remaining model POM and BC biases, future improvements are required related to (1) in-cloud scavenging and vertical transport in convective clouds and (2) emissions of anthropogenic and biomass burning aerosols.

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.

The presence of a large fraction of organic matter in primary sea spray aerosol (SSA) can strongly affect its cloud condensation nuclei activity and interactions with marine clouds. Global climate models require new parameterizations of the SSA composition in order to improve the representation of these processes. Existing proposals for such a parameterization use remotely sensed chlorophyll a concentrations as a proxy for the biogenic contribution to the aerosol. However, both observations and theoretical considerations suggest that existing relationships with chlorophyll a, derived from observations at only a few locations, may not be representative for all ocean regions. We introduce a novel framework for parameterizing the fractionation of marine organic matter into SSA based on a competitive Langmuir adsorption equilibrium at bubble surfaces. Marine organic matter is partitioned into classes with differing molecular weights, surface excesses, and Langmuir adsorption parameters. The classes include a lipid-like mixture associated with labile dissolved organic carbon (DOC), a polysaccharide-like mixture associated primarily with semilabile DOC, a protein-like mixture with concentrations intermediate between lipids and polysaccharides, a processed mixture associated with recalcitrant surface DOC, and a deep abyssal humic-like mixture. Box model calculations have been performed for several cases of organic adsorption to illustrate the underlying concepts. We then apply the framework to output from a global marine biogeochemistry model, by partitioning total dissolved organic carbon into several classes of macromolecules. Each class is represented by model compounds with physical and chemical properties based on existing laboratory data. This allows us to globally map the predicted organic mass fraction of the nascent submicron sea spray aerosol. Predicted relationships between chlorophyll a and organic fraction are similar to existing empirical parameterizations, but can vary between biologically productive and nonproductive regions, and seasonally within a given region. Major uncertainties include the bubble film thickness at bursting, and the variability of organic surfactant activity in the ocean, which is poorly constrained. In addition, polysaccharides may enter the aerosol more efficiently than Langmuir adsorption would suggest. Potential mechanisms for enrichment of polysaccharides in sea spray include the formation of marine colloidal particles that may be more efficiently swept up by rising bubbles, and cooperative adsorption of polysaccharides with proteins or lipids. These processes may make important contributions to the aerosol, but are not included here. This organic fractionation framework is an initial step towards a closer linking of ocean biogeochemistry and aerosol chemical composition in Earth system models. Future work should focus on improving constraints on model parameters through new laboratory experiments or through empirical fitting to observed relationships in the real ocean and atmosphere, as well as on atmospheric implications of the variable composition of organic matter in sea spray.