scholarly journals Source depopulation potential and surface-wave tomography using a crosscorrelation method in a scattering medium

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. SA51-SA61 ◽  
Author(s):  
Pierre Gouédard ◽  
Philippe Roux ◽  
Michel Campillo ◽  
Arie Verdel ◽  
Huajian Yao ◽  
...  

We use seismic prospecting data on a 40 × 40 regular grid of sources and receivers deployed on a 1 km × 1 km area to assess the feasibility and advantages of velocity analysis of the shallow subsurface by means of surface-wave tomography with Green’s functions estimated from crosscorrelation. In a first application we measure Rayleigh-wave dispersion curves in a 1D equivalent medium. The assumption that the medium is laterally homogeneous allows using a simple projection scheme and averaging of crosscorrelation functions over the whole network. Because averaging suppresses noise, this method yields better signal-to-noise ratio than traditional active-source approaches, and the improvement can be estimated a priori from acquisition parameters. We find that high-quality dispersion curves can be obtained even when we reduce the number of active sources used as input for the correlations. Such source depopulation can achieve significant reduction in the cost of active source acquisition. In a second application we compare Rayleigh-wave group velocity tomography from raw and reconstructed data. We can demonstrate that the crosscorrelation approach yields group velocity maps that are similar to active source maps. Scattering has an importance here as it may enhance the crosscorrelation performance. We quantify the scattering properties of the medium using mean free path measurements from coherent and incoherent parts of the signal. We conclude that for first-order velocity analysis of the shallow subsurface, the use of crosscorrelation offers a cost-effective alternative to methods that rely exclusively on active sources.

2014 ◽  
Vol 6 (1) ◽  
pp. 1-34 ◽  
Author(s):  
M. Abdetedal ◽  
Z. H. Shomali ◽  
M. R. Gheitanchi

Abstract. We applied seismic ambient noise surface wave tomography to estimate Rayleigh wave empirical Green's functions from cross-correlations to study crust and uppermost mantle structure beneath the Makran region in south-east Iran. We analysed 12 months of continuous data from January 2009 through January 2010 recorded at broadband seismic stations. We obtained group velocity of the fundamental mode Rayleigh-wave dispersion curves from empirical Green's functions between 10 and 50 s periods by multiple-filter analysis and inverted for Rayleigh wave group velocity maps. The final results demonstrate significant agreement with known geological and tectonic features. Our tomography maps display low-velocity anomaly with south-western north-eastern trend, comparable with volcanic arc settings of the Makran region, which may be attributable to the geometry of Arabian Plate subducting overriding lithosphere of the Lut block. At short periods (<20 s) there is a pattern of low to high velocity anomaly in northern Makran beneath the Sistan Suture Zone. These results are evidence that surface wave tomography based on cross correlations of long time-series of ambient noise yields higher resolution group speed maps in those area with low level of seismicity or those region with few documented large or moderate earthquake, compare to surface wave tomography based on traditional earthquake-based measurements.


2020 ◽  
Vol 110 (3) ◽  
pp. 1359-1371
Author(s):  
Lun Li ◽  
Yuanyuan V. Fu

ABSTRACT An understanding of mantle dynamics occurring beneath the Tibetan plateau requires a detailed image of its seismic velocity and anisotropic structure. Surface waves at long periods (&gt;50  s) could provide such critical information. Though Rayleigh-wave phase velocity maps have been constructed in the Tibetan regions using ambient-noise tomography (ANT) and regional earthquake surface-wave tomography, Love-wave phase velocity maps, especially those at longer periods (&gt;50  s), are rare. In this study, two-plane-wave teleseismic surface-wave tomography is applied to develop 2D Rayleigh-wave and Love-wave phase velocity maps at periods between 20 and 143 s across eastern and central Tibet and its surroundings using four temporary broadband seismic experiments. These phase velocity maps share similar patterns and show high consistency with those previously obtained from ANT at overlapping periods (20–50 s), whereas our phase velocity maps carry useful information at longer periods (50–143 s). Prominent slow velocity is imaged at periods of 20–143 s beneath the interior of the Tibetan plateau (i.e., the Songpan–Ganzi terrane, the Qiangtang terrane, and the Lhasa terrane), implying the existence of thick Tibetan crust along with warm and weak Tibetan lithosphere. In contrast, the dispersal of fast velocity anomalies coincides with mechanically strong, cold tectonic blocks, such as the Sichuan basin and the Qaidam basin. These phase velocity maps could be used to construct 3D shear-wave velocity and radial seismic anisotropy models of the crust and upper mantle down to 250 km across the eastern and central Tibetan plateau.


