Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic coherency”

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 37 records · Page 2

Evaluating multistation phase picking algorithm phase neural operator (PhaseNO) on local seismic networks

Reliable automatic phase picking is important for many seismic applications. With the development of machine learning approaches, many algorithms are proposed, evaluated and applied to different areas. Many of these algorithms are single station based, while recent proposed methods start to combine surrounding stations into consideration in the problem of phase picking. Among these algorithms, the phase neural operator (PhaseNO) shows promising results on regional data sets comparing to existing algorithms. But there are many use cases for the local seismic networks in our community, therefore in this paper we evaluate the performance of PhaseNO on four different local data sets and compare the results to PhaseNet and EQTransformer. We used both individual phase picking metrics as well as association metrics to illustrate the performance of PhaseNO. By manually reviewing the newly detected events, we find that the PhaseNO model outperforms the single station-based approaches in the local-scale use cases due to its consideration of coherent signals from multiple stations. We also explored PhaseNO’s behaviours when only using one station, as well as gradually increasing the number of stations in the seismic network to better understand its behaviour. Overall, using the off-the-shelf machine learning based phase pickers, PhaseNO demonstrated its good performance on local-scale seismic networks.

58 GEOSCIENCES↗

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗

Quantitative assessment of Distributed Acoustic Sensing at the Source Physics Experiment, Phase II

In this report, we assess the data recorded by a Distributed Acoustic Sensing (DAS) cable deployed during the Source Physics Experiment, Phase II (DAG) in comparison with the data recorded by nearby 4.5-Hz geophones. DAS is a novel recording method with unprecedented spatial resolution, but there are significant concerns around the data fidelity as the technology is ramped up to more common usage. Here we run a series of tests to quantify the similarity between DAS data and more conventional data and investigate cases where the higher spatial resolution of the DAS can provide new insights into the wavefield. These tests include 1D modeling with seismic refraction and bootstrap uncertainties, assessing the amplitude spectra with distance from the source, measuring the frequency dependent inter-station coherency, estimating time-dependent phase velocity with beamforming and semblance, and measuring the cross-correlation between the geophone and the particle velocity inferred from the DAS. In most cases, we find high similarity between the two datasets, but the higher spatial resolution of the DAS provides increased details and methods of estimating uncertainty.

47 OTHER INSTRUMENTATION↗

Quantitative assessment of Distributed Acoustic Sensing at the Source Physics Experiment (Phase II)

In this report, we assess the data recorded by a Distributed Acoustic Sensing (DAS) cable deployed during the Source Physics Experiment, Phase II (DAG) in comparison with the data recorded by nearby 4.5-Hz geophones. DAS is a novel recording method with unprecedented spatial resolution, but there are significant concerns around the data fidelity as the technology is ramped up to more common usage. Here we run a series of tests to quantify the similarity between DAS data and more conventional data and investigate cases where the higher spatial resolution of the DAS can provide new insights into the wavefield. These tests include 1D modeling with seismic refraction and bootstrap uncertainties, assessing the amplitude spectra with distance from the source, measuring the frequency dependent inter-station coherency, estimating time-dependent phase velocity with beamforming and semblance, and measuring the cross-correlation between the geophone and the particle velocity inferred from the DAS. In most cases, we find high similarity between the two datasets, but the higher spatial resolution of the DAS provides increased details and methods of estimating uncertainty.

58 GEOSCIENCES↗

The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves

Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.

58 GEOSCIENCES↗

Sub-Surface Windscreen for the Measurement of Outdoor Infrasound

A windscreen has been developed that features two advantages favorable for the measurement of outdoor infrasound. First, the sub-surface location, with the top of the windscreen flush with the ground surface, minimizes the mean velocity of the impinging wind. Secondly, the windscreen material (closed cell polyurethane foam) has a sufficiently low acoustic impedance (222 times that of air) and wall thickness (0.0127 m) to provide a transmission coefficient of nearly unity over the infrasonic frequency range (0-20 Hz). The windscreen, a tightly-sealed box having internal dimensions of 0.3048 x 0.3048 x 0.3556 m, contains a microphone, preamplifier, and a cable feed thru to an external power supply. Provisions are made for rain drainage and seismic isolation. A three-element array, configured as an equilateral triangle with 30.48 m spacing and operating continuously in the field, periodically receives highly coherent signals attributed to emissions from atmospheric turbulence. The time delays between infrasonic signals received at the microphones permit determination of the bearing and elevation of the sources, which correlate well with locations of pilot reports (PIREPS) within a 320 km radius about the array. The test results are interpreted to yield spectral information on infrasonic emissions from clear air turbulence.

