Explore the vast range of titles available.

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. Plain Language Summary The Energy Exascale Earth System Model (E3SM) is a relatively new fully coupled Earth system and climate model used in major international model simulation projects and mission‐defined efforts for the U.S. Department of Energy. This paper describes the first simulation of the model in its high‐resolution configuration. This higher‐resolution version is able to capture the most energetic motions in the ocean, which are poorly represented in standard resolution coupled climate models, as well as the largest of storms in the atmosphere. Evaluation of this simulation confirms the benefits of high resolution found by other models with a few notable exceptions. These discrepancies with other studies are interesting because they provide a richer understanding of how and why resolution affects model bias. Another key finding is that climate and aerosol sensitivity in E3SM is unaffected by resolution change. This affirms the usefulness of coarser‐resolution models for understanding global‐scale climate change. This study also confirms the benefits of increased resolution for studying fine‐scale features such as hurricanes and orographic precipitation. Finally, the high‐resolution version of E3SM is shown to compare favorably to its low‐resolution counterpart and to the models participating in Phase 5 of the Coupled Model Intercomparison Project. Key Points The high‐resolution E3SMv1 model was run for 50 years using 1950 forcing data according to the HighResMIP protocol Higher resolution and associated retuning improved bias relative to coarser versions of E3SMv1, particularly in ocean and sea ice metrics Aerosol and climate sensitivity were relatively unaffected by resolution change; resolution‐related tuning had a larger impact

We present the analysis of global sympagic primary production (PP) from 300 years of pre-industrial and historical simulations of the E3SMv1.1-BGC model. The model includes a novel, eight-element sea ice biogeochemical component, MPAS-Seaice zbgc, which is resolved in three spatial dimensions and uses a vertical transport scheme based on internal brine dynamics. Modeled ice algal chlorophyll-a concentrations and column-integrated values are broadly consistent with observations, though chl-a profile fractions indicate that upper ice communities of the Southern Ocean are underestimated. Simulations of polar integrated sea ice PP support the lower bound in published estimates for both polar regions with mean Arctic values of 7.5 and 15.5 TgC/a in the Southern Ocean. However, comparisons of the polar climate state with observations, using a maximal bound for ice algal growth rates, suggest that the Arctic lower bound is a significant underestimation driven by biases in ocean surface nitrate, and that correction of these biases supports as much as 60.7 TgC/a of net Arctic PP. Simulated Southern Ocean sympagic PP is predominantly light-limited, and regional patterns, particularly in the coastal high production band, are found to be negatively correlated with snow thickness.

The Energy Exascale Earth System Model (E3SM) is a new coupled Earth system model sponsored by the U.S Department of Energy. Here we present E3SM global simulations using active ocean and sea ice that are driven by the Coordinated Ocean‐ice Reference Experiments II (CORE‐II) interannual atmospheric forcing data set. The E3SM ocean and sea ice components are MPAS‐Ocean and MPAS‐Seaice, which use the Model for Prediction Across Scales (MPAS) framework and run on unstructured horizontal meshes. For this study, grid cells vary from 30 to 60 km for the low‐resolution mesh and 6 to 18 km at high resolution. The vertical grid is a structured z‐star coordinate and uses 60 and 80 layers for low and high resolution, respectively. The lower‐resolution simulation was run for five CORE cycles (310 years) with little drift in sea surface temperature (SST) or heat content. The meridional heat transport (MHT) is within observational range, while the meridional overturning circulation at 26.5°N is low compared to observations. The largest temperature biases occur in the Labrador Sea and western boundary currents (WBCs), and the mixed layer is deeper than observations at northern high latitudes in the winter months. In the Antarctic, maximum mixed layer depths (MLD) compare well with observations, but the spatial MLD pattern is shifted relative to observations. Sea ice extent, volume, and concentration agree well with observations. At high resolution, the sea surface height compares well with satellite observations in mean and variability. Key Points The Energy Exascale Earth System Model (E3SM) is a new climate model by the U.S. Department of Energy E3SM ocean and ice components use unstructured horizontal meshes for variable‐resolution simulations The 310‐year E3SM simulations agree well with observations in ocean currents and sea ice coverage

Isopycnal diffusivity due to stirring by mesoscale eddies in an idealized, wind-forced, eddying, midlatitude ocean basin is computed using Lagrangian, in Situ, Global, High-Performance Particle Tracking (LIGHT). Simulation is performed via LIGHT within the Model for Prediction across Scales Ocean (MPAS-O). Simulations are performed at 4-, 8-, 16-, and 32-km resolution, where the first Rossby radius of deformation (RRD) is approximately 30 km. Scalar and tensor diffusivities are estimated at each resolution based on 30 ensemble members using particle cluster statistics. Each ensemble member is composed of 303 665 particles distributed across five potential density surfaces. Diffusivity dependence upon model resolution, velocity spatial scale, and buoyancy surface is quantified and compared with mixing length theory. The spatial structure of diffusivity ranges over approximately two orders of magnitude with values of O (10 5 ) m 2 s −1 in the region of western boundary current separation to O (10 3 ) m 2 s −1 in the eastern region of the basin. Dominant mixing occurs at scales twice the size of the first RRD. Model resolution at scales finer than the RRD is necessary to obtain sufficient model fidelity at scales between one and four RRD to accurately represent mixing. Mixing length scaling with eddy kinetic energy and the Lagrangian time scale yield mixing efficiencies that typically range between 0.4 and 0.8. A reduced mixing length in the eastern region of the domain relative to the west suggests there are different mixing regimes outside the baroclinic jet region.

