scholarly journals High-precision relocation of seismic sequences above a dipping Moho: the case of the January–February 2014 seismic sequence on Cephalonia island (Greece)

Solid Earth ◽  
2015 ◽  
Vol 6 (1) ◽  
pp. 173-184 ◽  
Author(s):  
V. K. Karastathis ◽  
E. Mouzakiotis ◽  
A. Ganas ◽  
G. A. Papadopoulos

Abstract. Detailed velocity structure and Moho mapping is of crucial importance for a high precision relocation of seismicity occurring out of, or marginal to, the geometry of seismological networks. Usually the seismographic networks do not cover the boundaries of converging plates such as the Hellenic arc. The crustal thinning from the plate boundary towards the back-arc area creates significant errors in accurately locating the earthquake, especially when distant seismic phases are included in the analysis. The case of the Cephalonia (Ionian Sea, Greece) sequence of January–February 2014 provided an excellent example where the hypocentral precision was greatly affected by the crustal thinning from the plate boundary at the Ionian sea towards the Aegean sea. This effect was examined in detail by testing various velocity models of the region in order to determine an optimal model. Our tests resulted in the adoption of a velocity model that resembles the crustal thinning of the region. Then, a relocation procedure was performed in the Cephalonia sequence for the time period of 26 January to 15 May 2014 by applying probabilistic non-linear location algorithms. The high-precision relocation resulted in an improved spatial distribution of the seismicity with respect to the preliminary locations and provided a reliable basis to examine seismotectonic implications of the Cephalonia sequence.

2014 ◽  
Vol 6 (2) ◽  
pp. 2699-2733 ◽  
Author(s):  
V. K. Karastathis ◽  
E. Mouzakiotis ◽  
A. Ganas ◽  
G. A. Papadopoulos

Abstract. Detailed velocity structure and Moho mapping is of crucial importance for a high precision relocation of seismicity occurring out of, or marginally to, the geometry of seismological networks, such as at the boundary of converging plates. The crustal thinning from the plate boundary towards the back-arc area creates significant errors in accurately locating the earthquake, especially when distant seismic phases are included in the analysis. The case of the Cephalonia (Ionian Sea, Greece) sequence of January–February 2014 provided an excellent example where locations were greatly affected by the crustal thinning from the plate boundary at the Ionian sea towards the Aegean sea. This effect was examined in detail by testing various velocity models of the region in order to determine an optimal model. Our tests resulted in the adoption of a velocity model that resembles the crustal thinning of the region. Then, a relocation procedure was performed in the Cephalonia sequence for the time period from 26 January 2014 to 15 May 2014 by applying probabilistic non-linear location algorithms. The high-precision relocation resulted in an improved spatial distribution of the seismicity with respect to the preliminary locations and provided a reliable basis to examine seismotectonic implications of the Cephalonia sequence.


2021 ◽  
Author(s):  
Bryant Chow

