Layered lithospheric anisotropy beneath southeastern Tibet using harmonic decomposition of receiver functions

2021 ◽  
Author(s):  
Ashwani Kant Tiwari ◽  
Arun Singh ◽  
Dipankar Saikia ◽  
Chandrani Singh
Keyword(s):  
Lower Crust ◽  
Receiver Function ◽  
Seismic Station ◽  
Research Work ◽  
Receiver Functions ◽  
Ductile Deformation ◽  
Lithospheric Deformation ◽  
Indian Lithosphere ◽  
Harmonic Decomposition ◽  
Southeastern Tibet

<p>The present research work interrogates the depth-dependent lithospheric dipping and anisotropic fabrics that characterize major fault and suture zone rheology, essential to understanding the lithospheric deformation and geodynamic process beneath southeastern Tibet. The depth-dependent anisotropic trend has been investigated via harmonic stripping of receiver functions (RFs) at 70 stations of the Eastern Syntaxis experiment, operated between 2003-2004. First, 3683 good quality P-RFs are computed from 174 teleseismic events. All the events are of magnitude ≥5.5 and recorded in the epicentral distribution of 30° to 90°. After that, the harmonic stripping technique is performed at each seismic station to retrieve the first (k = 1) and second (k =2) degree harmonics from the receiver function dataset. Our study also characterizes the type (fast or slow) of the symmetric axis. The upper crustal (0-20 km) anisotropic orientations are orthogonal to the major faults and suture zones of the area and suggest the structure-induced anisotropy. However, the anisotropic orientations in the mid-to-lower crust and uppermost mantle orientations suggest the ductile deformation due to material flow towards the east. Comparison from depth-dependent lithospheric trend and fast polarization directions obtained from the core-refracted and direct-S phases suggest the decoupled crust and lithospheric mantle beneath the area.  The distinct anisotropic trends in the Namche Barwa Metamorphic Massif (NBMM) indicate the northward indentation of the Indian crust beneath the Lhasa block. However, the lower crust and uppermost anisotropic orientation suggest the fragmented Indian lithosphere beneath the area. Our results add new constraints in understanding the type of strain and its causes in the region.</p>

scholarly journals Lithospheric and sub-lithospheric deformation under the Borborema Province of NE Brazil from receiver function harmonic stripping

2019 ◽  
Author(s):  
Gaelle Lamarque ◽  
Jordi Julià
Keyword(s):  
Lithospheric Mantle ◽  
Seismic Anisotropy ◽  
Receiver Function ◽  
Seismic Station ◽  
Atlantic Coast ◽  
Shear Zones ◽  
Receiver Functions ◽  
Lithospheric Deformation ◽  
Borborema Province ◽  
Continental Interior

Abstract. The depth-dependent anisotropic structure of the lithosphere under the Borborema Province of northeast Brazil has been investigated through harmonic stripping of receiver functions developed at 39 stations in the region. This method retrieves the first (k = 1) and second (k = 2) degree harmonics of a receiver function dataset, which characterize seismic anisotropy beneath a seismic station. Anisotropic fabrics are in turn directly related to the deformation of the lithosphere from past and current tectonic processes. Our results reveal the presence of anisotropy within the crust and the lithospheric mantle throughout the entire Province, with the exception of a few stations in the continental interior that lack evidence for any anisotropic signatures. Most stations in the continental interior report consistent anisotropic orientations in the crust and lithospheric mantle, suggesting a dominant NE-SW pervasive deformation along lithospheric-scale shear zones developed during the Brasiliano-Pan African orogeny. The lack of anisotropy at a few stations along a NE-SW trend in the center on the Province is harder to explain, but might be related to heating of the lithosphere by an asthenospheric channel. Finally, several stations along the Atlantic coast reveal depth-dependent anisotropic orientations roughly (sub)perpendicular to the margin. These results suggest a more recent overprint, probably related to the presence of frozen anisotropy in the lithosphere due to stretching and rifting during the opening of the South Atlantic.

