P-wave tomographic model from local bulletin data for improved seismic location in and around Israel

2020 ◽  
Author(s):  
Lewis Schardong ◽  
Yochai Ben-Horin ◽  
Alon Ziv ◽  
Hillel Wust-Bloch ◽  
Yael Radzyner
Keyword(s):  
Eastern Mediterranean ◽  
A Priori ◽  
Warning System ◽  
Velocity Model ◽  
P Wave ◽  
Seismic Velocities ◽  
Fast Marching ◽  
Earthquake Early Warning System ◽  
Starting Model ◽  
Correction Terms

<p>For the past 40 years, the Geophysical Institute of Israel has been in charge of the recording, monitoring and relocating of local earthquakes. Due to the variety of data analysts and data sources, as well as several network upgrades, the resulting bulletin data has to be completed and homogenised, and station metadata needs to be tracked down, and sometimes corrected. For those reasons, as well as because of the lack of consensus on an accurate model for seismic velocities in the area, published source locations are often poorly constrained. We present a homogenised Israeli bulletin, including natural and man-made explosion data. We extract sets of seismic sources with location accuracy greater than 5 km (GT5), as well as GT0 explosions.</p><p>We select a set of events with the highest network coverage, comprising (1) natural earthquakes, (2) man-made quarry or mine blasts, (3) GT5 earthquakes or explosions, and (4) GT0 explosions. We relocate them altogether using the <em>BayesLoc</em> package, a Bayesian, hierarchical, multi-event locator which produces, after source relocation, event-, station- and phase-specific correction terms. We put different a priori constraints on the different categories of seismic events, allowing poorly constrained origin parameters to improve thanks to the more accurate GT locations. <em>BayesLoc</em> also produces traveltime correction terms that can be used to correct systematic errors in the dataset, as well as error estimates.</p><p>Eventually, we invert this homogenised local traveltime dataset in order to invert for a <em>P</em>-wave crustal velocity model of Israel and its surroundings. To do so, we use the <em>Fast Marching Tomography</em> package, which allows the representation of a wide variety of input structures (starting model and geometry of layer boundaries) and can take many different types of input data. We show preliminary inversion tests and results that are in good agreement with past local studies.</p><p>This crustal model of Israel is ultimately to be used as a starting model in a larger tomographic study of the Eastern Mediterranean and Middle East region, where the <em>Regional Seismic Travel Time</em> approach is to be expanded, in order to improve the CTBT’s capabilities in monitoring the regional seismicity. Eventually, such a velocity model could also be used to relocate the whole earthquake catalogue more accurately, and improve the Earthquake Early Warning System currently in development in Israel.</p>

Elastic reflection waveform inversion based on the decomposition of sensitivity kernels

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R235-R250 ◽  
Author(s):  
Zhiming Ren ◽  
Zhenchun Li ◽  
Bingluo Gu
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Background Model ◽  
Reverse Time ◽  
S Wave ◽  
Velocity Models ◽  
Time Migration ◽  
Starting Model

Full-waveform inversion (FWI) has the potential to obtain an accurate velocity model. Nevertheless, it depends strongly on the low-frequency data and the initial model. When the starting model is far from the real model, FWI tends to converge to a local minimum. Based on a scale separation of the model (into the background model and reflectivity model), reflection waveform inversion (RWI) can separate out the tomography term in the conventional FWI kernel and invert for the long-wavelength components of the velocity model by smearing the reflected wave residuals along the transmission (or “rabbit-ear”) paths. We have developed a new elastic RWI method to build the P- and S-wave velocity macromodels. Our method exploits a traveltime-based misfit function to highlight the contribution of tomography terms in the sensitivity kernels and a sensitivity kernel decomposition scheme based on the P- and S-wave separation to suppress the high-wavenumber artifacts caused by the crosstalk of different wave modes. Numerical examples reveal that the gradients of the background models become sufficiently smooth owing to the decomposition of sensitivity kernels and the traveltime-based misfit function. We implement our elastic RWI in an alternating way. At each loop, the reflectivity model is generated by elastic least-squares reverse time migration, and then the background model is updated using the separated traveltime kernels. Our RWI method has been successfully applied in synthetic and real reflection seismic data. Inversion results demonstrate that the proposed method can retrieve preferable low-wavenumber components of the P- and S-wave velocity models, which are reliable to serve as a starting model for conventional elastic FWI. Also, our method with a two-stage inversion workflow, first updating the P-wave velocity using the PP kernels and then updating the S-wave velocity using the PS kernels, is feasible and robust even when P- and S-wave velocities have different structures.

scholarly journals The scaling of PGA with IV2p and its potential for Earthquake Early Warning in Thessaly (Central Greece)

