Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic location”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗
NASA NTRS2023

Interchangeable Use of GNSS and Seismic Data for Rapid Earthquake Characterization: 2021 Chignik Earthquake, Alaska

Earthquake magnitude estimation using peak ground velocities (PGV) derived from 13 Global Navigation Satellite Systems (GNSS) data has shown promise for rapid 14 characterization of damaging earthquakes. Here we examine the feasibility of using 15 GNSS-derived velocity waveforms as interchangeable data for ground motion estimation 16 and other products that typically rely on strong-motion seismic records. Our study 17 compares PGVs derived from high-rate GNSS to those computed from high-rate seismic 18 records (strong-motion and velocity), at co- and closely-located stations. The recent 2021 Manuscript Click here to access/download;Manuscript;Manuscript_Final.docx 2 19 Mw 8.2 Chignik earthquake in Alaska that was recorded on co-located GNSS and strong20 motion sensors provides the perfect opportunity to compare the two data streams and 21 their application in rapid response. The Chignik velocity records appear almost identical 22 at co-located GNSS and strong-motion stations when observed at frequencies < 0.25 Hz. 23 GNSS and strong-motion derived velocity data are further employed to generate rapid 24 estimates of PGV-derived moment magnitudes for the earthquake. The moment 25 magnitude estimates from GNSS and joint GNSS/seismic data are within ~ ±0.4 26 magnitude units (Fang et al., 2020) of the final magnitude (Mw 8.2). ShakeMaps 27 generated for the 2021 Chignik earthquake using GNSS and seismic PGVs show notable 28 agreement between them, and show negligible shifts in PGV contours when co-/closely 29 located GNSS and seismic stations are substituted for one another. Therefore, we posit 30 that GNSS is a powerful alternative or addition to seismic data and vice versa.

Earthquake rapid response↗
NASA NTRS1985

Seismic analysis of the large 70-meter antenna. Part 2: General dynamic response and a seismic safety check

An extensive dynamic analysis for the new JPL 70-meter antenna structure is presented. Analytical procedures are based on the normal mode decomposition which include dumping and special forcing functions. The dynamic response can be obtained for any arbitrarily selected point on the structure. A new computer program for computing the time-dependent, resultant structural displacement, summing the effects of all participating modes, was developed also. Program compatibility with natural frequency analysis output was verified. The program was applied to the JPL 70-meter antenna structure and the dynamic response for several specially selected points was computed. Seismic analysis of structures, a special application of the general dynamic analysis, is based also on the normal modal decomposition. Strength specification of the antenna, with respect to the earthquake excitation, is done by using the common response spectra. The results indicated basically a safe design under an assumed 5% or more damping coefficient. However, for the antenna located at Goldstone, with more active seismic environment, this study strongly recommends and experimental program that determines the true damping coefficient for a more reliable safety check.

Kiedron, K.↗
DOE OSTI2021

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗
DOE OSTI2022

Crustal and uppermost mantle structure of Cape Verde from ambient noise tomography

SUMMARY We present a seismic ambient noise tomography of the Cape Verde archipelago, located in the Atlantic Ocean, approximately 600 km west of Senegal. We used 38 seismic broad-band stations that continuously recorded for 10 months, in order to construct the first 3-D model of Sv-wave velocities for the crust and uppermost mantle beneath the Cape Verde region. We started by computing phase cross-correlations for vertical component recordings using all possible inter-island station pairs. Next, a time–frequency phase-weighted stack was applied to obtain robust Rayleigh-wave group-velocity dispersion curves in the period band between 10 and 24 s. Group-velocity maps at different periods are obtained by inverting the dispersion curves. We then inverted the group-velocity maps to obtain the 3-D shear wave velocity structure of the crust and uppermost mantle beneath Cape Verde. The final 3-D model extends from 8 km down to 23 km and has a lateral resolution of about 50 km. The crust in the southwestern sector, encompassing Fogo, presents lower S-wave velocities that may be caused by the presence of melt pockets and/or hydrothermal fluids circulation. The uppermost mantle beneath the northwestern sector is characterized by higher S-wave velocities in agreement with previous results obtained from Ps and Sp receiver functions. Those high-velocity anomalies can reflect non-altered crust or remnants of magma chambers or solidified basaltic intrusions, which fed the volcanism in these islands. Our maps revealed the presence of crustal underplating across the entire archipelago, yet stronger beneath the groups Santo Antão—São Vicente—São Nicolau and Fogo—Santiago—Maio.

