Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volcanic”

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 253 records · Page 14

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗

IC Project: w19_hossfault “Modelling of stick-slip behavior in sheared granular fault gouge & nonlinear elasticity behavior in cracked solid”

Low frequency earthquakes, non-volcanic tremor, and acoustic emissions are examples of weak seismic signals that may help detect major seismic events, i.e., earthquakes. In a laboratory setting, in which stick-slip events are simulated, acoustic emissions are detected far from the stick-slip events. In both, field and laboratory scale cases, the acoustic emission signals are sourced or detected in a volume remote from the volume that spawns the earthquake. Therefore, it is relevant to establish the causal relationship between signals detected on passive, remote monitors and the dynamics of the elastic structures that launch important seismic events. The work conducted under this IC allocation allowed us to utilize a numerical model that let us follow this causality, i.e., examine and connect the dynamics in a granular system (fault gouge), to signals detected on passive remote monitors. It was demonstrated that stress chains are key in the dynamics of granular systems and that their evolution is the source of the acoustic emission.

58 GEOSCIENCES↗

Advancing Measurements and Understanding of the Rate and Structure of Atmospheric Warming

The rate and geographic pattern of atmospheric warming are key indicators of historical climate change and play a role in modulating future changes in Earth’s climate. The past and future evolution of atmospheric temperature is simulated using global climate models while observational estimates of past warming are derived from satellite microwave measurements. Individual climate models simulate widely varying rates of past and future atmospheric warming and most model simulations exhibit greater tropospheric (lowest ~10 km of atmosphere) temperature change than satellite observations between 1979 and 2020. This project examined intermodel differences in the pattern of atmospheric warming and how these differences influence climate feedbacks that can amplify or damp the rate of global surface warming. A key result is that climate model representation of the current climate can influence climate feedbacks and the simulation of future changes in climate. Another focus of this project was to analyze how models respond to different input datasets, such as volcanic aerosols or sea surface temperature (in atmosphere-only simulations). We find that different prescribed inputs can affect simulated changes in atmospheric temperature, even though “input uncertainty” is often unconsidered in model-observational comparisons. A final focus of this research project was to consider the influence of natural internal climate variability on satellite era changes in climate. Although the Earth is warming substantially due to anthropogenic emissions of greenhouse gases, the observed warming rate can be modulated by natural variations in the Earth’s climate. We find that natural climate variability has slowed the rate of tropical tropospheric warming, which explains model-satellite differences in the rate of tropospheric warming.

58 GEOSCIENCES↗

Acid and heat stable enzymes for economical fuel production from biomass

CinderBio’s vision is to harness three billion years of biology to transform many industrial practices into sustainable processes with superior performance. Environments such as volcanic hot springs have evolved extremophilic microbes and enzymes that are naturally heat and acid stable. CinderBio was founded based on technologies developed at Lawrence Berkeley National Laboratory. Our core technologies are the molecular tools and know how to efficiently express extremophilic genes in an extremophilic host such that they retain their natural heat and acid stability. During this project, we tested our enzymes at the Advanced Biofuels and Bioproducts Process Development Unit (ABPDU), automated our fermentation and downstream processing, and built a comprehensive technoeconomic assessment (TEA) model for our manufacturing costs.

09 BIOMASS FUELS↗

Shocked Transport: Experiments to Study Dynamic Particle Behavior under Varying Volume Fraction Conditions

Shocks interacting with densely-packed porous media is common in volcanic eruptions or in engineering applications such as ordinance. Large-scale, exploratory experiments have been performed by Frost and collaborators examining the behavior of a bed of particles subjected to a detonation wave. The detonation replaces the central explosive with a high-density, high-temperature gas that propagates outward. A shock wave travels through the particle bed, followed by a contact discontinuity between the shocked air and the explosive products. In the early times, the detonation wave causes many of the particles to collide and compact with each other. As the particles disperse, the experiments demonstrate a series of large-scale jets or fingers that form in the advancing particle front.

42 ENGINEERING↗

Discovering Hidden Geothermal Signatures using Unsupervised Machine Learning

Discovering hidden geothermal resources is a very challenging task. It requires the mining of large datasets, including various diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used Play Fairway Analysis (PFA) typically relies on subject-matter expertise to analyze site or regional data to estimate geothermal conditions and prospectivity. Here, we demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset of Southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. However, most of these systems are poorly characterized because of insufficient existing data and limited past explorative studies. This study aims to discover hidden patterns and relationships in the SWNM geothermal dataset to better understand regional hydrothermal conditions. This is achieved by applying an unsupervised machine learning algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden (latent) signatures characterizing datasets, (2) the optimal number of these signatures, (3) dominant data attributes associated with each signature, and (4) spatial distribution of the extracted signatures. Here, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, geothermal attributes at 44 locations in SWNM. NMFk successfully finds data patterns and identifies the spatial associations of hydrothermal signatures with the four physiographic provinces in SWNM (Colorado Plateau, Volcanic Field, Basin and Range, and the Rio Grande rift). The algorithm identified up to 5 hydrothermal signatures in the SWNM datasets that differentiate between low- and medium-temperature hydrothermal systems in different provinces. Also, the algorithm identifies two medium-temperature hydrothermal systems in SWNM that require further exploration for geothermal resource development. Based on our analyses, 12 of the attributes are important to identify medium-temperature hydrothermal systems, and the remaining six attributes are critical to characterize low-temperature hydrothermal systems. Based on the obtained results, we identify potential physiographic provinces for further exploration to characterize them as geothermal resources. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored areas.

