SEARCH · Engineering Papers

Results for “Response surface 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2019

The effect of plant physiological responses to rising CO 2 on global streamflow

River flow statistics are expected to change as a result of increasing atmospheric CO 2 but uncertainty in Earth system model projections is high. While this is partly driven by changing precipitation, with well-known Earth system model uncertainties, in this work we show that the influence of plant stomatal conductance feedbacks can cause equally large changes in regional flood extremes and even act as the main control on future low latitude streamflow. Over most tropical land masses, modern climate predictions suggest that plant physiological effects will boost streamflow, overwhelming opposing effects of soil drying driven by the effects of CO 2 on atmospheric radiation, warming and rainfall redistribution. The relatively unknown uncertainties in representing eco-physiological processes must therefore be better constrained in land-surface models. To this end, we identify a distinct plant physiological fingerprint on annual peak, low and mean discharge throughout the tropics and identify river basins where physiological responses dominate radiative responses to rising CO 2 in modern climate projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗

DOE OSTI2021

Machine learning enables national assessment of wind plant controls with implications for land use

Summary As deployment of wind energy continues to expand, computationally efficient tools for predicting wind plant performance over a wide range of layout designs, technology innovations, and spatial locations are increasingly important for policy and investment decisions. We demonstrate two approaches to training a surrogate model to predict annual energy production (AEP) of parameterized wind plant layouts: one using a Gaussian process (GP) and the other using a fully convolutional neural network (FCNN). We leverage the powerful FCNN architecture by encoding wind plant design parameters and output response surface as an image. The FCNN produces more accurate results than the GP with mean absolute errors equivalent to 1% and 1.9% of plant rated power, respectively, although the GP performs well under limited training data and provides useful uncertainty information. We also evaluate a surrogate model for wake steering, enabling a nationwide assessment of the impact of plant control strategies and plant layout decisions. Across two million locations, we find that wake steering strategies boost AEP with relative gains upwards of 3%. Gains are most pronounced at sites without a dominant wind direction and where layout optimization is less fruitful. Additionally, we perform a nationwide sensitivity analysis showing that wake steering can mitigate wake losses from higher density plant layouts. Our results suggest that regions which have not been previously viable for wind deployment due to moderate wind resources are especially well enhanced by wake steering strategies that could help overcome land constraints and inflexible layout options, potentially identifying new deployment opportunities.

17 WIND ENERGY↗

DOE OSTI2026

Wind Conditions for Solar Energy Facility Design Based on Field Measurements at Nevada Solar One

This study characterises near-surface wind conditions at the Nevada Solar One (NSO) concentrated solar power plant using two years of measurements from a 15-m mast located 30 m west of the solar arrays. The surrounding terrain is fairly flat within 3 km, but nearby solar photovoltaic (PV) and concentrated solar power (CSP) installations introduce strong directional wind inhomogeneity, which motivates a sector-based analysis that distinguishes undisturbed inflow from waked flow. The analysis uses data from three sonic anemometers to assess whether turbulence in the surface layer can be represented as Gaussian and whether standard surface-layer scaling remains valid below 10 m height. The roughness length and displacement height are estimated for twelve directional sectors under near-neutral wind conditions, and all three velocity components show departures from Gaussian behaviour. Spectral analysis indicates that the turbulence contains low-frequency fluctuations with amplitudes higher than those predicted by standard models such as the Kaimal model or the Simiu and Scanlan model. The study introduces a new turbulence spectral model based on the eddy-surface-layer framework to better represent near-surface wind turbulence. The analysis further quantifies the influence of the surrounding power plants on mean wind profiles, one-point spectra, and vertical coherence for all velocity components. The solar collector arrays located 30 m east of the mast significantly modify the turbulence structure, whereas PV installations situated a few hundred metres away have only a limited influence on the turbulence measured at the mast. These findings highlight the limitations of conventional turbulence models for near-surface wind over solar infrastructure and provide improved parameterisations to support the structural design. The implications for the structural response of parabolic troughs and PV modules to such turbulence remain uncertain and require further investigation.

14 SOLAR ENERGY↗

DOE OSTI2021

Doppler-Lidar Evaluation of HRRR-Model Skill at Simulating Summertime Wind Regimes in the Columbia River Basin during WFIP2