Carvalho, J. (ORCID:0000000310778994)↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m- 1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave’s traveltime, but also the path of a wave, which can be observed in deviations of the wave’s incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope’s regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2021

Seismically Detecting Nuclear Reactor Operations Using a Power Spectral Density (PSD) Misfit Detector

To explore the ability to indirectly detect and attribute various operations conducted at a nuclear reactor using waveform data, we investigated the seismic signals recorded near the High Flux Isotope Reactor (HFIR) located at Oak Ridge National Laboratory in Oak Ridge, Tennessee. Specifically, we processed seismic data collected from a single seismoacoustic station, WACO, near the HFIR facility, and employed a power spectral density misfit detector to identify signals of interest and associate the detections with operational events. Initial results suggest that this method provides a promising means of regularly detecting at least 19 unique operations. Furthermore, with additional station deployment and more comprehensive data logs, we anticipate that future analysis will offer an additional means to seismically monitor nuclear reactors (such as HFIR) health and performance more accurately.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Methodological Overview of Seismic Analysis for Nuclear Event Detection

Underground explosions generate potentially detectable signatures, including energy waves that travel through the Earth’s subsurface (i.e., seismic waves), low-frequency sound waves (i.e., infrasound and hydroacoustic waves), and radioactive gases and/or particles that might leak from the test cavity (if the event was nuclear). There can also be intelligence indicators of a test, such as observations of modified patterns of life and activity at a suspected test site. If all of these detectable signatures and intelligence indicators are present and self-consistent, then analysts have high confidence in classifying a signature generating event as an explosion. However, because only partial information about an event is likely to be available, determining whether an event was natural (e.g., an earthquake or landslide) or manmade (e.g., a chemical or nuclear explosion) is much more challenging. This primer describes how one category of event signatures—seismic signatures—can augment event analyses. While universities and government organizations have generated detailed technical descriptions of seismic analytic techniques, we seek to translate seismic event analysis for a broad, non-technical audience. When the geologic conditions near an event are well-characterized, seismic data can be used to calculate critical information, such as event location and depth, with relatively high accuracy. Moreover, specific features within seismic datasets can help determine whether an event was an explosion. However, a key challenge in seismic analysis is that geologic site conditions are often poorly characterized, complicating the ability to discern the true nature of the event. To overcome this challenge, geologists answer a series of questions (discussed in section 1) to guide seismic event analysis and determine the most probable nature of an event. As more information is gathered during each analytic step, confidence grows regarding the nature of the event. Section 2 addresses uncertainties in seismic analysis and the vital nature of high-fidelity geologic data for accurate seismic event analysis.

58 GEOSCIENCES↗
DOE OSTI2023

Pore Pressure Diffusion and Onset of Induced Seismicity

The Raton Basin of Southern Colorado and Northern New Mexico experienced a significant increase in seismicity over the last several decades, including two M5+ earthquakes. The increase in seismicity started soon after an increase in wastewater disposal associated with coal-bed methane production. Twenty-nine wastewater injection wells have operated in the Basin, with several high-rate injectors located within a kilometer of the most seismically active regions. We developed a pore pressure model to investigate the relationship between pore pressure diffusion and onset of induced seismicity, with an emphasis on the Trinidad fault zones. Our model indicates that sufficient pore pressure increase had accumulated on the northern extent of the Trinidad fault zone to induce seismicity at the time of early-stage seismic activity. The hydrogeologic architecture of the Basin played a key role in controlling pore pressure perturbation. The basal sedimentary interval with moderately low diffusivity restricted vertical pore pressure diffusion, creating a substantial lag time between wells becoming active and the expression of perturbation in the fault zones. The spatiotemporal relationship between early-stage seismicity and wastewater injection in conjunction with modeled pore pressures suggests that the early-stage seismicity in the Raton Basin was likely triggered by pore pressure diffusion.

