Lithosphere-asthenosphere Boundary Research Articles

Upper-mantle discontinuities beneath Australia from transdimensional Bayesian inversions using multimode surface waves and receiver functions

SUMMARYRadially anisotropic S-wave structures under the permanent seismic stations in Australia are reconstructed using multimode surface waves (SWs) and receiver functions (RFs) in a framework of the Bayesian inference. We have developed a fully nonlinear method of joint inversions incorporating P-RFs and multimode Rayleigh and Love waves, based on the transdimensional Hierarchical Bayesian formulation. The method allows us to estimate a probabilistic Earth model taking account of the complexity and uncertainty of Earth structure, by treating the model parameters and data errors as unknowns. The Parallel Tempering algorithm is employed for the effective parameter search based on the reversible-jump Markov Chain Monte Carlo method. The use of higher modes enables us to enhance the sensitivity to the depth below the continental asthenosphere. Synthetic experiments indicate the importance of higher mode SWs for the better recovery of radial anisotropy in the whole depth range of the upper mantle. The method is applied to five Global Seismographic Network stations in Australia. While the S-wave models in eastern Australia show shallow lithosphere–asthenosphere boundary (LAB) above 100 km depth, those in central and Western Australia exhibit both mid-lithosphere discontinuities (MLDs) and LAB. Also, seismic velocity jumps equivalent to the Lehmann discontinuity (L-D) are found in all seismic stations. The L-D under the Australian continents is found at around 200–300 km depth, depending on locations. Radial anisotropy in the depth range between LAB and L-D tends to show faster SH anomalies, which may indicate the effects of horizontal shear underneath the fast-moving Australian plate.

Seismic discontinuities in the lithospheric mantle at the Dead Sea Transform

SUMMARYThe Dead Sea Transform (DST) was formed in the mid-Cenozoic, about 18 Myr ago, as a result of the breakaway of the Arabian plate from the African plate. Higher resolution information about the sub-Moho structure is still sparse in this region. Here, we study seismic discontinuities in the mantle lithosphere in the region of the DST using a modified version of the P- and S-receiver function methods. We use open data from permanent and temporary seismic stations. The results are displayed in a number of depth profiles through the study area. The Moho is observed on both sides of the transform at nearly 40 km depth by S-to-p and in P-to-s converted signals. The lithosphere–asthenosphere boundary (LAB) on the eastern side of the DST is observed near 180–200 km depth, which is according to our knowledge the first LAB observation at that depth in this region. This observation could lead to the conclusion that the thickness of the Arabian lithosphere east of the DST is likely cratonic. In addition, we observe in the entire area a negative velocity gradient at 60–80 km depth, which was previously interpreted as LAB.

A Lithosphere‐Asthenosphere Boundary and Partial Melt Estimated Using Marine Magnetotelluric Data at the Central Middle Atlantic Ridge

AbstractThe differential motion between the lithosphere and the asthenosphere is aseismic, so the magnetotelluric (MT) method plays an important role in studying the depth and nature of the lithosphere‐asthenosphere boundary (LAB). In March 2016, we deployed 39 marine MT instruments across the Middle Atlantic Ridge (MAR), 2,000 km away from the African coast, to study the evolution of the LAB with ages out to 45 million years (My). The MT acquisition time was limited to about 60 days by battery life. After analyzing dimensionality and coast effects for the MT data, determinant data were inverted for two‐dimensional resistivity models along two profiles north and south of the Chain Fracture Zone (CFZ). The imaged thickness of the lithospheric lid (>100 Ωm) ranges from 20 to 80 km, generally thickening with age. In the north of CFZ, punctuated low‐resistivity anomalies (<1 Ωm), likely associated with potential partial melts, occur along its base. In the south of CFZ, the base of the resistive lid is demarcated by a low‐resistivity channel (<1 Ωm) most likely fed by deeper melts. Sensitivity analyses and structural recovery tests indicate the robustness of these features. Resistivity models are in good agreement with results of seismic data. These results imply that partial melt is persistent over geologic timescales and that the LAB is dynamic features fed by upward percolation of mantle melt. The melt fraction is about 1–7% based on the resistivity, temperature, pressure, and hydrous basalt models, which is consistent with petrophysical observations.

Sharpness of the Midlithospheric Discontinuities and Craton Evolution in North China