2021 ◽  
Vol 58 ◽  
pp. 177
Author(s):  
Ioannis Spingos ◽  
Filippos Vallianatos ◽  
George Kaviris
Keyword(s):  
Early Warning ◽  
Scaling Laws ◽  
Warning System ◽  
Seismic Signal ◽  
Earthquake Early Warning ◽  
P Wave ◽  
Scaling Relations ◽  
Initial Portion ◽  
Earthquake Early Warning System ◽  
Single Station

The main goal of an Earthquake Early Warning System (EEWS) is to estimate the expected peak ground motion of the destructive S-waves using the first few seconds of P-waves, thus becoming an operational tool for real-time seismic risk management in a short timescale. EEWSs are based on the use of scaling relations between parameters measured on the initial portion of the seismic signal, after the arrival of the first wave. Herein, using the abundant seismicity that followed the 3 March 2021 Mw=6.3 earthquake in Thessaly we propose scaling relations for PGA, from data recorded by local permanent stations, as a function of the integral of the squared velocity (IV2p). The IV2p parameter was estimated directly from the first few seconds-long signal window (tw) after the P-wave arrival. Scaling laws are extrapolated for both individual and across sites (i.e., between a near-source reference instrument and a station located close to a target). The latter approach is newly investigated, as local site effects could have a significant impact on recorded data. Considering that further study on the behavior of IV2p is necessary, there are indications that this parameter could be used in future on-site single‐station earthquake early warning operations for areas affected by earthquakes located in Thessaly, as itpresents significant stability.

scholarly journals A P-wave velocity model of the upper crust of the Sannio region (Southern Apennines, Italy)

1998 ◽  
Vol 41 (4) ◽  
Author(s):  
G. Iannaccone ◽  
L. Improta ◽  
P. Capuano ◽  
A. Zollo ◽  
G. Biella ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Velocity ◽  
Carbonate Platform ◽  
Velocity Model ◽  
P Wave ◽  
Seismic Velocities ◽  
Upper Crust ◽  
Near Surface ◽  
P Wave Velocity ◽  
Apulia Platform

This paper describes the results of a seismic refraction profile conducted in October 1992 in the Sannio region, Southern Italy, to obtain a detailed P-wave velocity model of the upper crust. The profile, 75 km long, extended parallel to the Apenninic chain in a region frequently damaged in historical time by strong earthquakes. Six shots were fired at five sites and recorded by a number of seismic stations ranging from 41 to 71 with a spacing of 1-2 km along the recording line. We used a two-dimensional raytracing technique to model travel times and amplitudes of first and second arrivals. The obtained P-wave velocity model has a shallow structure with strong lateral variations in the southern portion of the profile. Near surface sediments of the Tertiary age are characterized by seismic velocities in the 3.0-4.1 km/s range. In the northern part of the profile these deposits overlie a layer with a velocity of 4.8 km/s that has been interpreted as a Mesozoic sedimentary succession. A high velocity body, corresponding to the limestones of the Western Carbonate Platform with a velocity of 6 km/s, characterizes the southernmost part of the profile at shallow depths. At a depth of about 4 km the model becomes laterally homogeneous showing a continuous layer with a thickness in the 3-4 km range and a velocity of 6 km/s corresponding to the Meso-Cenozoic limestone succession of the Apulia Carbonate Platform. This platform appears to be layered, as indicated by an increase in seismic velocity from 6 to 6.7 km/s at depths in the 6-8 km range, that has been interpreted as a lithological transition from limestones to Triassic dolomites and anhydrites of the Burano formation. A lower P-wave velocity of about 5.0-5.5 km/s is hypothesized at the bottom of the Apulia Platform at depths ranging from 10 km down to 12.5 km; these low velocities could be related to Permo-Triassic siliciclastic deposits of the Verrucano sequence drilled at the bottom of the Apulia Platform in the Apulia Foreland.

scholarly journals Traces of the crustal units and the upper mantle structure in the southwestern part of the East European Craton

2014 ◽  
Vol 6 (1) ◽  
pp. 985-1021
Author(s):  
I. Janutyte ◽  
E. Kozlovskaya ◽  
M. Majdanski ◽  
P. H. Voss ◽  
M. Budraitis ◽  
...  
Keyword(s):  
Upper Mantle ◽  
Velocity Model ◽  
P Wave ◽  
Seismic Velocities ◽  
P Waves ◽  
Teleseismic Tomography ◽  
East European ◽  
Seismic Experiment ◽  
Passive Seismic ◽  
Synthetic Datasets

