Energy Exascale Earth System Model Research Articles

Simulation of ENSO Teleconnections to Precipitation Extremes over the United States in the High-Resolution Version of E3SM

Abstract We evaluate the simulated teleconnection of El Niño–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ño–associated extremes, reducing LR’s model bias. Both LR and HR produce weaker than observed increase in storm track activity during El Niño events there, but HR improves the ENSO-associated variability of moisture transport over SE-U.S. During El Niño, 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ña associated extremes is amplified in HR. Both models simulate stronger than observed moisture transport from the Pacific Ocean into the region during La Niña years. The amplified HR bias there is due to stronger orographically driven vertical updrafts that create stronger large-scale precipitation, despite weaker La Niña–induced storm track activity. Significance Statement New high-resolution Earth system models (ESMs) solve mathematical equations of fluid flow at much smaller spatial scales than prevalent ESMs, and thus are prohibitively expensive to compute. However, they can be useful for simulating accurate details of regional climate extremes that are driven by naturally occurring climate oscillations like El Niño–Southern Oscillation (ENSO). Here, we evaluate the simulation of ENSO-driven precipitation extremes over the United States in the high-resolution version of the U.S. Department of Energy’s new Energy Exascale Earth System Model version 1. We find that the high-resolution model improves upon its low-resolution counterpart over the southeastern United States by producing a better transport of moisture into the region from the Gulf of Mexico during El Niño. Over the U.S. Pacific Northwest, the high-resolution model simulates the atmospheric flow in more detail over the complex mountainous terrain. However, it also brings in more moisture from the Pacific Ocean just like the low-resolution model. This causes it to produce precipitation extremes during La Niña years there that are stronger than that observed in the real world.

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 1: assessing E3SM aerosol predictions using aircraft, ship, and surface measurements

Abstract. An Earth system model (ESM) aerosol–cloud diagnostics package is developed to facilitate the routine evaluation of aerosols, clouds, and aerosol–cloud interactions simulated by the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). The first version focuses on comparing simulated aerosol properties with aircraft, ship, and surface measurements, which are mostly measured in situ. The diagnostics currently cover six field campaigns in four geographical regions: eastern North Atlantic (ENA), central US (CUS), northeastern Pacific (NEP), and Southern Ocean (SO). These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the Atmospheric Radiation Measurement (ARM) program and other agencies. Various types of diagnostics and metrics are performed for aerosol number, size distribution, chemical composition, cloud condensation nuclei (CCN) concentration, and various meteorological quantities to assess how well E3SM represents observed aerosol properties across spatial scales. Overall, E3SM qualitatively reproduces the observed aerosol number concentration, size distribution, and chemical composition reasonably well, but it overestimates Aitken-mode aerosols and underestimates accumulation-mode aerosols over the CUS and ENA regions, suggesting that processes related to particle growth or coagulation might be too weak in the model. The current version of E3SM struggles to reproduce the new particle formation events frequently observed over both the CUS and ENA regions, indicating missing processes in current parameterizations. The diagnostics package is coded and organized in a way that can be extended to other field campaign datasets and adapted to higher-resolution model simulations.

Diurnal Rainfall Response to the Physiological and Radiative Effects of CO2 in Tropical Forests in the Energy Exascale Earth System Model v1