scholarly journals Lithospheric and sublithospheric deformation under the Borborema Province of northeastern Brazil from receiver function harmonic stripping

Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 893-905 ◽  
Author(s):  
Gaelle Lamarque ◽  
Jordi Julià
Keyword(s):  
Lithospheric Mantle ◽  
Seismic Anisotropy ◽  
Receiver Function ◽  
Seismic Station ◽  
Atlantic Coast ◽  
Shear Zones ◽  
Receiver Functions ◽  
Azimuthal Anisotropy ◽  
Northeastern Brazil ◽  
Borborema Province

Abstract. The depth-dependent anisotropic structure of the lithosphere under the Borborema Province in northeast Brazil has been investigated via harmonic stripping of receiver functions developed at 39 stations in the region. This method retrieves the first (k=1) and second (k=2) degree harmonics of a receiver function dataset, which characterize seismic anisotropy beneath a seismic station. Anisotropic fabrics are in turn directly related to the deformation of the lithosphere from past and current tectonic processes. Our results reveal the presence of anisotropy within the crust and the lithospheric mantle throughout the entire province. Most stations in the continental interior report consistent anisotropic orientations in the crust and lithospheric mantle, suggesting a dominant northeast–southwest pervasive deformation along lithospheric-scale shear zones developed during the Brasiliano–Pan-African orogeny. Several stations aligned along a northeast–southwest trend located above the (now aborted) Mesozoic Cariri–Potiguar rift display large uncertainties for the fast-axis direction. This non-azimuthal anisotropy may be related to a complex anisotropic fabric resulting from a combination of deformation along the ancient collision between Precambrian blocks, Mesozoic extension and thermomechanical erosion dragging by sublithospheric flow. Finally, several stations along the Atlantic coast reveal depth-dependent anisotropic orientations roughly (sub)perpendicular to the margin. These results suggest a more recent overprint, probably related to the presence of frozen anisotropy in the lithosphere due to stretching and rifting during the opening of the South Atlantic.

Crustal anisotropy beneath northeastern Tibetan Plateau from the harmonic decomposition of receiver functions

2019 ◽  
Vol 220 (3) ◽  
pp. 1585-1603
Author(s):  
Zhenxin Xie ◽  
Vadim Levin ◽  
Qingju Wu
Keyword(s):  
Tibetan Plateau ◽  
Seismic Anisotropy ◽  
Receiver Function ◽  
Flow Direction ◽  
Receiver Functions ◽  
Low Velocity ◽  
Northeastern Tibetan Plateau ◽  
Harmonic Decomposition ◽  
Anisotropic Layers ◽  
Qilian Orogen

SUMMARY A uniformly spaced linear transect through the northeastern Tibetan Plateau was constructed using 54 stations from ChinaArray Phase II. We used a set of colocated earthquakes to form receiver function beams that were then used to construct a 2-D image of main converting boundaries in our region and to investigate lateral changes in main impedance contrasts along the transect. The image revealed obvious mid-crustal low-velocity zones beneath the Qilian Orogen and the Alxa Block. We developed a new procedure that uses harmonically decomposed receiver functions to characterize seismic anisotropy, and that can determine both the orientations of symmetry axes and their type (fast or slow). We tested our technique on a number of synthetic models, and subsequently applied it to the data from the transect. We found that: (1) within the upper crust the orientations of slow symmetry axes are nearly orthogonal to the strike directions of faults, and thus anisotropy is likely caused by the shape preferred orientation of fluid-saturated cracks or fractures and (2) together with the low-velocity zones revealed from receiver functions stacks, anisotropic layers in the middle-to-lower crust could be explained by the crustal channel flow that was proposed for this region by previous studies. The shear within the boundary layers of crustal flow forms anisotropy with symmetry axes parallel to the flow direction.

Crustal Anisotropy beneath Northern Iran calculated by harmonic decomposition of Receiver Functions.

2020 ◽  
Author(s):  
Mohsen Azqandi ◽  
Mohammad Reza Abbassi ◽  
Meysam Mahmoodabadi ◽  
Ahmad Sadidkhouy
Keyword(s):  
Receiver Function ◽  
Receiver Functions ◽  
Upper Crust ◽  
Northern Iran ◽  
Crustal Anisotropy ◽  
Principal Stress Direction ◽  
South Central ◽  
Stress Direction ◽  
Harmonic Decomposition ◽  
Time Variations

