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.

Sub-grid topographic heterogeneity has large impacts on surface energy balance and land-atmosphere interactions. However, the impacts of representing sub-grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub-grid topographic representations on surface energy balance, turbulent heat flux and scalar (co-)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub-grid topographic representations in ELM were compared: (1) the default sub-grid structure (D), (2) the recently developed sub-grid topographic structure (T), and (3) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (1) the default plane-parallel radiative transfer scheme (PP) and (2) the parameterization scheme (TOP) that accounts for sub-grid topographic effects on solar radiation. A series of simulations with the three grid structures (D, T and 1KM) and two treatments of solar radiation (TOP and PP) were carried out in the Sierra Nevada, California. There are significant differences between TOP and PP in the 1-km simulated surface energy balance, but the differences in the mean values and standard deviations become small when aggregated to the grid-scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations than the D configuration, and better captures the sub-grid topographic effects on surface energy balance as well as surface boundary conditions. These results underline the importance of representing sub-grid topographic heterogeneities in LSMs and motivate future research to understand the sub-grid topographic effects on land-atmosphere interactions over mountain areas.