Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial variability”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modulating charge transport via 2 MeV He + irradiation in VO 2

Vanadium dioxide (VO 2 ) is of interest for adaptive electronic applications such as neuromorphic neuristor devices and variable emissivity or tunable thermal control materials, thanks to its key property—a metal–insulator transition (MIT) at 68 °C that is accompanied by a dramatic change in electrical and optical properties. To improve performance in these roles, it is critical to develop approaches to engineer transport properties and the MIT behavior. While many documented techniques exist to modulate the MIT and film resistivities via lattice strain and chemical doping, less is known about the effects of ion irradiation on the intrinsic properties of VO 2 , despite the ability to control the spatial distribution of irradiation beams and the prevalence of high energy ion implantation in the semiconductor industry. The impact of irradiation of different acceleration energies on the responses of VO 2 is of specific interest, as charged particle energy generally impacts both the resulting defect profile and corresponding transport behavior. Here, we demonstrate that 2 MeV He ions at equivalent calculated displacements per atom, in two different types of films, can create remarkable changes to the nature of charge transport in VO 2 , especially in the low-temperature insulating phase. Simulation of resulting changes in electrical conductivity reveals that He ion irradiation offers a strategy to increase both oscillation frequency and the signal transmission. These results provide insights into the intentional design of defect populations to modulate transport for neuromorphic VO 2 devices.

