Explore the vast range of titles available.

The interactions between boundary layer turbulence, including Langmuir turbulence, and submesoscale processes in the oceanic mixed layer are described using large-eddy simulations of the spindown of a temperature front in the presence of submesoscale eddies, winds, and waves. The simulations solve the surface-wave-averaged Boussinesq equations with Stokes drift wave forcing at a resolution that is sufficiently fine to capture small-scale Langmuir turbulence. A simulation without Stokes drift forcing is also performed for comparison. Spatial and spectral properties of temperature, velocity, and vorticity fields are described, and these fields are scale decomposed in order to examine multiscale fluxes of momentum and buoyancy. Buoyancy flux results indicate that Langmuir turbulence counters the restratifying effects of submesoscale eddies, leading to small-scale vertical transport and mixing that is 4 times greater than in the simulations without Stokes drift forcing. The observed fluxes are also shown to be in good agreement with results from an asymptotic analysis of the surface-wave-averaged, or Craik–Leibovich, equations. Regions of potential instability in the flow are identified using Richardson and Rossby numbers, and it is found that mixed gravitational/symmetric instabilities are nearly twice as prevalent when Langmuir turbulence is present, in contrast to simulations without Stokes drift forcing, which are dominated by symmetric instabilities. Mixed layer depth calculations based on potential vorticity and temperature show that the mixed layer is up to 2 times deeper in the presence of Langmuir turbulence. Differences between measures of the mixed layer depth based on potential vorticity and temperature are smaller in the simulations with Stokes drift forcing, indicating a reduced incidence of symmetric instabilities in the presence of Langmuir turbulence.

Ocean fronts are an important submesoscale feature, yet frontogenesis theory often neglects turbulence – even parameterized turbulence – leaving theory lacking in comparison with observations and models. A perturbation analysis is used to include the effects of eddy viscosity and diffusivity as a first-order correction to existing strain-induced inviscid, adiabatic frontogenesis theory. A modified solution is obtained by using potential vorticity and surface conditions to quantify turbulent fluxes. It is found that horizontal viscosity and diffusivity tend to be readily frontolytic – reducing frontal tendency to negative values under weakly non-conservative perturbations and opposing or reversing front sharpening, whereas vertical viscosity and diffusivity tend to only weaken frontogenesis by slowing the rate of sharpening of the front even under strong perturbations. During late frontogenesis, vertical diffusivity enhances the rate of frontogenesis, although perturbation theory may be inaccurate at this stage. Surface quasi-geostrophic theory – neglecting all injected interior potential vorticity – is able to describe the first-order correction to the along-front velocity and ageostrophic overturning circulation in most cases. Furthermore, local conditions near the front maximum are sufficient to reconstruct the modified solution of both these fields.

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

Wind-wave processes have generally been excluded from coupled Earth system models due to the high computational expense of spectral wave models, which resolve a frequency and direction spectrum of waves across space and time. Existing uniform-resolution wave modeling approaches used in Earth system models cannot appropriately represent wave climates from global to coastal ocean scales, largely because of tradeoffs between coastal resolution and computational costs. To resolve this challenge, we introduce a global unstructured mesh capability for the WAVEWATCH III (WW3) model that is suitable for coupling within the US Department of Energy's Energy Exascale Earth System Model (E3SM). The new unstructured WW3 global wave modeling approach can provide the accuracy of higher global resolutions in coastal areas at the relative cost of lower uniform global resolutions. This new capability enables simulation of waves at physically relevant scales as needed for coastal applications.

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