AbstractThe lithospheric layering and evolution of Archean cratons can be effectively constrained by seismic properties. We present fine‐scale lithospheric discontinuity images of the North China Craton (NCC) using the seismic daylight imaging (SDI) technique, which is also known as the autocorrelogram method. With the aid of large‐scale surface‐wave velocity models, the characters of autocorrelograms related to lithospheric discontinuities in the NCC are investigated. Spatial variation in discontinuity characteristics is revealed across the craton. In the modified/destroyed eastern NCC, both the midlithospheric discontinuity (MLD) and the lithosphere‐asthenosphere boundary (LAB) are rather sharp (~2 km). Various dynamic processes, including juvenile underplating, asthenospheric upwelling, and corner flow induced by paleo‐Pacific Plate subduction, may contribute to such sharp interfaces. In the stable western NCC, P wave reflections cannot easily reveal the LAB, suggesting a diffuse (thicker than 30 km) lithosphere‐asthenosphere transition (LAT). Multiple MLDs are also found that may have arisen during craton formation and multiple phases of rejuvenation of the cratonic lithosphere as evidenced by the petrological results. The use of higher‐frequency data enables us to distinguish among the lithospheric discontinuity styles and favor multiple origins of the MLDs and LAB in the cratonic lithosphere.

Lithosphere Mantle Density of the North China Craton

AbstractWe constrain the lithospheric mantle density of the North China Craton (NCC) at both in situ and standard temperature‐pressure (STP) conditions from gravity data. The lithosphere‐asthenosphere boundary (LAB) depth is constrained by our new thermal model, which is based on a new regional heat flow data set and a recent regional crustal model NCcrust. The new thermal model shows that the thermal lithosphere thickness is <120 km in most of the NCC, except for the northern and southern parts with the maximum depth of 170 km. The gravity calculations reveal a highly heterogeneous density structure of the lithospheric mantle with in situ and STP values of 3.22–3.29 and 3.32–3.40 g/cm3, respectively. Thick and reduced‐density cratonic‐type lithosphere is preserved mostly in the southern NCC. Most of the Eastern Block has a thin (90–140 km) and high‐density lithospheric mantle. Most of the Western Block has a high‐density lithospheric mantle and a thin (80–110 km) lithosphere typical of Phanerozoic regions, which suggests that the Archean lithosphere is no longer present there. We conclude that in almost the entire NCC the lithosphere has lost its cratonic characteristics by geodynamic processes that include, but are not limited to, the Paleozoic closure of the Paleo‐Asian Ocean in the north, the Mesozoic Yangtze Craton flat subduction in the south, the Mesozoic Pacific subduction in the east, the Cenozoic remote response to the Indian‐Eurasian collision in the west, and the Cenozoic extensional tectonics (possibly associated with the slab roll‐back) in the center.

Lithospheric Structure of Eastern Tibetan Plateau from Terrestrial and Satellite Gravity Data Modeling: Implication for Asthenospheric Underplating

Abstract The lithosphere of the eastern Tibetan plateau is underlain by a low-velocity zone at shallow depths which is interpreted as asthenospheric material in the upper-most mantle in various seismic tomography studies. The driving mechanism for the presence of asthenospheric material in the upper-most mantle is not well understood, and the spatial extent of the asthenospheric material is not well delineated. We use 2.5D gravity models to assess what drove the asthenospheric flow upwards in the past and determine the lateral extent of the asthenospheric material in the upper-most mantle. The models also allow us to determine the Indian slab configuration below the Tibetan plateau. The gravity models show that lithospheric thickness increases from ~120 km in the central and eastern parts of the plateau to ~150 km in the west, indicating that the lithosphere in the central and eastern parts of the plateau may have been delaminated. The ~30 km shallower Lithosphere-Asthenosphere Boundary in the central and eastern Tibetan plateau may indicate that asthenospheric flow could have been induced in the past by a combination of lithospheric delamination and a slab break-off event of the Greater Indian slab. The spatial extent of the asthenospheric material in the upper-most mantle beneath the Tibetan plateau is ~15,000 km2 (N−S length=500 km and thickness=30 km) between 85°E and 88°E, which could even extend east of 92°E. The Indian slab is dipping more steeply in the east. The slab dip along the Indian plate increases from ~10° in the west to ~18° in the central (~87°E) and ~25° in the eastern part (~91°E) of the plateau, indicating that the style of lithospheric deformation changes from underthrusting to slab roll-back from west to east.

Clinopyroxene and Garnet Mantle Cargo in Kimberlites as Probes of Dharwar Craton Architecture and Geotherms, with Implications for Post-1·1 Ga Lithosphere Thinning Events Beneath Southern India

