SEARCH · Engineering Papers
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.
Item purpose: SnT virtual poster showing LYNM work Internal peer review status: This poster has been peer reviewed HQ review status: this poster has been approved by HQ
Earthquake Early Warning (EEW) systems depend on the dense seismic network to make fast and accurate earthquake detections to issue useful early alerts. The recent development and implementation of the ShakeAlert system is relying on the west coast seismic network to detect and estimate the earthquake parameters in real-time. While working on making improvements on the current system, one potential approach is to include more seismic recordings from various sources, such as the low-cost smartphone seismic network, such as the MyShake network, or the Google Android smartphones, since these smartphone seismic networks have much more portable stations (phones) to potentially provide useful data to the system. This report summarizes the initial exploration of using a deep learning approach to combine the data from both traditional seismic stations and the smartphone data in California. Due to the heterogeneity nature of the data, we aggregate data from the phones (using simulation data), as well as that from traditional seismic stations to grid cells. By generating real-time triggering ratio grid cell maps, the designed deep learning algorithm can process the data from multiple sources and detect the earthquake faster than only using that from a traditional seismic network.
SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.
This project demonstrates the feasibility of a passive tracking capability for ground-based and airborne moving sources using multiple information streams from acoustic and seismic signals. These moving source types are characterized by their variability with time and distance (exemplified but not limited to the familiar doppler shift of a passing vehicle) and the signatures used for detection should have similar characteristics. We have demonstrated that time-varying sensor signatures (TVS) combined with time-stable sensor signatures provide a more complete understanding of moving source behavior. Importantly, we have collected a new, relevant LANL dataset of moving sources (small uncrewed aerial vehicle flights and ground vehicle movement), developed new signatures, and demonstrated the capability to identify, track and forecast the position of these moving sources.
Multimodal, curated data sets and nuisance event catalogs remain rare in the explosion monitoring community relative to curated seismic data sets. The source of this relative absence is the difficultly in deploying multimodal receivers that sense the seismic, acoustic, and other modalities from multiphysics sources. We provide such a data set in this study that delivers seismic, infrasound, and electromagnetic (magnetometer) sensor records collected over a two–week period, within 255 km of a 10 ton buried chemical explosion called DAG–4 that was located at 37.1146°, –116.0693° on 22 June 2019 21:06:19.88 UTC. This catalog includes 485 seismic, seismoacoustic, and infrasound–only events that an expert analyst manually built by reviewing waveforms from 29 seismic and infrasound sensors. Our data release includes waveforms from these 29 seismic, infrasound, and seismoacoustic stations and two magnetometer stations and their station metadata. We deliver these waveforms in NNSA KB Core CSS.w format (i4) with a corresponding wfdisc table that provides the header information. Here, we expect that this data set will provide a valuable, benchmark resource to develop signal processing algorithms and explosion monitoring methods against manual, human observations.
Abstract We performed the first 3D island‐scale tomography of the electrical conductivity of Stromboli volcano (Aeolian Islands, Italy) using 2D acquisition lines (37.2 km) and a total of 18,880 measurements and 2,402 unique electrode locations. This 3D data set was inverted using a Gauss‐Newton algorithm, parallel‐processing on an unstructured tetrahedral mesh containing 678,420 finite‐element nodes and 3,580,145 elements to account for the topography of the volcanic island. The tomogram exhibits a conductive body (10 −2 –1.0 S m −1 ) consistent with the location of CO 2 and temperature anomalies observed at the ground surface. It corresponds to the hydrothermal system with high electrical conductivity associated with alteration. In order to confirm this interpretation, a 2.5D large‐scale induced polarization tomography was performed crossing the volcano. The joint interpretation of the conductivity and normalized chargeability is done with a petrophysical model previously tested and verified at both shield‐ and strato‐volcanoes. This model implies that alteration (through the effect of the cation exchange capacity associated with clay minerals and zeolites) plays a strong role in both controlling the electrical conductivity and normalized chargeability at Stromboli volcano. A temperature tomogram, derived from the geoelectrical measurements, is consistent with surface temperature anomalies and the Very Long Period (VLP) seismicity related to the mild‐explosive activity. This survey displays at 600 m a.s.l. a lateral shift in the highest temperature location, also corresponding to the source of VLP seismicity. Structural boundaries have a major role in the hottest hydrothermal fluids rising below the active crater terrace of Stromboli volcano.
This dataset provides a set of polygons representing the zone of depletion for historic submarine landslides (also referred to as mass transport deposits) within four regions of the US Exclusive Economic Zone in the northern Gulf of Mexico. Landslides were digitized by geologists and spatial scientists at the National Energy Technology Laboratory by visually interpreting landslide boundaries from a seismic-derived, high resolution bathymetric hillshade provided by the Bureau of Ocean and Energy Management (BOEM, 2017). A portion of the landslide features are derived from other spatial sources including the Seismic Water Anomalies dataset by BOEM (2016) as well as from McAdoo et al., 2000 and Twichell et al, 2005. The scale that landslides depletion areas can be interpreted at is limited by the spatial resolution of the gridded bathymetry, which is 12.192 meters (BOEM, 2017). For each landslide feature in the dataset, geometry metrics were calculated including geodesic area (km2) and geodesic perimeter (km) using the North America Albers Equal Area Conic projected coordinate system. The same geometry metrics were calculated for the four inventory regions.
Although accurately classifying signals from earthquakes and explosions at local distance (<250 km) remains an important task for seismic network operations, the growing volume of available seismic data presents a challenge for analysts using traditional source discrimination techniques. In recent years, deep-learning models have proven effective at discriminating between low-magnitude earthquakes and explosions measured at local distances, but it is not clear how well these models are capable of generalizing across different geological settings. To address the issue of generalization between regions, we train deep-learning models (convolutional neural networks [CNNs]) on time–frequency representations (scalograms) of three-component earthquake and explosion signals from eight different regions in the continental United States. We explore scenarios where models are trained on data from all regions, individual regions, or all but one region. We find that although CNN models trained on individual regions do not necessarily generalize well across different settings, models trained on multiple regions that include diverse path coverage generalize to new regions, with station-level accuracy of up to 90% or more for data sets from unseen regions. In general, CNN-based discrimination models significantly outperform models based on uncorrected P/S ratio (measured in the 10–18 Hz frequency band), even when CNN models are tested on data from entirely unseen regions.
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.
Conference presentation for International Meeting for Applied Geoscience & Energy (IMAGE), Houston, TX, August 28 – September 1, 2023. As part of the North Dakota CarbonSAFE project, seismic, controlled-source electromagnetic, and microgravity baseline monitoring surveys were collected to assess the safe, permanent, commercial-scale geologic storage of CO 2 generated by the Milton R. Young coal-fired power plant. Data from five 2D seismic lines were assessed as an alternative to 3D seismic for CO 2 monitoring. 2D and 3D surveys were compared using wave equation-based (WEB)-AVO inversion. The encouraging WEB-AVO results demonstrated that 2D surveys could be included in CO 2 monitoring programs as a cost-effective method.
Results are presented for an analysis of all moonquake data obtained by the Apollo seismic stations during the period from November 1969 to May 1974 and a preliminary analysis of critical data obtained in the interval from May 1974 to May 1975. More accurate locations are found for previously located moonquakes, and additional sources are located. Consideration is given to the sources of natural seismic signals, lunar seismic activity, moonquake periodicities, tidal periodicities in moonquake activity, hypocentral locations and occurrence characteristics of deep and shallow moonquakes, lunar tidal control over moonquakes, lunar tectonism, the locations of moonquake belts, and the dynamics of the lunar interior. It is concluded that: (1) moonquakes are distributed in several major belts of global extent that coincide with regions of the youngest and most intense volcanic and tectonic activity; (2) lunar tides control both the small quakes occurring at great depth and the larger quakes occurring near the surface; (3) the moon has a much thicker lithosphere than earth; (4) a single tectonic mechanism may account for all lunar seismic activity; and (5) lunar tidal stresses are an efficient triggering mechanism for moonquakes.
Explore the source record for details and available documents.
We consider the problem of designing a sensor network to most efficiently locate the source of a seismic event. The left panel of Figure 1 shows 982 possible sensor locations in a region of interest. The objective is to infer the unknown location of an earthquake, indicated by the gold dot for illustration purposes. Assume that three sensors have already been placed for data collection. Travel times of the P-wave from the source to the three sensors are observed and used to infer the unknown source location via a Bayesian statistical model described subsequently. The right panel of Figure 1 shows samples from the resulting posterior distribution for source location along with a posterior estimate of this location.
The author has identified the following significant results. At least three seismically active faults were identified which had not been previously mapped. One of these passes near the proposed site of a hydroelectric project on the Susitna River. Evidence of the state's past deformational history was obtained, indicating that right lateral offset has occurred sequentially from the northern part of the state to the southern. An apparent fault passes near Fairbanks, and is presumably the source of much seismic activity in the area.
A systematic study has been conducted of oceanic intraplate seismicity in the 1964-1983 period occurring in lithosphere younger than 35 million years. Detailed findings are presented on the stress in young lithosphere, depth of seismicity, intraplate events in fracture zones, seismicity distributions, and regional concentrations of seismicity. It is concluded that there is no evidence of a general transition from tensional to compressive stresses at any particular age of the lithosphere. Extensional events are located mostly in the Indian Ocean and have tensional axes oriented at large angles to the spreading direction. Compressional axes of thrust faulting events show a weak preferred orientation parallel to the spreading direction. Normal faulting events occur at greater depths and temperatures than thrust faulting events. Differential thermal contraction along fracture zones is not a significant source of intraplate seismicity.
The Leo Brady Seismic Network (LBSN, originally the Sandia Seismic Network) was established in 1960 by Sandia National Laboratories to monitor underground nuclear tests (UGTs) at the Nevada National Security Site (NNSS, formerly named the Nevada Test Site). The LBSN has been in various configurations throughout its existence, but it has generally been comprised of four to six stations at regional distances (~150–400 km) from the NNSS with approximately evenly spaced azimuthal coverage. Between 1962 and the end of nuclear testing in 1992, the LBSN—and a sister network operated by Lawrence Livermore National Laboratories—was the most comprehensive United States source of regional seismic data of UGTs. Approximately 75% of all UGTs performed by the United States occurred in the predigital era. At that time, LBSN data were transmitted as frequency-modulated (FM) audio over telephone lines to a central location and recorded as analog waveforms on high-fidelity magnetic audio tapes. These tapes have been in dry temperature-stable storage for decades and contain the sole record of this irreplaceable data; full waveforms of LBSN-recorded UGTs from this era were not routinely digitized or otherwise published. We have developed a process to recover and calibrate data from these tapes. First, we play back and digitize the tapes as audio. Next, we demodulate the FM “audio” into individual waveforms. We then estimate the various instrument constants through careful measurement of “weight-lift” tests performed prior to each UGT on each instrument. Finally, these coefficients allow us to scale and shape the derived instrument response of the seismographs and compute poles and zeros. Finally, the result of this process is a digital record of the recorded seismic ground motion in a modern data format, stored in a searchable database. To date, we have digitized tapes from 592 UGTs.
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.
An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.