<p><b>Seismic tomography is a powerful tool for understanding Earth structure. In New Zealand, velocity models derived using ray-based tomography have been used extensively to characterize the complex plate boundary between the Australian and Pacific plates. Advances in computational capabilities now allow us to improve these velocity models using adjoint tomography, an imaging method which minimizes differences between observed and simulated seismic waveforms. We undertake the first application of adjoint tomography in New Zealand to improve a ray-based New Zealand velocity model containing the Hikurangi subduction zone and the North Island of New Zealand.</b></p> <p>In support of this work we deployed the Broadband East Coast Network (BEACON), a temporary seismic network aimed at improving coverage of the New Zealand permanent network, along the east coast of the North Island. We concurrently develop an automated, open-source workflow for full-waveform inversion using spectral element and adjoint methods. We employ this tool to assess a candidate velocity model’s suitability for adjoint tomography. Using a 3D ray-based traveltime tomography model of New Zealand, we generate synthetic seismic waveforms for more than 10 000 source–receiver pairs and evaluate waveform misfits. We subsequently perform synthetic checkerboard inversions with a realistic New Zealand source–receiver distribution. Reasonable systematic time shifts and satisfactory checkerboard resolution in synthetic inversions indicate that the candidate model is appropriate as an initial model for adjoint tomography. This assessment also demonstrates the relative ease of use and reliability of the automated tools.</p> <p>We then undertake a large-scale adjoint tomography inversion for the North Island of New Zealand using up to 1 800 unique source–receiver pairs to fit waveforms with periods 4–30 s, relating to minimum waveform sensitivities on the order of 5 km. Overall, 60 geographically well-distributed earthquakes and as many as 88 broadband station locations are included. Using a nonlinear optimization algorithm, we undertake 28 model updates of Vp and Vs over six distinct inversion legs which progressively increase resolution. The total inversion incurred a computational cost of approximately 500 000 CPU-hours. The overall time shift between observed and synthetic seismograms is reduced, and updated velocities show as much as ±30% change with respect to initial values. A formal resolution analysis using point spread tests highlights that velocity changes are strongly resolved onland and directly offshore, at depths above 30 km, with low-amplitude changes (> 1%) observed down to 100 km depth. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone.</p> <p>We interpret the updated velocity model in terms of New Zealand tectonics and geology, and observe good agreement with known basement terranes, and major structural elements such as faults, sedimentary basins, broad-scale subduction related features. We recover increased spatial heterogeneity in seismic velocities along the strike of the Hikurangi subduction zone with respect to the initial model. Below the East Coast, we interpret two localized high-velocity anomalies as previously unidentified subducted seamounts. We corroborate this interpretation with other work, and discuss the implications of deeply subducted seamounts on slip behavior along the Hikurangi margin. In the Cook Strait we observe a low-velocity zone that we interpret as a deep sedimentary basin. Strong velocity gradients bounding this low-velocity zone support hypotheses of a structural boundary here separating the North and South Islands of New Zealand. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate seismic velocities at depth to crustal magmatic activity below the Taupo Volcanic Zone.</p> <p>This new velocity model provides more accurate synthetic seismograms and additional constraints on enigmatic tectonic processes related to the North Island of New Zealand. Both the velocity model itself, and the underpinning methodological contributions, improve our ever-expanding understanding of the North Island of New Zealand, the Hikurangi subduction zone, and the broader Australian-Pacific plate boundary.</p>


Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 533-546 ◽  
Author(s):  
Robert G. Clapp ◽  
Biondo L. Biondi ◽  
Jon F. Claerbout

In areas of complex geology, prestack depth migration is often necessary if we are to produce an accurate image of the subsurface. Prestack depth migration requires an accurate interval velocity model. With few exceptions, the subsurface velocities are not known beforehand and should be estimated. When the velocity structure is complex, with significant lateral variations, reflection‐tomography methods are often an effective tool for improving the velocity estimate. Unfortunately, reflection tomography often converges slowly, to a model that is geologically unreasonable, or it does not converge at all. The large null space of reflection‐tomography problems often forces us to add a sparse parameterization of the model and/or regularization criteria to the estimation. Standard tomography schemes tend to create isotropic features in velocity models that are inconsistent with geology. These isotropic features result, in large part, from using symmetric regularization operators or from choosing a poor model parameterization. If we replace the symmetric operators with nonstationary operators that tend to spread information along structural dips, the tomography will produce velocity models that are geologically more reasonable. In addition, by forming the operators in helical 1D space and performing polynomial division, we apply the inverse of these space‐varying anisotropic operators. The inverse operators can be used as a preconditioner to a standard tomography problem, thereby significantly improving the speed of convergence compared with the typical regularized inversion problem. Results from 2D synthetic and 2D field data are shown. In each case, the velocity obtained improves the focusing of the migrated image.


1990 ◽  
Vol 80 (2) ◽  
pp. 395-410 ◽  
Author(s):  
Glenn D. Nelson ◽  
John E. Vidale