<p>This study concerns crustal anisotropy at 16 permanent seismic stations to investigate preferentially aligned cracks or structures and their relation to the stress-state in the South Central Alborz (northern Iran). We consider plunging anisotropy and dipping interfaces of multiple layers using harmonic functions to correct the arrival time variations of <em>Ps</em> phases from different back-azimuths.</p><p>The dominant fast orientation of integrated crustal anisotropy strikes NE, almost parallel to the stress direction in the upper crust. The magnitude of crustal anisotropy is found to be in range of 0.1 s to 0.5 s. In some stations, intracrustal interface is observed, for which we analyzed harmonic decomposition of receiver functions to consider anisotropy in the upper crust. Upper crustal anisotropy strikes NE, close to the principal stress direction, indicating that stress in the upper crust plays a major role in producing anisotropy and deformation. In a few stations, crustal anisotropy display different directions rather than NE, which maybe controlled by cracks and fractures of dominant faults.</p><p>Keywords: Anisotropy, Receiver function, harmonic decomposition, Northern Iran.</p>

Crust-Mantle Kinematics in and around the Hellenic Arc Elucidated by Local Shear Wave Splitting and Receiver Function Analyses

2021 ◽  
Author(s):  
Derya Keleş ◽  
Tuna Eken ◽  
Judith M. Confal ◽  
Tuncay Taymaz
Keyword(s):  
Seismic Anisotropy ◽  
Receiver Function ◽  
Receiver Functions ◽  
Data Sets ◽  
S Wave ◽  
Uppermost Mantle ◽  
Hellenic Arc ◽  
Wave Splitting ◽  
Harmonic Decomposition ◽  
Local Shear

<p>The fundamental knowledge on seismic anisotropy inferred from various data sets can enhance our understanding of its vertical resolution that is critical for a better interpretation of past and current dynamics and resultant crustal and mantle kinematics in the Hellenic Trench and its hinterland. To investigate the nature of deformation zones, we perform both local S-wave splitting (SWS) measurements and receiver functions (RFs) analysis. Our preliminary findings from the harmonic decomposition technique performed on radial and tangential RFs suggest relatively more substantial anisotropic signals in the lower crust and uppermost mantle with respect to upper and middle crustal structure in the region. Apparent anisotropic orientations obtained from RFs harmonic decomposition process show several consistencies with those discovered from local SWS measurements at selected stations. The actual anisotropic orientation for the structures, however, requires further modelling of the receiver functions obtained.</p>

scholarly journals Crustal and upper mantle seismic structure of the Svalbard Archipelago from the receiver function analysis

2015 ◽  
Vol 36 (2) ◽  
pp. 89-107 ◽  
Author(s):  
Monika Wilde−Piórko
Keyword(s):  
Upper Mantle ◽  
Barents Sea ◽  
Receiver Function ◽  
Function Analysis ◽  
Seismic Station ◽  
Receiver Functions ◽  
Seismic Velocities ◽  
Seismic Structure ◽  
Svalbard Archipelago ◽  
The North

Abstract Receiver function provides the signature of sharp seismic discontinuities and the information about the shear wave (S−wave) velocity distribution beneath the seismic station. This information is very valuable in areas where any or few reflection and/or refraction studies are available and global and/or regional models give only rough information about the seismic velocities. The data recorded by broadband seismic stations have been analysed to investigate the crustal and upper mantle structure of the Svalbard Archipelago. Svalbard Archipelago is a group of islands located in Arctic, at the north−western part of the Barents Sea continental platform, which is bordered to the west and to the north by passive continental margins. The new procedure of parameterization and selection of receiver functions (RFs) has been proposed. The back−azimuthal sections of RF show a strong variation for the HSPB and KBS stations. Significant amplitudes of transversal component of RF (T−RF) for the HSPB station indicate a shallow dipping layer towards the southwest. The structure of the crust beneath the SPITS array seems to be less heterogeneous, with very low amplitudes of converted phase comparing to the KBS and HSPB stations. Forward modelling by trial−and−error method shows a division of the crust into 3-4 layers beneath all stations and layering of the uppermost mantle beneath the SPITS array and the HSPB stations. The thickness of the mantle transition zone is larger for western part of archipelago and smaller for eastern part comparing to iasp91 model.