Abstract. The presented study is a part of the passive seismic experiment PASSEQ 2006–2008 which took place around the Trans-European Suture Zone (TESZ) from May 2006 to June 2008. The dataset of 4195 manually picked arrivals of teleseismic P waves of 101 earthquakes (EQs) recorded in the PASSEQ seismic stations deployed to the east of the TESZ was inverted using the non-linear teleseismic tomography algorithm TELINV. Two 3-D crustal models were used to estimate the crustal travel time (TT) corrections. As a result, we obtained a model of P wave velocity variations in the upper mantle beneath the TESZ and the EEC. In the study area beneath the craton we observed 5 to 6.5% higher and beneath the TESZ about 4% lower seismic velocities compared to the IASP91 velocity model. We found the seismic lithosphere-asthenosphere boundary (LAB) beneath the TESZ at a depth of about 180 km, while we observed no seismic LAB beneath the EEC. The inversion results obtained with the real and the synthetic datasets indicated a ramp shape of the LAB in the northern TESZ where we observed values of seismic velocities close to those of the craton down to about 150 km. The lithosphere thickness in the EEC increases going from the TESZ to the NE from about 180 km beneath Poland to 300 km or more beneath Lithuania. Moreover, in western Lithuania we possibly found an upper mantle dome. In our results the crustal units are not well resolved. There are no clear indications of the features in the upper mantle which could be related with the crustal units in the study area. On the other hand, at a depth of 120–150 km we possibly found a trace of a boundary of proposed palaeosubduction zone between the East Lithuanian Domain (EL) and the West Lithuanian Granulite Domain (WLG). Also, in our results we may have identified two anorogenic granitoid plutons.

scholarly journals Travel-time tomography imaging the Ecuadorian subduction, north of the Mw 7.8 Pedernales earthquake

2021 ◽  
Author(s):  
Alexandra Skrubej ◽  
Audrey Galve ◽  
Mireille Laigle ◽  
Andreas Rietbrock ◽  
Philippe Charvis ◽  
...  
Keyword(s):  
Oceanic Crust ◽  
A Priori ◽  
Structural Complexity ◽  
Velocity Model ◽  
P Wave ◽  
Ocean Bottom ◽  
Seismic Reflection Profile ◽  
Seismic Structure ◽  
Aseismic Slip ◽  
Velocity Models

<p>The Ecuadorian subduction regularly hosts large earthquakes. Among them, the Mw 8.8 1906 earthquake is the 7th biggest known event. Following the recent 2016 Mw 7.8 Pedernales earthquake, a large deployment of onshore/offshore seismological stations, in addition to the permanent seismological/geodetical network, revealed a complex slip behavior including the presence  of  seismic and aseismic slip.</p><p>During the geophysical experiment HIPER, in march 2020, 47 Ocean Bottom Seismometers (OBS), were densely deployed along a 93-km-long trench-perpendicular profile, recording airgun shots (4990 cu.inch.) performed by R/V Atalante to obtain a high-resolution P-wave velocity image. The profile was located north of the 2016 Pedernales rupture zone passing through an area experiencing aseismic slip and a region of contrasted geodetic interseismic coupling.    </p><p>We used the traveltime tomography code « tomo2d » (Korenaga et al., 2000) to invert first arrivals and reflected phases recorded by our OBS.  A joint 2D-seismic-reflection profile was acquired (abstract by L. Schenini) and provides details on the oceanic basement topography and on Vp velocities in shallow sedimentary layers.</p><p>Regarding the structural complexity in the region, we decided to start the inversion  using an a priori 2D velocity model. Several geophysical experiments have already been conducted offshore-onshore Ecuador (SISTEUR, 2000 ; SALIERI, 2001 and ESMERALDAS, 2005). Compilation of velocity models from tomographic images were used to build two a priori 1D Vp velocity models for both the Nazca oceanic crust and the forearc seismic structure. A 2D a priori Vp velocity model was built by merging the results of the two localized inversions using a selection of OBS on each side of the trench.</p><p>We obtain the crustal structure of the upper and subducting plates down to 20 km depth. Beneath the trench, a ~30-km-wide low-Vp anomaly is observed at lithospheric scale. This velocity is 10% lower than the typical Vp values observed for hydrated Pacific-type oceanic crust near the trench (Grevemeyer et al., 2018). Recorded PmP phases will allow us to further constrain the crustal thickness. While we observe PmP phases in areas of low-Vp, the Moho reflectivity weakens and even disappears from the coincident MCS line. This intriguing observation could highlight processes, such as the presence of fluids or serpentinization, that need to be identified and better understood.</p>

scholarly journals Remote Detection of the Electric Field Change Induced at the Seismic Wave Front from the Start of Fault Rupturing

