Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Resolution refinement”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Quantum state preparation with resolution refinement

We introduce a method called resolution refinement that allows one to bootstrap eigenstate preparation on a quantum computer. We first prepare an eigenstate of a low-resolution Hamiltonian using any method of choice. The eigenstate is then lifted to higher resolution and adiabatically evolved to produce the corresponding eigenstate of a higher-fidelity Hamiltonian. We give examples of resolution refinement applied to both single-particle basis states as well as a spatial lattice grid. For basis refinement, we compute few-body ground states of the Busch model for interacting particles in a harmonic trap in one dimension. For lattice refinement, we compute Hartree-Fock nuclear states for a central Woods-Saxon potential in three dimensions, and we compute bound states and continuum states in a multi-species Hubbard model of fermions in one dimension. In all cases, the method is efficient and requires an adiabatic evolution time that scales with the inverse of the energy gap times the square root of the system size. We show that this very favorable scaling arises from the fact that resolution refinement does not make large changes to the structure or energies of the low-energy eigenstates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Simulation of Compound Flooding Using River‐Ocean Two‐Way Coupled E3SM Ensemble on Variable‐Resolution Meshes

Abstract Coastal zone compound flooding (CF) can be caused by the interactive fluvial and oceanic processes, particularly when coastal backwater propagates upstream and interacts with high river discharge. The modeling of CF is limited in existing Earth System Models (ESMs) due to coarse mesh resolutions and one‐way coupled river‐ocean components. In this study, we present a novel multi‐scale coupling framework within the Energy Exascale Earth System Model (E3SM), integrating global atmosphere and land with interactively coupled river and ocean models using different meshes with refined resolutions near the coastline. To evaluate this framework, we conducted ensemble simulations of a CF event (Hurricane Irene in 2011) in a Mid‐Atlantic estuary. The results demonstrate that the novel E3SM configuration can reasonably reproduce river discharge and sea surface height variations. The two‐way river‐ocean coupling improves the representation of coastal backwater effects at the terrestrial‐aquatic interface that are caused by the combined actions of tide and storm surge during the CF event, thus providing a valuable modeling tool for better understanding the river‐estuary‐ocean dynamics in extreme events under climate change. Notably, our results show that the most significant CF impacts occur when the highest storm surge generated by a tropical cyclone meets with a moderate river discharge. This study highlights the state‐of‐the‐art advancements developed within E3SM for simulating multi‐scale coastal processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluating E3SM Global Storm‐Resolving Model Simulations of Deep Convection: Insights From DP‐SCREAM During TRACER

Global Storm-Resolving Models (GSRMs) are becoming increasingly vital for advancing climate modeling and improving the prediction of extreme weather events. Houston, a coastal region frequently affected by deep convective storms, offers an ideal setting to evaluate the ability of GSRMs to simulate deep convection. This study assesses the performance of the Doubly Periodic Simple Cloud-Resolving E3SM (Energy Exascale Earth System Model) Atmosphere Model (DP-SCREAM) using observations from the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign. DP-SCREAM effectively reproduces the diurnal cycles of clouds and precipitation, demonstrating much greater skill than the E3SM single column model. The DP-SCREAM is demonstrated to be applicable to coastal regions, partially due to the forcing data sets already capturing the influence of breezes. DP-SCREAM also replicates biases persistent in the global version of SCREAM: the underrepresentation of boundary layer shallow clouds, a lack of mid-level congestus clouds, and the popcorn convection, characterized by small and disorganized convective cells generating the strongest precipitation. To investigate these issues, two sensitivity experiments were conducted: increasing the mixing length and scaling up the buoyancy flux within the Simplified Higher Order Closure scheme. Increasing the mixing length improved mid-level congestus representation and reduced unrealistic early morning fog occurrence. Enhancing buoyancy flux only marginally improved the bias of underproduced big convective cells. In conclusion, an additional resolution sensitivity test at 0.5 km grid spacing demonstrated that a refined horizontal resolution alone is insufficient to resolve these biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗
DOE OSTI2026

Simulation of wind and solar energy generation over California with E3SM SCREAM regionally refined models at 3.25 km and 800 m resolutions