scholarly journals Reliable workflow for inversion of seismic receiver function and surface wave dispersion data: a “13 BB Star” case study

2019 ◽  
Vol 24 (1) ◽  
pp. 101-120
Author(s):  
Kajetan Chrapkiewicz ◽  
Monika Wilde-Piórko ◽  
Marcin Polkowski ◽  
Marek Grad
Keyword(s):  
Surface Wave ◽  
Inverse Problems ◽  
Receiver Function ◽  
Joint Inversion ◽  
Wave Dispersion ◽  
Receiver Functions ◽  
P Wave ◽  
Surface Wave Dispersion ◽  
Phase Velocity Dispersion ◽  
Wide Range

AbstractNon-linear inverse problems arising in seismology are usually addressed either by linearization or by Monte Carlo methods. Neither approach is flawless. The former needs an accurate starting model; the latter is computationally intensive. Both require careful tuning of inversion parameters. An additional challenge is posed by joint inversion of data of different sensitivities and noise levels such as receiver functions and surface wave dispersion curves. We propose a generic workflow that combines advantages of both methods by endowing the linearized approach with an ensemble of homogeneous starting models. It successfully addresses several fundamental issues inherent in a wide range of inverse problems, such as trapping by local minima, exploitation of a priori knowledge, choice of a model depth, proper weighting of data sets characterized by different uncertainties, and credibility of final models. Some of them are tackled with the aid of novel 1D checkerboard tests—an intuitive and feasible addition to the resolution matrix. We applied our workflow to study the south-western margin of the East European Craton. Rayleigh wave phase velocity dispersion and P-wave receiver function data were gathered in the passive seismic experiment “13 BB Star” (2013–2016) in the area of the crust recognized by previous borehole and refraction surveys. Final models of S-wave velocity down to 300 km depth beneath the array are characterized by proximity in the parameter space and very good data fit. The maximum value in the mantle is higher by 0.1–0.2 km/s than reported for other cratons.

scholarly journals The lithosphere-asthenosphere boundary observed with USArray receiver functions

2012 ◽  
Vol 4 (1) ◽  
pp. 1-31 ◽  
Author(s):  
P. Kumar ◽  
X. Yuan ◽  
R. Kind ◽  
J. Mechie
Keyword(s):  
Surface Wave ◽  
Receiver Function ◽  
The United States ◽  
Receiver Functions ◽  
Function Technique ◽  
Surface Wave Tomography ◽  
Entire Area ◽  
S Receiver Function ◽  
The Usa ◽  
Wave Tomography

Abstract. The dense deployment of seismic stations so far in the western half of the United States within the USArray project provides the opportunity to study in greater detail the structure of the lithosphere-asthenosphere system. We use the S receiver function technique for this purpose which has higher resolution than surface wave tomography, is sensitive to seismic discontinuities and has no problems with multiples like P receiver functions. Only two major discontinuities are observed in the entire area down to about 300 km depth. These are the crust-mantle boundary (Moho) and a negative boundary which we correlate with the lithosphere-asthenosphere boundary (LAB) since a low velocity zone is the classical definition of the seismic observation of the asthenosphere by Gutenberg (1926). Our S receiver function LAB is at a depth of 70–80 km in large parts of westernmost North America. East of the Rocky Mountains its depth is generally between 90 and 110 km. Regions with LAB depths down to about 140 km occur in a stretch from northern Texas over the Colorado Plateau to the Columbia Basalts. These observations agree well with tomography results in the westernmost USA and at the east coast. However, in the central cratonic part of the USA the tomography LAB is near 200 km depth. At this depth no discontinuity is seen in the S receiver functions. The negative signal near 100 km depth in the central part of the USA is interpreted by Yuan and Romanowicz (2010) or Lekic and Romanowicz (2011) as a recently discovered mid lithospheric discontinuity (MLD). A solution for the discrepancy between receiver function imaging and surface wave tomography is not yet obvious and requires more high resolution studies at other cratons before a general solution may be found. Our results agree well with petrophysical models of increased water content in the asthenosphere, which predict a sharp and shallow LAB also in continents (Mierdel et al., 2007).

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>

Sign in / Sign up

Export Citation Format

Share Document