Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “active source”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Insulation Activation and Contamination in Low-Power Experiments

Among the diverse advanced nuclear reactor concepts being developed are designs that expect to operate at much higher temperatures than conventional light-water reactors, and therefore require alternative materials for components capable of withstanding these extreme temperatures. Insulation is one such component where alternative high-temperature industrial composites are being considered. Low-power experiments can provide an opportunity for advancing our understanding of material behavior under irradiation; however, these experiments are not without occupational hazards. The focus of this paper is the potential for activation of insulating composites when exposed to neutron fluence during reactor experiments, and to discuss the subsequent contamination of the testing area. This paper will also discuss practical measures for preventing the creation and inhalation of activated particulate matter (i.e., dust), such as engineering and administrative controls and personal protective equipment. This paper evaluates the suitability of several insulation materials for use in reactor experiments based on their activation from neutron irradiation. The materials considered are Pyrogel XTE by Aspen Aerogel, Cerablanket and Kaowool by Morgan Advanced Materials, Thermo-12 Gold by Johns Manville, and Maxsil CF6-2000 by McAllister Mills. Although these materials are not yet widely used in the nuclear industry, they are common in other industries due to their thermal insulating properties. These materials produce dust while handling that can pose health and safety hazards, which could be further exacerbated by activation. The analysis presented in this paper utilizes a deterministic method developed to predict the activation source term of the provided insulation materials. The method calculates the atomic density of each element within the composites and simulates their exposure to a specified neutron flux density representative of low power experimental conditions. During the irradiation phase and throughout the extensive cooling period, neutron absorptions and nuclear decays are tracked for relevant isotopes. The residual activity within each composite is then assessed to identify the material with the lowest potential radiological hazard to workers. Dust will be created during installation and removal of insulating materials and may also be generated during the experiment. Contamination of surrounding surfaces due to dispersion is highly likely unless preventative measures are taken. The results of this analysis are used to estimate an internal dose from inhalation of activated dust, which is significantly lower than the allowable whole-body dose. However, the risk of internal exposure may increase if the dust becomes a transport vector for nuclear material. The goals of this paper are to enhance the understanding of the potential activation and contamination hazards from the use of these insulating materials in low-power nuclear experiments and to foster safer work environments through informed decision-making and strategic planning. This work ultimately contributes to a broader understanding of how new non-nuclear materials for advanced reactors impact the radiological safety of co-located workers, helping to ensure the responsible management of radioactive materials during installation and decommissioning activities.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2021

Assessing the Accuracy of a Borehole-Controlled Velocity Model of Yucca Flat, Nevada, Using Large-N Seismic Data

The Source Physics Experiment (SPE) aims to improve nuclear test monitoring capabilities through a series of chemical explosions in Yucca Flat within the Nevada National Security Site (NNSS) (Snelson et al., 2013). With geologic data from over 950 boreholes, Yucca Flat has excellent borehole control on stratigraphy. These data were used to create a Geologic Framework Model (GFM) of the basin to aid in the study of seismic wave propagation through the Earth. However, the seismic velocities of the geologic units in the GFM are less well-defined. With the acquisition of five new active source Large-N datasets, we can now quantitatively assess the accuracy of the seismic model previously controlled by borehole data alone.

58 GEOSCIENCES↗
DOE OSTI2022

Respiratory energy demands and scope for demand expansion and destruction