AbstractUnderstanding how the connection between rainfall and tropical forests will respond to increasing CO2 concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO2. The plant physiological and radiative impacts of CO2 on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1‐BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO2 levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO2 radiative and CO2 plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime Convective Available Potential Energy (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response.

Assessing the Potential Impact of River Chemistry on Arctic Coastal Production

The Arctic coastal margin receives a disproportionately large fraction of the global river discharge. The bio-geochemistry of the river water as it empties into the marine environment reflects inputs and processes that occur as the water travels from its headwaters. Climate-induced changes to Arctic vegetation and permafrost melt may impact river chemistry. Understanding the impact of river nutrients on coastal marine production, and how this may change in the future, are important for resource managers and community members who monitor and rely on coastal food resources. Using the Energy Exascale Earth System Model we explore the impact of timing and river nutrient concentrations on primary production in each coastal Arctic region and then assess how this influences secondary production and particle fluxes supporting the benthic food web. Our results indicate that while the concentration of Arctic river nitrogen can have a significant impact on annual average nitrogen and primary production in the coastal Arctic, with production increases of up to 20% in the river influenced interior Seas, the timing of the river nutrient inputs into the marine environment appears less important. Bloom timing and partitioning between small and large phytoplankton were minimally impacted by both river nutrient concentration and timing, suggesting that in general, coastal Arctic ecosystem dynamics will continue to be primarily driven by light availability, rather than nutrients. Under a doubling river nutrient scenario, the percentage increase in the POC flux to the benthos on river influenced Arctic coastal shelves was 2-4 times the percentage increase in primary production, suggesting changes to the river nutrient concentration has the potential to modify the Arctic food web structure and dynamics. Generally, the nutrient-induced changes to primary production were smaller than changes previously simulated in response to ice reduction and temperature increase. However, in the Laptev Sea, the production increase resulting from a doubling of river nutrients exceeded the production increase simulated with an atmospheric warming scenario. Dissolved organic carbon is presently poorly represented in the model so its impact on production is hard to simulate. Applying established relationships between modeled DOC, total DOC, and light absorption we illustrate that DOC could play a very important role in modulating production. Our findings highlight the importance of developing more realistic river nutrient and discharge forcing for Earth System Models such that their impact on the critical Arctic coastal domain can be more adequately resolved.

Description of historical and future projection simulations by the global coupled E3SMv1.0 model as used in CMIP6

Abstract. This paper documents the experimental setup and general features of the coupled historical and future climate simulations with the first version of the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv1.0). The future projected climate characteristics of E3SMv1.0 at the highest emission scenario (SSP5-8.5) designed in the Scenario Model Intercomparison Project (ScenarioMIP) and the SSP5-8.5 greenhouse gas (GHG) only forcing experiment are analyzed with a focus on regional responses of atmosphere, ocean, sea ice, and land. Due to its high equilibrium climate sensitivity (ECS of 5.3 K), E3SMv1.0 is one of the Coupled Model Intercomparison Project phase 6 (CMIP6) models with the largest surface warming by the end of the 21st century under the high-emission SSP5-8.5 scenario. The global mean precipitation change is highly correlated with the global temperature change, while the spatial pattern of the change in runoff is consistent with the precipitation changes. The oceanic mixed layer generally shoals throughout the global ocean. The annual mean Atlantic meridional overturning circulation (AMOC) is overly weak with a slower change from ∼ 11 to ∼ 6 Sv (Sverdrup) relative to other CMIP6 models. The sea ice, especially in the Northern Hemisphere, decreases rapidly with large seasonal variability. We detect a significant polar amplification in E3SMv1.0 from the atmosphere, ocean, and sea ice. Comparing the SSP5-8.5 all-forcing experiment with the GHG-only experiment, we find that the unmasking of the aerosol effects due to the decline of the aerosol loading in the future projection period causes transient accelerated warming in the all-forcing experiment in the first half of the 21st century. While the oceanic climate response is mainly controlled by the GHG forcing, the land runoff response is impacted primarily by forcings other than GHG over certain regions, e.g., southern North America, southern Africa, central Africa, and eastern Asia. However, the importance of the GHG forcing on the land runoff changes grows in the future climate projection period compared to the historical period.

MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes

Abstract. We present MPAS-Seaice, a sea-ice model which uses the Model for Prediction Across Scales (MPAS) framework and spherical centroidal Voronoi tessellation (SCVT) unstructured meshes. As well as SCVT meshes, MPAS-Seaice can run on the traditional quadrilateral grids used by sea-ice models such as CICE. The MPAS-Seaice velocity solver uses the elastic–viscous–plastic (EVP) rheology and the variational discretization of the internal stress divergence operator used by CICE, but adapted for the polygonal cells of MPAS meshes, or alternatively an integral (“finite-volume”) formulation of the stress divergence operator. An incremental remapping advection scheme is used for mass and tracer transport. We validate these formulations with idealized test cases, both planar and on the sphere. The variational scheme displays lower errors than the finite-volume formulation for the strain rate operator but higher errors for the stress divergence operator. The variational stress divergence operator displays increased errors around the pentagonal cells of a quasi-uniform mesh, which is ameliorated with an alternate formulation for the operator. MPAS-Seaice shares the sophisticated column physics and biogeochemistry of CICE and when used with quadrilateral meshes can reproduce the results of CICE. We have used global simulations with realistic forcing to validate MPAS-Seaice against similar simulations with CICE and against observations. We find very similar results compared to CICE, with differences explained by minor differences in implementation such as with interpolation between the primary and dual meshes at coastlines. We have assessed the computational performance of the model, which, because it is unstructured, runs with 70 % of the throughput of CICE for a comparison quadrilateral simulation. The SCVT meshes used by MPAS-Seaice allow removal of equatorial model cells and flexibility in domain decomposition, improving model performance. MPAS-Seaice is the current sea-ice component of the Energy Exascale Earth System Model (E3SM).

OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) sea spray organic aerosol emissions – implementation in a global climate model and impacts on clouds

Abstract. Sea spray aerosol is one of the major sources of atmospheric particulate matter globally. It has increasingly been recognized that organic matter derived from ocean biological precursors contributes significantly to the composition of submicron sea spray and may modify sea spray aerosol impacts on clouds and climate. This paper describes the implementation of the OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) parameterization for sea spray organic aerosol emissions in a global Earth system model, the Energy Exascale Earth System Model (E3SM). OCEANFILMS is a physically based model that links sea spray chemistry with ocean biogeochemistry using a Langmuir partitioning approach. We describe the implementation details of OCEANFILMS within E3SM, compare simulated aerosol fields with observations, and investigate impacts on simulated clouds and climate. Four sensitivity cases are tested, in which organic emissions either strictly add to or strictly replace sea salt emissions (in mass and number) and are either fully internally or fully externally mixed with sea salt. The simulation with internally mixed, added organics agrees reasonably well with observed seasonal cycles of organic matter in marine aerosol and has been selected as the default configuration of the E3SM. In this configuration, sea spray organic aerosol contributes an additional source of cloud condensation nuclei, adding up to 30 cm−3 to Southern Ocean boundary-layer cloud condensation nuclei concentrations (supersaturation = 0.1 %). The addition of this new aerosol source strengthens shortwave radiative cooling by clouds by −0.36 W m−2 in the global annual mean and contributes more than −3.5 W m−2 to summertime zonal mean cloud forcing in the Southern Ocean, with maximum zonal mean impacts of about −4 W m−2 around 50–60∘ S. This is consistent with a previous top-down, satellite-based empirical estimate of the radiative forcing by sea spray organic aerosol over the Southern Ocean. Through its mechanistic approach, OCEANFILMS offers a path towards improved understanding of the feedbacks between ocean biology, sea spray organic matter, and climate.

CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM)

Abstract. Numerical models used in weather and climate prediction take into account a comprehensive set of atmospheric processes (i.e., phenomena) such as the resolved and unresolved fluid dynamics, radiative transfer, cloud and aerosol life cycles, and mass or energy exchanges with the Earth's surface. In order to identify model deficiencies and improve predictive skills, it is important to obtain process-level understanding of the interactions between different processes. Conditional sampling and budget analysis are powerful tools for process-oriented model evaluation, but they often require tedious ad hoc coding and large amounts of instantaneous model output, resulting in inefficient use of human and computing resources. This paper presents an online diagnostic tool that addresses this challenge by monitoring model variables in a generic manner as they evolve within the time integration cycle. The tool is convenient to use. It allows users to select sampling conditions and specify monitored variables at run time. Both the evolving values of the model variables and their increments caused by different atmospheric processes can be monitored and archived. Online calculation of vertical integrals is also supported. Multiple sampling conditions can be monitored in a single simulation in combination with unconditional sampling. The paper explains in detail the design and implementation of the tool in the Energy Exascale Earth System Model (E3SM) version 1. The usage is demonstrated through three examples: a global budget analysis of dust aerosol mass concentration, a composite analysis of sea salt emission and its dependency on surface wind speed, and a conditionally sampled relative humidity budget. The tool is expected to be easily portable to closely related atmospheric models that use the same or similar data structures and time integration methods.

An evaluation of the E3SMv1 Arctic ocean and sea-ice regionally refined model

Abstract. The Energy Exascale Earth System Model (E3SM) is a state-of-the-science Earth system model (ESM) with the ability to focus horizontal resolution of its multiple components in specific areas. Regionally refined global ESMs are motivated by the need to explicitly resolve, rather than parameterize, relevant physics within the regions of refined resolution, while offering significant computational cost savings relative to the respective cost of configurations with high-resolution (HR) everywhere on the globe. In this paper, we document results from the first Arctic regionally refined E3SM configuration for the ocean and sea-ice components (E3SM-Arctic-OSI), while employing data-based atmosphere, land, and hydrology components. Our aim is an improved representation of the Arctic coupled ocean and sea-ice state, its variability and trends, and the exchanges of mass and property fluxes between the Arctic and the sub-Arctic. We find that E3SM-Arctic-OSI increases the realism of simulated Arctic ocean and sea-ice conditions compared to a similar low-resolution E3SM simulation without the Arctic regional refinement in ocean and sea-ice components (E3SM-LR-OSI). In particular, exchanges through the main Arctic gateways are greatly improved with respect to E3SM-LR-OSI. Other aspects, such as the Arctic freshwater content variability and sea-ice trends, are also satisfactorily simulated. Yet, other features, such as the upper-ocean stratification and the sea-ice thickness distribution, need further improvements, involving either more advanced parameterizations, model tuning, or additional grid refinements. Overall, E3SM-Arctic-OSI offers an improved representation of the Arctic system relative to E3SM-LR-OSI, at a fraction (15 %) of the computational cost of comparable global high-resolution configurations, while permitting exchanges with the lower-latitude oceans that cannot be directly accounted for in Arctic regional models.