Shams, Qamar A.↗

A1 moonquakes - Problems of determining their epicenters and mechanisms

Review of some of the problems associated with the determination of epicenters and mechanisms of moonquakes with the aid of the four seismic stations now operating on the moon. No decoupling of compressional, transverse, and surface waves takes place there, no clear first arrivals can be found, and coherence between vertical and horizontal components on a seismograph stationed on the moon is poor. New techniques for locating epicenters are presented.

Meissner, R.↗

Field measurement of penetrator seismic coupling in sediments and volcanic rocks

Field experiments were conducted to determine experimentally how well a seismometer installed using a penetrator would be coupled to the ground. A dry lake bed and a lava bed were chosen as test sites to represent geological environments of two widely different material properties. At each site, two half-scale penetrators were fired into the ground, a three-component geophone assembly was mounted to the aft end of each penetrator, and dummy penetrators were fired at various distances to generate seismic signals. The recorded signals were digitized, and cross-spectral analyses were performed to compare the observed signals in terms of power spectral density ratio, coherence and phase difference. The analyses indicate that seismometers deployed by penetrators will be as well coupled to the ground as are seismometers installed by conventional methods for the frequency range of interest in earthquake seismology, although some minor differences were observed at frequencies near the upper limit of the frequency band.

Nakamura, Y.↗

Deformation of Alaskan Volcanoes, Measured by Satellite Radar Inferometry

The purpose of this project was to determine the suitability of measuring active deformation of volcanoes in Alaska using Interferometric Synthetic Aperture Radar (INSAR) techniques. Work sponsored by this grant supported one graduate student (for almost 2 years) and one postdoc (for several months), and has resulted in two published peer-reviewed papers and a front-page article in EOS. An additional paper is in review and a fourth is in preparation. An additional paper in preparation was based in part on research supported by this grant and in part by a successor grant from NASA's Solid Earth Natural Hazards program. Over the course of this research, we documented measurable uplift of Trident volcano in the Katmai group, conducted a systematic study of the change in phase coherence over time on volcanic surfaces, and measured and modeled the spectacular 1.5 m deflation of Okmok caldera associated with its 1997 eruption. We also generated initial interferograms spanning the 1996 seismic swarm of Akutan volcano; however, during the period covered by this project we were not able to remove topography. That has been done under the subsequent funding and a paper is now in preparation. This report summarizes work done under two separate contracts because both were based on the same proposal to NASA's ADRO (Application Development and Research Opportunity) program. The first year was funded out of a grant from NASA Headquarters and the second and third years out of a grant through Goddard. The work, however, was a continuous three year effort.

Freymueller, Jeff↗

Imaging the Moon's Core with Seismology

Constraining the structure of the lunar core is necessary to improve our understanding of the present-day thermal structure of the interior and the history of a lunar dynamo, as well as the origin and thermal and compositional evolution of the Moon. We analyze Apollo deep moonquake seismograms using terrestrial array processing methods to search for the presence of reflected and converted energy from the lunar core. Although moonquake fault parameters are not constrained, we first explore a suite of theoretical focal spheres to verify that fault planes exist that can produce favorable core reflection amplitudes relative to direct up-going energy at the Apollo stations. Beginning with stacks of event seismograms from the known distribution of deep moonquake clusters, we apply a polarization filter to account for the effects of seismic scattering that (a) partitions energy away from expected components of ground motion, and (b) obscures all but the main P- and S-wave arrivals. The filtered traces are then shifted to the predicted arrival time of a core phase (e.g. PcP) and stacked to enhance subtle arrivals associated with the Moon s core. This combination of filtering and array processing is well suited for detecting deep lunar seismic reflections, since we do not expect scattered wave energy from near surface (or deeper) structure recorded at varying epicentral distances and stations from varying moonquakes at varying depths to stack coherently. Our results indicate the presence of a solid inner and fluid outer core, overlain by a partial-melt-containing boundary layer (Table 1). These layers are consistently observed among stacks from four classes of reflections: P-to-P, S-to-P, P-to-S, and S-to-S, and are consistent with current indirect geophysical estimates of core and deep mantle properties, including mass, moment of inertia, lunar laser ranging, and electromagnetic induction. Future refinements are expected following the successful launch of the GRAIL lunar orbiter and SELENE 2 lunar lander missions.