Photosynthesis is the primary energy input to plants, but most plant metabolic processes are powered much or all of the time by ATP and NAD(P)H that come from respiratory oxidation of photosynthetically produced carbohydrates, that is, “dark respiration” (Amthor, 1994). The processes of growth, nutrient uptake and assimilation, active transport, and maintenance are thus all clients of respiration and compete for their share of a respiratory energy budget that is capped by photosynthate income. The economic principle of opportunity cost applies to these clients’ competing demands: spending respiratory energy on process A means missing the benefit of investing that energy in process B (Mahmoudabadi et al., 2019). This principle is key to assessing prospects for crop improvement by metabolic engineering. As de Lorenzo (2015) put it: “metabolism…frames and ultimately resolves whether a given genetic program (existing…or engineered) can be deployed or not.” Foreign or reconfigured native processes bolted on to a crop-plant’s metabolic chassis must compete for respiratory energy with native ones without crashing the energy economy. A proviso on plant carbon budgets is that photosynthetic carbon fixation (“source activity”) can in certain cases increase to meet increased carbon demand (“sink activity”), that is, budget envelopes are not always fixed (Smith et al., 2018). However, as there is some consensus that productivity is most often limited or co-limited by carbon supply (Ainsworth and Long, 2005; Körner, 2015; Sonnewald and Fernie, 2018), we make this our basal assumption in the analyses below.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Template and target-site recognition by human LINE-1 in retrotransposition

The long interspersed element-1 (LINE-1, hereafter L1) retrotransposon has generated nearly one-third of the human genome and serves as an active source of genetic diversity and human disease. L1 spreads through a mechanism termed target-primed reverse transcription, in which the encoded enzyme (ORF2p) nicks the target DNA to prime reverse transcription of its own or non-self RNAs. Here we purified full-length L1 ORF2p and biochemically reconstituted robust target-primed reverse transcription with template RNA and target-site DNA. We report cryo-electron microscopy structures of the complete human L1 ORF2p bound to structured template RNAs and initiating cDNA synthesis. The template polyadenosine tract is recognized in a sequence-specific manner by five distinct domains. Among them, an RNA-binding domain bends the template backbone to allow engagement of an RNA hairpin stem with the L1 ORF2p C-terminal segment. Moreover, structure and biochemical reconstitutions demonstrate an unexpected target-site requirement: L1 ORF2p relies on upstream single-stranded DNA to position the adjacent duplex in the endonuclease active site for nicking of the longer DNA strand, with a single nick generating a staggered DNA break. Our research provides insights into the mechanism of ongoing transposition in the human genome and informs the engineering of retrotransposon proteins for gene therapy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Data Analysis of the 2020 Central Idaho Mainshock-Aftershock Sequence

In an effort to inform the Senior Seismic Hazard Analysis Committee for the Idaho National Laboratory, we provide an improved aftershock catalog related to the March 31, 2020, Mw6.5 Stanley, Idaho earthquake from picks related to a temporary network of two real-time and 15 non-telemetered seismometers within the epicentral area. From the permanent and temporary (XP) real-time network, the USGS cataloged 1,946 aftershocks between April 1, 2020 and October 31, 2020. To improve aftershock location and magnitudes, we manually picked arrival times of P and S waves from off-line stations in the XP temporary network, generated a new crustal velocity model, and independently relocated each event using the HypoDD double-difference earthquake algorithm. We created our new velocity model from existing broadband and active source seismic campaign data that were acquired near the epicentral region prior to the 2020 earthquake. We compare arrival time differences, epicentral locations and depths between aftershocks recorded with the two catalogs. We find the addition of local stations provides tighter aftershock clustering that suggests an improved aftershock locations. To detect lower magnitude events, we employed deep learning. Our method solves common problems associated with detecting many events that have a low signal-to-noise ratio. From the machine learning database, we detected more than 74,000 aftershocks. Based on the number of identified earthquakes and Gutenberg-Richter relationships derived from the USGS catalog, we estimate that we have reduced the completion magnitude for the Stanley earthquake sequence to below M1 using this machine learning approach. We located each aftershock with our new velocity model. Our new velocity model and picks suggests aftershocks occurred mostly at shallower depths than assessed in the USGS catalog. These aftershocks align along two linear trends that suggest the activation of two unnamed primary faults.

58 GEOSCIENCES↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗
DOE OSTI2023

High Spatial Melt Rate Variability Near the Totten Glacier Grounding Zone Explained by New Bathymetry Inversion