Abstract We present a new method for locating earthquakes in a region with arbitrarily complex three-dimensional velocity structure, called QUAKE3D. Our method searches a gridded volume and finds the global minimum travel-time residual location within the volume. Any minimization criterion may be employed. The L1 criterion, which minimizes the sum of the absolute values of travel-time residuals, is especially useful when the station coverage is sparse and is more robust than the L2 criterion (which minimizes the RMS sum) employed by most earthquake location programs. On a UNIX workstation with 8 Mbytes memory, travel-time grids of size 150 by 150 by 50 are reasonably employed, with the actual geographic coverage dependent on the grid spacing. Location precision is finer than the grid spacing. Earthquake recordings at six stations in Bear Valley are located as an example, using various layered and laterally varying velocity models. Locations with QUAKE3D are nearly identical to HYPOINVERSE locations when the same flat-layered velocity model is used. For the examples presented, the computation time per event is approximately 4 times slower than HYPOINVERSE, but the computation time for QUAKE3D is dependent only on the grid size and number of stations, and independent of the velocity model complexity. Using QUAKE3D with a laterally varying velocity model results in locations that are physically more plausible and statistically more precise. Compared to flat-layered solutions, the earthquakes are more closely aligned with the surface fault trace, are more uniform in depth distribution, and the event and station travel-time residuals are much smaller. Hypocentral error bars computed by QUAKE3D are more realistic in that the trade-off of depth versus origin time is implicit in our error estimation, but ignored by HYPOINVERSE.


2017 ◽  
Vol 54 (2) ◽  
pp. 163-172 ◽  
Author(s):  
Shutian Ma ◽  
Pascal Audet

Models of the seismic velocity structure of the crust in the seismically active northern Canadian Cordillera remain poorly constrained, despite their importance in the accurate location and characterization of regional earthquakes. On 29 August 2014, a moderate earthquake with magnitude 5.0, which generated high-quality Rayleigh wave data, occurred in the Northwest Territories, Canada, ∼100 km to the east of the Cordilleran Deformation Front. We carefully selected 23 seismic stations that recorded the Rayleigh waves and divided them into 13 groups according to the azimuth angle between the earthquake and the stations; these groups mostly sample the Cordillera. In each group, we measured Rayleigh wave group velocity dispersion, which we inverted for one-dimensional shear-wave velocity models of the crust. We thus obtained 13 models that consistently show low seismic velocities with respect to reference models, with a slow upper and lower crust surrounding a relatively fast mid crustal layer. The average of the 13 models is consistent with receiver function data in the central portion of the Cordillera. Finally, we compared earthquake locations determined by the Geological Survey of Canada using a simple homogenous crust over a mantle half space with those estimated using the new crustal velocity model, and show that estimates can differ by as much as 10 km.


2021 ◽  
Author(s):  
Gregor Rajh ◽  
Josip Stipčević ◽  
Mladen Živčić ◽  
Marijan Herak ◽  
Andrej Gosar

&lt;p&gt;The investigated area of the NW Dinarides is bordered by the Adriatic foreland, the Southern Alps, and the Pannonian basin at the NE corner of the Adriatic Sea. Its complex crustal structure is the result of interactions among different tectonic units. Despite numerous seismic studies taking place in this region, there still exists a need for a detailed, smaller scale study focusing mainly on the brittle part of the Earth's crust. Therefore, we decided to investigate the velocity structure of the crust using concepts of local earthquake tomography (LET) and minimum 1-D velocity model. Here, we present the results of the 1-D velocity modeling and the catalogue of the relocated seismicity. A minimum 1-D velocity model is computed by simultaneous inversion for hypocentral and velocity parameters together with seismic station corrections and represents the best fit to the observed arrival times.&lt;/p&gt;&lt;p&gt;We used 15,579 routinely picked P wave arrival times from 631 well-located earthquakes that occurred in Slovenia and in its immediate surroundings (mainly NW Croatia). Various initial 1-D velocity models, differing in velocity and layering, were used as input for velocity inversion in the VELEST program. We also varied several inversion parameters during the inversion runs. Most of the computed 1-D velocity models converged to a stable solution in the depth range between 0 and 25 km. We evaluated the inversion results using rigorous testing procedures and selected two best performing velocity models. Each of these models will be used independently as the initial model in the simultaneous hypocenter-velocity inversion for a 3-D velocity structure in LET. Based on the results of the 1-D velocity modeling, seismicity distribution, and tectonics, we divided the study area into three parts, redefined the earthquake-station geometry, and performed the inversion for each part separately. This way, we gained a better insight into the shallow velocity structure of each subregion and were able to demonstrate the differences among them.&lt;/p&gt;&lt;p&gt;Besides general structural implications and a potential to improve the results of LET, the new 1-D velocity models along with station corrections can also be used in fast routine earthquake location and to detect systematic travel time errors in seismological bulletins, as already shown by some studies using similar methods.&lt;/p&gt;