This study presents wind and solar power generation estimates derived from the US Department of Energy’s Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) Regionally Refined Models (RRM) over California at 3.25 km and 800 m horizontal resolutions, using the Python wrapper for the System Advisor Model (PySAM). The resulting wind and solar generation estimates are compared to monthly capacity factors reported to the Energy Information Administration (EIA), High-Resolution Rapid Refresh (HRRR, 3 km resolution) forecast model, and E3SM North American regionally refined model (NARRM, 25 km resolution). We systematically assess the impacts of generation modeling assumptions, meteorological models, and horizontal resolution. Results show that resolution plays a dominant role for wind energy: increasing from 25 to 3.25 km brings qualitative and quantitative improvements, most notably by resolving the phase error in the seasonal cycle found in coarser simulations. However, further refinement to 800 m offers minimal gains. SCREAM performs better than HRRR for solar power generation in single- and dual-axis tracking systems, likely due to more accurate surface radiation. The sensitivity of PySAM to system configuration, particularly for axis-tracking modeling in photovoltaics, is also highlighted. Overall, SCREAM-RRM shows strong potential for high-resolution energy assessments, with future progress depending on more in situ observations and clearer quantification of uncertainties in generation modeling.

Geosciences↗
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Atmospheric-river-induced precipitation in California as simulated by the regionally refined Simple Convective Resolving E3SM Atmosphere Model (SCREAM) Version 0

Abstract. Using the regionally refined mesh (RRM) configuration of the US Department of Energy's Simple Cloud-Resolving Energy Exascale Earth System Model (E3SM) Atmosphere Model (SCREAM), we simulate and evaluate four meteorologically distinct atmospheric river events over California. We test five different RRM configurations, each differing in terms of the areal extent of the refined mesh and the resolution (ranging from 800 m to 3.25 km). We find that SCREAM RRM generally has a good representation of the AR-generated precipitation in CA, even for the control simulation which has a very small 3 km refined patch, and is able to capture the fine-scale regional distributions that are controlled largely by the fine-scale topography of the state. It is found that SCREAM generally has a wet bias over topography, most prominently over the Sierra Nevada mountain range, with a corresponding dry bias on the lee side. We find that refining the resolution beyond 3 km (specifically 1.6 km and 800 m) has virtually no benefit towards reducing systematic precipitation biases but that improvements can be found when increasing the areal extent of the upstream refined mesh. However, these improvements are relatively modest and only realized if the size of the refined mesh is expanded to the scale where employing RRM no longer achieves the substantial cost benefit it was intended for.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗
DOE OSTI2025

A staged deep learning approach to spatial refinement in 3D temporal atmospheric transport

High-resolution spatiotemporal simulations effectively capture the complexities of atmospheric plume dispersion in complex terrain. However, their high computational cost makes them impractical for applications requiring rapid responses or iterative processes, such as optimization, uncertainty quantification, or inverse modeling. To address this challenge, this work introduces the Dual-Stage Temporal Three-dimensional UNet Super-resolution (DST3D-UNet-SR) model, a highly efficient deep learning model for plume dispersion predictions. DST3D-UNet-SR is composed of two sequential modules: the temporal module (TM), which predicts the transient evolution of a plume in complex terrain from low-resolution temporal data, and the spatial refinement module (SRM), which subsequently enhances the spatial resolution of the TM predictions. We train DST3D-UNet-SR using a comprehensive dataset derived from high-resolution large eddy simulations (LES) of plume transport. We propose the DST3D-UNet-SR model to significantly accelerate LES of three-dimensional (3D) plume dispersion by three orders of magnitude. Additionally, the model demonstrates the ability to dynamically adapt to evolving conditions through the incorporation of new observational data, substantially improving prediction accuracy in high-concentration regions near the source.

3D temporal sequences↗
DOE OSTI2025

Impacts of Resolution on Heavy‐Precipitating Storms in Climate Model Hindcasts

The present study investigates the impact of horizontal resolutions on heavy‐precipitating storms using the Energy Exascale Earth System Model version2 (E3SMv2) at low (∼100 km, LR) and high (∼25 km, HR) resolutions through short‐range hindcasts. The short‐range hindcast approach ensures a faithful comparison of model resolution in simulating the same storm events under a controlled large‐scale environment. Using a phenomenon‐based framework, we attribute precipitation to specific storm types: tropical cyclones (TCs), extratropical cyclones, atmospheric rivers, and mesoscale convective systems (MCSs). Our findings show that E3SM hindcasts with both HR and LR configurations significantly underestimate storm‐associated precipitation intensity but overestimate precipitation from other sources. Furthermore, both HR and LR hindcasts face significant challenges in accurately simulating extreme precipitation events, particularly over MCS hotspots. Nevertheless, HR simulations capture more detailed and intense precipitation patterns with an improved representation of storm dynamics. HR hindcasts produce 16% more storm precipitation compared to LR. For precipitation extremes, HR simulates a 33% higher 99th percentile precipitation magnitudes compared to LR, and most of the increment comes from these four heavy‐precipitating storm types. The increase in precipitation mainly comes from stratiform precipitation rather than convective precipitation. The improvement in HR simulations varies across different storm types with TC showing the largest improvement. The phenomenon‐based approach provides important insights into precipitation simulations especially for extremes. Our results emphasize the need for further refinement in high‐resolution models to improve the accuracy of precipitation predictions, which is crucial for better understanding and mitigating climate change impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data-informed grid refinement to improve traveltime accuracy in the regional seismic traveltime (RSTT) model