2021 ◽  
Author(s):  
Maryam Rezaei ◽  
Taghi Shirzad

&lt;p&gt;Classical surface wave tomography based on waveform scattering through seismic data has an important role in studying the structure of the Earth&amp;#8216;s crust and upper mantle on regional and global scale. The shallow crustal velocity structure is studied using earthquake waveforms in Khorasan/E-NE Iran. For this purpose, 522 local recorded waveforms with M&amp;#8805;4, which occurred between 53&amp;#176;-63&amp;#176;E and 30&amp;#176;-42&amp;#176;N, were selected. Therefore, all available vertical components of waveforms recorded at the stations in the Iranian Seismological Center (IrSC), the International Institute of Earthquake Engineering and Seismology (IIEES), and the IRIS global network collected in the period between January 2006 to October 2020. Then, some data selection criteria were applied for each waveform, including (i)SNR&gt;4, (ii) the gap time less than 2 s within the expected signal window (1.5-4.5 km/s), (iii) epicentral distance greater than 20 km. The multiple-filter analysis technique was then applied by the computer program in seismology Hermann and Ammon (2013) to measure Rayleigh wave dispersion curves in the period range of 3-50 s. Finally, Rayleigh wave 2D horizontal group velocity maps are calculated by the fast marching surface wave tomography method. Our tomographic results indicate some local low velocity anomalies appeared in the W-NW of the study area where it connected with the East-Alborz tectonic structure. Also, the Doruneh Fault System clearly separates Kopeh-Dagh tectonic zone and Central Iran micro-plateau. However, a high velocity anomaly appears in Kopeh-Dagh tectonic zone at periods larger than 16 s and 30 s.&lt;/p&gt;


2020 ◽  
Author(s):  
Petr Kolínský ◽  
Tena Belinić ◽  
Josip Stipčević ◽  
Irene Bianchi ◽  
Florian Fuchs ◽  
...  

&lt;p&gt;The Alpine-Dinarides are a complex orogenic system, with its tectonic evolution controlled by the ongoing convergence between Eurasian and African plates with the Adriatic microplate wedged between them. Our study focuses on the upper mantle of the wider Alpine-Dinarides region, and we present surface-wave tomography of two overlapping subregions, interpreting the seismic velocity features in the context of regional geodynamics.&lt;/p&gt;&lt;p&gt;In the first part, we use records of 151 teleseismic earthquakes (2010-2018) at 98 stations distributed across the wider Dinarides region. Surface-wave phase velocities are measured in the range of 30 &amp;#8211; 160 s by the two-station method at pairs of stations aligned along the great circle paths with the epicenters. We apply several data-quality tests before the dispersion curves are measured. We use Rayleigh waves recorded on both radial and vertical components. Only the dispersions measured coherently at both components are used for the tomography. In total, we reach the number of 9000 phase velocity measurements for the period of 50 s. Tomographic results including resolution estimates are provided for various frequencies; the local dispersion curves are inverted for depths from the surface down to 300 km. Results are shown as maps for various depths and as cross-sections along several profiles of shear-wave velocities in the whole region.&lt;/p&gt;&lt;p&gt;The other study focuses on the Alps. The AlpArray seismic network stretches hundreds of kilometers in width and more than thousand kilometers in length. It is distributed over the greater Alpine region (Europe) and consists of around 250 temporary and around 400 permanent broadband stations with interstation distances around 40 km. The earthquakes are selected between years 2016-2019. The methodology differs from the Dinarides case in a sense, that while before we used many earthquakes and less stations pairs (due to sparser station coverage), for the Alps, we use less earthquakes (32) and many more stations pairs (tens of thousands) making use of the dense station coverage of the AlpArray network.&lt;/p&gt;&lt;p&gt;Results of the depth inversion of the local dispersion measurements for the Alps are compared with local surface-wave phase-velocity measurement obtained from the (sub)array approach.&lt;/p&gt;


Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. EN13-EN26
Author(s):  
Ilaria Barone ◽  
Emanuel Kästle ◽  
Claudio Strobbia ◽  
Giorgio Cassiani

