Lithospheric structure of the Arabian Shield from joint inversion of P- and S-wave receiver functions and dispersion velocities

2015 ◽  
Vol 65 (2) ◽  
pp. 239-255 ◽  
Author(s):  
Abdullah M. Al-Amri
Keyword(s):  
Green's Functions ◽  
Velocity Structure ◽  
Green’S Functions ◽  
Joint Inversion ◽  
Receiver Functions ◽  
Gulf Of Aqaba ◽  
Arabian Shield ◽  
S Wave ◽  
Velocity Models ◽  
Seismic Networks

Abstract New velocity models of lithospheric thickness and velocity structure have been developed for the Arabian Shield by three tasks: 1) Computing P-Wave Receiver Functions (PRFs) and S-Wave Receiver Functions (SRFs) for all the broadband stations within the Saudi seismic networks. The number of receiver function waveforms depends on the recording time window and quality of the broadband station. 2) Computing ambient noise correlation Green’s functions for all available station pairs within the Saudi seismic networks to image the shear velocity in the crust and uppermost mantle beneath the Arabian Peninsula. Together they provided hundreds of additional, unique paths exclusively sampling the region of interest. Both phase and group velocities for all the resulting empirical Green’s functions have been measured and to be used in the joint inversion. 3) Jointly inverted the PRFs and SRFs obtained in task 1 with dispersion velocities measured on the Green’s functions obtained in task 2 and with fundamental-mode, Rayleigh-wave, group and phase velocities borrowed from the tomographic studies to precisely determine 1D crustal velocity structure and upper mantle. The analysis of the PRFs revealed values of 25-45 km for crustal thickness, with the thin crust next to the Red Sea and Gulf of Aqaba and the thicker crust under the platform, and Vp/Vs ratios in the 1.70-1.80 range, suggesting a range of compositions (felsic to mafic) for the shield’s crust. The migrated SRFs suggest lithospheric thicknesses in the 80-100 km range for portions of the shield close to the Red Sea and Gulf of Aqaba and near the Arabian Gulf. Generally, the novelty of the velocity models developed under this paper has consisted in the addition of SRF data to extend the velocity models down to lithospheric and sub-lithospheric depths.

Mapping crustal structure of the Nechako–Chilcotin plateau using teleseismic receiver function analysis,

2014 ◽  
Vol 51 (4) ◽  
pp. 407-417 ◽  
Author(s):  
H.S. Kim ◽  
J.F. Cassidy ◽  
S.E. Dosso ◽  
H. Kao
Keyword(s):  
Wave Velocity ◽  
Crustal Structure ◽  
Velocity Structure ◽  
Receiver Function ◽  
Function Analysis ◽  
Crustal Thickness ◽  
Receiver Functions ◽  
S Wave ◽  
Low Velocity ◽  
Velocity Models

This paper presents results of a passive-source seismic mapping study in the Nechako–Chilcotin plateau of central British Columbia, with the ultimate goal of contributing to assessments of hydrocarbon and mineral potential of the region. For the present study, an array of nine seismic stations was deployed in 2006–2007 to sample a wide area of the Nechako–Chilcotin plateau. The specific goal was to map the thickness of the sediments and volcanic cover, and the overall crustal thickness and structural geometry beneath the study area. This study utilizes recordings of about 40 distant earthquakes from 2006 to 2008 to calculate receiver functions, and constructs S-wave velocity models for each station using the Neighbourhood Algorithm inversion. The surface sediments are found to range in thickness from about 0.8 to 2.7 km, and the underlying volcanic layer from 1.8 to 4.7 km. Both sediments and volcanic cover are thickest in the central portion of the study area. The crustal thickness ranges from 22 to 36 km, with an average crustal thickness of about 30–34 km. A consistent feature observed in this study is a low-velocity zone at the base of the crust. This study complements other recent studies in this area, including active-source seismic studies and magnetotelluric measurements, by providing site-specific images of the crustal structure down to the Moho and detailed constraints on the S-wave velocity structure.

scholarly journals Joint Body‐ and Surface‐Wave Tomography of Yucca Flat, Nevada, Using a Novel Seismic Source

2019 ◽  
Vol 109 (5) ◽  
pp. 1922-1934 ◽  
Author(s):  
Liam D. Toney ◽  
Robert E. Abbott ◽  
Leiph A. Preston ◽  
David G. Tang ◽  
Tori Finlay ◽  
...  
Keyword(s):  
Surface Wave ◽  
Phase Velocity ◽  
Wave Velocity ◽  
Velocity Structure ◽  
Joint Inversion ◽  
P Wave ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Velocity Models ◽  
Phase Velocity Dispersion

Abstract In preparation for the next phase of the Source Physics Experiments, we acquired an active‐source seismic dataset along two transects totaling more than 30 km in length at Yucca Flat, Nevada, on the Nevada National Security Site. Yucca Flat is a sedimentary basin which has hosted more than 650 underground nuclear tests (UGTs). The survey source was a novel 13,000 kg modified industrial pile driver. This weight drop source proved to be broadband and repeatable, richer in low frequencies (1–3 Hz) than traditional vibrator sources and capable of producing peak particle velocities similar to those produced by a 50 kg explosive charge. In this study, we performed a joint inversion of P‐wave refraction travel times and Rayleigh‐wave phase‐velocity dispersion curves for the P‐ and S‐wave velocity structure of Yucca Flat. Phase‐velocity surface‐wave dispersion measurements were obtained via the refraction microtremor method on 1 km arrays, with 80% overlap. Our P‐wave velocity models verify and expand the current understanding of Yucca Flat’s subsurface geometry and bulk properties such as depth to Paleozoic basement and shallow alluvium velocity. Areas of disagreement between this study and the current geologic model of Yucca Flat (derived from borehole studies) generally correlate with areas of widely spaced borehole control points. This provides an opportunity to update the existing model, which is used for modeling groundwater flow and radionuclide transport. Scattering caused by UGT‐related high‐contrast velocity anomalies substantially reduced the number and frequency bandwidth of usable dispersion picks. The S‐wave velocity models presented in this study agree with existing basin‐wide studies of Yucca Flat, but are compromised by diminished surface‐wave coherence as a product of this scattering. As nuclear nonproliferation monitoring moves from teleseismic to regional or even local distances, such high‐frequency (>5  Hz) scattering could prove challenging when attempting to discriminate events in areas of previous testing.