The regional seismic traveltime (RSTT) model predicts traveltimes of regional seismic phases accounting for 3-D structure of the crust and the upper mantle on a global scale. Previous versions of the RSTT model have been implemented using nodes separated by ∼1° spacing across the globe. A regional-scale study using regional Pn and Pg traveltimes across Israel and the Middle East demonstrated that data driven, systematic grid refinement reduces traveltime residuals and enhances resolution of smaller tectonic features in regions having dense ray coverage. High density Pn ray coverage in the western US, Europe, Middle East and East Asia can likewise provide the resolution that allows systematic global grid refinement of the RSTT model. In this study, we use a large number of Pn ray paths originating from events located with an epicentral location uncertainty of 25 km (GT25) or better. We conduct targeted grid refinements at 1.0°, 0.5°, 0.25° and 0.125° on a global scale, producing a refined RSTT model that yields a 21.6 per cent reduction in median event location error in Europe and the Middle East, when compared with the original global RSTT model presented in Begnaud et al. The new model also resolves finer tectonic structures in regions with high Pn ray density.

58 GEOSCIENCES↗
DOE OSTI2026

Cooperative Agreement To Analyze variabiLity, change and predictabilitY in the earth SysTem (CATALYST)

CATALYST proposes to perform foundational coordinated research in a team-oriented collaborative effort aimed at advancing a robust understanding of modes of Earth system variability and change using models, observations and process studies. The proposed research will address the DOE/BER mission by exploring the limits to predictability, identifying fundamental underlying mechanisms, quantifying interactions among modes of variability, and discovering tipping points in the Earth system to understand the current and future impacts of these phenomena on regional and global climate. Four fundamental gaps are identified in our knowledge of the Earth system: 1) What are the limits to predictability on various timescales? 2) What are the interactions among modes of Earth system variability? 3) How may modes of Earth system variability change in response to changes in external forcing, and what are the tipping points involved with those changes? 4) How are high impact events connected to modes of Earth system variability and how may they change in the future? Related to those gaps in our knowledge, we formulate four research objectives to address those gaps using a combination of Earth system models (ESMs) and machine learning (ML) methods. Research Objective 1 (RO1) addresses the first gap above and proposes to understand modes of variability and their limits of predictability on subseasonal to decadal timescales using ESMs and ML. Research Objective 2 (RO2) addresses the second gap and proposes to use a hierarchy of models to understand relevant processes and feedbacks related to how modes of variability interact with each other. Research Objective 3 (RO3) is designed to study the third gap and proposes to examine the role of external forcings in changes of modes of Earth system variability and their interactions, and the likelihood and predictability of tipping points and irreversible changes. Research Objective 4 (RO4) will address the fourth gap and proposes to use high resolution ESMs, regionally refined models (RRMs), and ML methods to investigate the relationships between high impact events (e.g. flash droughts and precipitation extremes, atmospheric rivers (ARs), tropical cyclones (TCs), storm surge/sea level rise), the synoptic systems that produce them, and their changes related to modes of Earth system variability. The research will involve the use of the Community Earth System Model (CESM), Energy Exascale Earth System Model (E3SM), CMIP multi-model data sets, a hierarchy of simpler models, and numerous observational data sets. In the course of the proposed research, CATALYST will contribute to metrics and diagnostics that will be integrated in Coordinated Model Evaluation Capabilities (CMEC), particularly with regards to the Quasi-biennial Oscillation (QBO) and its interactions with the Madden-Julian Oscillation (MJO), high atmospheric pressure blocking, and new precipitation metrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Resolution requirements for high-speed droplet impact along a rigid wall

