OVERALL PERFORMANCE MEASURES
2nd QUARTER METRIC COMPLETED:
Demonstrating and Evaluating Kilometer-Scale Regional Refinement Modeling Capability in Simulated Water Cycle Extremes in the U.S.
Product Definition
Increasing the horizontal grid resolution in atmosphere models has been shown to lead generally to improvements in modeling extreme events such as extreme precipitation produced by thunderstorms (e.g., Pope and Stratton 2002, Roekner et al. 2006, Wehner et al. 2010, Jong et al. 2023). At higher resolutions, atmosphere models can produce stronger vertical motions, for example, by better resolving the temperature and moisture gradients associated with fronts and tropical cyclones and representing the orographic uplift associated with mountains, and subsequently the higher condensation rates needed to generate intense precipitation associated with storms. However, atmosphere models simulate the precipitation processes using physics parameterizations such as cloud microphysics and convection with resolution- and time step-dependent behaviors, which complicate the relationship between model resolution and the simulated mean and extreme precipitation (e.g., Kopparla et al. 2013, Yang et al. 2014, Wehner et al. 2021).
Global kilometer-scale models that have emerged in the last decade (Stevens et al. 2019) hold great promise to further improve modeling of extreme weather events by explicitly resolving deep convection, a major mechanism for generating extreme precipitation, and potentially improving the simulation of mesoscale and large-scale atmospheric environments for the storms. Recently, the Energy Exascale Earth System Model (E3SM) project funded by the U.S. Department of Energy developed a global kilometer-scale model called Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) (Caldwell et al. 2021). Implementation of SCREAM, including innovative algorithmic and software engineering advancements to run on the Frontier Exascale computer, has demonstrated a groundbreaking performance of computational throughput > 1 simulated year per day (SYPD) (Taylor et al. 2023). However, running decadal-to-century-scale simulations using global kilometer-scale models remains computationally challenging, especially for quantifying the statistics of extreme events and their future changes (which typically requires large ensemble simulations). To circumvent such computational challenges, SCREAM has a regional refinement capability that allows kilometer-scale modeling to be performed in regions of interest within the context of global modeling at coarser resolutions outside those regions. This capability was first demonstrated in a study that used SCREAM with regional refinement over the central-eastern U.S. for simulating the strong winds and convection associated with a derecho event (Liu et al. 2023).
Building on the demonstration of Liu et al. (2023), this report documents the configurations and evaluation of SCREAM using regional refinement for kilometer-scale simulations of cold-season and warm-season storm events that produce extreme precipitation across diverse geographical regions in the U.S. Three examples are included to demonstrate the regional refinement capability of SCREAM to simulate (1) atmospheric rivers that made landfall in the U.S. Pacific Northwest Puget Sound basin that produced flooding, (2) mesoscale convective systems and isolated deep convection in the central-eastern U.S. that produced heavy precipitation in the region, and (3) a blizzard in the northeastern U.S. that produced snow and ice that damaged energy infrastructure. The results provide clear evidence that using regional refinement at 3.25-km grid spacing, the SCREAM model provides a scientifically robust and computationally efficient capability for simulating a variety of extreme weather events, including atmospheric rivers, mesoscale convective systems, and winter storms, across different geographies and regions of the contiguous U,S. Such capability is critical for understanding water cycle extremes such as flooding caused by heavy precipitation, rain-on-snow events, and compound events, as well as ice storms that cause disruptions and power outages.
Product Documentation
This report documents the modeling of water cycle extremes in the U.S. using SCREAM with regional refinement capability for kilometer-scale modeling. This capability is demonstrated using three examples of simulating cold-season and warm-season storms in the U.S. underlined by different processes. SCREAM is configured for a global domain at 25-km grid spacing with regional refinement at 3.25-km grid spacing covering part of the U.S. and the surrounding ocean where the storms developed. As SCREAM does not include any parameterization of deep convection, which is needed at 25-km grid spacing, nudging was applied to constrain the SCREAM simulations with the European Center for Medium-Range Forecast Reanalysis version 5 (ERA5) (Hersbach et al. 2020) outside the refined region. Initialized using atmospheric conditions from ERA5, the simulations, each covering 5-50 days, were compared with observations to demonstrate the model skill in simulating important water cycle aspects of the storm events.
1st QUARTER METRIC COMPLETED:
Developing and demonstrating a Kilometer-Scale Land Simulation Capability for Modeling the Terrestrial Water Cycle over the Contiguous U.S.
Product Definition
Kilometer-scale (k-scale) modeling allows for explicit modeling of physical processes that are poorly represented in climate models with typical resolutions of 25-100 km, thus providing opportunities to significantly improve the accuracy of climate simulations (Slingo et al. 2022). The recent advancements in computing power are making k-scale regional and global simulations using Land Surface Models (LSMs) and Earth System Models (ESMs) increasingly feasible (Condon et al. 2020, Caldwell et al. 2021). Groundwater is a vital human water resource that provides 20-30% of global freshwater withdrawals (Döll 2009). Previous k-scale LSM studies for specific watersheds or basins have demonstrated the impacts of fine scale-structures on terrestrial hydrologic processes including groundwater dynamics (Maxwell and Kollet 2008, Fan et al. 2013).
However, the parameters of LSMs within ESMs being run at the k scale are typically derived from coarse-resolution data sets or outdated data sets. Consequently, k-scale modeling may not accurately represent fine-scale land surface heterogeneity unless high-resolution land surface parameters at the kilometer or finer scales are used. Additionally, LSMs will have to be recalibrated at k scale to accurately simulate terrestrial processes (Ruiz‐Vásquez et al. 2023) and the recalibration is expected to have a significantly large computational cost at such high spatial resolution.
In this report, we develop a first-of-a-kind k-scale global land simulation capability and demonstrate this for modeling the terrestrial water cycle over the contiguous U.S. in ESMs. To enable this capability, we first develop a new set of global land surface parameters at 1-km resolution by using the newest high‑resolution data sources for multiple years by combining the latest and most accurate available global data sets to provide input data for ESMs, including the U.S. Department of Energy’s flagship Energy Exascale Earth System Model (E3SM) (Leung et al. 2020). Second, we produce an initial 5-year simulation using E3SM Land Model version 2 (ELMv2) over the contiguous United States (CONUS) at 1-km resolution using the newly developed land surface parameters. A spatial scaling analysis was performed to underscore the value of the high-resolution land surface parameters, and AI/ML methods were used to identify the most important land surface parameters and climate conditions that drive the spatial variability at k scale and hence, spatial information loss as resolutions were coarsened in ELMv2 simulations. Using the newly developed 1-km input data, the k-scale simulations show significant spatial heterogeneity of soil moisture, latent heat, emitted longwave radiation, and absorbed shortwave radiation, which is expected to have important effect on modeling of land-atmosphere interactions. Third, we demonstrate significant improvements in the ELM-simulated water table depth over CONUS by calibrating ELMv2’s subsurface drainage parameterization for k-scale modeling.
Product Documentation
This report documents a k-scale land simulation capability, which includes the development of global k-scale land surface parameters for LSMs (2.1), the application of the newly developed surface parameters in k-scale ELMv2 simulation over CONUS and analysis that demonstrates the importance of the k-scale surface parameters data set (2.2), and the calibration of ELMv2’s subsurface hydrologic processes to improve the prediction of groundwater at a 1-km resolution over CONUS (2.3).