AbstractThe Wajrakarur Kimberlite Field (WKF) on the Eastern Dharwar Craton in southern India hosts several occurrences of Mesoproterozoic kimberlites, lamproites and ultramafic lamprophyres, for which mantle-derived xenoliths are rare and only poorly preserved. The general paucity of mantle cargo has hampered the investigation of the nature and evolution of the continental lithospheric mantle (CLM) beneath cratonic southern India. We present a comprehensive study of the major and trace element compositions of clinopyroxene and garnet xenocrysts recovered from heavy mineral concentrates for three c.1·1 Ga old WKF kimberlite pipes (P7, P9, P10), with the goal to improve our understanding of the cratonic mantle architecture and its evolution beneath southern India. The pressure-temperature conditions recorded by peridotitic clinopyroxene xenocrysts, estimated using single-pyroxene thermobarometry, suggest a relatively moderate cratonic mantle geotherm of 40 mW/m2 at 1·1 Ga. Reconstruction of the vertical distribution of clinopyroxene and garnet xenocrysts, combined with some rare mantle xenoliths data, reveals a compositionally layered CLM structure. Two main lithological horizons are identified and denoted as layer A (∼80–145 km depth) and layer B (∼160–190 km depth). Layer A is dominated by depleted lherzolite with subordinate amounts of pyroxenite, whereas layer B comprises mainly refertilised and Ti-metasomatized peridotite. Harzburgite occurs as a minor lithology in both layers. Eclogite stringers occur within the lower portion of layer A and at the bottom of layer B near the lithosphere–asthenosphere boundary at 1·1 Ga. Refertilisation of layer B is marked by garnet compositions with enrichment in Ca, Ti, Fe, Zr and LREE, although Y is depleted compared to garnet in layer A. Garnet trace element systematics such as Zr/Hf and Ti/Eu indicate that both kimberlitic and carbonatitic melts have interacted with and compositionally overprinted layer B. Progressive changes in the REE systematics of garnet grains with depth record an upward percolation of a continuously evolving metasomatic agent. The intervening zone between layers A and B at ∼145–160 km depth is characterized by a general paucity of garnet. This ‘garnet-paucity’ zone and an overlying type II clinopyroxene-bearing zone (∼115–145 km) appear to be rich in hydrous mineral assemblages of the MARID- or PIC kind. The composite horizon between ∼115–160 km depth may represent the product of intensive melt/rock interaction by which former garnet was largely reacted out and new metasomatic phases such as type II clinopyroxene and phlogopite plus amphibole were introduced. By analogy with better-studied cratons, this ‘metasomatic horizon’ may be a petrological manifestation of a former mid-lithospheric discontinuity at 1·1 Ga. Importantly, the depth interval of the present-day lithosphere–asthenosphere boundary beneath Peninsular India as detected in seismic surveys coincides with this heavily overprinted metasomatic horizon, which suggests that post-1·1 Ga delamination of cratonic mantle lithosphere progressed all the way to mid-lithospheric depth. This finding implies that strongly overprinted metasomatic layers, such as the ‘garnet-paucity’ zone beneath the Dharwar craton, present structural zones of weakness that aid lithosphere detachment and foundering in response to plate tectonic stresses.

Receiver Function Velocity Analysis Technique and Its Application to Remove Multiples

AbstractResults of P wave receiver function (RF) analysis (e.g., H‐κ stacking, common conversion point (CCP) stacking, and migration) depend on the velocity models employed. Converted phases of crustal and upper mantle discontinuities are often interfered by multiple reflections from shallower structures, adding further complexity to the interpretation of P wave RF images. In this study, we propose a receiver function velocity analysis technique (RFVAT) that can constrain the velocity structure above discontinuities and remove multiple reflections arising from target discontinuities. By scanning the amplitude spectra of the Ps and PpPs phases, we can obtain precise arrival times of these phases for interfaces above 250 km (errors < 0.13 s) and derive the average velocities and thicknesses of the layers between interfaces in the crust and upper mantle. Synthetic experiments are used to demonstrate that crustal and upper mantle discontinuities can be coherently identified and that multiple reflections can be effectively removed in the CCP stacking images of P wave RFs using the RFVAT even with low signal‐to‐noise ratio data. By applying the proposed RFVAT to real P wave RF data from the northeastern part of the North China Craton, we achieve an improvement in CCP stacking image of lithospheric discontinuities, especially for the lithosphere‐asthenosphere boundary (LAB). Our P wave RF images reveal distinct lateral variations of the LAB depth. The lithospheric thickness extends to ~120 km beneath the central orogenic belt and decreases sharply to ~80 km at the west edge of the eastern block, consistent with previous S wave RF studies.