58 GEOSCIENCES↗
NASA NTRS1977

Lunar seismicity and tectonics

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.

Lammlein, D. R.↗
DOE OSTI2022

On component isolation of conceptual advanced reactors

Implementation of component isolation in nuclear industry is challenging due to gaps in research and the lack of specific guidelines. In this study, parameters affecting component-level isolation of advanced reactor vessels are identified based on a representative numerical model with explicit consideration of nonlinear soil-structure interaction (SSI). The objective of this study is to evaluate the effectiveness of, and to identify potential limitations of using conventional friction pendulum and lead-rubber bearings to seismically isolate vessels. It is found that slender vessels or components are particularly vulnerable to rotational accelerations at the isolation interface, which are caused by rotation of the sub-structure and by excitation of higher modes in the horizontal direction of the seismically isolated system. Component isolation is found to be more effective for relatively stiffer vessels and at sites with stiff soil properties. Considering that conventional isolators are deficient in resisting axial tension, it is observed that the optimum location for supporting a component to achieve seismic isolation, is at a cross-sectional plane passing through the center of mass of the vessel. These findings are corroborated by numerous simulations of the response of seismically isolated reactor vessels at different nuclear sites subject to a variety of ground motions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Emulation of seismic-phase traveltimes with machine learning

SUMMARY We present a machine learning (ML) method for emulating seismic-phase traveltimes that are computed using a global-scale 3-D earth model and physics-based ray tracing. Accurate traveltime predictions based on 3-D earth models are known to reduce the bias of event location estimates, increase our ability to assign phase labels to seismic detections and associate detections to events. However, practical use of 3-D models is challenged by slow computational speed and the unwieldiness of pre-computed lookup tables that are often large and have prescribed computational grids. In this work, we train a ML emulator using pre-computed traveltimes, resulting in a compact and computationally fast way to approximate traveltimes that are based on a 3-D earth model. Our model is trained using approximately 850 million P-wave traveltimes that are based on the global LLNL-G3D-JPS model, which was developed for more accurate event location. The training-set consists of traveltimes between 10 393 global seismic stations and randomly sampled event locations that provide a prescribed, distance-dependent geographic sample density for each station. Prediction accuracy is dependent on event-station distance and whether the station was included in the training set. For stations included in the training set the mean absolute deviation (MAD) of the difference between traveltimes computed using ray tracing through the 3-D model and the ML emulator for local, regional, and teleseismic distances are 0.090, 0.125 and 0.121 s, respectively. For tested station locations not included in the training set, MAD values for the three distance ranges increase to 0.173, 0.219 and 0.210 s, respectively. Empirical traveltime residuals for a global reference data are indistinguishable when ML emulation or the 3-D model is used to compute traveltimes. This result holds regardless of whether the recording station is used in ML training or not.

58 GEOSCIENCES↗
DOE OSTI2023

Fourier-DeepONet: Fourier-enhanced deep operator networks for full waveform inversion with improved accuracy, generalizability, and robustness

In this article, full waveform inversion (FWI) infers the subsurface structure information from seismic waveform data by solving a non-convex optimization problem. Data-driven FWI has been increasingly studied with various neural network architectures to improve accuracy and computational efficiency. Nevertheless, the applicability of pre-trained neural networks is severely restricted by potential discrepancies between the source function used in the field survey and the one utilized during training. Here, we develop a Fourier-enhanced deep operator network (Fourier-DeepONet) for FWI with the generalization of seismic sources, including the frequencies and locations of sources. Specifically, we employ the Fourier neural operator as the decoder of DeepONet, and we utilize source parameters as one input of Fourier-DeepONet, facilitating the resolution of FWI with variable sources. To test Fourier-DeepONet, we develop three new and realistic FWI benchmark datasets (FWI-F, FWI-L, and FWI-FL) with varying source frequencies, locations, or both. Our experiments demonstrate that compared with existing data-driven FWI methods, Fourier-DeepONet obtains more accurate predictions of subsurface structures in a wide range of source parameters. Moreover, the proposed Fourier-DeepONet exhibits superior robustness when handling data with Gaussian noise or missing traces and sources with Gaussian noise, paving the way for more reliable and accurate subsurface imaging across diverse real conditions.