Gurrola, Rebeca M. [Texas A & M Univ., College Sta↗
DOE OSTI2021

Correlative analysis of structure and chemistry of Li x FePO 4 platelets using 4D-STEM and X-ray ptychography

Lithium iron phosphate (Li x FePO 4 ), a cathode material used in rechargeable Li-ion batteries, phase separates upon de/lithiation under equilibrium. The interfacial structure and chemistry within these cathode materials affects Li-ion transport, and therefore battery performance. Here, correlative imaging of Li x FePO 4 was performed using four-dimensional scanning transmission electron microscopy (4D-STEM), scanning transmission X-ray microscopy (STXM), and X-ray ptychography in order to analyze the local structure and chemistry of the same particle set. Over 50,000 diffraction patterns from 10 particles provided measurements of both structure and chemistry at a nanoscale spatial resolution (16.6–49.5 nm) over wide (several micron) fields-of-view with statistical robustness. Li x FePO 4 particles at varying stages of delithiation were measured to examine the evolution of structure and chemistry as a function of delithiation. In lithiated and delithiated particles, local variations were observed in the degree of lithiation even while local lattice structures remained comparatively constant, and calculation of linear coefficients of chemical expansion suggest pinning of the lattice structures in these populations. Partially delithiated particles displayed broadly core–shell-like structures, however, with highly variable behavior both locally and per individual particle that exhibited distinctive intermediate regions at the interface between phases, and pockets within the lithiated core that correspond to FePO 4 in structure and chemistry. The results provide insight into the Li x FePO 4 system, subtleties in the scope and applicability of Vegard’s law (linear lattice parameter-composition behavior) under local versus global measurements, and demonstrate a powerful new combination of experimental and analytical modalities for bridging the crucial gap between local and statistical characterization.

25 ENERGY STORAGE↗
DOE OSTI2024

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG 1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG 1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~200 anti-GCGR IgG 1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2–10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Disparities in the air quality monitoring stations and PM₂.₅ in Chicago’s air quality landscape

Fine particulate matter (PM₂.₅) poses significant public and environmental health risks in urban areas. Chicago’s dense industry and traffic create variable air quality, yet monitoring is unevenly distributed, resulting in undersampling of air quality data in some city areas. This study applied a hybrid approach using GIS-based kernel density mapping, interpolation modeling (IDW, Spline, Kriging) of USEPA monitoring data, multi-scale temporal trend analyses (hourly to annual), and ESDA. Accordingly, the density surface showed that monitors are concentrated in the affluent north, northwest, and southwest sides of Chicago (up to ~ 0.07 stations per sq mile), while the south and southeast regions, with predominantly minority communities, have virtually no coverage. Overall, citywide coverage is minimal (~ 4–5 monitors total; ~0.02 per sq mile; ≈1 per 600,000 residents). Temporal analyses showed that the city’s mean annual PM₂.₅ (~ 10.8 µg/m³) exceeds USEPA/WHO standards (9 µg/m³), with summer means (~ 17.1 µg/m³) significantly higher than other seasons. Diurnally, a clear pattern was observed, with PM₂.₅ concentrations peaking overnight (00:00–03:00) and during the morning rush hours, and dipping during midday to late afternoon. Spatial distribution of PM₂.₅ identified hotspots near O’Hare Airport, the downtown Loop area, and south-side neighborhoods, contrasting with lower concentrations on the north side, revealing Chicago’s socioeconomic divides and resulting environmental inequities. The findings underscore the need for expanded monitoring and targeted interventions in under-monitored, high-pollution communities to advance equitable community health.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessing the Effectiveness of Large-Scale Environmental Restoration: Challenges and Opportunities

A recent National Academies consensus report addresses monitoring and assessment of cumulative effects of large-scale and multiple restoration projects within the context of long-term environmental change. Fines and penalties from the Deepwater Horizon oil spill in the Gulf of Mexico (GoM) have supported hundreds of restoration projects at spatial scales not often possible in the past. Here, in this paper, the report committee members and staff provide personal reflections from our time working on the study. We found that gaps in data collection, issues with data accessibility, and a lack of synthesis and analysis are hindering the ability to answer a basic question: What are the impacts of these many restoration efforts on improving ecosystem health and productivity in the GoM at the regional and Gulf wide scale? Restoration efforts are occurring in environments where many trends are changing and exhibiting higher variability than in the past, suggesting that previously successful restoration practices may no longer be adequate to compensate for the effects of environmental changes and variability. Our proposed approach to these challenges includes employing emerging monitoring technologies; using conceptual models; devising an adaptive management framework; rethinking restoration outcome goals; assessing cumulative effects; and undertaking rigorous synthesis and analysis of existing information on long-term environmental trends and restoration efforts. Restoration scientists and practitioners working in the GoM have an unprecedented opportunity to demonstrate large-scale environmental recovery if advances in monitoring, synthesis, assessment, and action are taken quickly. We are cautiously optimistic that, with mid-course adjustments, continued progress toward large-scale environmental recovery is possible.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulating Sea-Ice Deformation in Viscous-Plastic Sea-Ice Models With CD-Grids

Linear kinematic features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. In viscous-plastic (VP) sea-ice models, the simulation of LKFs depends on several factors such as the grid resolution, the numerical solver convergence, and the placement of the variables on the mesh. In this study, we compare two recently proposed discretization with a CD-grid placement with respect to their ability to reproduce LKFs. The first (CD1) is based on a nonconforming finite element discretization, whereas the second (CD2) uses a conforming subgrid discretization. To analyze their resolution properties, we evaluate runs from different models (e.g., FESOM, MPAS) on a benchmark problem using quadrilateral, hexagonal and triangular meshes. Our findings show that the CD1 setup simulates more deformation structure than the CD2 setup. This highlights the importance of the type of spatial discretization for the simulation of LKFs. Due to the higher number of degrees of freedom, both CD-grids resolve more LKFs than traditional A, B, and C-grids at fixed mesh level. This is an advantage of the CD-grid approach, as high spatial mesh resolution is needed in VP sea-ice models to simulate LKFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

reV (The Renewable Energy Potential Model - Open Source) [SWR-21-59, SWR-20-20 and SWR-17-34]

The Renewable Energy Potential (reV) model is a platform for the detailed assessment of renewable energy resources and their geospatial intersection with grid infrastructure and land use characteristics. The reV model currently supports photovoltaic (PV), concentrating solar power (CSP), and land-based wind turbine technologies. Modules in the reV framework function at different spatial and temporal resolutions, allowing for the assessment of resource potential, technical potential, and supply curves at varying levels of detail. The platform runs on the National Renewable Energy Laboratory’s (NREL’s) high-performance computing system, providing scalable and efficient performance from a single location up to a continent, for a single year or decades of time-series resource data. Coupled with NREL’s System Advisor Model (SAM), reV supports resource assessments from 5-minute to hourly temporal resolutions and supports the analysis of long-term (i.e., year-on-year) variability of renewable generation (e.g., interannual variability and exceedance probabilities).

Maclaurin, Galen↗
DOE OSTI2021

Wastewater Disposal Has Not Significantly Altered the Regional Stress State in Southern Kansas

Wastewater disposal is primarily responsible for the increased seismicity rate since ~2013 in southern Kansas. Previous work that used shear-wave splitting (SWS) in southern Kansas interpreted an ~90° temporal rotation in the fast polarization direction and attributed it to increased pore pressures resulting from fluid injection. However, this interpreted rotation coincided with a change in the stations used to make the SWS measurements. We investigate the temporal variability of fast azimuths in southern Kansas by making SWS measurements on earthquake families with similar source–receiver paths recorded on a stable local seismic network. We select high-quality SWS measurements by investigating the stability of results across 65 different frequency bands between 0.5 and 15 Hz. We find that the fast polarization direction in southern Kansas is relatively constant with an average east-northeast (~N79°E) orientation between 2014 and 2017. Our fast polarization measurements are primarily a reflection of the maximum principal horizontal stress direction (SHmax). We observe a slight spatial change in SHmax to the northeast (~N55°E) near the Nemaha ridge in Oklahoma. However, we do not observe any significant temporal rotation of SHmax or variation in delay time (i.e., crack density) in southern Kansas, contrary to the earlier study. The previously interpreted ~90° rotation may either be a reflection of a very local stress change or a misinterpretation of SWS results potentially due to the use of inconsistent source–receiver paths. Our SWS measurements cover the period of peak wastewater disposal and seismicity rates and suggest an absence of significant temporal rotations in the local anisotropy and stress orientations associated with wastewater disposal.

Geochemistry & Geophysics↗
DOE OSTI2025

Microstructure-sensitive mechanical behavior of an additively manufactured psuedoelastic shape memory alloy

The additive manufacturing of shape memory alloys into complex geometries enables fabrication of advanced functional systems across a variety of fields and domains. This work presents results focused on the mechanical behavior of additively manufactured shape memory pseudoelastic NiTi. The deformation induced solid state phase transformation from austenite to martensite allows this system to accommodate large recoverable strains. This deformation behavior is fundamentally driven by crystal-scale transformation physics. Laser powder bed fusion processing reveals that the resulting microstructure, both grain morphology and crystallographic texture, is strongly dependent on the manufacturing processing history. Exhaustive mechanical testing demonstrates that these microstructural factors strongly impact both tensile and cyclic stress–strain behavior. Cyclic dissipative behavior, however, is similar across all tested microstructures following an initial transient period. Remarkably, analysis of spatial strain fields during tensile loading reveals two distinctly different localization “modes”. The first is initiation of localized deformation bands which continuously propagate through the tensile bar during loading. In the second mode localization is observed but lacks propagation; instead additional localization cites nucleate during subsequent loading. The latter phenomena is suspected to be driven by grain-scale deformation physics as the localized band morphologies coincide with grain morphologies. These phenomena strongly impact the resulting aggregate stress–strain behavior. Hence, manufacturers and designers of psuedoelastic functional components must at the very least consider the potential variability in properties when considering additive manufacturing processing. More ideally the process–structure–property relations can be used to further tailor and optimize final functional performance.

Additive manufacturing↗
DOE OSTI2022

A geo-visual analysis for exploring the socioeconomic benefits of the heating electrification using geothermal energy

In parallel to population growth and climate change, the rapid pace of urbanization worldwide has led to an enormous increase in energy demand and costs in urban areas. The subsequent energy burden has become an increasing concern for many households in the U.S. Previous studies have revealed that geothermal resources can effectively lower the electricity demand and carbon emissions in large cities. In this paper, we focus on the socioeconomic impacts of geothermal energy on urban systems by presenting an interactive visual analytics dashboard. The dashboard allows urban planners to spatially examine geothermal energy's practical benefits on energy affordability, urban livability, and resilience across the U.S. We compiled a list of socioeconomic metrics by integrating the simulation results from multiple geothermal and building models with multi-domain urban datasets (socioeconomic, demographic, and electricity utility). These metrics are created to characterize the benefits of the heating electrification of buildings using Geothermal Heat Pumps (GHP) for lowering the energy burden of middle-and low-income households nationwide. The visual dashboard employs a combination of multivariate, glyph-based, and geospatial visualization to reveal the variability and patterns in our metrics. We present a pilot study to demonstrate the GHPs' potential as a renewable and affordable solution for increasing the economic and energy grid resilience in U.S cities.

Xu, Haowen↗
DOE OSTI2023

Biogeophysical Effects of Land-Use and Land-Cover Change Not Detectable in Warmest Month

Land-use and land-cover changes (hereafter simply “land use”) alter climates biogeophysically by affecting surface fluxes of energy and water. Yet, near-surface temperature responses to land use across observational versus model-based studies and spatial-temporal scales can be inconsistent. Here we assess the prevalence of the historical land use signal of daily maximum temperatures averaged over the warmest month of the year (t LU ) using regularized optimal fingerprinting for detection and attribution. We use observations from the Climatic Research Unit and Berkeley Earth alongside historical simulations with and without land use from phase 6 of the Coupled Model Intercomparison Project to reconstruct an experiment representing the effects of land use on climate. To assess the signal of land use at spatially resolved continental and global scales, we aggregate all input data across reference regions and continents, respectively. At both scales, land use does not comprise a significantly detectable set of forcings for two of four Earth system models and their multimodel mean. Furthermore, using a principal component analysis, we find that t LU is mostly composed of the nonlocal effects of land use rather than its local effects. These findings show that, at scales relevant for climate attribution, uncertainties in Earth system model representations of land use are too high relative to the effects of internal variability to confidently assess land use.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Active layer depth and soil properties impact specific leaf area variation and ecosystem productivity in a boreal forest

Specific leaf area (SLA, leaf area per unit dry mass) is a key canopy structural characteristic, a measure of photosynthetic capacity, and an important input into many terrestrial process models. Although many studies have examined SLA variation, relatively few data exist from high latitude, climate-sensitive permafrost regions. We measured SLA and soil and topographic properties across a boreal forest permafrost transition, in which dominant tree species changed as permafrost deepened from 54 to >150 cm over 75 m hillslope transects in Caribou-Poker Creeks Research Watershed, Alaska. We characterized both linear and threshold relationships between topographic and edaphic variables and SLA and developed a conceptual model of these relationships. We found that the depth of the soil active layer above permafrost was significantly and positively correlated with SLA for both coniferous and deciduous boreal tree species. Intraspecific SLA variation was associated with a fivefold increase in net primary production, suggesting that changes in active layer depth due to permafrost thaw could strongly influence ecosystem productivity. While this is an exploratory study to begin understanding SLA variation in a non-contiguous permafrost system, our results indicate the need for more extensive evaluation across larger spatial domains. These empirical relationships and associated uncertainty can be incorporated into ecosystem models that use dynamic traits, improving our ability to predict ecosystem-level carbon cycling responses to ongoing climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Data from: “Bald Cypress (Taxodium distichum) Knees Are Methane Sources Controlled by Geomorphology, Climate, and Hydrologic Extremes”

This dataset is associated with the manuscript “Bald Cypress (Taxodium distichum) Knees Are Methane Sources Controlled by Geomorphology, Climate, and Hydrologic Extremes”. Bald cypress “knees” (aboveground woody roots) have been shown to contribute to wetland methane (CH4) efflux, with large variation within and between studies. To explain this variation, we investigated spatial (i.e., across knee surface, within sites, between sites) and temporal dynamics of CH4 fluxes from knees. Methane fluxes were collected from September 2022 to August 2024 at three locations in western Kentucky, USA, within the Mississippi Alluvial Valley: a main channel (semi-permanently flooded), side channel (seasonally flooded), and reservoir edge (artificially flooded). Knee CH4 fluxes (“Ross_et_al_Knee_Flux_Data.csv”) were measured from multiple heights on knees (20, 40, and 60 cm) of various sizes (knee straight height ranged from 24 to 93 cm) using a LiCOR LI-7810 CH4/CO2/H2O Trace Gas Analyzer. The dataset also includes environmental variables collected with each knee measurement, including water level adjusted for knee-to-knee elevational differences, subsurface and air temperature, and humidity. Soil CH4 fluxes (“Ross_et_al_Soil_Flux_Data.csv”) were also collected adjacent to knees (starting in April 2023) when water levels didn’t overtop soil collars, using a LiCOR Smart Chamber and calculated in SoilFluxPro software. The soil flux dataset includes associated variables collected by the Smart Chamber. Three separate files (“*_Water_Level.csv”) are included for water level and subsurface temperature data collected at each site using HOBO U20L barometric pressure loggers. Each file type (knee flux, soil flux, water level) has an associated data dictionary (“*_dd.csv”). For specifics on methodology used and calculations, see the associated manuscript. The R script includes code used for figures and analyses reported in the manuscript.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advancing spatiotemporal forecasts of CO 2 plume migration using deep learning networks with transfer learning and interpretation analysis

Accurate and timely forecasts of CO 2 plume distribution throughout the injection and post-injection phases are crucial for detecting plume migration, assessing leakage risks, and supporting operational decisions in geologic carbon storage (GCS). Current convolutional neural network-based approaches primarily focus on spatial information and overlook temporal dependencies in plume distributions, thus limiting their ability to capture dynamic movement effects and provide accurate predictions of plume migration. In this work, we propose two deep learning models, Auto-Encoder (AE)-LSTM and Encoder-Decoder (ED)-ConvLSTM, each uniquely designed to capture both spatial and temporal features. We apply the proposed methods to forecast the dynamic distribution of CO 2 plumes based on 108 reservoir simulations over a 30-year injection and a 30-year post-injection period. The results indicate that the ED-ConvLSTM model outperforms the AE-LSTM model in accurately predicting the spatiotemporal dynamics of CO 2 plume migration, achieving R 2 values above 0.99. To provide a deeper understanding of these model predictions, we employ a gradient-based explanation method on the trained models. This approach provides insights into the influence of input variables on plume migration forecasts and uncovers the underlying prediction mechanisms of the proposed models. Furthermore, we introduce a transfer learning technique, enabling fast and accurate plume migration forecasting in the post-injection phase by leveraging the trained model during the injection phase. This reduces the necessity for extensive data collection or re-training. In conclusion, the methods proposed in our work enhances the performance and interpretability of CO 2 plume migration forecasts, thereby facilitating informed decision-making throughout the entire lifecycle of GCS applications.

58 GEOSCIENCES↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2024

A Comprehensive Stress Drop Map From Trench to Depth in the Northern Chilean Subduction Zone

We compute stress drops for earthquakes in Northern Chile recorded between 2007 and 2021. By applying two analysis techniques, (a) the spectral ratio (SR) method and (b) the spectral decomposition (SDC) method, a stress drop map for the subduction zone consisting of 51,510 stress drop values is produced. We build an extended set of empirical Green’s functions (EGF) for the SR method by systematic template matching. Outputs are used to compare with results from the SDC approach, where we apply cell-wise obtained global EGF's to compensate for the structural heterogeneity of the subduction zone. We find a good consistency of results of the two methods. The increased spatial coverage and quantity of stress drop estimates from the SDC method facilitate a consistent stress drop mapping of the different seismotectonic domains. Albeit only small differences of median stress drop, strike-perpendicular depth sections clearly reveal systematic variations, with earthquakes at different seismotectonic locations exhibiting distinct values. In particular, interface seismicity is characterized by the lowest observed median value, whereas upper plate earthquakes show noticeably higher stress drop values. Intermediate depth earthquakes show comparatively high average stress drop and a rather strong depth-dependent increase of median stress drop. Additionally, we observe spatio-temporal variability of stress drops related to the occurrence of the two megathrust earthquakes in the study region. The presented study is the first coherent large scale 3D stress drop mapping of the Northern Chilean subduction zone. It provides an important component for further detailed analysis of the physics of earthquake ruptures.

58 GEOSCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