2021 ◽  
Author(s):  
Yueyang Xia ◽  
Dirk Klaeschen ◽  
Heidrun Kopp ◽  
Michael Schnabel

Abstract. Accurate subsurface velocity models are crucial for geological interpretations based on seismic depth images. Seismic reflection tomography is an effective iterative method to update and refine a preliminary velocity model for depth imaging. Based on residual move-out analysis of reflectors in common image point gathers an update of the velocity is estimated by a ray-based tomography. To stabilize the tomography, several preconditioning strategies exist. Most critical is the estimation of the depth error to account for the residual move-out of the reflector in the common image point gathers. Because the depth errors for many closely spaced image gathers must be picked, manual picking is extremely time-consuming, human biased, and not reproducible. Data-driven picking algorithms based on coherence or semblance analysis are widely used for hyperbolic or linear events. However, for complex-shaped depth events, pure data-driven picking is difficult. To overcome this, the warping method named Non-Rigid Matching is used to estimate a depth error displacement field. Warping is used, e.g., to merge photographic images or to match two seismic images from time-lapse data. By calculating the displacements between an offset to its neighbouring offset in the common image point domain, a locally smooth-shaped displacement field is defined for each data sample. Depending on the complexity of the subsurface, sample tracking through the displacement field along predefined horizons or on a simple regular grid yields discrete depth error values for the tomography. The application to a multi-channel seismic line across the Sunda subduction zone offshore Lombok island, Indonesia, illustrates the approach and documents the advantages of the method to estimate a detailed velocity structure in a complex tectonic regime. By incorporating the warping scheme into the reflection tomography, we demonstrate an increase in the velocity resolution and precision by improving the data-driven accuracy of depth error picks with arbitrary shapes. This approach offers the possibility to use the full capacities of tomography and further leads to more accurate interpretations of complex geological structures.


2020 ◽  
Author(s):  
Jonas Preine ◽  
Benjamin Schwarz ◽  
Alexander Bauer ◽  
Dirk Gajewski ◽  
Christian Hübscher

&lt;p&gt;The active seismic method is a standard tool for studying the Earth&amp;#8217;s lithosphere. On scales from centimetres to kilometres, academic research is generally interested in highly complex geological targets such as volcanic edifices, crustal faults or salt environments. In order to properly image these structures, large and expensive multichannel acquisitions with a high offset-to-target depth ratio are required. In practice, however, these are often hardly affordable for academic institutions, with the result that reflections often only poorly illuminate laterally variable structures, which in turn compromises imaging and interpretation. As in common practice, most of the processing and interpretational steps are tailored to the reflected wavefield, faint diffracted contributions are typically considered as an unwanted by-product.&lt;/p&gt;&lt;p&gt;In recent works, however, it has been shown that diffractions possess unique properties which bear the potential to overcome the aforementioned limitations. Wave diffraction occurs at geodynamically important features like faults, pinch-outs, erosional surfaces or other small-scale scattering objects and encodes sub-wavelength information on the scattering geometry. Since diffracted waves do not obey Snell&amp;#8217;s Law, they provide superior illumination compared to reflected waves. Moreover, due to their passive-source like radiation, they encode their full multichannel response in prominent data subsets like the zero-offset section. In order to explore what can be learned from the faint diffracted wavefield, we use academic seismic data from the Santorini-Amorgos Tectonic Zone (SATZ) in the Southern Aegean Sea. This is an area well known for its local complexity, indicated by the occurrence of extended fault systems and volcanic edifices as well as a complex acoustic basement. As the available seismic data in this region were acquired using a relatively short streamer, the SATZ represents a classical example for the need of innovative methods for seismic processing and interpretation.&lt;/p&gt;&lt;p&gt;By means of a robust and computationally efficient scheme for the extraction of diffractions that models and adaptively subtracts the reflected wavefield from the data, we reveal a rich diffracted wavefield from zero-offset data. On the one hand, we use the diffraction-only sections for analysing the small-scale structural complexity and demonstrate that the geological interpretation can benefit from these observations. On the other hand, we use the diffractions to estimate insightful wavefront attributes in the zero-offset domain. Based on these attributes, we perform wavefront tomography to obtain depth-velocity models. Compared to depth-velocity models derived from the reflected contributions, the diffraction-based velocity model fits the data significantly better. After refining this velocity model, we perform prestack depth migration and obtain highly valuable depth converted seismic sections. Concluding our results, we strongly encourage the incorporation of diffractions in standard processing and interpretational schemes.&lt;/p&gt;


