Imaging the Tectonic Grain of the Northern Cordillera Orogen Using Transportable Array Receiver Functions

2020 ◽  
Vol 91 (6) ◽  
pp. 3086-3105 ◽  
Author(s):  
Vera Schulte-Pelkum ◽  
Jonathan Saul Caine ◽  
James V. Jones ◽  
Thorsten W. Becker
Keyword(s):  
Receiver Function ◽  
Shear Zones ◽  
Receiver Functions ◽  
Denali Fault ◽  
Structural Intensity ◽  
Aleutian Arc ◽  
Flat Slab Subduction ◽  
Geographic Gradient ◽  
Orogenic Belts ◽  
Low Amplitude

Abstract Azimuthal variations in receiver function conversions can image lithospheric structural contrasts and anisotropic fabrics that together compose tectonic grain. We apply this method to data from EarthScope Transportable Array in Alaska and additional stations across the northern Cordillera. The best-resolved quantities are the strike and depth of dipping fabric contrasts or interfaces. We find a strong geographic gradient in such anomalies, with large amplitudes extending inboard from the present-day subduction margin, the Aleutian arc, and an influence of flat-slab subduction of the Yakutat microplate north of the Denali fault. An east–west band across interior Alaska shows low-amplitude crustal anomalies. Anomaly amplitudes correlate with structural intensity (density of aligned geological elements), but are the highest in areas of strong Cenozoic deformation, raising the question of an influence of current stress state. Imaged subsurface strikes show alignment with surface structures. We see concentric strikes around arc volcanoes implying dipping magmatic structures and fabric into the middle crust. Regions with present-day weaker deformation show lower anomaly amplitudes but structurally aligned strikes, suggesting pre-Cenozoic fabrics may have been overprinted or otherwise modified. We observe general coherence of the signal across the brittle-plastic transition. Imaged crustal fabrics are aligned with major faults and shear zones, whereas intrafault blocks show imaged strikes both parallel to and at high angles to major block-bounding faults. High-angle strikes are subparallel to neotectonic deformation, seismicity, fault lineaments, and prominent metallogenic belts, possibly due to overprinting and/or co-evolution with fault-parallel fabrics. We suggest that the underlying tectonic grain in the northern Cordillera is broadly distributed rather than strongly localized. Receiver functions thus reveal key information about the nature and continuity of tectonic fabrics at depth and can provide unique insights into the deformation history and distribution of regional strain in complex orogenic belts.

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.

Mid-crustal magma reservoirs at Cleveland and Akutan Volcano imaged through novel receiver function analyses

2020 ◽  
Author(s):  
Helen Janiszewski ◽  
Lara Wagner ◽  
Diana Roman
Keyword(s):  
Receiver Function ◽  
Imaging Techniques ◽  
Seismic Source ◽  
Scattered Wave ◽  
Receiver Functions ◽  
Aleutian Arc ◽  
Source Studies ◽  
End Members ◽  
Future Work ◽  
Lower Crustal

<p>Processes related to magma formation, transport, emplacement, and eruption at volcanoes are linked by structures that transect the entire crust, but imaging the mid- to lower-crustal portions of these magmatic systems has been a longstanding challenge. Tomography, local seismic source studies, geodetic, and geochemical constraints are typically most sensitive to shallow storage and/or have insufficient resolution at these depths. Scattered wave seismic imaging techniques, particularly receiver function analyses, provide a promising pathway towards imaging the mid- to deep-crustal magmatic structure beneath volcanoes with only a modest number of broadband seismic instruments (N < 10). Using seismic data from two recently-active volcanoes in Alaska’s Aleutian arc, Akutan and Cleveland, we demonstrate the feasibility of seismically imaging crustal magmatic structure with only three and seven local broadband seismometers at each volcano, respectively. The two volcanoes have significantly differing eruptive histories: Akutan last erupted in 1992 and has since experienced only experienced a shallow dike intrusion in 1996, whereas Cleveland is one of the most frequently-erupting volcanoes in the Aleutian arc. Both also have significantly different depths-to-slab, with Cleveland representing one of the global shallow end members at ~ 70 km depth, and a more globally-average depth of 85 km at Akutan. Receiver functions reveal different underlying crustal magmatic structures, with a mid-crustal sill-like structure that has a well-defined top and base beneath Akutan, and a thicker and deeper magmatic region with less abrupt boundaries beneath Cleveland. Future work using similar approaches will enable an unprecedented comparative examination of magmatic systems beneath sparsely instrumented volcanoes globally.</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 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).