Large iron isotope variation in the eastern Pacific mantle as a consequence of ancient low-degree melt metasomatism

Studies of mid-ocean ridge basalts (MORB) have revealed a heterogeneous asthenospheric mantle in chemical elements and radiogenic isotopes. Here we report that MORB mantle is also heterogeneous in Fe isotopes through studying the glass samples from seamounts flanking the northern East Pacific Rise between 5° and 15°N. These samples show large Fe isotope variation with δ56Fe values (+0.03‰ to +0.36‰) exceeding the known range of MORB (+0.05‰ to +0.17‰). Such highly varied δ56Fe values cannot be well explained by seafloor alteration, fractional crystallization or partial melting processes, but instead require a source mantle significantly heterogeneous in Fe isotope compositions. Importantly, the δ56Fe values of these basalts correlate significantly with major and trace elements and Sr-Nd-Pb-Hf radiogenic isotopes, reflecting melting-induced mixing of a two-component mantle with the enriched component having heavy Fe isotope compositions dispersed as physically distinct domains in the depleted mantle matrix. The major and trace element characteristics of the enriched mantle component, as inferred from these basalts, are consistent with a low-degree melting origin. Such low-degree melts with heavy Fe isotope compositions most likely formed at sites such as the lithosphere-asthenosphere boundary beneath ocean basins, which can metasomatize the overlying oceanic lithosphere by crystallizing dikes/veins of garnet pyroxenite lithologies. Recycling of these dikes/veins with isotopically heavier Fe can readily contribute to the Fe isotope heterogeneity in the MORB mantle. However, the extremely high primitive δ56Fe values of the two alkali basalts (up to 0.34‰) require an enriched source component with unusually high δ56Fe values. We suggest that partial melts from the recycled dikes/veins of garnet pyroxenite lithologies can react with the ambient peridotitic mantle and generate a secondary garnet pyroxenite with heavier Fe isotope compositions than, but similar radiogenic isotope compositions as its precursor. Melting-induced mixing between these garnet pyroxenites (recycled and newly formed) and depleted mantle matrix can readily explain the compositional variations in elements, radiogenic isotopes and Fe isotopes observed in these seamount lavas. These new data and correlated variations offer a new dimension for understanding the origin of mantle chemical and isotopic heterogeneity in the context of chemical differentiation of the Earth.

Metasomatic control of hydrogen contents in the layered cratonic mantle lithosphere sampled by Lac de Gras xenoliths in the central Slave craton, Canada

Whether hydrogen incorporated in nominally anhydrous mantle minerals plays a role in the strength and longevity of the thick cratonic lithosphere is a matter of debate. In particular, the percolation of hydrogen-bearing melts and fluids could potentially add hydrogen to the mantle lithosphere, weaken its olivines (the dominant mineral in mantle peridotite), and cause delamination of the lithosphere's base. The influence of metasomatism on hydrogen contents of cratonic mantle minerals can be tested in mantle xenoliths from the Slave Craton (Canada) because they show extensive evidence for metasomatism of a layered cratonic mantle. Minerals from mantle xenoliths from the Diavik mine in the Lac de Gras kimberlite area located at the center of the Archean Slave craton were analyzed by FTIR for hydrogen contents. The 18 peridotites, two pyroxenites, one websterite and one wehrlite span an equilibration pressure range from 3.1 to 6.6 GPa and include samples from the shallow (≤145 km), oxidized ultra-depleted layer; the deeper (∼145–180 km), reduced less depleted layer; and an ultra-deep (≥180 km) layer near the base of the lithosphere. Olivine, orthopyroxene, clinopyroxene and garnet from peridotites contain 30–145, 110–225, 105–285, 2–105 ppm H2O, respectively. Within each deep and ultra-deep layer, correlations of hydrogen contents in minerals and tracers of metasomatism (for example light over heavy rare-earth-element ratio (LREE/HREE), high-field-strength-element (HFSE) content with equilibration pressure) can be explained by a chromatographic process occurring during the percolation of kimberlite-like melts through garnet peridotite. The hydrogen content of peridotite minerals is controlled by the compositions of the evolving melt and of the minerals and by mineral/melt partition coefficients. At the beginning of the process, clinopyroxene scavenges most of the hydrogen and garnet most of the HFSE. As the melt evolves and becomes enriched in hydrogen and LREE, olivine and garnet start to incorporate hydrogen and pyroxenes become enriched in LREE. The hydrogen content of peridotite increases with decreasing depth, overall (e.g., from 75 to 138 ppm H2O in the deep peridotites). Effective viscosity calculated using olivine hydrogen content for the deepest xenoliths near the lithosphere-asthenosphere boundary overlaps with estimates of asthenospheric viscosities. These xenoliths cannot be representative of the overall cratonic root because the lack of viscosity contrast would have caused basal erosion of lithosphere. Instead, metasomatism must be confined in narrow zones channeling kimberlite melts through the lithosphere and from where xenoliths are preferentially sampled. Such localized metasomatism by hydrogen-bearing melts therefore does not necessarily result in delamination of the cratonic root.

