Mantle heterogeneities and their significance: results from Lithoprobe seismic reflection and refraction – wide-angle reflection studiesThis article is one of a series of papers published in this Special Issue on the theme Lithoprobe — parameters, processes, and the evolution of a continent.Lithoprobe Contribution 1486.

2010 ◽  
Vol 47 (4) ◽  
pp. 409-443 ◽  
Author(s):  
Ron M. Clowes ◽  
Don J. White ◽  
Zoltan Hajnal

Within Lithoprobe’s 10 transects, data from more than 20 000 km of multichannel seismic (MCS) reflection profiling and 12 refraction – wide-angle reflection (R/WAR) surveys were acquired. While the main results related to crustal structure, the data also indicated substantial heterogeneity in the lithospheric mantle. Images of fossilized subduction zones from the Eocene to the Neoarchean demonstrate that current plate tectonic processes have been active for more than 2.6 Ga. The Canadian Cordillera has a thin (50–60 km) lithosphere that is likely receiving some dynamic support from the asthenosphere below. Vestiges of the last stage of accretionary tectonic processes that formed the Archean Superior craton are indicated by an unusual anisotropic high velocity layer that may represent relic oceanic lithosphere. Within the Paleoproterozoic Trans-Hudson Orogen, a restricted region of upper mantle P-wave velocity anisotropy is identified with the continental collision between the bounding Hearne and Superior cratons. In the Archean Hearne and Wyoming provinces, two dipping structures within the sub-crustal lithosphere are interpreted as subduction features related to the assembly of the two cratons. Finite-difference modeling of long-offset data (over 1300 km) reveals fine-scale heterogeneities within a layer between 90 and 150 km in the continental lithosphere, perhaps formed through lateral flow or deformation within the upper mantle. Based on Lithoprobe data, heterogeneities within the lithospheric mantle are reasonably common. They have a wide range of seismic signatures, include many different types and show differing scales. Nevertheless, their extent in the lithospheric mantle is considerably less than in the crust.

2021 ◽  
Author(s):  
Banafsheh Vahdati ◽  
Seyed Ahmad Mazaheri

<p>Mashhad granitoid complex is part of the northern slope of the Binalood Structural Zone (BSZ), Northeast of Iran, which is composed of granitoids and metamorphic rocks. This research presents new petrological and geochemical whole-rock major and trace elements analyses in order to determine the origin of granitoid rocks from Mashhad area. Field and petrographic observations indicate that these granitoid rocks have a wide range of lithological compositions and they are categorized into intermediate to felsic intrusive rocks (SiO<sub>2</sub>: 57.62-74.39 Wt.%). Qartzdiorite, tonalite, granodiorite and monzogranite are common granitoids with intrusive pegmatite and aplitic dikes and veins intruding them. Based on geochemical analyses, the granitoid rocks are calc-alkaline in nature and they are mostly peraluminous. On geochemical variation diagrams (major and minor oxides versus silica) Na<sub>2</sub>O and K<sub>2</sub>O show a positive correlation with silica while Al<sub>2</sub>O<sub>3</sub>, TiO<sub>2</sub>, CaO, Fe<sub>2</sub>O<sub>3</sub>, and MgO show a negative trend. Therefore fractional crystallization played a considerable role in the evolution of Mashhad granitoids. Based on the spider diagrams, there are enrichments in LILE and depletion in HFSE. Low degrees of melting or crustal contamination may be responsible for LILE enrichment. Elements such as Pb, Sm, Dy and Rb are enriched, while Ba, Sr, Nd, Zr, P, Ti and Yb (in monzogranites) are all depleted. LREE enrichment and HREE depletion are observed in all samples on the Chondrite-normalized REE diagram. Similar trends may be evidence for the granitoids to have the same origin. Besides, LREE enrichment relative to HREE in some samples can indicate the presence of garnet in their source rock. Negative anomalies of Eu and Yb are observed in monzogranites. Our results show that Mashhad granitoid rocks are orogenic related and tectonic discrimination diagrams mostly indicate its syn-to-post collisional tectonic setting. No negative Nb anomaly compared with MORB seems to be an indication of non-subduction zone related magma formation. According to the theory of thrust tectonics of the Binalood region, the oceanic lithosphere of the Palo-Tethys has subducted under the Turan microplate. Since the Mashhad granitoid outcrops are settled on the Iranian plate, this is far from common belief that these granitoid rocks are related to the subduction zones and the continental arcs. The western Mashhad granitoids show more mafic characteristics and are possibly crystallized from a magma with sedimentary and igneous origin. Thus, Western granitoid outcrops in Mashhad are probably hybrid type and other granitoid rocks, S and SE Mashhad are S-type. Evidences suggest that these continental collision granitoid rocks are associated with the late stages of the collision between the Iranian and the Turan microplates during the Paleo-Tethys Ocean closure which occurred in the Late Triassic.</p>