42 ENGINEERING↗
DOE OSTI2025

Microseismicity Modulation Due To Changes in Geothermal Production at San Emidio, Nevada, USA

Brief cessations of geothermal production can induce seismicity, a phenomenon that has drawn increasing attention in recent years. Such observations are rare, and the underlying mechanism requires careful analysis. In April 2022, a dense seismic and hydrologic monitoring system was deployed at the San Emidio geothermal field, Nevada, to accompany a planned power plant shutdown. Using the dense seismic array data, we detected and located ∼1,800 microseismic events (MSEs) and developed a high-resolution tomographic P-wave velocity model. We observed substantially increased microseismicity during shutdown. Most MSEs occurred on pre-existing normal faults, which are contained within extremely low-velocity zones that are likely damaged, fluid-filled, and hydraulically connected to nearby production wells. Hydrologic data show rapid fluid pressure increases of <60 kPa following the shutdown. We suggest that the cessation of production rapidly increased fluid pressures along pre-existing fault zones, activating critically stressed fault patches and fractures and producing microseismicity.

Guo, Hao [University of Wisconsin-Madison, WI (Uni↗
DOE OSTI2026

Prairie State Generating Company 2D Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, two-dimensional (2D) seismic data was acquired at both sites. This report summarizes the findings from two phases of 2D acquisition and seismic interpretation at the Prairie State Generating Company site near Marissa, Illinois. The seismic data indicates the presence of three faults that completely transect the storage and confining units. The initial 2021 2D seismic acquisition revealed the presence of a feature two miles east of the Lively Grove #1 characterization well. The feature is a polyphase fault with a component of strike-slip motion, forming a small positive flower structure. The second phase of 2D seismic acquisition in 2022 constrained the maximum extent of the fault, indicating it has a relatively limited length. Other 2022 seismic lines revealed the presence of two other faults that transect the storage and confining units in the southwest and northeast portions of the 2D seismic acquisition area. The 2021 and 2022 2D seismic surveys have identified three specific locations within the project area that may have an elevated risk of out of zone CO 2 migration due to faulting. Further technical work will be needed to quantify this risk and determine how and if this risk will impact the placement of CO 2 injection wells across the project area.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1981

Lunar seismology - The internal structure of the moon

It is pointed out that seismology has provided the most detailed information concerning the structure and state of the earth's interior. Beginning in 1969, seismometers were landed on the moon by the Apollo missions, providing the first opportunity to attempt similar studies on another planetary body. In September 1977 the operation of these instruments was terminated. A description is presented of the internal structure of the moon, as determined from the obtained lunar seismic data. The analysis of the lunar data is approached in a systematic fashion, using appropriate techniques to minimize the number of necessary assumptions, extract the maximum amount of structural information, and determine its reliability. The completed lunar seismic network consists of four stations located at the landing sites of Apollo missions 12, 14, 15, and 16. Attention is given to crustal structure, the structure of the lunar mantle, the attenuating region, and the core.

Goins, N. R.↗
NASA NTRS1993

Characterization of the Marquez Dome buried impact crater using gravity and magnetic data

The buried impact crater, Marquez Dome, located in Leon County in east central Texas, is an approximately 15 km diameter structure whose central uplift is now partially exposed due to headward erosion of the post-impact cover. The central uplift is approximately 3 km in diameter and the rocks within it have been uplifted more than 1200 m above their regional level. The crater rim remains buried and previous attempts to determine its location have had to rely on seismic reflection data and geologic well logs. These attempts have been somewhat successful in mapping the extent of the disturbed zone around Marquez Dome, but more limited in their ability to image the shallow buried rim. In an attempt to define accurately the whole Marquez Dome structure and assist in the selection of drilling sites, a geophysical investigation involving gravity and magnetic data over the central uplift and the surrounding area has been undertaken.

Wong, A. M.↗