Lithospheric architecture of a Phanerozoic orogen from magnetotellurics: AusLAMP in the Tasmanides, southeast Australia

We present a resistivity model of the southern Tasmanides of southeastern Australia using Australian Lithospheric Architecture Magnetotelluric Project (AusLAMP) data. Modelled lower crustal conductivity anomalies resemble concentric geometries revealed in the upper crust by potential field and passive seismic data. These geometries are a key part of the crustal architecture predicted by the Lachlan Orocline model for the evolution of the southern Tasmanides, in which the Proterozoic Selwyn Block drives oroclinal rotation against the eastern Gondwana margin during the Silurian period. For the first time, we image these structures in three dimensions (3D) and show they persist below the Moho. These include a lower crustal conductor largely following the northern Selwyn Block margin.Spatial association between lower crustal conductors and both Paleozoic to Cenozoic mafic to intermediate alkaline volcanism and gold deposits suggests a genetic association i.e. fluid flow into the lower crust resulting in the deposition of conductive phases such as hydrogen, iron, sulphides and/or graphite.The 3D model resolves a different pattern of conductors in the lithospheric mantle, including northeast trending anomalies in the northern part of the model. Three of these conductors correspond to Cenozoic leucitite volcanoes along the Cosgrove mantle hotspot track which likely map the metasomatised mantle source region of these volcanoes. The northeasterly alignment of the conductors correlates with variations in the lithosphere-asthenosphere boundary (LAB) and the direction of Australian plate movement, and may be related to movement of an irregular LAB topography over the asthenosphere.By revealing the tectonic architecture of a Phanerozoic orogen and the overprint of more recent tectono-magmatic events, our resistivity model enhances our understanding of the lithospheric architecture and geodynamic processes in southeast Australia, demonstrating the ability of magnetotelluric data to image geological processes over time.

Accumulation of Stress and Strain due to an Infinite Strike-Slip Fault in an Elastic Layer Overlying a Viscoelastic Half Space of Standard Linear Solid (SLS)

Viscoelastic behaviour of materials at the lithosphere-asthenosphere boundary is observed in this paper. The observed post-seismic and aseismic deformation help us to understand the rheological properties of the mantle and asthenosphere. A quasi-static model having homogeneous, isotropic, elastic material overlying viscoelastic material of a standard linear solid (SLS) is considered. A long strike-slip fault of finite width inclined to the free surface due to a sudden movement has been studied in this work. Analytical solutions for displacement, stresses and strains are obtained before and after fault movement using a technique involving the use of Green’s functions and integral transforms, assuming that tectonic forces maintain a shear strain far away from the fault. The effect of aseismic fault movement across it is found to depend on distance, dimension, relative position and other characteristics of the fault. Also, the effect of different inclinations, elastic layer thickness and slip magnitude have been studied. The study of such earthquake fault dynamical model helps us to understand the mechanism of the lithosphere-asthenosphere boundary.

Phlogopite-Olivine Nephelinites Erupted During Early Stage Rifting, North Tanzanian Divergence