2021 ◽  
pp. M56-2019-44
Author(s):  
Philip T. Leat ◽  
Aidan J. Ross ◽  
Sally A. Gibson

AbstractAbundant mantle-derived ultramafic xenoliths occur in Cenozoic (7.7-1.5 Ma) mafic alkaline volcanic rocks along the former active margin of West Antarctica, that extends from the northern Antarctic Peninsula to Jones Mountains. The xenoliths are restricted to post-subduction volcanic rocks that were emplaced in fore-arc or back-arc positions relative to the Mesozoic-Cenozoic Antarctic Peninsula volcanic arc. The xenoliths are spinel-bearing, include harzburgites, lherzolites, wehrlites and pyroxenites, and provide the only direct evidence of the composition of the lithospheric mantle underlying most of the margin. The harzburgites may be residues of melt extraction from the upper mantle (in a mid-ocean ridge type setting), that accreted to form oceanic lithosphere, which was then subsequently tectonically emplaced along the active Gondwana margin. An exposed highly-depleted dunite-serpentinite upper mantle complex on Gibbs Island, South Shetland Islands, supports this interpretation. In contrast, pyroxenites, wehrlites and lherzolites reflect percolation of mafic alkaline melts through the lithospheric mantle. Volatile and incompatible trace element compositions imply that these interacting melts were related to the post-subduction magmatism which hosts the xenoliths. The scattered distribution of such magmatism and the history of accretion suggest that the dominant composition of sub-Antarctic Peninsula lithospheric mantle is likely to be harzburgitic.


Geosphere ◽  
2019 ◽  
Vol 15 (6) ◽  
pp. 2018-2042 ◽  
Author(s):  
Melissa V. Bernardino ◽  
Craig H. Jones ◽  
William Levandowski ◽  
Ian Bastow ◽  
Thomas J. Owens ◽  
...  

Abstract The Isabella anomaly, a prominent upper-mantle high-speed P-wave anomaly located within the southern Great Valley and southwestern foothills of the Sierra Nevada, has been interpreted either as foundering sub-Sierran lithosphere or as remnant oceanic lithosphere. We used Vp/Vs anisotropy tomography to distinguish among the probable origins of the Isabella anomaly. S waveforms were rotated into the Sierran SKSFast and SKSSlow directions determined from SKS-splitting studies. Teleseismic P-, SFast-, SSlow-, SKSFast-, and SKSSlow-wave arrival times were then inverted to obtain three-dimensional (3-D) perturbations in Vp, Vp/VsMean, and percent azimuthal anisotropy using three surface wave 3-D starting models and one one-dimensional (1-D) model. We observed the highest Vp/Vs anomalies associated with slower velocities in regions marked by young volcanism, with the largest of these anomalies being the Mono anomaly under the Long Valley region, which extends to depths of at least 75 km. Peak Vp/Vs perturbations of +4% were found at 40 km depth. The low velocities and high Vp/Vs values of this anomaly could be related to partial melt. The high wave speeds of the Isabella anomaly coincide with low Vp/Vs values with peak perturbations of −2%, yet they do not covary spatially. The P-wave inversion imaged the Isabella anomaly as a unimodal eastward-plunging body. However, the volume of that Isabella anomaly contains three separate bodies as defined by varying Vp/Vs values. High speeds, regionally average Vp/Vs values (higher than the other two anomalies), and lower anisotropy characterize the core of the Isabella anomaly. The western and shallowest part has high wave speeds and lower Vp/Vs values than the surrounding mantle. The eastern and deepest part of the anomaly also contains high speeds and lower Vp/Vs values but exhibits higher anisotropy. We considered combinations of varying temperature, Mg content (melt depletion), or modal garnet to reproduce our observations. Our results suggest that the displaced garnet-rich mafic root of the Mesozoic Sierra Nevada batholith is found in the core of the Isabella anomaly. If remnant oceanic lithosphere exists within the Isabella anomaly, it most likely resides in the shallow, westernmost feature. Within the Sierra Nevada, the highest upper-mantle anisotropy is largely contained within the central portion of the range and the adjacent Great Valley. Anisotropy along the Sierra crest is shallow and confined to the lithosphere between 20 and 40 km depth. Directly below, there is a zone of low anisotropy (from 170 to 220 km depth), low velocities, and high Vp/Vs values. These features suggest the presence of vertically upwelling asthenosphere and consequent horizontal flow at shallower depths. High anisotropy beneath the adjacent western foothills and Great Valley is found at ∼120 km depth and could represent localized mantle deformation produced as asthenosphere filled in a slab gap.