Complex-terrain locations often have repeatable near-surface wind patterns, such as synoptic gap flows and local thermally forced flows. An example is the Columbia River Valley in east-central Oregon-Washington, a significant wind-energy-generation region and the site of the Second Wind-Forecast Improvement Project (WFIP2). Data from three Doppler lidars deployed during WFIP2 define and characterize summertime wind regimes and their large-scale contexts, and provide insight into NWP model errors by examining differences in the ability of a model [NOAA’s High-Resolution Rapid-Refresh (HRRR-version1)] to forecast wind-speed profiles for different regimes. Seven regimes were identified based on daily time series of the lidar-measured rotor-layer winds, which then suggested two broad categories. First, in three regimes the primary dynamic forcing was the large-scale pressure gradient. Second, in two regimes the dominant forcing was the diurnal heating-cooling cycle (regional sea-breeze-type dynamics), including the marine intrusion previously described, which generates strong nocturnal winds over the region. The other two included a hybrid regime and a non-conforming regime. For the large-scale pressure-gradient regimes, HRRR had wind-speed biases of ~1 m s -1 and RMSEs of 2-3 m s -1 . Errors were much larger for the thermally forced regimes, owing to the premature demise of the strong nocturnal flow in HRRR. Thus, the more dominant the role of surface heating in generating the flow, the larger the errors. Major errors could result from surface heating of the atmosphere, boundary-layer responses to that heating, and associated terrain interactions. Finally, measurement/modeling research programs should be aimed at determining which modeled processes produce the largest errors, so those processes can be improved and errors reduced.

17 WIND ENERGY↗

DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

DOE OSTI2021

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗

DOE OSTI2021

Updated respiration routines alter spatio-temporal patterns of carbon cycling in a global land surface model

Abstract We updated the routines used to estimate leaf maintenance respiration (MR) in the Energy Land Model (ELM) using a comprehensive global respiration data base. The updated algorithm includes a temperature acclimating base rate, an updated instantaneous temperature response, and new plant functional type specific parameters. The updated MR algorithm resulted in a very large increase in global MR of 16.1 Pg (38%), but the signal was not geographically uniform. The increase was concentrated in the tropics and humid warm-temperate forests. The increase in MR led to large but proportionally smaller decreases in global net primary production (19%) and in average global leaf area index (15%). The effect on global gross primary production (GPP) was a more modest 5.7 Pg (4%). A detailed site level analysis also demonstrated a wide range of effects the updated algorithm can have on the seasonal cycle of GPP. Output from the updated and old models did not differ markedly in how closely they matched a suite of benchmarks. Given the substantial impact on the land surface carbon cycle, a neutral influence on model benchmarks, and better alignment with empirical evidence, an MR algorithm similar to the one presented here should be adopted into ELM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Biogeochemical Equilibrium Responses to Maximal Productivity in High Nutrient Low Chlorophyll Regions

Whether ocean iron fertilization in high nutrient low chlorophyll (HNLC) regions is a practical geoengineering strategy to combat climate warming has been debated because the fate of absorbed carbon by the ocean remains unclear. Here, we use an optimized oceanic biogeochemical inverse model to explore the equilibrium responses to increased productivity in major HNLC regions. We find that an overall downward shift of nutrients and carbon from the surface and intermediate to the deep waters with maximal productivity in the HNLC Southern Ocean (SO), North Pacific (NP) and eastern equatorial Pacific (EP). The efficiency of soft tissue pump is increased by 37.8%, 2.7% and 5.4% in the SO, EP and NP simulation while atmospheric CO 2 concentration is drawn down by 45.1, 3.4 and 6.7 ppm, respectively. Meanwhile, global mean oxygen is decreased by 21.1%, 1.3% and 2.4% in the SO, EP and NP simulation, accompanied by an expansion of oxygen minimum zones by 400%, 26% and 63%, respectively. Our study indicates that iron fertilization has the potential to be an effective geoengineering approach to mitigate climate warming but also has a strong impact on the development of ocean hypoxia.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

E3SM‐GCAM: A Synchronously Coupled Human Component in the E3SM Earth System Model Enables Novel Human‐Earth Feedback Research

Abstract Modeling human‐environment feedbacks is critical for assessing the effectiveness of climate change mitigation and adaptation strategies under a changing climate. The Energy Exascale Earth System Model (E3SM) now includes a human component, with the Global Change Analysis Model (GCAM) at its core, that is synchronously coupled with the land and atmosphere components through the E3SM coupling software. Terrestrial productivity is passed from E3SM to GCAM to make climate‐responsive land use and CO 2 emission projections for the next 5‐year period, which are interpolated and passed to E3SM annually. Key variables affected by the incorporation of these feedbacks include land use/cover change, crop prices, terrestrial carbon, local surface temperature, and climate extremes. Regional differences are more pronounced than global differences because the effects are driven primarily by differences in land use. This novel system enables a new type of scenario development and provides a powerful modeling framework that facilitates the addition of other feedbacks between these models. This system has the potential to explore how human responses to climate change impacts in a variety of sectors, including heating/cooling energy demand, water management, and energy production, may alter emissions trajectories and Earth system changes. Plain Language Summary Earth system models help us understand how humans are changing the climate. Currently, these models do not include human systems, so predetermined greenhouse gas, aerosol, and land use change data are input to these models. These data do not reflect human responses to changes projected by Earth system models. We have added a human component to an Earth system model to represent human responses to environmental change and calculate corresponding greenhouse gas and land use change data instead of using predetermined data. Including this human component changes projections of land use, land carbon storage, and regional climate. Key Points We have incorporated a novel, advanceable human component in an Earth system model to simulate human‐Earth feedbacks Including terrestrial productivity feedbacks from the Earth to the human systems affects land change, crop prices, carbon, and climate Regional effects of including terrestrial productivity feedbacks are greater than global effects because land change is the main driver