Surface wave tomography (SWT) is a powerful and well-established technique to retrieve 3D shear-wave (S-wave) velocity models at the regional scale from earthquakes and seismic noise measurements. We have applied SWT to 3D active-source data, in which higher modes and heterogeneous spatial sampling make phase extraction challenging. First, synthetic traveltimes calculated on a dense, regular-spaced station array are used to test the performance of three different tomography algorithms (linearized inversion, Markov chain Monte Carlo [MCMC], and eikonal tomography). The tests suggest that the lowest misfit to the input model is achieved with the MCMC algorithm, at the cost of a much longer computational time. Then, real phases were extracted from a 3D exploration data set at different frequencies. This operation included an automated procedure to isolate the fundamental mode from higher order modes, phase unwrapping in two dimensions, and the estimation of the zero-offset phase. These phases are used to compute traveltimes between each source-receiver couple, which are input into the previously tested tomography algorithms. The resulting phase-velocity maps show good correspondence, highlighting the same geologic structures for all three methods. Finally, individual dispersion curves obtained by the superposition of phase-velocity maps at different frequencies are depth inverted to retrieve a 3D S-wave velocity model.


2020 ◽  
Vol 224 (3) ◽  
pp. 1724-1741
Author(s):  
Jeremy M Gosselin ◽  
Pascal Audet ◽  
Andrew J Schaeffer ◽  
Fiona A Darbyshire ◽  
Clément Estève

SUMMARY Surface wave tomography is a valuable tool for constraining azimuthal anisotropy at regional scales. However, sparse and uneven coverage of dispersion measurements make meaningful uncertainty estimation challenging, especially when applying subjective model regularization. This paper considers azimuthal anisotropy constrained by measurements of surface wave dispersion data within a Bayesian trans-dimensional (trans-d) tomographic inversion. A recently proposed alternative model parametrization for trans-d inversion is implemented in order to produce more realistic models than previous studies considering trans-d surface wave tomography. The reversible-jump Markov chain Monte Carlo sampling technique is used to numerically estimate the posterior probability density of the model parameters. Isotropic and azimuthally anisotropic components of surface wave group velocity maps (and their associated uncertainties) are estimated while avoiding model regularization and allowing model complexity to be determined by the data information content. Furthermore, data errors are treated as unknown, and solved for within the inversion. The inversion method is applied to measurements of surface wave dispersion from regional earthquakes recorded over northern Cascadia and Haida Gwaii, a region of complex active tectonics but highly heterogeneous station coverage. Results for isotropic group velocity are consistent with previous studies that considered the southern part of the study region over Cascadia. Azimuthal anisotropic fast-axis directions are generally margin-parallel between Vancouver Island and Haida Gwaii, with a small change in direction and magnitude along the margin which may be attributed to the changing tectonic regime (from subduction to transform tectonics). Estimated errors on the dispersion data (solved for within the inversion) reveal a correlation between surface wave period and the dependence of data errors on travel path length. This paper demonstrates the value of considering azimuthal anisotropy within Bayesian tomographic inversions. Furthermore, this work provides structural context for future studies of tectonic structure and dynamics of northern Cascadia and Haida Gwaii, with the aim of improving our understanding of seismic and tsunami hazards.


2020 ◽  
Vol 224 (2) ◽  
pp. 1287-1300
Author(s):  
Małgorzata Chmiel ◽  
Philippe Roux ◽  
Marc Wathelet ◽  
Thomas Bardainne

SUMMARY We propose a new surface wave tomography approach that benefits from densely sampled active-source arrays and brings together elements from active-source seismic-wave interferometry, full waveform inversion and dense-array processing. In analogy with optical interferometry, seismic Michelson interferometer (SMI) uses seismic interference patterns given by the data-based diffraction kernels in an iterative inversion scheme to image a medium. SMI requires no traveltime measurements and no spatial regularization, and it accounts for bent rays. Furthermore, the method does not need computation of complex synthetic models, as it works as a data-driven inversion technique that makes it computationally very fast. In an automatic way, it provides high-resolution phase-velocity maps and their error estimation. SMI can complete traditional surface wave tomography studies, as its use can be easily extended from land active seismic data to the virtual source gathers of ambient-noise-based studies with dense arrays.


Sign in / Sign up

Export Citation Format

Share Document