Author(s):  
R.S. Alekseev, ◽  
◽  
Yu.L. Rebetsky ◽  

The Himalayan-Tibetan orogen is one of the active orogens on Earth. The processes caused by the collision of two continents have attracted attention of many researchers, and over the past decades, a large amount of geological and geophysical data has accumulated, on which models of the evolution of the region are based. The paper presents a model of the evolution of the Tibet plateau and the adjacent mountain chains, which complies with the modern concepts of the structure of the crust. The reference parameters of this model are the data on the values of stresses and on the patterns of the spatial distribution of principal stresses obtained in our own tectonophysical studies in region, as well as in other intracontinental orogens and in subduction zones between lithospheric plates. The basic assumptions of the model are the factors of the long stage of the Indian plate underthrusting beneath the Eurasian continent, metamorphic processes in the submerged slab (oceanic lithosphere) and in the continental lithosphere above it, combination of absolute horizontal movements of the Eurasian and Indian plates, small-scale convection in the upper mantle and vertical movements of matter, both in the continental lithosphere itself and in the upper mantle.


2020 ◽  
Vol 105 (8) ◽  
pp. 1191-1203
Author(s):  
Yanfei Zhang ◽  
Xuran Liang ◽  
Chao Wang ◽  
Zhenmin Jin ◽  
Lüyun Zhu ◽  
...  

Abstract Sedimentary diapirs can be relaminated to the base of the lithosphere during slab subduction, where they can interact with the ambient lithospheric mantle to form variably metasomatized zones. Here, high-pressure experiments in sediment-harzburgite systems were conducted at 1.5–2.5 GPa and 800–1300 °C to investigate the interaction between relaminated sediment diapirs and lithospheric mantle. Two end-member processes of mixed experiments and layered (reaction) experiments were explored. In the first end-member, sediment and harzburgite powders were mixed to a homogeneous proportion (1:3), whereas in the second, the two powders were juxtaposed as separate layers. In the first series of experiments, the run products were mainly composed of olivine + orthopyroxene + clinopyroxene + phlogopite in subsolidus experiments, while the phase assemblages were then replaced by olivine + orthopyroxene + melt (or trace phlogopite) in supersolidus experiments. Basaltic and foiditic melts were observed in all supersolidus mixed experiments (~44–52 wt% SiO2 at 1.5 GPa, ~35–43 wt% SiO2 at 2.5 GPa). In the phlogopite-rich experiment (PC431, 1.5 GPa and 1100 °C), the formed melts had low alkali contents (~<2 wt%) and K2O/Na2O ratios (~0.4–1.1). In contrast, the quenched melt in phlogopite-free/poor experiments showed relatively higher alkali contents (~4–8 wt%) and K2O/Na2O ratios (~2–5). Therefore, the stability of phlogopite could control the bulk K2O and K2O/Na2O ratios of magmas derived from the sediment-metasomatized lithospheric mantle. In layered experiments, a reaction zone dominated by clinopyroxene + amphibole (or orthopyroxene) was formed because of the reaction between harzburgite and bottom sediment-derived melts (~62.5–67 wt% SiO2). The total alkali contents and K2O/Na2O ratios of the formed melts were about 6–8 wt% and 1.5–3, respectively. Experimentally formed melts from both mixed and reaction experiments were rich in large ion lithosphile elements and displayed similar patterns with natural potassium-rich arc lavas from oceanic subduction zones (i.e., Mexican, Sunda, Central American, and Aleutian). The experimental results demonstrated that bulk sediment diapirs, in addition to sediment melt, may be another possible mechanism to transfer material from a subducting slab to an upper mantle wedge or lithospheric mantle. On the other hand, the breakdown of phlogopite may play an important role in the mantle source that produces potassium-rich arc lavas in subduction zones.


Author(s):  
Brandon VanderBeek ◽  
Manuele Faccenda

