This work documents the first version of the U.S. Department of Energy (DOE) new EnergyExascale 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 2 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 Key Points: • This work documents E3SMv1, the first version of the U.S. DOE Energy Exascale Earth System Model • The performance of E3SMv1 is documented with a set of standard CMIP6 DECK and historical simulations comprising nearly 3,000 years • E3SMv1 has a high equilibrium climate sensitivity (5.3 K) and strong aerosol-related effective radiative forcing (-1.65 W/m 2 ) Correspondence to: Chris Golaz, golaz1@llnl.gov Citation: Golaz, J.-C., Caldwell, P. M., Van Roekel, L. P., Petersen, M. R., Tang, Q., Wolfe, J. D., et al. (2019). The DOE E3SM coupled model version 1: Overview and evaluation at standard resolution. 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 2 ) and high equilibrium climate sensitivity (ECS = 5.3 K). Plain Language Summary The U.S. Department of Energy funded the development of a new state-of-the-art Earth system model for research and applications relevant to its mission. The Energy Exascale Earth System Model version 1 (E3SMv1) consists of five interacting components for the global atmosphere, land surface, ocean, sea ice, and rivers. Three of these components (ocean, sea ice, and river) are new and have not been coupled into an Earth system model previously. The atmosphere and land surface components were created by extending existing components part of the Community Earth System Model, Version 1. E3SMv1's capabilities are demonstrated by performing a set of standardized simulation experiments described by...

The Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25°, with the meteorology nudged toward the Year Of Tropical Convection analysis, and cloud simulators and the collocated A‐Train satellite observations are used to explore the resolution dependence of aerosol‐cloud interactions. The higher‐resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those of droplet nucleation and precipitation parameterizations. The annual average AIF in the Northern Hemisphere midlatitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m−2 (−30%) compared to the 2° model, leading to a 0.26 W m−2 reduction (−15%) in the global annual average AIF.

In the Amazon basin, particles containing mixed sodium salts are routinely observed and are attributed to marine aerosols transported from the Atlantic Ocean. Using chemical imaging analysis, we show that, during the wet season, fungal spores emitted by the forest biosphere contribute at least 30% (by number) to sodium salt particles in the central Amazon basin. Hydration experiments indicate that sodium content in fungal spores governs their growth factors. Modeling results suggest that fungal spores account for ~69% (31–95%) of the total sodium mass during the wet season and that their fractional contribution increases during nighttime. Contrary to common assumptions that sodium-containing aerosols originate primarily from marine sources, our results suggest that locally-emitted fungal spores contribute substantially to the number and mass of coarse particles containing sodium. Hence, their role in cloud formation and contribution to salt cycles and the terrestrial ecosystem in the Amazon basin warrant further consideration.

With increasing computational capabilities, cumulus parameterizations that are adaptable to the smaller grid spacing and temporal interval for high‐resolution climate model simulations are needed. In this study, we propose a method to improve the resolution adaptability of the Zhang‐McFarlane (ZM) scheme, by implementing spatial and temporal averaging to the CAPE tendency. This method allows for a more consistent application of the quasi‐equilibrium (QE) hypothesis at high spatial and temporal resolutions. The resolution adaptability of the original ZM scheme, the scheme with spatial averaging, and the scheme with spatiotemporal averaging at 4–32 km grid spacings are assessed using the Weather Research and Forecasting (WRF) model by comparing to cloud resolving model (CRM) simulation results coarse‐grained to these same grid spacings. We show the original ZM scheme has poor resolution adaptability, with spatiotemporally averaged subgrid convective transport and convective precipitation increasing significantly as the resolution increases. The spatial averaging method improves the resolution adaptability of the ZM scheme and better conserves total transport and total precipitation. Temporal averaging further improves the resolution adaptability of the scheme. With better constrained (although smoothed) convective transport and precipitation, both the spatial distribution and time series of total precipitation at 4 and 8 km grid spacings are improved with the averaging methods. The results could help develop resolution adaptability for other cumulus parameterizations that are based on the QE assumption.

Climate models with variable-resolution grids offer a computationally less expensive way to provide more detailed information and increased accuracy by resolving processes that cannot be adequately represented by a coarser grid. This study uses the Model for Prediction Across ScalesAtmosphere (MPAS-A), consisting of a nonhydrostatic dynamical core and a subset of Weather Research and Forecasting (WRF) model physics, to investigate the potential benefits of using tropical channel refinement. The simulations are performed with an idealized aquaplanet configuration using 30 and 240 km global grid spacing, and two variable-resolution grids spanning the same grid spacing range; one with a narrow (208S-208N) and one with a wide (308S-308N) tropical channel refinement. Increasing resolution in the tropics impacts both the tropical and extratropical circulation. Compared to the 30 km global grid, both refined channel simulations exhibit slightly stronger updrafts inside the Hadley cell resulting in more resolved precipitation. Using a wider tropical refinement leads to a closer correspondence with the global high-resolution grid. While different grid spacings produce similar cloud size distributions that are consistent with observations, the dependence of precipitation rate on cloud size varies among simulations. The refined channel simulations show improved tropical and extratropical precipitation relative to the global coarse simulation. All simulations show a single precipitation peak centered on the equator. Although the results show that tropical refinement is an effective method for avoiding artifacts due to grid resolution sensitivities seen in earlier studies that only refined a portion of the tropics, some biases remain well inside of the refinement region.