This paper introduces Omega, the Ocean Model for E3SM Global Applications. Omega is a new ocean model designed to run efficiently on high performance computing (HPC) platforms, including exascale heterogeneous architectures with accelerators, such as Graphics Processing Units (GPUs). Omega is written in C++ and uses the Kokkos performance portability library. These were chosen because they are well-supported and will help future-proof Omega for upcoming HPC architectures. Omega will eventually replace the Model for Prediction Across Scales-Ocean (MPAS-Ocean) in the US Department of Energy's (DOE's) Energy Exascale Earth System Model (E3SM). Omega runs on unstructured horizontal meshes with variable-resolution capability and implements the same horizontal discretization as MPAS-Ocean. This work documents the design and performance of Omega Version 0.1.0 (Omega-V0), which solves the shallow water equations with passive tracers and is the first step towards the full primitive equation ocean model. On Central Processing Units (CPUs), Omega-V0 is 1.4 times faster than MPAS-Ocean with the same configuration. Omega-V0 is more efficient on GPUs than CPUs on a per-watt basis – by a factor of 5.3 on Frontier and 3.6 on Aurora, two of the world's fastest exascale computers.

The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS), ice‐shelf basal melting and iceberg calving, are generally poorly represented in current Earth System Models (ESMs). Here we document the cryosphere configuration of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) v1.2. This includes simulating Antarctic ice‐shelf basal melting, which has been implemented through simulating the ocean circulation within static Antarctic ice‐shelf cavities, allowing for the ability to calculate ice‐shelf basal melt rates from the associated heat and freshwater fluxes. In addition, we added the capability to prescribe forcing from iceberg melt, allowing for realistic representation of the other dominant mass loss process from the AIS. In standard resolution simulations (using a noneddying ocean) under preindustrial climate forcing, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a transition to a high basal melt regime under the Filchner‐Ronne Ice Shelf (FRIS), presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy‐induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic ice‐shelf basal melt rates across the continent that are generally within the range inferred from observations. The accurate representation of ice‐shelf basal melting within a coupled ESM is an important step toward reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea‐level rise. Plain Language Summary The future of the Antarctic Ice Sheet (AIS) has the potential to have broad impacts on global climate, perhaps most notably in contributing to sea‐level rise. The current generation of Earth System Models (ESMs) do not accurately represent the two primary means in which ice is lost from the AIS, through melting at the base of ice shelves floating on the ocean and the calving of icebergs. This limits our ability to make climate projections that incorporate the impacts of the AIS in a changing climate. Here, we demonstrate a novel capability to model one of those processes, ice‐shelf basal melting, in an ESM. We demonstrate the ability to simulate ice‐shelf basal melt rates across many Antarctic ice shelves that are in line with present day observations. We also find that, for certain ice shelves, modeled ice‐shelf basal melting can experience a rapid transition to high melting far above present‐day estimates, and this simulated high melting can be mitigated through improved ocean physics. Key Points Capabilities have been added to an Earth System Model to model realistic Antarctic ice‐shelf basal melt fluxes and prescribe iceberg forcing Simulated basal melt rates have a strong sensitivity to the ocean mesoscale eddy parameterization For one choice of the mesoscale eddy parameterization, the Filchner‐Ronne Ice Shelf transitions to a high melt regime