<p>Despite the well-established anisotropic nature of Earth’s upper mantle, the influence of elastic anisotropy on teleseismic tomographic images remains largely ignored. In subduction zones, unmodeled anisotropic heterogeneity can lead to substantial isotropic velocity artefacts that may be misinterpreted as compositional heterogeneities (e.g. Bezada et al., 2016). Recent studies have demonstrated the possibility of inverting P-wave delay times for the strength and orientation of seismic anisotropy assuming a hexagonal symmetry system (e.g. Huang et al., 2015; Munzarová et al., 2018). However, the ability of P-wave delay times to constrain complex anisotropic patterns, such as those expected in subduction settings, remains unclear as the aforementioned methods are tested using ideal self-consistent data (i.e. data produced using the assumptions built into the tomography algorithm) generated from simplified synthetic models. Here, we test anisotropic P-wave imaging methods on data generated from geodynamic simulations of subduction. Micromechanical models of polymineralic aggregates advected through the simulated flow field are used to create an elastic model with up to 21 independent coefficients. We then model the teleseismic wavefield through this fully anisotropic model using SPECFEM3D coupled with AxiSEM. P-wave delay times across a synthetic seismic array are measured using conventional cross-correlation techniques and inverted for isotropic velocity and the strength and orientation of anisotropy using travel-time tomography methods. We propose and validate approximate analytic finite-frequency sensitivity kernels for the simplified anisotropic parameters. Our results demonstrate that P-wave delays can reliably recover horizontal and vertical changes in the azimuth of anisotropy. However, substantial isotropic artefacts remain in the solution when only inverting for azimuthal anisotropy parameters. These isotropic artefacts are largely removed when inverting for the dip as well as the azimuth of the anisotropic symmetry axis. Due to the relative nature of P-wave delay times, these data generally fail to reconstruct anisotropic structure that is spatially uniform over large scales. To overcome this limitation, we propose a joint inversion of SKS splitting intensity with P-wave delay times. Preliminary results demonstrate that this approach improves the recovery of the magnitude and azimuth of anisotropy. We conclude that teleseismic P-wave travel-times are a useful observable for probing the 3D distribution of upper mantle anisotropy and that anisotropic inversions should be explored to better understand the nature of isotropic velocity anomalies in subduction settings.</p>


1992 ◽  
Vol 29 (7) ◽  
pp. 1530-1548 ◽  
Author(s):  
B. C. Zelt ◽  
R. M. Ellis ◽  
R. M. Clowes ◽  
E. R. Kanasewich ◽  
I. Asudeh ◽  
...  

As part of the Lithoprobe Southern Cordillera transect, seismic refraction data were recorded along a 330 km long strike profile in the Intermontane belt. An iterative combination of two-dimensional traveltime inversion and amplitude forward modelling was used to interpret crust and upper mantle P-wave velocity structure. This region is characterized by (i) a thin near-surface layer with large variations in velocity between 2.8 and 5.4 km/s, and low-velocity regions that correlate well with surface expressions of Tertiary sedimentary and volcanic rocks; (ii) an upper and middle crust with low average velocity gradient, possibly a weak low-velocity zone, and lateral velocity variations between 6.0 and 6.4 km/s; (iii) a distinctive lower crust characterized by significantly higher average velocities relative to midcrustal values beginning at 23 km depth, approximately 8 km thick with average velocities of 6.5 and 6.7 km/s at top and base; (iv) a depth to Moho, as defined by wide-angle reflections, that averages 33 km with variations up to 2 km; and (v) a Moho transition zone of depth extent 1–3 km, below which lies the upper mantle with velocities decreasing from 7.9 km/s in the south to 7.7 km/s in the north. Where the refraction line obliquely crosses a Lithoprobe deep seismic-reflection profile, good agreement is obtained between the interpreted reflection section and the derived velocity structure model. In particular, depths to wide-angle reflectors in the upper crust agree with depths to prominent reflection events, and Moho depths agree within 1 km. From this comparison, the upper and middle crust probably comprise the upper part of the Quesnellia terrane. The lower crust from the refraction interpretation does not show the division into two components, parautochthonous and cratonic North America, that is inferred from the reflection data, indicating that their physical properties are not significantly different within the resolution of the refraction data. Based on these interpretations, the lower lithosphere of Quesnellia is absent and presumably was recycled in the mantle. At a depth of ~ 16 km below the Moho, an upper mantle reflector may represent the base of the present lithosphere.


Geophysics ◽  
1989 ◽  
Vol 54 (12) ◽  
pp. 1564-1574 ◽  
Author(s):  
B. S. Byun ◽  
D. Corrigan ◽  
J. E. Gaiser