58 GEOSCIENCES↗

Multicontinuum Flow Models for Assessing Two-Phase Flow in Containment Science

We present a new pre-processor tool written in Python that creates multicontinuum meshes for PFLOTRAN to simulate two-phase flow and transport in both the fracture and matrix continua. We discuss the multicontinuum modeling approach to simulate potentially mobile water and gas in the fractured volcanic tuffs at Aqueduct Mesa, at the Nevada National Security Site.

58 GEOSCIENCES↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

LL18-Signal Analysis-NDD2Ad (Final Report)

Under favorable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially-distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection, and scaling with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning, and possible application to studies of repeating events in induced, volcanic, and glacial seismicity. Our results place limits on the widely-held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗

Combined Cycle Integrated Thermal Energy Storage “CiTES” (Final Scientific/Technical Report)

The Phase I of this project confirmed the technical feasibility of a Combined Cycle integrated Thermal Energy Storage “CiTES” system, calculated the key performance parameters like power efficiency and costs, and proved its commercial value with full-year simulations for several US electricity markets with high degree of variable renewable generation and volatile hourly electricity prices. The core element of this project is the Electro Thermal Energy Storage (ETES) technology from Siemens Gamesa Renewable Energy GmbH, using thermally stable and inexpensive volcanic rocks as storage material and air as heat transfer medium. This technology is backed by more than 10 years of experience and a 440MMBTU (130MWh-th) pilot plant in Hamburg, Germany, which is in operation since 2019. The integration of this thermal storage in an existing combined cycle power plant (CCPP) is typical power plant technology without any major technology risks. It allows the storage of inexpensive renewable energy during times of surplus renewable generation and the discharge of this energy in times of high energy demand when the fossil plant is in operation. This supplements the fossil power generation with CO2-emission-free energy. The secondary effect of the CiTES system is that a small part of the stored thermal energy is used to keep the heat recovery steam generator (HRSG) and steam turbine (ST) of the combined cycle power plant in hot and ready-to-start condition. This enables the plant to start rapidly when fossil generation is required to satisfy demand as soon variable generation drops off in the evenings or during cloud cover and calm wind periods. Without pre-warming of the HRSG and ST, the CCPP would need several hours for a cold or warm start, burn a lot of gas and release high NOx emissions during start and wouldn’t be able to use the short times of high energy prices in an efficient or economical manner. The economic parameters of CiTES were determined by a full year “8760” simulation using a data set calculation for each of the hours of the year, and historical electricity and gas prices. For consistency, the simulations were focused on the pre-COVID year 2019. The financially most attractive markets were in the Energy Reliability Council of Texas (ERCOT) region, which allowed substantial value generation with arbitrage (charge with cheap energy during renewable surplus times and discharge when energy is needed and expensive). The improvement of flexibility with the CiTES system by pre-heating and warm-keeping of the CCPP allowed for additional power generation during short time periods when demand is high but renewable generation is down; when the hourly energy prices are highest in these markets. The simulations are based on 2019 data, when ERCOT had 27GW of installed photovoltaic (PV) and wind generation. They showed that the created revenue with the prototypically sized CiTES system of 1,000MMBTU (300MWh-th) falls a little bit short of what is expected from a commercially viable investment. The system has optimization opportunities for cost reduction and increased effectiveness which will be realized during a potential Phase II Pre-FEED study following this project. Furthermore, it is safe to assume that a lot of renewable generation capacity will be added all over the US in the coming years. As an example, ERCOT is predicting to more than double its renewable generation from 27GW in 2019 to a forecasted 63GW in 2023. This will increase the amount of renewable overproduction exponentially. This rapid increase of local overproduction and the need to curtail renewable generation is well documented by the California ISO (www.CAISO.com / managing oversupply). However, the simulations also revealed a weakness in the structure of the electricity markets in the US. More specifically, when electricity prices are very low and approaching negative levels, the owners of Variable Renewable Energy (VRE) will curtail a part of their facility to stabilize the price by reducing supply. This results in a situation in which storage facilities, which are integrated in existing fossil assets and don’t have the behind-the-meter benefit of a VRE, won’t be able to purchase low cost – otherwise curtailed – renewable energy off the grid. A special tariff, which motivates VRE owners to sell otherwise curtailed renewable energy to storage facilities (Hydrogen, thermal, pumped hydro, etc.) can solve this issue. The implementation of such a regulating tariff by Independent System Operators, thus avoiding renewable curtailment, is a pre condition for successful commercialization for renewable energy storage technologies. With this advancement of design and technology and improvements in the market environment, it can be expected that the Combined Cycle integrated Thermal Energy Storage proves itself as an important innovation to keep highly efficient, natural gas-based power generation economically successful and relevant for the power industry in the United States of America.