Totten Glacier is a fast-moving East Antarctic outlet with the potential for significant future sea-level contributions. We deployed four autonomous phase-sensitive radars on its ice shelf to monitor ice-ocean interactions near its grounding zone and made active source seismic observations to constrain gravity-derived bathymetry models. We observe an asymmetry in basal melting with mean melt rates along the grounding zone differing by up to 20 m/a. Our new bathymetry model reveals that this melt rate asymmetry coincides with an asymmetry in water column thickness and that the low-melting ice-shelf portion is shielded from the main cavity circulation. A 2-year record yields year-to-year melt rate variability of 7–9 m/a with no seasonal cycle. Our results highlight the key role of bathymetry near grounding lines for accurate modeling of ice-shelf melt, and the importance of sustained multi-year monitoring, especially at ice-shelf cavities where the dominant melt rate drivers vary primarily inter-annually.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Moisture-induced autonomous surface potential oscillations for energy harvesting

A variety of autonomous oscillations in nature such as heartbeats and some biochemical reactions have been widely studied and utilized for applications in the fields of bioscience and engineering. Here, we report a unique phenomenon of moisture-induced electrical potential oscillations on polymers, poly([2-(methacryloyloxy)ethyl] dimethyl-(3-sulfopropyl) ammonium hydroxide-co-acrylic acid), during the diffusion of water molecules. Chemical reactions are modeled by kinetic simulations while system dynamic equations and the stability matrix are analyzed to show the chaotic nature of the system which oscillates with hidden attractors to induce the autonomous surface potential oscillation. Using moisture in the ambient environment as the activation source, this self-excited chemoelectrical reaction could have broad influences and usages in surface-reaction based devices and systems. As a proof-of-concept demonstration, an energy harvester is constructed and achieved the continuous energy production for more than 15,000 seconds with an energy density of 16.8 mJ/cm 2 . A 2-Volts output voltage has been produced to power a liquid crystal display toward practical applications with five energy harvesters connected in series.

42 ENGINEERING↗
DOE OSTI2022

Direct neutrino-mass measurement with sub-electronvolt sensitivity

Since the discovery of neutrino oscillations, we know that neutrinos have non-zero mass. However, the absolute neutrino-mass scale remains unknown. Here we report the upper limits on effective electron anti-neutrino mass, m ν , from the second physics run of the Karlsruhe Tritium Neutrino experiment. In this experiment, m ν is probed via a high-precision measurement of the tritium β-decay spectrum close to its endpoint. This method is independent of any cosmological model and does not rely on assumptions whether the neutrino is a Dirac or Majorana particle. By increasing the source activity and reducing the background with respect to the first physics campaign, we reached a sensitivity on m ν –2 at a 90% confidence level (CL). The best fit to the spectral data yields m$^2_ν$ = (0.26 ± 0.34) eV 2 c –4 , resulting in an upper limit of m ν < 0.9 eV c –2 at 90% CL. By combining this result with the first neutrino-mass campaign, we find an upper limit of m ν < 0.8 eV c –2 at 90% CL.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Physics informed neural network can retrieve rate and state friction parameters from acoustic monitoring of laboratory stick-slip experiments

