Seismo-Acoustic Shockwave Isolation for Low-Yield Local Explosions

Author(s):  
Joshua Dickey ◽  
Michael Pasyanos ◽  
Richard Martin ◽  
Raúl Peña

<p>Seismic and acoustic recordings have long been used for the forensic analysis of various natural and anthropogenic events, especially in the realm of nuclear treaty monitoring. More recently, multi-phenomenological analysis has been applied to these signals with great success, providing unique constraints for studying a broad range of source events, including man-made noise, earthquakes and explosions. In particular, the fusion of seismic and infrasonic data has proven valuable for the analysis of explosive yield, significantly improving on the yield estimates obtained from either seismic or acoustic analysis alone.</p><p>Unfortunately, the seismo-acoustic analysis of local explosions is complicated by the fact that the two phenomena are potentially co-dependent. Large seismic waves displace the earth like a piston, potentially inducing acoustic waves into the atmosphere as they pass. Similarly, large acoustic waves can couple into the earth, inducing ground motion along their path. This co-dependence can be problematic, particularly when the passing acoustic shockwave couples into the earth coincident with a seismic phase arrival, thereby corrupting the signal.</p><p>To address this problem, we present a method for isolating the shockwave response of a seismic sensor, such that any underlying seismic phase arrivals can be recovered. This is accomplished by employing the adaptive noise cancellation model, where a co-located infrasound sensor is used as a reference measurement for the shockwave. In this model, the adaptive filter learns the transform between the relative atmospheric pressure (as recorded by the infrasound sensor), and the resulting ground motion (as recorded by the seismometer). In this way, the filtered infrasound recording approximates the seismic shockwave response, and can be subtracted from the seismograph to recover the phase arrivals.</p><p>The experimental data comes from a set of three low-yield near-surface chemical explosions conducted by LLNL as part of a field experiment, known as FE2. The explosions were recorded at eight stations, located at varying distances from the source (between 64m and 2km), with each station consisting of a co-located three-component seismic velocity transducer and differential infrasound sensor. The adaptive technique is demonstrated for recovering seismic arrivals in both the vertical and horizontal channels across all eight stations, and evaluated using leave-one-out cross-validation across the three explosions.</p>

Author(s):  
Athanasius Cipta ◽  
Phil Cummins ◽  
Masyhur Irsyam ◽  
Sri Hidayati

We use earthquake ground motion modelling via Ground Motion Prediction Equations (GMPEs) and numerical simulation of seismic waves to consider the effects of site amplification and basin resonance in Jakarta, the capital city of Indonesia. While spectral accelerations at short periods are sensitive to near-surface conditions (i.e., Vs30), our results suggest that, for basins as deep as Jakarta’s, available GMPEs cannot be relied upon to accurately estimate the effect of basin depth on ground motions at long periods (>1 s). Amplitudes at such long periods are influenced by entrapment of seismic waves in the basin, resulting in longer duration of strong ground motion, and interference between incoming and reflected waves as well as focusing at basin edges may amplify seismic waves. In order to simulate such phenomena in detail, a basin model derived from a previous study is used as a computational domain for deterministic earthquake scenario modeling in a 2-dimensional cross-section. A Mw 9.0 megathrust, a Mw 6.5 crustal thrust and a Mw 7.0 instraslab earthquake are chosen as scenario events that pose credible threats to Jakarta, and the interactions with the basin of seismic waves generated by these events were simulated. The highest PGV amplifications are recorded at sites near the middle of the basin and near its southern edge, with maximum amplifications of PGV in the horizontal component of 200% for the crustal, 600% for the megathrust and 335% for the deep intraslab earthquake scenario, respectively. We find that the levels of ground motion response spectral acceleration fall below those of the 2012 Indonesian building Codes's design response spectrum for short periods (< 1 s), but closely approach or may even exceed these levels for longer periods.


Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. Q37-Q48 ◽  
Author(s):  
Joachim Place ◽  
Deyan Draganov ◽  
Alireza Malehmir ◽  
Christopher Juhlin ◽  
Chris Wijns