In this study, we analyze the realism with which tropical cyclones (TCs) are simulated in the fully coupled low‐ and high‐resolution Energy Exascale Earth System Model (E3SM) version 1, with a focus on the latter. Compared to the low‐resolution (grid spacing of ∼1°), the representation of TCs improves considerably in the high‐resolution configuration (grid spacing of ∼0.25°). Significant improvements are found in the global TC frequency, TC lifetime maximum intensities, and the relative distribution of TCs among the different basins. However, at both resolutions, spurious TC activity is found in some basins, notably in the subtropical regions. Contrasting the simulated large‐scale TC environment with observations reveals that the model environment is unrealistically conducive for TC development in those regions. Further analysis indicates that these biases are likely related to those in thermodynamic potential intensity, caused by systematic SST biases, and vertical wind shear in the coupled model. TC‐ocean interaction is also examined in the high‐resolution configuration of the model. The salient features of the ocean's response to TC‐induced mixing and the ocean's impact on TC intensification are well‐reproduced. Finally, an evaluation of the influence of El Niño Southern Oscillation (ENSO) on TCs in the high‐resolution configuration of the model reveals that the ENSO‐TC relationship in the model has the right sign and is significant for the North Atlantic and Northwest Pacific, albeit weaker than in observations. In summary, the high‐resolution configuration of the E3SM model simulates TC activity reasonably and hence could be a useful tool for TC‐related research. Plain Language Summary Considering the wide range of spatiotemporal scales governing tropical cyclones (TCs), high‐resolution numerical models that can resolve both small‐scale and large‐scale processes are needed to simulate TCs realistically. In this study, we describe the simulation of TCs in version 1 of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) model with a focus on results from the high‐resolution configuration (spatial resolution of 0.25° in the atmosphere). When compared to the simulation at the low or standard resolution configuration (spatial resolution of 1° in the atmosphere), the salient features of TCs, such as their frequency distribution among the various basins and their strength, improve considerably in the high‐resolution configuration, underlining the value of high‐resolution modeling. Further, the high‐resolution configuration of the model realistically simulates the interaction between TCs and the ocean beneath them, which is a fundamental aspect of their development. However, our analysis also reveals that systematic errors in the simulated climate can cause misrepresentation of the model TC environment, leading to spurious TC activity in certain regions and a weaker control exerted by El Niño‐Southern Oscillation. Overall, the high‐resolution E3SM well represents the global TC climatology, supporting it as a valuable tool for studying TCs. Key Points Under a variety of metrics, simulated TCs are significantly improved in high‐resolution versus low‐resolution E3SM configurations Spurious TC activity in the model can be well‐explained by errors in the simulation of sea surface temperatures and the vertical wind shear TC‐induced oceanic cooling and the ocean's negative feedback on TCs are well‐simulated in the high‐resolution configuration of the model

Current submesoscale restratification parameterizations, which help set mixed layer depth in global climate models, depend on a simplistic scaling of frontal width shown to be unreliable in several circumstances. Observations and theory indicate that frontogenesis is common, but stable frontal widths arise in the presence of turbulence and instabilities that participate in keeping fronts at the scale observed, the arrested scale. Here we propose a new scaling law for arrested frontal width as a function of turbulent fluxes via the turbulent thermal wind (TTW) balance. A variety of large-eddy simulations (LES) of strain-induced fronts and TTW-induced filaments are used to evaluate this scaling. Frontal width given by boundary layer parameters drawn from observations in the General Ocean Turbulence Model (GOTM) are found qualitatively consistent with the observed range in regions of active submesoscales. The new arrested front scaling is used to modify the mixed layer eddy restratification parameterization commonly used in coarse-resolution climate models. Results in CESM-POP2 reveal the climate model’s sensitivity to the parameterization update and changes in model biases. A comprehensive multimodel study is in planning for further testing.

We evaluate the simulated teleconnection of El Ni˜no—Southern Oscillation (ENSO) to winter season precipitation extremes over the United States in a long (98 years) 1950 control high-resolution version (HR; 25-km nominal atmosphere model horizontal resolution) of the U.S. Department of Energy’s (DOE) Energy Exascale Earth System Model version 1 (E3SMv1). The model bias and spatial pattern of ENSO teleconnections to mean and extreme precipitation in HR overall are similar to the low-resolution model’s (LR; 110 km) historical simulation (four-member ensemble, 1925–59). However, over the southeastern United States (SE-U.S.), HR produces stronger El Ni˜no—associated extremes, reducing LR’s model bias. Both LR and HR produce weaker than observed increase in storm track activity during El Ni˜no events there, but HR improves the ENSO-associated variability of moisture transport over SE-U.S. During El Ni˜no, stronger vertical velocities in HR produce stronger large-scale precipitation, causing larger latent heating of the troposphere that pulls in more moisture from the Gulf of Mexico into the SE-U.S. This positive feedback also contributes to the stronger mean and extreme precipitation response in HR. Over the Pacific Northwest, LR’s bias of stronger than observed La Ni˜na associated extremes is amplified in HR. Both models simulate stronger than observed moisture transport from the Pacific Ocean into the region during La Ni˜na years. The amplified HR bias there is due to stronger orographically driven vertical updrafts that create stronger large-scale precipitation, despite weaker La Ni˜na—induced storm track activity.