Various machine learning (ML) and deep learning (DL) techniques have been recently applied to the forecasting of laboratory earthquakes from friction experiments. The magnitude and timing of shear failures in stick-slip cycles are predicted using features extracted from the recorded ultrasonic or acoustic emission (AE) signals. In addition, the Rate and State Friction (RSF) constitutive laws are extensively used to model the frictional behavior of faults. In this work, we use data from shear experiments coupled with passive acoustic (variance, kurtosis, and AE rate) interleaved with active source ultrasonic monitoring (transmitted wave amplitude) to develop physics-informed neural network (PINN) models incorporating the RSF law and AE rate generation equation with wave amplitude serving as a proxy for friction state variable. This PINN framework allows learning RSF parameters from stick-slip experiments rather than measuring them through a series of velocity step experiments. We observe that when the stick-slip cycles are irregular, the PINN models outperform the data-driven DL models. Transfer learning (TL) PINN models are also developed by pre-training on data collected at one normal stress level followed by forecasting shear failures and retrieving RSF parameters at other stress levels (i.e., with different recurrence intervals) after retraining on a limited amount of new data. Our findings suggest that TL models perform better compared to standalone models. Both standalone and TL PINN-estimated RSF parameters and their ground truth values show excellent agreements thus demonstrating that RSF parameters can be retrieved from laboratory stick-slip experiments using the corresponding acoustic data and that the transmitted wave amplitude provides a good representation of the evolving frictional state during stick-slips.

58 GEOSCIENCES↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗
DOE OSTI2021

Phase-weighted slant stacking for surface wave dispersion measurement

SUMMARY Surface wave retrieval from ambient noise records using seismic interferometry techniques has been widely used for multiscale shear wave velocity (Vs) imaging. One key step during Vs imaging is the generation of dispersion spectra and the extraction of a reliable dispersion curve from the retrieved surface waves. However, the sparse array geometry usually affects the ability for high-frequency (>1 Hz) seismic signals’ acquisition. Dispersion measurements are degraded by array response due to sparse sampling and often present smeared dispersion spectra with sidelobe artefacts. Previous studies usually focus on interferograms’ domain (e.g. cross-correlation function) and attempt to enhance coherent signals before dispersion measurement. We propose an alternative technique to explicitly deblur dispersion spectra through use of a phase-weighted slant-stacking algorithm. Numerical examples demonstrate the strength of the proposed technique to attenuate array responses as well as incoherent noise. Three different field examples prove the flexibility and superiority of the proposed technique: the first data set consists of ambient noise records acquired using a nodal seismometer array; the second data set utilizes distributed acoustic sensing (DAS) and a marine fibre-optic cable to acquire a similar ambient noise data set; the last data set is a vibrator-based active-source surface wave data. The enhanced dispersion measurements provide cleaner and higher-resolution spectra without distortions which will assist both human interpreters as well as ML algorithms in efficiently picking curves for subsequent Vs inversion.

58 GEOSCIENCES↗
DOE OSTI2023

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗
DOE OSTI2021

The feasibility of MT tipper data to monitor CO2 storage sites

Monitoring carbon storage sites using geophysical techniques is a critical component to the success and safety of storage programs. Currently, the primary methods of monitoring such sites are seismic and, to a much lesser extent, electromagnetics using active sources. The cost of such methods, especially seismic, can be prohibitively expensive. Natural source lectromagnetics, or magnetotellurics (MT), represents a low-cost, and underutilized method with the potential to aid CO2 monitoring efforts. Specifically, the tipper of MT data gives insight to the dimensionality of the subsurface and is able to detect the expansion front of a CO2 plume in a saline reservoir. We analyze the feasibility of using the tipper to monitor two shallow CO2 plumes and conclude that the tipper may be a suitable method for long-term monitoring.<br>

Kohnke, Colton↗
DOE OSTI2024

Utah FORGE 3-2535: Joint Imaging of Fracture Growth and Estimation of Fracture Properties During EGS Development - 2024 Annual Workshop Presentation

This is a presentation on the Joint Electromagnetic/Seismic/InSAR Imaging of Spatial-Temporal Fracture Growth and Estimation of Physical Fracture Properties During EGS Resource Development by Lawrence Berkeley National Laboratory, presented by David Alumbaugh. This is a video presentation on developing a set of technologies and workflow to image induced fracture generation and growth for an Enhanced Geothermal System (EGS). Using a combination of passive seismic and active source borehole EM and InSAR technology it is anticipated imaging and fracture generation and growth monitoring will be actualized. The workflow will provide a template for imaging induced fracture systems for future EGS systems. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 13, 2024.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