Author(s):  
Г.Д. Етирмишли ◽  
С.Э. Казымова

При изучении скоростной структуры земной коры Азербайджана по сейсмологическим данным ис- пользовались в основном два подхода. Первый состоит в уточнении модели среды на основании наблюда- емых отклонений времен пробега волн от землетрясений относительно стандартного годографа. Второй основан на использовании разности времен пробега от источников до станции для групп близко располо- женных событий. Одномерные скоростные модели Р и S-волн были построены на основе одномерных моделей, пред- ложенных в работе Гасанова А.Г. Построение модели, расчет станционных поправок и перелокация со- бытий производились в программе Velest. Исследуемый объем до глубины 60 км был разбит на мелкие слои толщиной 2 км в интервале глубин 010 км и толщиной 510 км в интервале глубин 1060 км. В ходе исследования рассматривались сейсмологические данные о параметрах локальных землетрясений и вре- менах прихода P и S-волн зарегистрированных сетью телеметрических станций за период 20042011 гг. Анализировались данные 28-ми сейсмических станций Азербайджана, охватывающие всю исследуемую территорию. Для расчета трехмерного скоростного поля использовалась программа TomotetraFD. В этой програм- ме реализован классический сейсмотомографический метод для случая, когда источники и приемники находятся внутри исследуемого региона. Two approaches were used for investigation of crust velocity structure of Azerbaijan on the basis of seismological data. The first one consists in medium model adjustment on the basis of observed deviation of travel times of waves from earthquakes relative to standard hodograph. The second is based on difference in travel times from source to station for a group of close located stations. One dimensional velocity models of P- and S-waves were constructed on the basis of one dimensional models proposed by A.G.Gasanov. Model design, calculation of stations corrections and relocationing of events 74 Геология и геофизика Юга России, 1, 2012 were performed in Velest program. Investigating volume to depth of 60 km was divided in small layers of 2 km thickness in 0-10 km interval and 5-10 km in 10-60 km interval. Seismological data about parameters of local earthquakes registered by network of telemetric stations in 2004-2011 and arrival times of P- and S-waves were used. Data of 28 seismic stations of Azerbaijan covering all the investigating territory were analyzed. Three dimensional velocity field was calculated by means of TomotetraFD program. Classical seismotomographical method for the case when sources and receivers are located within investigating region is realized in the program.


Author(s):  
Mungunsuren Dashdondog ◽  
Odonbaatar Chimed ◽  
Anne Meltzer ◽  
Nomin-Erdene Erdenetsogt ◽  
Josh Stachnik

One dimensional (1D) velocity models are still widely used for computing earthquake locations at seismological centers. The location accuracy of an earthquake strongly depends on the velocity model used to compute the location. In the past, the local velocity model developed for the Hangay region was lacking precision due to insufficient data. Within the framework of the “Intracontinental Deformation and Surface Uplift- Geodynamic Evolution of the Hangay Dome, Mongolia, Central Asia” project [15], 72 seismic Broadband stations network were deployed in the Hangay Dome. This gives us an opportunity to estimate the crustal velocity structure of the South Hangay region using recorded local earthquake data. For this purpose, available velocity models for the South Hangay region have been re-evaluated.  By simultaneous invertion P- and S-wave arrival times using VELEST algorithm, we estimated minimum 1D velocity models, station corrections, hypocentre locations, and origin times for the south Hangay region. Consequently, 1D crustal velocity model is proposed for the South Hangay region. This new model is expected to improve the accuracy of the routine hypocenter determination and as initial reference models for seismic tomography study.


Sign in / Sign up

Export Citation Format

Share Document