Weber, Renee C.↗

Geologic Provinces Beneath the Greenland Ice Sheet Constrained by Geophysical Data Synthesis

Present understanding of Greenland's subglacial geology is derived mostly from interpolation of geologic mapping of its ice-free margins and unconstrained by geophysical data. Here we refine the extent of its geologic provinces by synthesizing geophysical constraints on subglacial geology from seismic, gravity, magnetic and topographic data. North of 72°N, no province clearly extends across the whole island, leaving three distinct subglacial regions yet to be reconciled with margin geology. Geophysically coherent anomalies and apparent province boundaries are adjacent to the onset of faster ice flow at both Petermann Glacier and the Northeast Greenland Ice Stream. Separately, based on their subaerial expression, dozens of unusually long, straight and sub-parallel subglacial valleys cross Greenland's interior and are not yet resolved by current syntheses of its subglacial topography.

Joseph A. MacGregor↗

Seismicity as a guide to global tectonics and earthquake prediction.

From seismicity studies, evidence is presented for several aspects of plate-tectonic theory, including ideas of sea-floor spreading, transform faulting and underthrusting of the lithosphere in island arcs. Recent advances in seismic instrumentation, the use of computers in earthquake location, and the installation of local networks of instruments are shown to have vastly increased the data available for seismicity studies. It is pointed out that most of the world's earthquakes are located in very narrow zones along active plate margins and are intimately related to global processes in an extremely coherent manner. Important areas of uncertainty calling for further studies are also pointed out.

Sykes, L. R.↗

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↗

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗

Nowcasting Earthquakes in Southern California With Machine Learning: Bursts, Swarms, and Aftershocks May Be Related to Levels of Regional Tectonic Stress

Abstract Seismic bursts in Southern California are sequences of small earthquakes strongly clustered in space and time and include seismic swarms and aftershock sequences. A readily observable property of these events, the radius of gyration ( R G ), allows us to connect the bursts to the temporal occurrence of the largest M ≥ 7 earthquakes in California since 1984. In the Southern California earthquake catalog, we identify hundreds of these potentially coherent space‐time structures in a region defined by a circle of radius 600 km around Los Angeles. We compute R G for each cluster then filter them to identify those bursts with large numbers of events closely clustered in space, which we call “compact” bursts. Our basic assumption is that these compact bursts reflect the dynamics associated with large earthquakes. Once we have filtered the burst catalog, we apply an exponential moving average to construct a time series for the Southern California region. We observe that the R G of these bursts systematically decreases prior to large earthquakes, in a process that we might term “radial localization.” The R G then rapidly increases during an aftershock sequence, and a new cycle of “radial localization” then begins. These time series display cycles of recharge and discharge reminiscent of seismic stress accumulation and release in the elastic rebound process. The complex burst dynamics we observe are evidently a property of the region as a whole, rather than being associated with individual faults. This new method allows us to improve earthquake nowcasting, which is a technique to evaluate the current state of hazard in a seismically active region.

58 GEOSCIENCES↗

Analysis of Temperature Effects on GS-13 Short Period Sensors Utilizing Sandia's Seismic Sensor Temperature Testbed

Sandia National Laboratories has fabricated a Seismic Sensor Temperature Testbed (SSTT) suitable for exposing sensors under test to a range of reasonably stable temperatures while a co-located reference sensor is maintained at room temperature. The testbed has proven sufficiently quiet to allow recording of high-coherence signals from regional earthquakes in the passband of the seismometers, allowing a direct comparison of signals between the sensors under test and the reference sensor.

47 OTHER INSTRUMENTATION↗

Nested-crater model of lunar ringed basins

We propose a model for the origin of impact-basin rings whereby the main topographic rim of a basin approximates the limit of excavation and inner rings approximate the rims of craters formed inside the transient crater by some perturbation in the cratering process. The cause of this complexity in transient cavities may be the presence of discontinuities in the target material. The second inward ring may have formed at the seismic discontinuity about 20 km deep in the lunar crust, and the third, innermost ring of a few large basins at the crust-mantle interface about 60 km deep. Slumping increased the original diameters of many rings and split some initially coherent rings into subsidiary or partial rings. Deformation outside the transient crater produced external arcs. This model differs from prevalent hypotheses of ring formation whereby an inner ring approximates the transient crater rim and major faulting of the flank produced the outer ring structures.

Wilhelms, D. E.↗