The North Tanzanian Divergence (NTD, eastern branch of the East African Rift) corresponds to an early stage of continental breakup. In the southern NTD, two quaternary volcanoes of the Manyara-Balangida rift (Labait, Kwaraha) have erupted primary nephelinite lavas (Mg# = 79–57) that allow characterization of their deep mantle source and the alkaline magmas that percolated through the lithosphere during rift initiation. Nephelinites are olivine- and clinopyroxene-rich, and contain up to 4 vol% magmatic phlogopite that crystallized as a liquidus phase with olivine and clinopyroxene. The presence of hydrous mineral (phlogopite) phenocrysts in Kwaraha and Labait lavas strongly suggests that the alkaline melts were H2O-bearing at the time of phlogopite crystallization (1.57–2.12 wt% H2O in phlogopite), demonstrating that phlogopite may have influenced the partitioning of water between the silicate melt and anhydrous silicate minerals (<1 ppm wt H2O clinopyroxene, 1–6 ppm wt H2O in olivine). Geochemical modeling indicates that the nephelinite magmas resulted from a low degree of partial melting (0.2–1%) of a carbonate-rich (0.3%) garnet peridotite containing ∼2 vol% phlogopite. We estimate the depth of partial melting based on primary melt compositions and empirical relations, and suggest that melting occurred at depths of 110–130 km (4 GPa) for craton-edge lavas (Kwaraha volcano) and 150 km (5 GPa) for on-craton lavas (Labait volcano), close to or below the lithosphere-asthenosphere boundary in agreement with the presence of deep refractory mantle xenoliths in Labait lavas. The depth of melting becomes gradually deeper toward the southern NTD: highly alkaline magmas in the north (Engaruka-Natron Basin) are sourced from amphibole- and CO2-rich peridotite at 75–90 km depth, whereas magmatism in the south (south Manyara Basin) is sourced from deep phlogopite- and CO2-rich garnet-peridotite beneath the Tanzania craton (e.g., at the on-craton Labait volcano). Percolation of deep asthenospheric CO2-rich alkaline magmas during their ascent may have produced strong heterogeneities in the thick sub-continental lithospheric mantle by inducing metasomatism and phlogopite crystallization in glimmerite lithologies.

Post-critical SsPmp and its applications to virtual deep seismic sounding (VDSS)—3: back-projection imaging of the crust–mantle boundary in a heterogeneous lithosphere, theory and application

SUMMARY Virtual deep seismic sounding (VDSS) uses the arrival time of post-critical SsPmp relative to the direct S wave to infer Moho depth at the Pmp reflection point. Due to the large offset between the virtual source and the receiver, SsPmp is more sensitive to lateral variations of structures than near-vertical phases such as Ps, which is used to construct conventional P receiver functions. However, the way post-critical SsPmp is affected by lateral variations in lithospheric structure is not well understood, and previous studies largely assumed a 1-D structure when analysing SsPmp waveforms. Here we present synthetic tests with various 2-D models to show that lateral variations in lithospheric structures, from the lithosphere–asthenosphere boundary (LAB) to sedimentary basins, profoundly affect traveltime, phase and amplitude of post-critical SsPmp, and that a 1-D approximation is usually inappropriate when analysing 2-D data. Despite these strong effects we show, with synthetic examples and the ChinArray data from the Ordos Block in northern China, that a simple ray-theory-based back-projection method can retrieve the geometry of the crust–mantle boundary (CMB) given array observations in cases with moderate lateral variations in the CMB and/or the LAB. The success of our back-projection method indicates that ray-theory approximations are sufficient in modelling SsPmp traveltimes in the presence of moderate lateral heterogeneity. In contrast, we show that the ray theory is generally insufficient in modelling SsPmp phase shifts in a strongly heterogeneous lithosphere due to non-planar downgoing P waves incident at the CMB. Nonetheless, our results demonstrate the feasibility of direct imaging of the CMB with post-critical SsPmp even in the presence of 2-D variations of lithospheric structure.

Thermal and decompression history of the Lanzo Massif, northern Italy: Implications for the thermal structure near the lithosphere–asthenosphere boundary

The lithosphere–asthenosphere boundary (LAB) is a zone where thermal, mechanical, and material interactions take place between the conductive mantle and the underlying convective mantle, and it plays an important role in plate tectonics. In this paper we focus on the thermal aspects of the LAB zone, based on a petrological study of a large peridotite complex that experienced exhumation in the solid state, for the most part. The complex of interest is the Lanzo Massif in the western Alps, northern Italy, where we were able to clarify its thermal and exhumation history and estimate the original thermal structure before exhumation. We examined plagioclase-bearing lherzolites collected from 16 localities, covering the entire massif. All the constituent minerals show compositional heterogeneities on the grain scale. The patterns of Ca, Al, and Cr zoning in pyroxene, the fluorescence-corrected Ca zoning of olivine adjacent to pyroxenes, the Cr and Al zoning of spinel adjacent to plagioclase, and the Ca and Na zoning of plagioclase suggest that an early, nearly isothermal decompression of the Lanzo Massif partly crossing the dry solidus was followed by monotonous cooling through the plagioclase-facies peridotite field. Various deformation microstructures allowed us to specify the timing of the deformation in the framework of the decompression history by carefully observing their relationships with the compositional zoning of the minerals. We show that the deformation took place mainly when effective cooling had started following a period of nearly isothermal exhumation. By applying several geothermometers and evaluating the compositional zoning of the minerals, we were able to quantify the spatial variations in the thermal and decompression history of the Lanzo Massif and constrain the timescales of decompression and cooling. All the estimated temperatures decrease from the southern body towards the northern body. The grain-scale patterns of zoning indicate that the temperatures recorded by the cores of orthopyroxene (1000–1200 °C) indicate a long period of residence in the mantle, whereas those recorded by the rims of pyroxenes and the cores and rims of olivine (600–1100 °C) represent closure temperatures at various times during the decompression. All the closure temperatures decrease from south to north, while the temperature differences between the cores and rims of orthopyroxene increase. This suggests that the cooler and probably shallower northern body cooled at a relatively slow rate than the hotter and probably deeper southern body. The decrease in temperature of ~60 K from south to north, calculated from orthopyroxene cores, may represent the geotherm near the LAB zone. A thermal gradient of ~10 K/km is indicated, which is significantly greater than that estimated for deep subcontinental lithosphere in a steady thermal state. Such a high geotherm might have been caused by thermal perturbation from the underlying hotter asthenospheric mantle.

Inversion of Longer‐Period OBS Waveforms for P Structures in the Oceanic Lithosphere and Asthenosphere

AbstractWe performed waveform inversion of the P waveforms recorded by our BroadBand Ocean Bottom Seismometers (BBOBSs) deployed in the Northwestern Pacific. Consequently, the depth profile of the P velocity of the oceanic upper mantle, which has not been well resolved by previous surface wave or receiver function analyses, was revealed. We considered the azimuthal anisotropy in the lithosphere, which significantly improved the variance reduction from 34% to 44%. The resulting P model exhibited higher and lower velocities in the lithosphere and asthenosphere, respectively. The velocity contrast was found to be second/third of that observed in the previous S models; however the obtained model appeared to have some trade‐off with the VS structures in the vicinity of the source. We compared our P model with the previous S model obtained using our BBOBSs and obtained the VP/VS model. The resulting VP/VS model has two notable features. First, the lithosphere is characterized by a rapid increase in the VP/VS values with depth, which implies chemical stratification. Second, the VP/VS values in the vicinity of the lithosphere‐asthenosphere boundary (LAB) are larger than the synthetic values of any major upper mantle mineral predicted by considering the anharmonic effects, which suggests the effects of anelasticity or melt.

Delineation of Average 1-D Shear Velocity Structure below North India by Surface Wave Dispersion Study

Abstract During 2014-16, a semi-permanent network of four 3-component broadband seismographs was operational in the Rajasthan craton and Aravalli mobile belt in the NW Indian shield. The reliable and accurate broadband data from 16 selected regional Indian earthquakes of Mw5.5-7.8 from this seismic network enabled to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of northern India, covering the region below north India between the entire Himalaya (from the Pakistan Himalaya in west to the Burmese arc in the east) and Rajasthan (including Aravalli mobile belt). First, Rayleigh- (at 7-87s) and Love- (at 7-82s) wave group velocity dispersion curves were measured and then these curves were inverted to estimate the crustal and upper mantle structure below north India. It is observed that group velocities are of variable nature within the region. This could be attributed to the complex crust-mantle structure in the study region resulted from the magmatism episodes associated with the Proterozoic collision, 65 Ma Deccan volcanism and the Himalayan collision. The best model in the study region reveals a two-layered crust, with a 15-km thick upper-crust (UC) of average shear velocity (Vs) of 3.12 km/s and a 25-km thick lower-crust(LC) of average Vs of 3.44 km/sec. The modeling detects a drop in Vs (~1-2%) at 79-120 km depths, underlying north India, representing the probable seismic lithosphere-asthenosphere boundary (LAB) at 79 km depth. A geothermal gradient extrapolated from the surface heat flow (~74 mW/m2) shows that such a gradient would intercept CO2-bearing mantle peridotite solidus at 100 km depth, and thus could signal the presence of small amounts of partially melted magma below 100 km depth. Therefore, this 1-2% drop in Vs could be attributed to the presence of carbonatite melts in the upper mantle related to magmatic episode of 65 Ma Deccan plume activity as also suggested by existing geological and seismological evidence.

Diversity in the peninsular Indian lithosphere revealed from ambient noise and earthquake tomography

Significant diversity in the Indian lithosphere is seen in 3-D shear wave velocity images inverted from the dispersion of group velocities measured from ambient noise and earthquake waveforms along 21,600 paths. Velocity models constructed at 1° intervals to a depth of 200 km suggest a bipolar nature of the Indian lithosphere. We observe a two layer lithospheric mantle beneath the eastern Peninsular India, comprising the Archean cratons of eastern Dharwar, Bastar, Singhbhum, Chotanagpur, Bundelkhand and the Proterozoic Vindyhan basin. The intra-lithospheric mantle boundary is at a depth of ~90 km where Vs increases from 4.5 km/s to >4.7 km/s. The velocity increase continues to a depth of 140 to 180 km followed by its reversal with a minimum Vs of 4.3–4.4 km/s, possibly representing the lithosphere-asthenosphere boundary. Geologically, these regions correlate with the diamond bearing kimberlite fields. In contrast, the Precambrian terrains of western Dharwar, eastern Deccan Volcanic Province (DVP), southern Granulite terrane and the Marwar block in western India are characterized by an almost uniform lithosphere mantle with Vs of 4.4–4.55 km/s. Presence of lower than expected velocity in the lithosphere could be due to its fertile composition, such as clinopyroxene-rich cumulate or metasomatised peridotite. The imprint of Deccan volcanism is observed as the thinned lithosphere (100–120 km) restricted primarily to the westernmost part of DVP, defining the regional extent of Reunion plume-Indian lithosphere interaction. The seismic velocity image suggests evidence for segmentation of Indian lithosphere beneath the Ganga basin at 75° E and 84° E coinciding the Delhi-Haridwar and the Monghyr-Saharsa ridges. Subduction of such a heterogeneous lithosphere could have contributed to the varying style of deformation and subduction pattern in the Himalayan-Tibet region.

The influence of spreading rate and permeability on melt focusing beneath mid-ocean ridges

At mid-ocean ridges, oceanic crust is emplaced in a narrow neovolcanic region on the seafloor, whereas basaltic melt that forms this oceanic crust is generated in a wide region beneath as suggested by a few geophysical surveys. The combined observations suggest that melt generated in a wide region at depths has to be transported horizontally to a small region at the surface. We present results from a suite of two-phase models applied to the mid-ocean ridges, varying half-spreading rate and intrinsic mantle permeability using new openly available models, with the goal of understanding melt focusing beneath mid-ocean ridges and its relevance to the lithosphere-asthenosphere boundary (LAB). Three distinct melt focusing mechanisms are recognized in these models: 1) melting pressure focusing, 2) decompaction layers and 3) ridge suction, of which the first two play dominant roles in focusing melt. All three of these mechanisms exist in the fundamental two phase flow formulation but the manifestation depends largely on the choice of rheological model. The models also show that regardless of spreading rates, the amount of melt and melt transport patterns are sensitive to changes in intrinsic permeability, K0. In these models, the LAB is delineated by the melt-rich decompaction layers, which are essentially defined by the temperature dependent rheological and freezing boundaries. Geophysical observations place the LAB at a steeper incline as compared to the gentler profile suggested by most of our models. The models suggest that one way to reconcile this discrepancy is to have stronger melting pressure focusing mechanism as it is the only mechanism in these models that can focus melt before reaching the typical model thermal LAB. The apparent lack of observable decompaction layers in the geophysical observations hints at the possibility that melting pressure focusing could be significant. These models help improve our understanding of melt focusing beneath mid-ocean ridges and could provide new constraints for mantle rheology and permeability.