E3SM↗

DOE OSTI2026

Indirect climate impacts of the Hunga eruption

Injection of sulfur and water vapour by the Hunga volcanic eruption significantly altered chemical composition and radiative budget of the stratosphere. Yet, whether the eruption could also affect surface climate, especially via indirect pathways, remains poorly understood. Here we investigate these effects using large ensembles of simulations with the CESM2(WACCM6) Earth system model, incorporating interactive chemistry and aerosols in both coupled ocean and atmosphere-only configurations. We find some statistically significant extratropical regional climate responses to the eruption driven by circulation changes; these are partially linked to the modulation of El Niño–Southern Oscillation, and its associated teleconnections, and to perturbations of the stratospheric polar vortex in both hemispheres. The stratospheric anomalies affect surface climate through modulating the North Atlantic Oscillation in the Northern Hemisphere (up to three boreal winters following the eruption) and the Southern Annular Mode in the Southern Hemisphere in late 2023. The latter is partly related to a concurrent reduction in Antarctic ozone, as increased stratospheric aerosols and water vapor reach the polar vortex. Our study suggests that the eruption could have had a non-negligible influence on regional surface climate, and discusses the mechanisms via which such an influence could occur. However, the results also highlight that this forcing is relatively weak compared to interannual variability, and is subject to model uncertainties in the representation of key processes. More research is thus needed before definitive statements on the role of the eruption in contributing to surface climate and weather events in the following years are made.

Geosciences↗

DOE OSTI2024

Improved subseasonal-to-seasonal precipitation prediction of climate models with nudging approach for better initialization of Tibetan Plateau-Rocky Mountain Circumglobal wave train and land surface conditions

Abstract Reliable subseasonal-to-seasonal (S2S) precipitation prediction is highly desired due to the great socioeconomical implications, yet it remains one of the most challenging topics in the weather/climate prediction research area. As part of the Impact of Initialized Land Temperature and Snowpack on Sub-seasonal to Seasonal Prediction (LS4P) project of the Global Energy and Water Exchanges (GEWEX) program, twenty-one climate models follow the LS4P protocol to quantify the impact of the Tibetan Plateau (TP) land surface temperature/subsurface temperature (LST/SUBT) springtime anomalies on the global summertime precipitation. We find that nudging towards reanalysis winds is crucial for climate models to generate atmosphere and land surface initial conditions close to observations, which is necessary for meaningful S2S applications. Simulations with nudged initial conditions can better capture the summer precipitation responses to the imposed TP LST/SUBT spring anomalies at hotspot regions all over the world. Further analyses show that the enhanced S2S prediction skill is largely attributable to the substantially improved initialization of the Tibetan Plateau-Rocky Mountain Circumglobal (TRC) wave train pattern in the atmosphere. This study highlights the important role that initial condition plays in the S2S prediction and suggests that data assimilation technique (e.g., nudging) should be adopted to initialize climate models to improve their S2S prediction.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2022

A directional microcrack damage theory for brittle solids based on continuous hyperplasticity

We present a novel microcrack-damage theory for brittle solids under compression. Instead of using internal variables like zeroth, second or fourth rank damage tensors, the state of material damage is represented by an internal function that encapsulates the information of direction, density and size of microcracks. Just like other internal variables, the evolution of this state function must obey the second law of thermodynamics for arbitrary loading paths. This is done by casting the model in the framework of continuous hyperplasticity and enforcing a non-negative dissipation rate functional. The proposed framework offers predictions on the continuous evolution of microcrack density and the induced material anisotropy along with the macroscopic stress-strain curves. The use of continuous damage function grants the model significantly enhanced resolution in characterizing the direction-dependent response of cracked solids compared to classical models that are based on damage tensors. Two scenarios are considered in developing the theory, one assumes frictionless cracks and the other incorporates friction between crack surfaces. The results highlight that inelasticity, pressure dependence, and loading-unloading hysteresis exhibited by brittle solids are natural consequences of frictional microcracks. The proposed theory offers a generic and versatile framework to upscale micromechanical processes operating at individual crack scale to explain the macroscopic behavior of cracked solids.

Materials Science↗

DOE OSTI2026

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