Emissions‐driven (prognostic CO2) simulations are essential for representing two‐way carbon‐climate feedback in Earth System Models. We present an emissions‐driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land‐atmosphere emission‐hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions‐driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land‐atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO2 concentrations by 11–23 ppm yet stays within the 40‐ppm spread CMIP6 emission‐driven models and retains physical climate properties comparable to the reference runs. The CO2 biases are partly attributed to underrepresented oceanic CO2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin‐up assessment, atmospheric CO2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter‐model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions‐driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon‐budget metrics for mitigation scenarios. Plain Language Summary Understanding the impact of carbon dioxide (CO2) emissions on climate is vital for predicting future changes and crafting effective policies. Earth System Models (ESMs) are essential tools for simulating Earth's climate and assessing various influencing factors. In this study, we extended the Energy Exascale Earth System Model (E3SM)'s capabilities so that CO2 levels are calculated directly from human and natural emissions instead of being prescribed as a single global value. This extension allows for a more realistic representation of CO2 exchange between the atmosphere and land. We conducted historical simulations from 1850 to 2014 using this new development and compared results with observations and other models. Our model slightly overestimates atmospheric CO2 levels compared to measurements but is comparable to other models in capturing key climate features. To help other researchers build and test similar “emission‐driven” models, we created a step‐by‐step evaluation framework that checks CO2 behavior, climate variables, and land‐atmosphere interactions. Our work advances E3SM modeling by accurately representing how CO2 emissions affect Earth's systems. This enhancement lays the groundwork for modeling interactions between human‐Earth interactions, thereby enabling future studies that can inform mitigation and adaption. Key Points Implemented emissions‐driven land–atmosphere biogeochemistry in E3SMv2.1 (BGCLNDATM_progCO2), enabling prognostic CO2 simulations Established a structured evaluation protocol ensuring scientific rigor and facilitating inter‐model comparisons of model performance Emissions‐driven BGCLNDATM_progCO2 simulations maintain a physical climate similar to reference runs with prescribed CO2 concentrations

Over the past two decades, sea level trends have increased in the western tropical Pacific Ocean with rates that are approximately three times the global average. A general circulation model is used to show that the high rates are caused by a gradual intensification of Pacific trade winds since the early 1990s. The modeled sea level change captures the spatial trend pattern in satellite altimeter sea surface heights and the temporal trend shift in tide gauge observations. In addition to the sea level response, the model is used to show how other aspects of the ocean circulation have increased appreciably in amplitude as a consequence of the trade wind intensification, including tropical surface currents, the shallow meridional over‐turning circulation, the Equatorial Undercurrent, and the Indonesian Throughflow. These results highlight an ongoing shift in the state of the tropical Pacific Ocean that will continue as long as the trade wind trend persists. Key Points Trade wind intensification drives regional sea level change The intensification leads to broad changes in the Pacific circulation The modeling study focuses attention on the need for wind product validation

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. Plain Language Summary The U.S. Department of Energy recently released version two of its Energy Exascale Earth System Model (E3SM). E3SMv2 experienced a significant evolution in many of its model components (most notably the atmosphere and sea ice models), and its supporting software infrastructure. In this work, we document the computational performance of E3SMv2 and analyze its ability to reproduce the observed climate. To accomplish this, we utilize the standard Diagnosis and Evaluation and Characterization of Klima experiments augmented with historical simulations for the period 1850–2015. We find that E3SMv2 is nearly twice as fast as its predecessor and more accurately reproduces the observed climate in a number of metrics, most notably clouds and precipitation. We also find that the model's simulated response to increasing carbon dioxide (the equilibrium climate sensitivity) is much more realistic. Unfortunately, E3SMv2 underestimates the global mean surface temperature compared to observations during the second half of historical period. Using sensitivity experiments, where forcing agents (carbon dioxide, aerosols) are selectively disabled in the model, we determine that correcting this problem would require a strong reduction in the impact of aerosols. Key Points E3SMv2 is nearly twice as fast as E3SMv1 with a simulated climate that is improved in many metrics (e.g., precipitation and clouds) Climate sensitivity is substantially lower with a more plausible equilibrium climate sensitivity of 4.0 K (compared to an unlikely value of 5.3 K in E3SMv1) E3SMv2 underestimates the warming in the late historical period due to excessive aerosol‐related forcing

A climatologically forced high-resolution model is used to examine variability of subtropical mode water (STMW) in the northwestern Pacific Ocean. Despite the use of annually repeating atmospheric forcing, significant interannual to decadal variability is evident in the volume, temperature, and age of STMW formed in the region. This long time-scale variability is intrinsic to the ocean. The formation and characteristics of STMW are comparable to those observed in nature. STMW is found to be cooler, denser, and shallower in the east than in the west, but time variations in these properties are generally correlated across the full water mass. Formation is found to occur south of the Kuroshio Extension, and after formation STMW is advected westward, as shown by the transport streamfunction. The ideal age and chlorofluorocarbon tracers are used to analyze the life cycle of STMW. Over the full model run, the average age of STMW is found to be 4.1 yr, but there is strong geographical variation in this, from an average age of 3.0 yr in the east to 4.9 yr in the west. This is further evidence that STMW is formed in the east and travels to the west. This is qualitatively confirmed through simulated dye experiments known as transit-time distributions. Changes in STMW formation are correlated with a large meander in the path of the Kuroshio south of Japan. In the model, the large meander inhibits STMW formation just south of Japan, but the export of water with low potential vorticity leads to formation of STMW in the east and an overall increase in volume. This is correlated with an increase in the outcrop area of STMW. Mixed layer depth, on the other hand, is found to be uncorrelated with the volume of STMW.

This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km.A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.