Exhumation of crust exposes rocks to weathering agents that weaken the rocks’ mechanical strength. Weakened rocks will have lower seismic velocity than intact rocks and can therefore be mapped using seismic methods. However, if the rocks are heavily weathered, they will attenuate controlled-source seismic waves to such a degree that the recorded wavefield would become dominated by ambient noise and/or surface waves. Therefore, we have examined the structure of differential weathering by first-break traveltime tomography over a seismic profile extending approximately 3.5 km and acquired at a mining site in Zambia using explosive sources and a source based on the swept-impact seismic technique (SIST). Seismic interferometry has been tested for the retrieval of supervirtual first arrivals masked by uncorrelated noise. However, use of crosscorrelation in the retrieval process makes the method vulnerable to changes in the source signal (explosives and SIST). Thus, we have developed a crosscoherence-based seismic-interferometry method to tackle this shortcoming. We investigate the method’s efficiency in retrieving first arrivals and, simultaneously, correctly handling variations in the source signal. Our results illustrate the superiority of the crosscoherence- over crosscorrelation-based method for retrieval of the first arrivals, especially in alleviating spurious ringyness and in terms of the signal-to-noise ratio. These benefits are observable in the greater penetration depth and the improved resolution of the tomography sections. The tomographic images indicate isolated bodies of higher velocities, which may be interpreted as fresh rocks embedded into a heavily weathered regolith, providing a conspicuous example of differential weathering. Our study advances the potential of seismic methods for providing better images of the near surface (the critical zone).


2021 ◽  
Author(s):  
Hossam Marzouk ◽  
Tarek Arafa-Hamed ◽  
Michael Becken ◽  
Mohamed Abdel Zaher ◽  
Matthew Comeau

<p>We present electrical resistivity models of the crust and upper mantle estimated from 2D inversions of broadband magnetotellurics (MT) data acquired from two profiles in the western desert of Egypt, which can contribute to the understanding of the structural setup of this region. The first profile data are collected from 14 stations along a 250 km profile, in EW direction profile runs along latitude ~25.5°N from Kharga oasis to Dakhla oasis. The second profile comprises 19 stations measured along a 130 km profile in NS direction centered at longitude 28°E and crossing the Farafra. The acquisition for both profiles continued for 1 to 3 days at each station, which allowed for the calculation of impedances for periods from 0.01 sec up to  4096 sec at some sites. The wide frequency band corresponds to a maximal skin depths of up to 150 km that can provide penetration to the top of the asthenosphere. The inversion models display high-conductivity sediments cover at the near surface (<1-2 km), which can be associated with the Nubian aquifer. Along the EW-profile from Kaharge to Dhakla, the crustal basement is overly highly resistive and homogeneous und underlain by a more conductive lithospheric mantle below depths of 30-40 km. Along the N-S profile across Farafra, only the southern portion exhibits a highly resistive crust, whereas beneath Farafra northwards, moderate crustal conductivities are encountered. A comparison has been made between the resultant resistivity models with the 1° tessellated updated crust and lithospheric model of the Earth (LITHO1.0) which was developed by <em>Pasyanos, 2014</em> on the basis of seismic velocity data. The obtained results show a remarkable consistency between the resistivity models and the calculated crustal boundaries. Especially at the Kharga-Dakhla profile a clear matching can be noticed at the upper and lower boundaries of a characteristic anomaly with the Moho and LAB boundaries respectively.</p>


2019 ◽  
Vol 7 (1) ◽  
pp. T141-T154 ◽  
Author(s):  
Md. Iftekhar Alam

Seismic imaging of the shallow subsurface (approximately 5 m) can be very challenging when reflections are absent and the data are dominated by ground roll. I analyzed the transmission coda to produce fine-scale, interpretable vertical and horizontal component seismic velocity ([Formula: see text] and [Formula: see text]) models using full-waveform inversion (FWI). Application of FWI is tested through imaging two buried targets. The first target is a pair of well-documented utility pipes with known diameters (0.8 m) and burial depths (approximately 1.5 m). The second target is a poorly documented former location of the pipe(s), which is now a backfilled void. Data are acquired along a 23 m 2D profile using a static array with single-component vertical and horizontal geophones. Our results indicate considerable velocity updates in the [Formula: see text] and [Formula: see text] models across the pipes and backfill. The pipes appear as negative velocity updates in the final inverted [Formula: see text] and [Formula: see text] models, whereas the backfilled area represents negative and positive velocity updates in the [Formula: see text] and [Formula: see text] models, respectively. Variations of the polarities in the inverted models ([Formula: see text] and [Formula: see text]) across the backfill can be indicative of the medium, which respond differently to the vertical and horizontal component seismic waves. The attenuation models show a general decreasing trend with increasing depth. Therefore, simultaneous applications of vertical ([Formula: see text]) and horizontal ([Formula: see text]) component seismic velocity modeling can be an effective tool to understand the subsurface medium in near-surface characterization.


2020 ◽  
Author(s):  
Shahar Shani-Kadmiel ◽  
Gil Averbuch ◽  
Pieter Smets ◽  
Jelle Assink ◽  
Läslo Evers