New insights into the structural elements of the upper mantle beneath the contiguous United States fromS-to-Pconverted seismic waves

SUMMARYThe S-receiver function (SRF) technique is an effective tool to study seismic discontinuities in the upper mantle such as the mid-lithospheric discontinuity (MLD) and the lithosphere–asthenosphere boundary (LAB). This technique uses deconvolution and aligns traces along the maximum of the deconvolved SV signal. Both of these steps lead to acausal signals, which may cause interference with real signals from below the Moho. Here we go back to the origin of the SRF method and process S-to-P converted waves using S-onset times as the reference time and waveform summation without any filter like deconvolution or bandpass. We apply this ‘causal’ SRF (C-SRF) method to data of the USArray and obtain partially different results in comparison with previous studies using the traditional acausal SRF method. The new method does not confirm the existence of an MLD beneath large regions of the cratonic US. The shallow LAB in the western US is, however, confirmed with the new method. The elimination of the MLD signal below much of the cratonic US reveals lower amplitude but highly significant phases that previously had been overwhelmed by the apparent MLD signals. Along the northern part of the area with data coverage we see relics of Archean or younger northwest directed low-angle subduction below the entire Superior Craton. In the cratonic part of the US we see indications of the cratonic LAB near 200 km depth. In the Gulf Coast of the southern US, we image relics of southeast directed shallow subduction, likely of mid-Palaeozoic age.