Wolf, Thorsten↗

Bayesian Seismoacoustic Source Location: Acoustic Approach

Seismic waves and infrasound are key technologies in the International Monitoring System (IMS) to monitor explosive events in the solid Earth and atmosphere. Energetic man-made or natural events (e.g., chemical/nuclear explosions, volcanic eruptions, and earthquakes) near the Earth’s surface produce both ground motion and atmospheric pressure disturbances which propagate as seismic waves and infrasound, respectively. Seismic waves have been generally used to detect and identify underground and near-surface events (Myers, et. al., 2007), and infrasound are sensitive to events near the surface or in the atmosphere (Modrak et. al., 2010). Due to their different sensitivities to events, they can complement to each other to improve the event detection and discrimination. The framework of joint seismoacoustic event location has recently reviewed by theoretical research (Koch and Arrowsmith, 2019). Although the early applications showed promising results to improve the accuracy of event location, their application were still limited to a small set of events selected to prove the concepts, and practical capability of the method for operational purpose is not fully evaluated with data. Our final goal is to apply the method of seismoacoustic event location to a larger set of events and evaluate its applicability for operational event location in practice. To that end, we focus on developing and verifying acoustic source location method in this study.

58 GEOSCIENCES↗

FY2022 Progress on Imbibition Testing in Containment Science

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

42 ENGINEERING↗

Capturing Carbonation: Understanding Kinetic Complexities through a New Era of Electron Microscopy

Cryogenic plasma focused ion beam (PFIB) electron microscopy analysis is applied to visualizing ex situ (surface industrial) and in situ (subsurface geologic) carbonation products, to advance understanding of carbonation kinetics. Ex situ carbonation is investigated using NIST fly ash standard #2689 exposed to aqueous sodium bicarbonate solutions for brief periods of time. In situ carbonation pathways are investigated using volcanic flood basalt samples from Schaef et al. (2010) exposed to aqueous CO 2 solutions by them. The fly ash reaction products at room temperature show small amounts of incipient carbonation, with calcite apparently forming via surface nucleation. Reaction products at 75° C show beginning stages of an iron carbonate phase, e.g., siderite or ankerite, common phases in subsurface carbon sequestration environments. This may suggest an alternative to calcite in carbonation low calcium-bearing fly ashes. Flood basalt carbonation reactions show distinct zonation with high calcium and calcium-magnesium bearing zones alternating with high iron-bearing zones. The calcium-magnesium zones are notable with occurrence of localized pore space. Oscillatory zoning in carbonate minerals is distinctly associated with far-from-equilibrium conditions where local chemical environments fluctuate via a coupling of reaction with transport. The high porosity zones may reflect a precursor phase (e.g., aragonite) with higher molar volume that then “ripens” to the high-Mg calcite phase-plus-porosity. These observations reveal that carbonation can proceed with evolving local chemical environments, formation and disappearance of metastable phases, and evolving reactive surface areas. Together this work shows that future application of cryo-PFIB in carbonation studies would provide advanced understanding of kinetic mechanisms for optimizing industrial-scale and commercial-scale applications.

36 MATERIALS SCIENCE↗

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗

Distributed Acoustic Sensing as a Monitoring Tool at LANL [Slides]

DAS records a variety of signals, including (1) earthquakes and volcanic eruptions, (2) explosions (DAG, mining, etc.), (3) acoustic waves, (4) ocean (waves, whales), (5) vehicle, ship and pedestrian traffic, (6) glaciology and landslides, and (7) oil and gas. Main questions include: (1) what is the extent of the signal content? (2) how does DAS compare to traditional sensors? (3) how do interrogators, fiber types, and fiber depth burials impact measurements? (4) special applications offshore to overcome the lack of instrumentation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The geology, geochemistry, and geographic distribution of southern Idaho obsidian

The Snake River Plain Volcanic Province of Southern Idaho is host to an abundance of pyroclastic obsidian deposits that constitute the principal lithic raw material utilized by Native Americans of the region throughout the Precontact Period. Archaeologists have long looked to obsidian as a source of information on past patterns of raw material selection, landscape use, and mobility through geochemical analysis of obsidian artifacts and lithic raw materials. Formation of the Snake River Plain through cataclysmic “super-eruptions” of the Yellowstone hotspot, however, presents one of the most complex and poorly understood obsidian source-scapes of North America. To clarify the geographic distribution of Southern Idaho obsidian source groups, the INL CRMO has compiled a comprehensive reference library of geologic obsidian from 138 source locales distributed across the southern half of the state. In this report, we define 28 geochemically distinct obsidian source groups that occur at these deposits through X-ray Fluorescence spectrometry and contextualize the geographic distribution of each source group relative to the history of silicic magmatism on the Snake River Plain. This research facilitates ongoing investigations at the INL CRMO into patterns of obsidian procurement, conveyance, and exchange in Southern Idaho throughout the Precontact Period.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