2011 ◽  
Vol 2011 ◽  
pp. 1-11 ◽  
Author(s):  
Yukio Fujinawa ◽  
Kozo Takahashi ◽  
Yoichi Noda ◽  
Hiroshi Iitaka ◽  
Shinobu Yazaki
Keyword(s):  
Electric Field ◽  
Seismic Waves ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Warning System ◽  
Wave Mode ◽  
P Wave ◽  
Remote Detection ◽  
Electrokinetic Effect ◽  
Earthquake Early Warning System

Seismic waves are generally observed through the measurement of undulating elastic ground motion. We report the remote detection of the Earth's electric field variations almost simultaneously with the start of fault rupturing at about 100 km from the fault region using a special electric measurement. The rare but repeated detection indicates that the phenomenon is real. The characteristic time of diffusion is almost instantaneous, that is, less than 1 second to travel 100 km, more than ten times faster than ordinary seismic P wave propagation. We suggest that the measured electric field changes are produced by the electrokinetic effect through increased pore water pressure of the seismic pulse. It is also suggested that the long range propagation is due to the surface wave mode confined near the interface of the different conductivity. The length scale of the finite strength of the electric field is 16 km, 160 km for electric conductivity of 0.01, 0.001, Sm−1, respectively. This phenomenon suggests a new seismic sensing method and a new earthquake early warning system providing more seconds of lead time.

scholarly journals Development of Earthquake Early Warning System Using ADXL335 Accelerometer

2018 ◽  
Author(s):  
Yoga Priyana ◽  
Folkes E. Laumal ◽  
Emir E. Husni
Keyword(s):  
Early Warning ◽  
Early Warning System ◽  
Seismic Waves ◽  
Human Life ◽  
Warning System ◽  
Earthquake Early Warning ◽  
P Wave ◽  
Earthquake Early Warning System ◽  
Wave Data ◽  
Seismic Sensors

Indonesia is an archipelago located at three earthquake belts. This condition cause an earthquake can occur anytime and threaten human life. A quick and accurate early warning system by using the seismic wave data processing is required so, the number of victims affected by the earthquake can be shortened. Here, ADXL335 accelerometers are used as seismic sensors with an Arduino minimum system. The results show that when the first earthquake’s vibration occurs, P wave data detected by the ADXL335 sensor is successfully buffered, calibrated, transmitted and displayed on the server. When there are errors on the transmission, server will request for retransmission. The alarm of the earthquake early warning system will be activated if there are at least three sensors from different locations successfully transmit P wave data with the same scale. This is needed to prevent fake seismic waves.

scholarly journals Traces of the crustal units and the upper-mantle structure in the southwestern part of the East European Craton

Solid Earth ◽  
2014 ◽  
Vol 5 (2) ◽  
pp. 821-836 ◽  
Author(s):  
I. Janutyte ◽  
E. Kozlovskaya ◽  
M. Majdanski ◽  
P. H. Voss ◽  
M. Budraitis ◽  
...  
Keyword(s):  
Upper Mantle ◽  
Synthetic Data ◽  
Velocity Model ◽  
East European Craton ◽  
P Wave ◽  
Seismic Velocities ◽  
P Waves ◽  
Data Set ◽  
Teleseismic Tomography ◽  
East European

Abstract. The presented study is a part of the passive seismic experiment PASSEQ 2006–2008, which took place around the Trans-European Suture Zone (TESZ) from May 2006 to June 2008. The data set of 4195 manually picked arrivals of teleseismic P waves of 101 earthquakes (EQs) recorded in the seismic stations deployed to the east of the TESZ was inverted using the non-linear teleseismic tomography algorithm TELINV. Two 3-D crustal models were used to estimate the crustal travel time (TT) corrections. As a result, we obtain a model of P-wave velocity variations in the upper mantle beneath the TESZ and the East European Craton (EEC). In the study area beneath the craton, we observe up to 3% higher and beneath the TESZ about 2–3% lower seismic velocities compared to the IASP91 velocity model. We find the seismic lithosphere–asthenosphere boundary (LAB) beneath the TESZ at a depth of about 180 km, while we observe no seismic LAB beneath the EEC. The inversion results obtained with the real and the synthetic data sets indicate a ramp shape of the LAB in the northern TESZ, where we observe values of seismic velocities close to those of the craton down to about 150 km. The lithosphere thickness in the EEC increases going from the TESZ to the NE from about 180 km beneath Poland to 300 km or more beneath Lithuania. Moreover, in western Lithuania we find an indication of an upper-mantle dome. In our results, the crustal units are not well resolved. There are no clear indications of the features in the upper mantle which could be related to the crustal units in the study area. On the other hand, at a depth of 120–150 km we indicate a trace of a boundary of proposed palaeosubduction zone between the East Lithuanian Domain (EL) and the West Lithuanian Granulite Domain (WLG). Also, in our results, we may have identified two anorogenic granitoid plutons.

Sign in / Sign up

Export Citation Format

Share Document