<p>The amplitude of ground motions caused by earthquakes and subsurface explosions generally decreases with distance from the epicenter. However, in the near-source region, other factors, e.g., near surface geology, topography, and the source radiation pattern, may significantly vary the amplitude of ground motions. Although source location and magnitude (or yield), can be rapidly determined using distant seismic stations, without a dense seismological network in the epicentral region, the ability to resolve such variations is limited.</p><p>Besides seismic waves, earthquakes and subsurface explosions generate infrasound, i.e., inaudible acoustic waves in the atmosphere. The mechanical ground motions from such sources, including the effects from the above mentioned factors, are encapsulated by the acoustic pressure perturbations over the source region. Due to the low frequency nature of infrasound and facilitated by waveguides in the atmosphere, such perturbations propagate over long ranges with limited attenuation and are detected at ground-based stations. In this work we demonstrate a method for resolving ground motions and the source mechanism from remotely detected infrasound. This is illustrated for the 2010 Mw 7.0 Port-au-Prince, Haiti earthquake, and the 6th and largest nuclear test conducted by the Democratic People's Republic of Korea in 2017.</p><p>Such observations are made possible by: (1) An advanced array processing technique that enables the detection of coherent wavefronts, even when amplitudes are below the noise level, and (2) A backprojection technique that maps infrasound detections in time to their origin on the Earth's surface.</p><p>Infrasound measurements are conducted globally for the verification of the Comprehensive Nuclear-Test-Ban Treaty and together with regional infrasound networks allow for an unprecedented global coverage. This makes infrasound as an earthquake disaster mitigation technique feasible for the first time and contributes to the Treaty's verification capacity.</p>


2021 ◽  
Vol 13 (14) ◽  
pp. 2684
Author(s):  
Eldert Fokker ◽  
Elmer Ruigrok ◽  
Rhys Hawkins ◽  
Jeannot Trampert

Previous studies examining the relationship between the groundwater table and seismic velocities have been guided by empirical relationships only. Here, we develop a physics-based model relating fluctuations in groundwater table and pore pressure with seismic velocity variations through changes in effective stress. This model justifies the use of seismic velocity variations for monitoring of the pore pressure. Using a subset of the Groningen seismic network, near-surface velocity changes are estimated over a four-year period, using passive image interferometry. The same velocity changes are predicted by applying the newly derived theory to pressure-head recordings. It is demonstrated that the theory provides a close match of the observed seismic velocity changes.


2014 ◽  
Vol 501-504 ◽  
pp. 1447-1452
Author(s):  
Yan Yan Yu ◽  
Qi Fang Liu

Seismic response of the Shidian basin to moderate scenario earthquake is investigated considering 3D basin model incorporated with real topography by using the spectral-element method and parallel computing technique. The wave propagation process, the generation of surface wave, and the impact of soil deposits velocity to the basin-induced surface wave are studied in this paper. The results show that the amplification behavior of the basin is the interactions of basin geometry and low velocity soil deposits. First, locally small hollows in the basin are apt to trap seismic waves and produce much stronger ground motion, basin edge and areas with deep sediments are also characterized with large amplification. Then, basin with softer soil deposits produces stronger surface waves with lower propagation velocity and higher mode.


2021 ◽  
pp. 875529302110513
Author(s):  
Eleftheria Efthymiou ◽  
Alfredo Camara

The definition of the spatial variability of the ground motion (SVGM) is a complex and multi-parametric problem. Its effect on the seismic response of cable-stayed bridges is important, yet not entirely understood to date. This work examines the effect of the SVGM on the seismic response of cable-stayed bridges by means of the time delay of the ground motion at different supports, the loss of coherency of the seismic waves, and the incidence angle of the seismic waves. The focus herein is the effect of the SVGM on cable-stayed bridges with various configurations in terms of their length and of design parameters such as the pylon shape and the pylon–cable system configuration. The aim of this article is to provide general conclusions that are applicable to a wide range of canonical cable-stayed bridges and to contribute to the ongoing effort to interpret and predict the effect of the SVGM in long structures. This work shows that the effect of the SVGM on the seismic response of cable-stayed bridges varies depending on the pylon shape, height, and section dimensions; on the cable-system configuration; and on the response quantity of interest. Furthermore, the earthquake incidence angle defines whether the SVGM is important to the seismic response of the cable-stayed bridges. It is also confirmed that the SVGM excites vibration modes of the bridges that do not contribute to their seismic response when identical support motion is considered.


Sign in / Sign up

Export Citation Format

Share Document