Shallow velocity structure and poisson's ratio at the Tarzana, California, strong-motion accelerometer site

1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee
Keyword(s):  
Urban Areas ◽  
Velocity Structure ◽  
Strong Motion ◽  
P Wave ◽  
S Wave ◽  
Near Surface ◽  
Velocity Models ◽  
Wave Velocities ◽  
Poisson's Ratios ◽  
Poisson’S Ratios

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.

Joint receiver function and gravity inversion for new constraints on the Ivrea body structure: a 2D high-resolution view along the Val Sesia profile (N. Italy).

2021 ◽  
Author(s):  
Matteo Scarponi ◽  
György Hetényi ◽  
Jaroslava Plomerová ◽  
Stefano Solarino
Keyword(s):  
High Resolution ◽  
Velocity Structure ◽  
Receiver Function ◽  
Joint Inversion ◽  
Gravity Data ◽  
Rock Physics ◽  
Receiver Functions ◽  
Gravity Inversion ◽  
European Alps ◽  
Horizontal Distance

<p>We present results from a joint inversion study of new seismic and gravity data to constrain a 2D high-resolution image of one of the most prominent geophysical anomalies of the European Alps: the Ivrea geophysical body (IGB). Our work exploits both new data and multidisciplinary a priori constraints, to better resolve the shallow crustal structure in the Ivrea-Verbano zone (IVZ), where the IGB is known to reach anomalously shallow depths and partially outcrop at the surface.</p><p>A variety of previous studies, ranging from gravity surveys to vintage refraction seismics and recent local earthquake tomographies (Solarino et al. 2018, Diehl et al. 2009), provide comprehensive but spatially sparse information on the IGB structure, which we aim at investigating at higher resolution, along a linear profile crossing the IVZ. To this purpose, we deployed 10 broadband seismic stations (MOBNET pool, IG CAS Prague), 5 km spaced along a linear West-East profile, along Val Sesia and crossing Lago Maggiore. This network operated for 27 months and allowed us to produce a new database of ca. 1000 seismic high-quality receiver functions (RFs). In addition, we collected new gravity data in the IVZ, with a data coverage of 1 gravity point every 1-2 km along the seismic profile. The newly collected data was used to set up an inversion scheme, in which RFs and gravity anomalies are jointly used to constrain the shape and the physical property contrasts across the IGB interface.</p><p>We model the IGB as a single interface between far-field constraints, whose geometry is defined by the coordinates of four nodes which may vary in space, and  density and V<sub>S</sub> shear-wave velocity contrasts associated with the interface itself, varying independently. A Markov chain Monte Carlo (MCMC) sampling method with Metropolis-Hastings selection rule was implemented to efficiently explore the model space, directing the search towards better fitting areas.</p><p>For each model, we perform ray-tracing and RFs migration using the actual velocity structure both for migration and computation of synthetic RFs, to be compared with the observations via cross-correlation of the migration images. Similarly, forward gravity modelling for a 2D density distribution is implemented and the synthetic gravity anomaly is compared with the observations along the profile. The joint inversion performance is the product of these two misfits.</p><p>The inversion results show that the IGB reaches the shallowest depths in the western part of the profile, preferentially locating the IGB interface between 3 and 7 km depth over a horizontal distance of ca. 20 km (between Boccioleto and Civiasco, longitudes 8.1 and 8.3). Within this segment, the shallowest point reaches up to 1 km below sea level. The found density and velocity contrasts are in agreement with rock physics properties of various units observed in the field and characterized in earlier studies.</p>

Non‐Double‐Couple Moment Tensors of Earthquakes Calculated Using Empirical Green’s Functions

2019 ◽  
Vol 91 (1) ◽  
pp. 390-398
Author(s):  
Václav Vavryčuk ◽  
Petra Adamová
Keyword(s):  
Green's Functions ◽  
Green’S Functions ◽  
Joint Inversion ◽  
Velocity Model ◽  
Moment Tensors ◽  
Empirical Green’S Functions ◽  
Fracturing Process ◽  
Double Couple ◽  
Focal Zone ◽  
Empirical Green's Functions

Abstract We present a joint inversion for empirical Green’s functions (EGFs) and high‐resolution non‐double‐couple (non‐DC) moment tensors. First, the EGFs are constructed using known moment tensors of earthquakes occurring in a small focal zone. Second, the estimated EGFs are applied to refine the original moment tensors used for constructing the EGFs. Because the EGFs describe the velocity model better than the standard GFs, the refined moment tensors are more accurate. The method is applied to real observations of earthquakes of the 2008 swarm in West Bohemia, Czech Republic, where tiny details in fracturing in the focal zone are revealed. Refined moment tensors indicate fault closing caused by compaction of fault gouge during fracturing process related to fault weakening by fluids in the focal zone. The application of the proposed inversion can improve moment tensors reported in existing local, regional, or global catalogs for areas with a concentrated seismicity.

Sign in / Sign up

Export Citation Format

Share Document