Petrologic and thermal constraints on the origin of leucogranites in collisional orogens

2004 ◽  
Vol 95 (1-2) ◽  
pp. 73-85 ◽  
Author(s):  
Peter I. Nabelek ◽  
Mian Liu
Keyword(s):  
Partial Melting ◽  
Shear Zones ◽  
Source Rocks ◽  
Radioactive Isotopes ◽  
The Usa ◽  
Appalachian Orogen ◽  
Shear Heating ◽  
Heating Model ◽  
Orogenic Belts ◽  
Collisional Orogens

ABSTRACTLeucogranites are typical products of collisional orogenies. They are found in orogenic terranes of different ages, including the Proterozoic Trans-Hudson orogen, as exemplified in the Black Hills, South Dakota, and the Appalachian orogen in Maine, both in the USA, and the ongoing Himalayan orogen. Characteristics of these collisional leucogranites show that they were derived from predominantly pelitic sources at the veining stages of deformation and metamorphism in upper plates of thickened crusts. Once generated, the leucogranite magmas ascended as dykes and were emplaced within shallower parts of their source sequences. In these orogenic belts, there was a strong connection between deformation, metamorphism and granite generation. However, the heat sources needed for partial melting of the source rocks remain controversial. Lack of evidence for significant intrusion of mafic magmas necessary to cause melting of upper plate source rocks suggests that leucogranite generation in collisional orogens is mainly a crustal process.The present authors evaluate five types of thermal models which have previously been proposed for generating leucogranites in collisional orogens. The first, a thickened crust with exponentially decaying distribution of heat-producing radioactive isotopes with depth, has been shown to be insufficient for heating the upper crust to melting conditions. Four other models capable of raising the crustal temperatures sufficiently to initiate partial melting of metapelites in thickened crust include: (1) thick sequences of sedimentary rocks with high amounts of internal radioactive heat production; (2) decompression melting; (3) thinning of mantle lithosphere; and (4) shear-heating. The authors show that, for reasonable boundary conditions, shear-heating along crustal-scale shear zones is the most viable process to induce melting in upper plates of collisional orogens where pelitic source lithologies are usually located. The shear-heating model directly links partial melting to the deformation and metamorphism that typically precede leucogranite genera

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>

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.

A Discussion on natural strain and geological structure - Archaean deformation patterns in southern Africa

1976 ◽  
Vol 283 (1312) ◽  
pp. 313-331 ◽  
Keyword(s):  
Shear Zones ◽  
Geological Structure ◽  
Mobile Belt ◽  
Deformation Pattern ◽  
Strain Measurements ◽  
Regional Deformation ◽  
Greenstone Belts ◽  
Extension Direction ◽  
Orogenic Belts ◽  
Limpopo Belt

Strain measurements have been made to help quantify the intensity of deformation and amount of displacement across Archaean greenstone belts in Rhodesia and Botswana and across the gneisses of the Limpopo mobile belt. The area has been divided into three domains based on the orientation of the finite strain fabric and the orientation of the maximum extension direction in associated shear zones. The domains are considered to have different movement patterns and to be similar to small orogenic belts. Early deformation within the greenstone belts accompanied the intrusion of the diaipric granites, but there was also bulk translation and rotation of greenstone belt and gneiss leading to imbrication of the stratigraphic pile and the formation of large nappes of overturned rock. This was followed by regional phases of deformation which affected all the greenstone belts and the gneisses of the Limpopo belt. Detailed strain measurements show a variation in amount of shortening during this phase, from under 30 % across the Shabani-Bellingwe belt in central Rhodesia, to over 60 % across the Tati and Matsitama belts in northern Botswana. Many local variations in intensity of deformation occur within large ductile shear zones and deviations from plane strain may be partly due to such rotational deformation. The regional deformation pattern suggests that there was movement of the Rhodesian craton approximately 200 km to the southwest relative to the gneisses of the Limpopo belt, producing a dominantly flattening deformation in the southwest of Rhodesia, but dominantly simple shear with a nearly horizontal sinistral movement, in the southeast.

Sign in / Sign up

Export Citation Format

Share Document