The mechanical and thermal loads that are produced on high-speed projectiles from a droplet are key in the projectile development.To accurately predict the pressure produced along the wall surface during impact, it is important to correctly model the air/water interface of the droplet. Novel modeling approaches such as the Phase-Field enable control of the (numerical) interface thickness. However, the interplay between the grid resolution (including potential adaptive mesh refinement) and the interface thickness is unclear. In this work, the resolution requirements for numerical simulations of a two-dimensional cylindrical water droplet impacting a rigid wall at Mach four are computationally investigated. Numerical simulations are performed using a second-order accurate scheme with adaptive mesh refinement and a Phase-Field interface model. The choice of numerical interface thickness and resolution parameters affects the maximum pressure and location of the maximum pressure along the wall. The overall maximum pressure does not converge even at high resolutions. The number of levels and thickness (number of cells across the interface) play a role in obtaining a well behaved solution.

Burrell, E.↗
DOE OSTI2026

High-Resolution South American Wind Resource Data Downscaled with Generative Machine Learning Conditioned on Near-Surface Observations

High-resolution historical wind data was developed for the entirety of South America using the innovative Super-Resolution for Renewable Resource Data (sup3r) machine learning framework. The publicly available Sup3rWind South America dataset represents a significant advancement in wind resource data generation, leveraging generative machine learning conditioned on near-surface observations from the Meteorological Assimilation Data Ingest System (MADIS) to efficiently and accurately downscale coarse reanalysis data from the European Centre for Medium-Range Weather Forecasts (ERA5). This approach produces fine-scale, spatially and temporally coherent wind and meteorological fields hundreds of times more computationally efficient than traditional numerical weather modeling methods, enabling access to high-fidelity wind information across both continental and offshore regions. Sup3rWind South America builds on the earlier Sup3rWind Ukraine dataset through improvements in model architecture and outputs conditioned on near-surface observation inputs. As with the Ukraine data release, this dataset includes wind speed, wind direction, temperature, relative humidity, and pressure at a horizontal resolution of ~2 km, representing a 15x spatial enhancement relative to the 31 km ERA5 grid. Wind speed and direction are provided at 5-minute resolution, a 12x temporal refinement compared to the hourly ERA5 data, while temperature, relative humidity, and pressure remain at hourly resolution. The data covers all years from 2005 to 2024. Before downscaling, ERA5 inputs were bias-corrected using long-term monthly means and a limited number of quality-controlled observations to align large-scale statistics with regional conditions. The resulting dataset is the first publicly available high-resolution timeseries wind record that provides full spatial coverage of South America. Model validation demonstrates strong agreement with observations across several statistical metrics, consistent with other state-of-the-art high-resolution wind resource datasets. The potential applications of Sup3rWind South America span renewable energy resource assessment, energy system modeling, and grid resilience analysis. The 20-year record and high spatial and temporal resolution support accurate estimation of long-term energy yield and the economic feasibility of potential wind development sites. Continuous coverage across both continental and offshore regions enables comprehensive site prospecting within exclusive economic zones. The 2 km, 5-minute resolution data provide the spatial and temporal variability required for power system simulation, operational planning, and regional risk assessments.

17 WIND ENERGY↗
DOE OSTI2024

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗
DOE OSTI2024

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗
DOE OSTI2025

Improved Understanding of Multicentury Greenland Ice Sheet Response to Strong Warming in the Coupled CESM2‐CISM2 With Regional Grid Refinement

The simulation of ice sheet‐climate interactions, such as surface mass balance fluxes, is sensitive to model grid resolution. Here we simulate the multi‐century evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO 2 increases by 1% yr -1 until quadrupling the pre‐industrial level and then is held fixed). A variable‐resolution (VR) grid with 1/4° regional refinement over the broader Arctic and 1° resolution elsewhere is applied to the atmosphere and land components, and the results are compared with conventional 1° lat‐lon grid simulations to investigate the impact of grid refinement. Compared with the 1° runs, the VR run features a slower rate of surface melt, especially over the western and northern GrIS, where the ice surface slopes gently toward the periphery. This difference pattern originates primarily from higher snow albedo and, thus, weaker albedo feedback in the VR run. The VR grid better captures the CISM ice sheet topography by reducing elevation discrepancies between CAM and CISM and is, therefore, less reliant on the downscaling algorithm, which is known to underestimate albedo gradients. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% less than that of the 1° runs, respectively.

Earth System Model↗
DOE OSTI2026

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