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

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.

This study provides an overview of the coupled high‐resolution Version 1 of the Energy Exascale Earth System Model (E3SMv1) and documents the characteristics of a 50‐year‐long high‐resolution control simulation with time‐invariant 1950 forcings following the HighResMIP protocol. In terms of global root‐mean‐squared error metrics, this high‐resolution simulation is generally superior to results from the low‐resolution configuration of E3SMv1 (due to resolution, tuning changes, and possibly initialization procedure) and compares favorably to models in the CMIP5 ensemble. Ocean and sea ice simulation is particularly improved, due to better resolution of bathymetry, the ability to capture more variability and extremes in winds and currents, and the ability to resolve mesoscale ocean eddies. The largest improvement in this regard is an ice‐free Labrador Sea, which is a major problem at low resolution. Interestingly, several features found to improve with resolution in previous studies are insensitive to resolution or even degrade in E3SMv1. Most notable in this regard are warm bias and associated stratocumulus deficiency in eastern subtropical oceans and lack of improvement in El Niño. Another major finding of this study is that resolution increase had negligible impact on climate sensitivity (measured by net feedback determined through uniform +4K prescribed sea surface temperature increase) and aerosol sensitivity. Cloud response to resolution increase consisted of very minor decrease at all levels. Large‐scale patterns of precipitation bias were also relatively unaffected by grid spacing.

This paper documents the biogeochemistry configuration of the Energy Exascale Earth System Model (E3SM), E3SMv1.1‐BGC. The model simulates historical carbon cycle dynamics, including carbon losses predicted in response to land use and land cover change, and the responses of the carbon cycle to changes in climate. In addition, we introduce several innovations in the treatment of soil nutrient limitation mechanisms, including explicit dependence on phosphorus availability. The suite of simulations described here includes E3SM contributions to the Coupled Climate‐Carbon Cycle Model Intercomparison Project and other projects, as well as simulations to explore the impacts of structural uncertainty in representations of nitrogen and phosphorus limitation. We describe the model spin‐up and evaluation procedures, provide an overview of results from the simulation campaign, and highlight key features of the simulations. Cumulative warming over the twentieth century is similar to observations, with a midcentury cold bias offset by stronger warming in recent decades. Ocean biomass production and carbon uptake are underpredicted, likely due to biases in ocean transport leading to widespread anoxia and undersupply of nutrients to surface waters. The inclusion of nutrient limitations in the land biogeochemistry results in weaker carbon fertilization and carbon‐climate feedbacks than exhibited by other Earth System Models that exclude those limitations. Finally, we compare with an alternative representation of terrestrial biogeochemistry, which differs in structure and in initialization of soil phosphorus. While both configurations agree well with observational benchmarks, they differ significantly in their distribution of carbon among different pools and in the strength of nutrient limitations.

Man-made reservoirs play a key role in the terrestrial water system. They alter water fluxes at the land surface and impact surface water storage through water management regulations for diverse purposes such as irrigation, municipal water supply, hydropower generation, and flood control. Although most developed countries have established sophisticated observing systems for many variables in the land surface water cycle, long-term and consistent records of reservoir storage are much more limited and not always shared. Furthermore, most land surface hydrological models do not represent the effects of water management activities. Here, the contribution of reservoirs to seasonal water storage variations is investigated using a large-scale water management model to simulate the effects of reservoir management at basin and continental scales. The model was run from 1948 to 2010 at a spatial resolution of 0.25° latitude–longitude. A total of 166 of the largest reservoirs in the world with a total capacity of about 3900 km3 (nearly 60% of the globally integrated reservoir capacity) were simulated. The global reservoir storage time series reflects the massive expansion of global reservoir capacity; over 30 000 reservoirs have been constructed during the past half century, with a mean absolute interannual storage variation of 89 km3. The results indicate that the average reservoir-induced seasonal storage variation is nearly 700 km3 or about 10% of the global reservoir storage. For some river basins, such as the Yellow River, seasonal reservoir storage variations can be as large as 72% of combined snow water equivalent and soil moisture storage.

Hydrological drought is not only caused by natural hydroclimate variability but can also be directly altered by human interventions including reservoir operation, irrigation, groundwater exploitation, etc. Understanding and forecasting of hydrological drought in the Anthropocene are grand challenges due to complicated interactions among climate, hydrology and humans. In this paper, five decades (1961–2010) of naturalized and observed streamflow datasets are used to investigate hydrological drought characteristics in a heavily managed river basin, the Yellow River basin in north China. Human interventions decrease the correlation between hydrological and meteorological droughts, and make the hydrological drought respond to longer timescales of meteorological drought. Due to large water consumptions in the middle and lower reaches, there are 118–262 % increases in the hydrological drought frequency, up to 8-fold increases in the drought severity, 21–99 % increases in the drought duration and the drought onset is earlier. The non-stationarity due to anthropogenic climate change and human water use basically decreases the correlation between meteorological and hydrological droughts and reduces the effect of human interventions on hydrological drought frequency while increasing the effect on drought duration and severity. A set of 29-year (1982–2010) hindcasts from an established seasonal hydrological forecasting system are used to assess the forecast skill of hydrological drought. In the naturalized condition, the climate-model-based approach outperforms the climatology method in predicting the 2001 severe hydrological drought event. Based on the 29-year hindcasts, the former method has a Brier skill score of 11–26 % against the latter for the probabilistic hydrological drought forecasting. In the Anthropocene, the skill for both approaches increases due to the dominant influence of human interventions that have been implicitly incorporated by the hydrological post-processing, while the difference between the two predictions decreases. This suggests that human interventions can outweigh the climate variability for the hydrological drought forecasting in the Anthropocene, and the predictability for human interventions needs more attention.