Better calibration of cloud parameterizations and subgrid effects increases the fidelity of the E3SM Atmosphere Model version 1

Abstract. Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

AbstractSub‐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: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in 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 with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and 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 mountainous areas.

Building a machine learning surrogate model for wildfire activities within a global Earth system model

Abstract. Wildfire is an important ecosystem process, influencing land biogeophysical and biogeochemical dynamics and atmospheric composition. Fire-driven loss of vegetation cover, for example, directly modifies the surface energy budget as a consequence of changing albedo, surface roughness, and partitioning of sensible and latent heat fluxes. Carbon dioxide and methane emitted by fires contribute to a positive atmospheric forcing, whereas emissions of carbonaceous aerosols may contribute to surface cooling. Process-based modeling of wildfires in Earth system land models is challenging due to limited understanding of human, climate, and ecosystem controls on fire counts, fire size, and burned area. Integration of mechanistic wildfire models within Earth system models requires careful parameter calibration, which is computationally expensive and subject to equifinality. To explore alternative approaches, we present a deep neural network (DNN) scheme that surrogates the process-based wildfire model with the Energy Exascale Earth System Model (E3SM) interface. The DNN wildfire model accurately simulates observed burned area with over 90 % higher accuracy with a large reduction in parameterization time compared with the current process-based wildfire model. The surrogate wildfire model successfully captured the observed monthly regional burned area during validation period 2011 to 2015 (coefficient of determination, R2=0.93). Since the DNN wildfire model has the same input and output requirements as the E3SM process-based wildfire model, our results demonstrate the applicability of machine learning for high accuracy and efficient large-scale land model development and predictions.

Relative importance of high-latitude local and long-range-transported dust for Arctic ice-nucleating particles and impacts on Arctic mixed-phase clouds

Abstract. Dust particles, serving as ice-nucleating particles (INPs), may impact the Arctic surface energy budget and regional climate by modulating the mixed-phase cloud properties and lifetime. In addition to long-range transport from low-latitude deserts, dust particles in the Arctic can originate from local sources. However, the importance of high-latitude dust (HLD) as a source of Arctic INPs (compared to low-latitude dust, LLD) and its effects on Arctic mixed-phase clouds are overlooked. In this study, we evaluate the contribution to Arctic dust loading and INP population from HLD and six LLD source regions by implementing a source-tagging technique for dust aerosols in version 1 of the US Department of Energy's Energy Exascale Earth System Model (E3SMv1). Our results show that HLD is responsible for 30.7 % of the total dust burden in the Arctic, whereas LLD from Asia and North Africa contributes 44.2 % and 24.2 %, respectively. Due to its limited vertical transport as a result of stable boundary layers, HLD contributes more in the lower troposphere, especially in boreal summer and autumn when the HLD emissions are stronger. LLD from North Africa and East Asia dominates the dust loading in the upper troposphere with peak contributions in boreal spring and winter. The modeled INP concentrations show better agreement with both ground and aircraft INP measurements in the Arctic when including HLD INPs. The HLD INPs are found to induce a net cooling effect (−0.24 W m−2 above 60∘ N) on the Arctic surface downwelling radiative flux by changing the cloud phase of the Arctic mixed-phase clouds. The magnitude of this cooling is larger than that induced by North African and East Asian dust (0.08 and −0.06 W m−2, respectively), mainly due to different seasonalities of HLD and LLD. Uncertainties of this study are discussed, which highlights the importance of further constraining the HLD emissions.

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

AbstractEarth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m−3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

Assessing Two Approaches for Enhancing the Range of Simulated Scales in the E3SMv1 and the Impact on the Character of Hourly US Precipitation

AbstractImproving the representation of precipitation in Earth system models is essential for understanding and projecting water cycle changes across scales. Progress has been hampered by persistent deficiencies in representing precipitation frequency, intensity, and timing in current models. Here, we analyze simulated US precipitation in the low‐resolution (LR) configuration of the Energy Exascale Earth System Model (E3SMv1) and assess the effect of two approaches to enhance the range of explicitly resolved scales: high‐resolution (HR) and multiscale modeling framework (MMF), which incur similar computational expense. Both E3SMv1‐MMF and E3SMv1‐HR capture more intense and less frequent precipitation on hourly and daily timescales relative to E3SMv1‐LR. E3SMv1‐HR improves the intensity over the Eastern and Northwestern US during winter, while E3SMv1‐MMF improves the intensity over the Eastern US and summer diurnal timing over the Central US. These results indicate that both methods may be needed to improve simulations of different storm types, seasons, and regions.

