AI-Assisted Parameter Tuning Will Speed Development and Clarify Uncertainty in E3SM
Focal Area: This idea is best aligned with predictive modeling. It targets Earth System Model improvement and uncertainty quantification
54 ENVIRONMENTAL SCIENCES↗
SEARCH · Search NASA
Results for “Earth Systems Model”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales
Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.
54 ENVIRONMENTAL SCIENCES↗
Arctic Ocean Hydroacoustics
As the Arctic warms and loses its perennial ice cover, it is becoming more attractive for a variety of human uses. Hydroacoustic monitoring of this activity will grow in importance over the coming years and decades. Changes to the physical environment affect acoustic propagation and noise, with ramifications for our ability to detect and locate events and activities of interest. In this report, we use two long-term data sets from the Beaufort Sea, western Arctic Ocean, to determine how acoustic propagation conditions and seismic source detections are impacted by the changing environment: 1) oceanographic observations from ice-tethered profilers, and 2) passive acoustic recordings from a hydrophone. We find that changes to Beaufort Sea thermohaline stratification is stabilizing a subsurface duct, leading to more focused acoustic energy arrivals. Detections of catalogued submarine earthquakes show geographic differences in signal strength between seismic and acoustic waves. Signal strength increases with earthquake magnitude, but relationships to other source and path factors are less clear. Ambient noise also has clear seasonal patterns in the Arctic, with relatively low noise in the spring, higher noise near 1 Hz in summer, and higher noise near 10 Hz in winter. Climate change is expected to modify these seasonal noise patterns, impacting event detection. Future work will further investigate the mechanisms of ice effects on sound and couple acoustic modeling to an Earth System Model.
54 ENVIRONMENTAL SCIENCES↗
SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice
Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.
54 ENVIRONMENTAL SCIENCES↗
Demeter calibration and uncertainty analysis
Demeter is a community spatial downscaling model that disaggregates land use and land cover changes projected by integrated human-Earth system models. Demeter has not been intensively calibrated, and we still lack a good knowledge about its sensitivity to key parameters and the parameter uncertainties. We used long-term global satellite-based land cover records to calibrate key Demeter parameters. The results identified the optimal parameter values and showed that the parameterization substantially improved the model’s performance. The parameters of intensification ratio and selection threshold were the most sensitive and needed to be carefully tuned, especially for regional applications. Further, small parameter uncertainties after calibration can be inflated when propagated into future scenarios, suggesting that users should consider the parameterization equifinality to better account for the uncertainties in the Demeter downscaled products. Our study provides a key reference for Demeter users, and ultimately contribute to reducing the uncertainties in Earth system model simulations.
54 ENVIRONMENTAL SCIENCES↗
Projecting Future Energy Production from Operating Wind Farms in North America. Part II: Statistical Downscaling
Abstract Capacity factors (CFs) derived from daily expected power at 22 operating wind farms in different regions of North America are used as predictands to train statistical downscaling algorithms using output from ERA5. The statistical downscaling models are then used to make CF projections for a suite of CMIP6 Earth System Models (ESMs). Downscaling is performed using a hybrid statistical approach that employs synoptic types derived using k -means clustering applied to sea level pressure fields with variance corrections applied as a function of the pressure gradient intensity. ESMs exhibit marked variability in terms of the skill with which the frequency of synoptic types and pressure gradients are reproduced relative to ERA5, and that differential skill is used to infer differential credibility in the associated CF projections. Projections of median annual mean CF [P50(CF)] in each 20-yr period from 1980 to 2099 show evidence of declines at most wind farms except in parts of the southern Great Plains, although the magnitude of the changes is strongly dependent on the ESM. For example, P50(CF) in 2080–99 deviate from those in 1980–99 by from −3.1 to +0.2 percentage points in the Northeast. The largest-magnitude declines in P50(CF) ranging from −3.9 to −2 percentage points are projected for the southern West Coast. CF trends exhibit marked seasonality and are strongly linked to changes in the relative intensity of future synoptic patterns, with much less impact from shifts in the occurrence of synoptic types over time. Internal climate modes continue to play a significant role in inducing interannual variability in wind power production, even under high radiative forcing scenarios. Significance Statement We describe how future climate changes may affect wind resources and wind power generation. Near-term changes in projected wind power electricity generation potential at operating wind farms over North America are small, but by the end of the current century electricity production is projected to decrease in many areas but may increase in parts of the southern Great Plains. The amount of change in projected wind power production is a strong function of the Earth system model that is downscaled and also depends on the continued presence of internally forced climate variability. An additional dependence on the amount of greenhouse gas–induced global warming indicates the transition of the energy sector to low-carbon sources may assist in maintaining the abundant U.S. wind resource.
Meteorology & Atmospheric Sciences↗
Examining the Ice-Nucleating Particles from the North Slope of Alaska (ExINP-NSA) (Final Campaign Report)
The Examining the Ice-Nucleating Particles from the North Slope of Alaska (ExINP-NSA) campaign was conducted at the National Oceanic and Atmospheric Administration’s (NOAA) Barrow Atmospheric Baseline Observatory (71.3230° N, 156.6114° W, “BRW” hereafter), next to the Atmospheric Radiation Measurement (ARM) user facility NSA site and ~ 6 km northeast of the town of Utqiaġvik. The location of the NSA site is shown in Figure 1. Our observing period began in October 2021 and continued until May 2024. This campaign was funded by the U.S. Department of Energy (DOE) Office of Science Early Career Research Program through grant number DE-SC001879. This grant contains funding for three ARM field campaigns. The previous field campaigns were performed at ARM’s Southern Great Plains atmospheric observatory in Oklahoma (SGP; 36.6073° N, 97.4876° W) and the ARM’s Eastern North Atlantic atmospheric observatory on Graciosa Island, Azores (ENA; 39.0916° N, 28.0257° W). The ExINP-NSA campaign aims included: • Determining the number concentration of ice-nucleating particles (n INPs ) active at temperatures spanning the range of heterogeneous freezing processes (from ≈−30 °C to 0 °C) using a combination of online and offline measurements, • Determining whether local meteorological conditions and/or synoptic scale air mass transport impact INP abundance and/or ice nucleation efficiency, • Examining if the physicochemical properties of INPs relate to aerosol chemistry, and • Assessing if there is a similarity in INP properties across three ARM sites Multi-seasonal datasets of INP abundance in the NSA region were delivered from this campaign. This campaign also allowed researchers to perform a comprehensive analysis of atmospheric INPs based on long- term ground-based measurements in the Alaskan Arctic. Our data and results from the ExINP-NSA campaign will help refine current earth system models. One of the stated goals of ARM is to advance aerosol-cloud ice interaction, which will be a direct result of this campaign. Current earth system models poorly represent INPs, and the 15-minute time resolution data over several seasons generated during this campaign will provide an invaluable resource, especially combined with the datasets generated during two previous ARM ExINP campaigns. These datasets will allow for a greater understanding of ice nucleation processes as they may (or may not) relate to local meteorological processes and aerosol chemistry and will eventually help further the understanding of the Earth’s atmospheric processes and energy balance.
54 ENVIRONMENTAL SCIENCES↗
CMIP7 Data Request: atmosphere priorities and opportunities
This paper presents a comprehensive overview of the Coupled Model Intercomparison Project Phase 7 (CMIP7) request for data unlocking key research avenues in atmospheric science and provides justification for the resources needed to produce this data. Topics within the CMIP7 Atmosphere Theme centre around processes and feedbacks in atmospheric science such as clouds, aerosols and atmospheric chemistry, atmospheric circulation, temperature variability and extremes, radiative forcings, and Earth system model evaluation. These topics are summarised in this paper as scientific “opportunities” which will be realised through CMIP7 experiments and Earth system model outputs. These opportunities were submitted by a thematic group of atmospheric science community representatives combined with an extended consultation process. The production of these variables will close key gaps and uncertainties identified during previous rounds of CMIP, and will be broadly used by scientific, policy, governmental, industry, and other communities that rely on climate model projections for research and decision making, including supporting the 7th Intergovernmental Panel on Climate Change Assessment Report (AR7). As an author group, we also reflect on the process used to collate this data request and make recommendations to future CMIP governance on implementing a consultation on this scale in the future.
58 GEOSCIENCES↗
Late Holocene glacier and climate fluctuations in the Mackenzie and Selwyn mountain ranges, northwestern Canada
Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10 Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10 Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 °C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km 3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10 Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.
54 ENVIRONMENTAL SCIENCES↗
Impacts of Topography-Based Subgrid Scheme and Downscaling of Atmospheric Forcing on Modeling Land Surface Processes in the Conterminous US
The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.
54 ENVIRONMENTAL SCIENCES↗
Simulating Marine Ecosystem Dynamics and Biogeochemical Cycling With Multiple Plankton Functional Types
Current representations of marine ecosystems in Earth System Models are greatly simplified, neglecting key interactions between dynamic food webs, biogeochemistry, and climate change. We use the Marine Biogeochemistry Library code base within the Community Earth System Model 2.2.2 to create an expanded ecosystem model with eight phytoplankton groups and four zooplankton size classes (MARBL-8P4Z). Incorporating more specific plankton types and size classes has the potential to capture a wider range of possible behaviors of the ecosystem, its complex interactions with biogeochemistry, and its feedback to climate change. It also permits stronger observational constraints, including in situ group-specific biomass and various observational estimates of plankton community composition. MARBL-8P4Z broadly captures observed global-scale patterns in biomass and community composition for both phytoplankton and zooplankton, with a good performance in simulating broad biogeochemistry fields. The model shows comparable spatial patterns and magnitudes to the observed picophytoplankton biomass (Prochlorococcus, Synechococcus, picoeukaryotes), and captures the seasonal cycle of mesozooplankton biomass. Picophytoplankton groups and microzooplankton dominate biomass and production in oligotrophic, subtropical regions, while nano-phytoplankton, diatoms and the larger zooplankton groups prevail at higher latitudes and within upwelling zones. The model simulates reasonable energy transfer efficiency through the food web, with tight linkages between the phytoplankton community composition, zooplankton grazing, and carbon export, with the potential to link to fisheries models. Thus, MARBL-8P4Z has the potential to account for key climate-driven ecological shifts in the plankton that will modify ocean biogeochemistry in the future.
54 ENVIRONMENTAL SCIENCES↗
Scaling Ultrahigh-Resolution E3SM Land Model for Leadership-Class Supercomputers
This paper presents advancements in scaling the ultrahigh-resolution E3SM Land Model (uELM) for deployment on leadership-class supercomputers, addressing the increased demand for km-scale Earth system modeling. By focusing on km-scale ELM simulations, we enhance predictive capabilities for climate interactions, facilitating improved responses to climate change impacts on energy systems, agriculture, and water resources. Our approach leverages innovative software architecture optimizations, sophisticated data handling techniques, and advanced parallel processing, achieving strong scalability on two leadership supercomputers (2400 nodes (105,600 cores) on Summit, and 1200 nodes (76,800 cores) on Frontier). Results from extensive scalability assessments on the Summit and Frontier also demonstrate outstanding I/O performance (close to 400 GB/s write throughput) and the model's ability to efficiently handle increasing computational demands. This study not only establishes uELM's capability for high-resolution simulations over vast geographical domains, but also sets a foundation for future Earth system modeling breakthroughs.
Wang, Dali [ORNL] (ORCID:0000000168065108)↗
End-to-End Differentiable Modeling and Management of the Environment
Focal Area: (2) Data acquisition and assimilation enabled by machine learning, AI, and advanced methods including experimental/network design/optimization. We emphasize the importance of leveraging optimization techniques from AI/machine learning (ML) to solve challenging problems in Earth system modeling. Science Challenge: Automatic differentiation has had a transformative effect on ML by allowing the calculation of gradients of arbitrary functions in an incredibly large class of models. We can potentially realize similar improvements in parameter estimation and control for Earth system models (ESMs) by reimplementing them in computational frameworks from ML. Practitioners working with large (>10 7 parameters) models in ML and AI can obtain good predictive performance in a range of spatiotemporal tasks by making use of optimization via stochastic gradient descent and incorporating prior knowledge at multiple levels. We propose writing ESMs in open-source computational frameworks such as Torch, Tensorflow, and JAX to greatly expand the scope of environmental forecasting and management challenges, which can be addressed by leveraging automatic differentiation and gradient descent-like algorithms. We do not call for a wholesale replacement of physical models with data-driven surrogates, but rather advocate for interleaving physical and empirical equations in a manner that is most faithful to the extent of our scientific knowledge and observational data. Central to this topic is the merging of differentiable physical simulations with differentiable optimization layers, which are now both beginning to come to the forefront.
54 ENVIRONMENTAL SCIENCES↗
How Can NASA Science Benefit Solar Energy Development and Assessment?
This is a keynote presentation for a workshop highlighting the use of remote sensing and modeling products derived with NASA missions and Earth system modeling for solar energy system development. We broadly introduce NASA's fleet of satellite producing publicly available data products that would be useful for the information needed solar energy development. We also introduce the work in Earth system science modeling that provides additional global gridded data products that are also relevant. We identify that one issue is obtaining and using this data is the level of expertise needed to find it. NASA has developed a new Earth Action program to address this gap and bring the data closer to the user. Now under the Energy Resources element, the NASA LaRC's POWER (Prediction of Worldwide Renewable Energy) project has developed a web services platform that provides key information needed for solar energy development as analysis-ready. This means that the data for particular locations are immediately useful for help provide the information needed by engineers and architects to plan for solar systems. Thus, the presentation shows why NASA products would be relevant and then shows examples of finding the needed data through the POWER web portal.
remote sensing↗
Reversal of Increasing Tropical Ocean Hypoxia Trends With Sustained Climate Warming
Dissolved oxygen (O 2 ) is essential for the survival of marine animals. Climate change impacts on future oxygen distributions could modify species biogeography, trophic interactions, biodiversity, and biogeochemistry. The Coupled Model Intercomparison Project Phase 5 models predict a decreasing trend in marine O 2 over the 21st century. Here we show that this increasing hypoxia trend reverses in the tropics after 2100 in the Community Earth System Model forced by atmospheric CO 2 from the Representative Concentration Pathway 8.5 and Extended Concentration Pathway 8.5. In tropical intermediate waters between 200 and 1,000 m, the model predicts a steady decline of O 2 and an expansion of oxygen minimum zones (OMZs) during the 21st century. By 2150, however, the trend reverses with oxygen concentration increasing and OMZ volume shrinking through 2300. A novel five-box model approach in conjunction with output from the full Earth system model is used to separate the contributions of biological and physical processes to the trends in tropical oxygen. The tropical O 2 recovery is caused mainly by reductions in tropical biological export, coupled with a modest increase in ventilation after 2200. The time-evolving oxygen distribution impacts marine nitrogen cycling, with potentially important climate feedbacks.
54 ENVIRONMENTAL SCIENCES↗
Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers
Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.
Flood Simulation↗
Subseasonal Earth System Prediction with CESM2
Abstract Prediction systems to enable Earth system predictability research on the subseasonal time scale have been developed with the Community Earth System Model, version 2 (CESM2) using two configurations that differ in their atmospheric components. One system uses the Community Atmosphere Model, version 6 (CAM6) with its top near 40 km, referred to as CESM2(CAM6). The other employs the Whole Atmosphere Community Climate Model, version 6 (WACCM6) whose top extends to ∼140 km, and it includes fully interactive tropospheric and stratospheric chemistry [CESM2(WACCM6)]. Both systems are utilized to carry out subseasonal reforecasts for the 1999–2020 period following the Subseasonal Experiment’s (SubX) protocol. Subseasonal prediction skill from both systems is compared to those of the National Oceanic and Atmospheric Administration CFSv2 and European Centre for Medium-Range Weather Forecasts (ECMWF) operational models. CESM2(CAM6) and CESM2(WACCM6) show very similar subseasonal prediction skill of 2-m temperature, precipitation, the Madden–Julian oscillation, and North Atlantic Oscillation to its previous version and to the NOAA CFSv2 model. Overall, skill of CESM2(CAM6) and CESM2(WACCM6) is a little lower than that of the ECMWF system. In addition to typical output provided by subseasonal prediction systems, CESM2 reforecasts provide comprehensive datasets for predictability research of multiple Earth system components, including three-dimensional output for many variables, and output specific to the mesosphere and lower-thermosphere (MLT) region from CESM2(WACCM6). It is shown that sudden stratosphere warming events, and the associated variability in the MLT, can be predicted ∼10 days in advance. Weekly real-time forecasts and reforecasts with CESM2(CAM6) and CESM2(WACCM6) are freely available. Significance Statement We describe here the design and prediction skill of two subseasonal prediction systems based on two configurations of the Community Earth System Model, version 2 (CESM2): CESM2 with the Community Atmosphere Model, version 6 [CESM2(CAM6)] and CESM 2 with Whole Atmosphere Community Climate Model, version 6 [CESM2(WACCM6)] as its atmospheric component. These two systems provide a foundation for community-model based subseasonal prediction research. The CESM2(WACCM6) system provides a novel capability to explore the predictability of the stratosphere, mesosphere, and lower thermosphere. Both CESM2(CAM6) and CESM2(WACCM6) demonstrate subseasonal surface prediction skill comparable to that of the NOAA CFSv2 model, and a little lower than that of the ECMWF forecasting system. CESM2 reforecasts provide a comprehensive dataset for predictability research of multiple aspects of the Earth system, including the whole atmosphere up to 140 km, land, and sea ice. Weekly real-time forecasts, reforecasts, and models are publicly available.
54 ENVIRONMENTAL SCIENCES↗
Assessing Simulations of Forest Hurricane Disturbance and Recovery in Puerto Rico by ELM-FATES Using Field Measurements
In the past three decades, Puerto Rico (PR) experienced five hurricanes that met or exceeded category three, and they caused severe forest structural damage and elevated tree mortality. To improve our mechanistic understanding of hurricane impacts on tropical forests and assess hurricane-affected forest dynamics in Earth system models, we use in situ forest measurements at the Bisley Experimental Watersheds in Northeast PR to evaluate the Functionally Assembled Terrestrial Ecosystem Simulator coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES). The observations show that before Hurricane Hugo, 77.3% of the aboveground biomass (AGB) is from the shade-tolerant plant function type (PFT). The Hugo-induced mortality rates are over ~50%, and they induce a ~39% AGB reduction, which recovers to a level like the pre-Hugo condition in 2014, following a second, lower intensity hurricane, Georges. We perform numerical experiments that simulate damage from Hugo and Georges on the forests, including defoliation, sapwood and structural biomass damage, and hurricane-induced mortality. ELM-FATES can reasonably represent coexistence between the two PFTs–light-demanding and shade-tolerant–for both the pre-Hugo and post-Hugo conditions. The model represents a reasonable size distribution of mid-and large-sized trees although it underestimates AGB, likely due to the overestimated nonhurricane mortality. ELM-FATES temporarily stimulated leaf biomass and diameter increment after Georges, an effect that should be tested with observations of future hurricane defoliation events. This research indicates that addressing model-data mismatches in tree mortality and understory dynamics are essential to simulation of more extreme hurricane effects under climate change.
58 GEOSCIENCES↗