A new velocity analysis technique is presented for analyzing moveout of signals on multichannel surface seismic or VSP data. An approximate, skewed hyperbolic moveout formula is derived for horizontally layered, transversely isotropic media. This formula involves three measurement parameters: the average vertical velocity and horizontal and skew moveout velocities. By extending Dix‐type hyperbolic moveout analysis, we obtain improved coherence over large source‐geophone offsets for more accurate moveout correction. Compared with the stacking velocity obtained by simple hyperbolic analysis methods, the three velocity parameters estimated by this technique contain more physically meaningful geologic information regarding the anisotropy and/or velocity heterogeneity of the subsurface. Synthetic P‐wave model experiments demonstrate that the skewed hyperbolic moveout formula yields an excellent fit to time‐distance curves over a wide range of ray angles. Consequently, the measurement parameters are shown to reflect adequately the characteristics of velocity dependence on ray angle, i.e., velocity anisotropy. The technique is then applied to two field offset VSP data sets to measure and analyze the velocity parameters. The results show that the apparent anisotropy, defined as the ratio between the horizontal moveout velocity and average vertical velocity, correlates reasonably well with lithology. Highly anisotropic shale and chalk exhibit higher horizontal‐to‐vertical velocity ratios and sandstones show lower ratios.


Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1292-1309 ◽  
Author(s):  
Ilya Tsvankin

Although orthorhombic (or orthotropic) symmetry is believed to be common for fractured reservoirs, the difficulties in dealing with nine independent elastic constants have precluded this model from being used in seismology. A notation introduced in this work is designed to help make seismic inversion and processing for orthorhombic media more practical by simplifying the description of a wide range of seismic signatures. Taking advantage of the fact that the Christoffel equation has the same form in the symmetry planes of orthorhombic and transversely isotropic (TI) media, we can replace the stiffness coefficients by two vertical (P and S) velocities and seven dimensionless parameters that represent an extension of Thomsen's anisotropy coefficients to orthorhombic models. By design, this notation provides a uniform description of anisotropic media with both orthorhombic and TI symmetry. The dimensionless anisotropic parameters introduced here preserve all attractive features of Thomsen notation in treating wave propagation and performing 2-D processing in the symmetry planes of orthorhombic media. The new notation has proved useful in describing seismic signatures outside the symmetry planes as well, especially for P‐waves. Linearization of P‐wave phase velocity in the anisotropic coefficients leads to a concise weak‐anisotropy approximation that provides good accuracy even for models with pronounced polar and azimuthal velocity variations. This approximation can be used efficiently to build analytic solutions for various seismic signatures. One of the most important advantages of the new notation is the reduction in the number of parameters responsible for P‐wave velocities and traveltimes. All kinematic signatures of P‐waves in orthorhombic media depend on just the vertical velocity [Formula: see text] and five anisotropic parameters, with [Formula: see text] serving as a scaling coefficient in homogeneous media. This conclusion, which holds even for orthorhombic models with strong velocity anisotropy, provides an analytic basis for application of P‐wave traveltime inversion and data processing algorithms in orthorhombic media.


2021 ◽  
Author(s):  
Michael J. Schmidtke ◽  
Ruth Keppler ◽  
Jacek Kossak-Glowczewski ◽  
Nikolaus Froitzheim ◽  
Michael Stipp

Abstract. Subduction and exhumation are key processes in the formation of orogenic systems across the world, for example, in the European Alps. For geophysical investigations of these orogens, it is essential to understand the petrophysical properties of the rocks involved. These are the result of a complex interaction of mineral composition and rock fabric including mineral textures (i.e. crystallographic preferred orientations). In this study we present texture-derived elastic anisotropy data for a representative set of different lithologies involved in the Alpine orogeny. Rock samples were collected in the Lago di Cignana area in Valtournenche, in the Italian Northwestern Alps. At this locality a wide range of units of continental and oceanic origin with varying paleogeographic affiliations and tectono-metamorphic histories are accessible. Their mineral textures were determined by time-of-flight neutron diffraction at the Frank Laboratory of Neutron Physics at the JINR in Dubna, Russia. From these data the elastic properties of the samples were calculated. The data set includes representative lithologies from a subduction-exhumation-setting. In subducted lithologies originating from the oceanic crust, the elastic anisotropies range from 1.4 to 5.0 % with average P-wave velocities of 7.01–8.24 km/s and VP / VS-ratios of 1.71–1.76. In the metasediments of the former accretionary prism the elastic anisotropies range from 4.7 to 8.2 %. This tectonic setting displays average P-wave velocities of 6.47–7.23 km/s and VP / VS-ratios of 1.60–1.76. Continental crust which is incorporated in the collisional orogen shows elastic anisotropies ranging from 1.8 to 2.8 % with average P-wave velocities of 6.42–6.51 km/s and VP / VS-ratios of 1.56–1.60. Our results suggest that mafic and felsic rocks in subduction zones at depth may be discriminated by a combination of seismic signatures: lower anisotropy and higher VP / VS ratio for mafic rocks, higher anisotropy and lower VP / VS ratio for felsic rocks and metasediments.


Sign in / Sign up

Export Citation Format

Share Document