Reconstructing High Resolution ESM Data Through a Novel Fast Super Resolution Convolutional Neural Network (FSRCNN)

AbstractWe present the first application of a fast super resolution convolutional neural network (FSRCNN) based approach for downscaling earth system model (ESM) simulations. Unlike other SR approaches, FSRCNN uses the same input feature dimensions as the low resolution input. This allows it to have smaller convolution layers, avoiding over‐smoothing, and reducing computational costs. We adapt the FSRCNN to improve reconstruction on ESM data, we term the FSRCNN‐ESM. We use high‐resolution (∼0.25°) monthly averaged model output of five surface variables over North America from the US Department of Energy's Energy Exascale Earth System Model's control simulation. These high‐resolution and corresponding coarsened low‐resolution (∼1°) pairs of images are used to train the FSRCNN‐ESM and evaluate its use as a downscaling approach. We find that FSRCNN‐ESM outperforms FSRCNN and other super‐resolution methods in reconstructing high resolution images producing finer spatial scale features with better accuracy for surface temperature, surface radiative fluxes, and precipitation.

The Madden–Julian Oscillation in the Energy Exascale Earth System Model Version 1

AbstractThe present study examines the characteristics of the Madden–Julian Oscillation (MJO) events represented in the Energy Exascale Earth System Model version 1 (E3SMv1), DOE’s new Earth system model. The coupled E3SMv1 realistically simulates the eastward propagation of precipitation and moist static energy (MSE) anomalies associated with the MJO. As in observations, horizontal moisture advection and longwave radiative feedback are found to be the dominant processes in E3SMv1 that lead to the eastward movement and maintenance of the MJO MSE anomalies, respectively. Modulation of the diurnal cycle of precipitation in the Maritime Continent region by the MJO is also well represented in the model despite systematic biases in the magnitude and phase of the precipitation diurnal cycle. On the MJO impact over the midlatitude, E3SMv1 reasonably captures the pattern of the MJO teleconnections across the North Pacific and North America, with improvement in the performance in a high‐resolution version, despite the magnitude being a bit weaker than the observed. Regarding the interannual variability of the MJO, the El Niño‐Southern Oscillation (ENSO) modulation of the zonal extent of MJO’s eastward propagation, as well as associated changes in the mean state moisture gradient in the tropical west Pacific, are well reproduced in the model. However, MJO in E3SMv1 exhibits no sensitivity to the Quasi‐Biennial Oscillation (QBO), with the MJO propagation characteristics being almost identical between easterly QBO and westerly QBO years. Processes that have been suggested as critical to MJO simulation are also examined by utilizing recently developed process‐oriented diagnostics.

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

AbstractIn global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

Incorporating Microtopography in a Land Surface Model and Quantifying the Effect on the Carbon Cycle

AbstractNorthern peatlands are a terrestrial carbon store, with an annual sink of 0.1 Pg C yr−1 and a total storage estimate of 547 Pg C. Northern peatlands are also major contributors of atmospheric methane. Most land surface models do not accurately represent peatland carbon emissions, partly because they do not represent the hydrologic cycle and/or microtopography adequately. Interactions between water table depth and microtopography in peatlands influence decomposition and modulate CO2 and CH4 fluxes. A modified version of the land surface component of the Energy Exascale Earth System Model, was recently created to represent the microtopography and hydrology of a raised dome bog in northern Minnesota, USA. In this study, three microtopographic parameters are analyzed in the modified version: hummock height, hummock‐hollow spacing, and percent hollow. Terrestrial laser scanning observations are used to set uncertainty bounds for these parameters. Our model experiment results suggest that carbon‐related quantities of interest (QOI) were typically the most sensitive to hummock height, and those QOI (especially net ecosystem exchange, NEE) were sensitive to interactions between parameters. Furthermore, NEE was most relatively influenced by microtopographic parameters in the model, varying by 35%. We found that increasing hummock height resulted in more C being stored in plant tissue and less in soil organic matter. This coincided with decreases in Sphagnum and increases in Picea and shrub net primary production. These results suggest that future studies may consider extending prognostic capabilities of carbon cycling by incorporating hummock hollow microtopography into earth system models.