Lithosphere-asthenosphere Boundary Research Articles

Crust and upper mantle structure beneath Southwest China and its implications for Mesozoic multistage gold deposits

The Youjiang Basin (YB) on the southwestern margin of the Yangtze Craton (YC) contains >200 Triassic or Cretaceous gold deposits, but deep dynamic processes contributing to large-scale gold accumulation remain unknown. Based on a south–north trending linear broadband seismic array, we obtained the crust and mantle structure beneath the YB and southwestern YC by P-wave and S-wave receiver functions and ambient noise tomography. Our results show a remarkable lithospheric boundary beneath the northern margin of the YB, which separates the YB with felsic thin crust and thin mantle lithosphere from the YC with both a thick crust and thick mantle lithosphere. A new structural image shows that the southern end of the YC Moho dips southward and overlaps with the YB Moho. A 5-km-thick high-velocity layer is also found at a depth of 10 km beneath the northern YB. Combining this information with previous tectonic and geological results, we reveal the process of Mesozoic lithospheric evolution in the southwestern South China Block. We infer that both the pre-Mesozoic weak crust and mantle lithosphere and the intense intracontinental orogenesis in the Early-Middle Triassic contributed to the formation of first-stage Au deposits in the YB. Lithospheric delamination in the Early Cretaceous thinned the lithosphere of YB rather than YC, forming a steep lithosphere–asthenosphere boundary beneath the transition zone between YB and YC, which favored upwelling of Au-rich asthenosphere below the thick cratonic root to the bottom of the crust in the northern YB and contributed to the second-stage Au deposits in the YB.

The lithospheric structure beneath the northeastern Tibetan Plateau and western North China Craton from S receiver function imaging and its implications for lateral growth of the plateau

The mantle lithospheric structures beneath the northeastern Tibetan Plateau (NETP) and adjacent areas provide critical information about the role of the mantle lithosphere in the lateral growth of the plateau. In this study, we construct a comprehensive image of the lithospheric structure across the NETP and western North China Craton (NCC) using the S-wave receiver function method based on the dense seismic array of ChinArray II. We detect one or multiple negative conversions at the depth ranging from 80 to 200 km across most of the study region, which is associated with the seismic lithosphere-asthenosphere boundary (LAB) or mid-lithospheric discontinuity (MLD) within the lithospheric mantle. The main features include that the Beishan Block, West Qinling Orogen, and Qilian Orogen present a relatively shallow LAB at a depth range of 120–140 km, while the Longxi Basin and Alashan Block have the similar lithospheric structure with a deep LAB at depth > 150 km and some MLDs within the lithospheric mantle. We find that the Qilian Orogen in the western NETP overlying a thin mantle lithosphere has been significantly elevated, whereas the Longxi Basin in the eastern NETP with a thick mantle lithosphere lacks both high elevations and significant crustal thickening. It implies that the thick mantle lithosphere beneath the eastern NETP and western NCC likely prevents the plateau from invading the Longxi Basin, while the thin mantle lithosphere across the western NETP and Beishan Block favors outward growth of the Qilian Orogen. Our seismic images provide new insights into the lithospheric heterogeneities and their roles in the lateral growth of the Tibetan Plateau.

Recycling process and proto-kimberlite melt metasomatism in the lithosphere-asthenosphere boundary beneath the Amazonian Craton recorded by garnet xenocrysts and mantle xenoliths from the Carolina kimberlite

Here we present new data on the major and trace element compositions of silicate and oxide minerals from mantle xenoliths brought to the surface by the Carolina kimberlite, Pimenta Bueno Kimberlitic Field, which is located on the southwestern border of the Amazonian Craton. We also present Sr-Nd isotopic data of garnet xenocrysts and whole-rocks from the Carolina kimberlite. Mantle xenoliths are mainly clinopyroxenites and garnetites. Some of the clinopyroxenites were classified as GPP–PP–PKP (garnet-phlogopite peridotite, phlogopite-peridotite, phlogopite-K-richterite peridotite) suites, and two clinopyroxenites (eclogites) and two garnetites are relicts of an ancient subducted slab. Temperature and pressure estimates yield 855–1102 °C and 3.6–7.0 GPa, respectively. Clinopyroxenes are enriched in light rare earth elements (LREE) (LaN/YbN = 5–62; CeN/SmN = 1–3; where N = primitive mantle normalized values), they have high Ca/Al ratios (10–410), low to medium Ti/Eu ratios (742–2840), and low Zr/Hf ratios (13–26), which suggest they were formed by metasomatic reactions with CO2-rich silicate melts. Phlogopite with high TiO2 (>2.0 wt.%), Al2O3 (>12.0 wt.%), and FeOt (5.0–13.0 wt.%) resemble those found in the groundmass of kimberlites, lamproites and lamprophyres. Conversely, phlogopite with low TiO2 (<1.0 wt.%) and lower Al2O3 (<12.0 wt.%) are similar to those present in GPP-PP-PKP, and in MARID (mica-amphibole-rutile-ilmenite-diopside) and PIC (phlogopite-ilmenite-clinopyorxene) xenoliths. The GPP-PP-PKP suite of xenoliths, together with the clinopyroxene and phlogopite major and trace element signatures suggests that an intense proto-kimberlite melt metasomatism occurred in the deep cratonic lithosphere beneath the Amazonian Craton. The Sr-Nd isotopic ratios of pyrope xenocrysts (G3, G9 and G11) from the Carolina kimberlite are characterized by high 143Nd/144Nd (0.51287–0.51371) and εNd (+4.55 to +20.85) accompanied with enriched 87Sr/86Sr (0.70405–0.71098). These results suggest interaction with a proto-kimberlite melt compositionally similar with worldwide kimberlites. Based on Sr-Nd whole-rock compositions, the Carolina kimberlite has affinity with Group 1 kimberlites. The Sm-Nd isochron age calculated with selected eclogitic garnets yielded an age of 291.9 ± 5.4 Ma (2 σ), which represents the cooling age after the proto-kimberlite melt metasomatism. Therefore, we propose that the lithospheric mantle beneath the Amazonian Craton records the Paleozoic subduction with the attachment of an eclogitic slab into the cratonic mantle (garnetites and eclogites); with a later metasomatic event caused by proto-kimberlite melts shortly before the Carolina kimberlite erupted.

Evidence for a lithospheric step and pervasive lithospheric thinning beneath southern New England, northeastern USA

AbstractIn this study, we use data from the SEISConn seismic experiment to calculate Sp receiver functions in order to characterize the geometry of upper-mantle structure beneath southern New England (northeastern United States). We image robust negative-velocity-gradient discontinuities beneath southern New England that we interpret as corresponding to the lithosphere-asthenosphere boundary (LAB) and identify a well-defined step of 15 km in LAB depth at a longitude of 73°W, which we interpret to be the boundary between Laurentian and Appalachian lithosphere, although the offset may be larger if the putative LAB phase is reinterpreted to be a mid-lithospheric discontinuity. We infer that the lithosphere throughout the region is substantially thinner than elsewhere in the continental interior, consistent with regional tomographic studies and previously published Sp receiver function results. The presence of thinned lithosphere suggests that the low-velocity Northern Appalachian Anomaly (NAA) in the upper mantle may extend as far south as coastal Connecticut. The presence of regionally thinned lithosphere and a step in lithospheric thickness suggests that inherited structure may be preserved in present-day lithosphere, even in the presence of more recent dynamic processes associated with the NAA.

Off-rift Axis Channelized Melt and Lithospheric Metasomatism along Mid-ocean Ridges—A Case Study from Iceland on the Limits of Melt Channelling

AbstractMagmatism in Iceland is classically explained by the interaction of the mid-Atlantic ridge with the Iceland plume. The growth of Iceland through time is the result of volcanic activity at the rift axis. However, Holocene volcanism (0–11.5 ka) is not restricted to the rift zone (RZ) but also occurs off-axis, specifically in the western Snæfellsnes Volcanic Belt (SVB) and in the Southern Flank Zone (SFZ). The rift and off-axis postglacial volcanic zones are separated by a gap of ~60–80 km. While the volcanic activity of the SFZ seems correlated with the actual location of the Iceland plume, the plume relation to the SVB is uncertain. To address the origin and relationships between rift and off-rift magmas, we present new data from two transects perpendicular to the Reykjanes rift. The lavas in the SVB are characterized by transitional to alkaline compositions, with elevated incompatible trace element content. In contrast, the RZ volcanic rocks have tholeiitic compositions with trace element signatures slightly more enriched than Mid-Ocean Ridge Basalt (MORB). Rift and off-rift Iceland lavas are all characterized by positive Ba and Nb anomalies, particularly in alkaline rocks. Tholeiitic and alkaline lavas show distinct differentiation sequences, with the main difference being the delayed crystallization of plagioclase in the fractionating assemblage of alkaline magmas. We apply these sequences to calculate primary magma compositions, which are then used to constrain melting conditions. Geochemical modelling indicates that Iceland rift and off-rift magmas can be produced from a peridotitic mantle source if lithospheric processes are involved. We demonstrate that recycled crust in the form of pyroxenite is not required to generate Snæfellsnes alkaline lavas. The low solidus temperature and high productivity of pyroxenite favour early and more extensive melting producing primary magmas that are not sufficiently enriched in incompatible trace elements to explain the compositional variation of Snæfellsnes magmas. An alternative mechanism to involve pyroxenite in the source of Snæfellsnes lavas relates to the hybridization of recycled oceanic crust with peridotite, but such reacted pyroxenite requires specific compositions to reproduce the Ba and Nb anomalies. As an alternative, we suggest that Snæfellsnes alkaline lavas are the result of channelized low-degree melts produced on the periphery of the melting column at distances exceeding 65 km from the ridge axis. These melts accumulate and percolate into the lithosphere producing amphibole ± phlogopite-bearing hydrous cumulates. Incongruent melting of these cumulates via renewed magmatic activity and melt-peridotite reaction can reproduce the alkaline compositions observed in the SVB, including the Ba and Nb anomalies. Numerical simulations of melt extraction below mid-ocean ridges suggest that low-degree melts produced as far as ~65 km from the central ridge axis rise vertically to the base of the lithosphere and are then focussed towards the ridge axis in decompaction channels. We propose that these melts interact with hydrous cumulates previously formed during the development of decompaction channels at the lithosphere–asthenosphere boundary and acquire specific Ba and Nb anomalies. The mixing of these distal enriched melts with more depleted melts extracted from the central part of the melting regime explains the composition of RZ lavas. Alkaline lavas observed in the SFZ show numerous analogies with the Snæfellsnes magmas, suggesting that similar lithospheric processes control their chemistry. The high thermal regime in Iceland and thick lithosphere explains the difference between Icelandic tholeiite and typical MORB. Our results highlight the importance of mantle dynamics below mid-ocean ridges and lithospheric interaction to produce off-axis magmatism with enriched alkaline signatures.

Discovery of distinct lithosphere-asthenosphere boundary and the Gutenberg discontinuity in the Atlantic Ocean.

The plate tectonic theory requires a rigid lithosphere floating over a weak asthenosphere, separated by the lithosphere-asthenosphere boundary, which has been sometimes interpreted as the Gutenberg discontinuity. Using a deep seismic reflection technique, we report the presence of two continuous reflections covering 27 Ma to 58 Ma oceanic lithosphere in the Atlantic Ocean. We find that the upper reflection deepens with age and follows the ~1250°C isotherm, whereas the deeper reflection lies at a constant depth of ~75 km. We suggest that the upper reflection represents the thermally controlled lithosphere-asthenosphere boundary, whereas the lower reflection is the Gutenberg discontinuity, a frozen-in dehydration boundary separating the dry mantle melting region above from the hydrated mantle below formed at the ridge axis. We also find that thermal mantle anomalies rejuvenate the lithosphere, uplift the lithosphere-asthenosphere boundary, and destroy the Gutenberg discontinuity.

Geoscience Poster G16: Potential field interpretation and modelling in the deep-water Otway Basin

Poster G16 The NW–SE trending Otway Basin is a well-established hydrocarbon province with onshore and shallow-water offshore discoveries. However, the deep-water offshore is underexplored due to the sparsity of data in the past. New 2D seismic, gravity and magnetic data, acquired during the 2020 Otway Basin Seismic Program, were combined with legacy seismic data to improve the understanding of this part of the basin. The new data have been processed and integrated with existing open-file data to derive Moho and lithosphere–asthenosphere boundary surfaces. Additional wavelength filtering of the newly acquired magnetic and gravity data, using constraints from seismic data, provides valuable information on the geometry and the spatial extent of igneous rocks in the deep-water basin. While extrusive igneous rocks and related intrusive rocks in the Otway Basin can be inferred from seismic reflection data, the new magnetic and gravity data provide additional constraints that help to map their distribution across the basin. To access the poster click the link on the right. To read the full paper click here

Three-dimensional seismic anisotropy in the Pacific upper mantle from inversion of a surface-wave dispersion data set

SUMMARYWe present a new, 3-D model of seismic velocity and anisotropy in the Pacific upper mantle, PAC13E. We invert a data set of single-station surface-wave phase-anomaly measurements sensitive only to Pacific structure for the full set of 13 anisotropic parameters that describe surface-wave anisotropy. Realistic scaling relationships for surface-wave azimuthal anisotropy are calculated from petrological information about the oceanic upper mantle and are used to help constrain the model. The strong age dependence in the oceanic velocities associated with plate cooling is also used as a priori information to constrain the model. We find strong radial anisotropy with vSH &amp;gt; vSV in the upper mantle; the signal peaks at depths of 100–160 km. We observe an age dependence in the depth of peak anisotropy and the thickness of the anisotropic layer, which both increase with seafloor age, but see little age dependence in the depth to the top of the radially anisotropic layer. We also find strong azimuthal anisotropy, which typically peaks in the asthenosphere. The azimuthal anisotropy at asthenospheric depths aligns better with absolute-plate-motion directions while the anisotropy within the lithosphere aligns better with palaeospreading directions. The relative strengths of radial and azimuthal anisotropy are consistent with A-type olivine fabric. Our findings are generally consistent with an explanation in which corner flow at the ridge leads to the development and freezing-in of anisotropy in the lithosphere, and shear between the lithosphere and underlying asthenosphere leads to anisotropy beneath the plate. We also observe large regions within the Pacific basin where the orientation of anisotropy and the absolute-plate-motion direction differ; this disagreement suggests the presence of shear in the asthenosphere that is not aligned with absolute-plate-motion directions. Azimuthal-anisotropy orientation rotates with depth; the depth of the maximum vertical gradient in the fast-axis orientation tends to be age dependent and agrees well with a thermally controlled lithosphere–asthenosphere boundary. We observe that azimuthal-anisotropy strength at shallow depths depends on half-spreading rate, with higher spreading rates associated with stronger anisotropy. Our model implies that corner flow is more efficient at aligning olivine to form lattice-preferred orientation anisotropy fabrics in the asthenosphere when the spreading rate at the ridge is higher.

Potential field interpretation and modelling in the deep-water Otway Basin

The NW–SE trending Otway Basin is a well-established hydrocarbon province with onshore and shallow-water offshore discoveries. However, the deep-water offshore is underexplored due to the sparsity of data in the past. New 2D seismic, gravity and magnetic data, acquired during the 2020 Otway Basin Seismic Program, were combined with legacy seismic data to improve the understanding of this part of the basin. The new data have been processed and integrated with existing open-file data to derive Moho and lithosphere–asthenosphere boundary surfaces. Additional wavelength filtering of the newly acquired magnetic and gravity data, using constraints from seismic data, provides valuable information on the geometry and the spatial extent of igneous rocks in the deep-water basin. While extrusive igneous rocks and related intrusive rocks in the Otway Basin can be inferred from seismic reflection data, the new magnetic and gravity data provide additional constraints that help to map their distribution across the basin.

Thermal Modification of the Northwest Indian Shield: As Evidenced by Integrated Geopotential Modelling

Abstract The architecture of the continental lithospheric mantle is vital for understanding the wide-ranging processes that govern the creation of Earth’s continents. The present investigation delineates the crustal and lithospheric boundaries underneath northwest India, which provides essential constraints as an imprint of tectonomagmatic events under the distinct tectonic setting. This study integrates three data sets (gravity, geoid, and topography data) and thermal parameters in a three-dimensional inversion algorithm to northwest India. Inversion results are constrained by available passive seismic crustal thickness, which we are using as a priori information. Our modelling results show crust down to a depth of 33 to 39 km in the Saurashtra and Kachchh region, which increases gradually to about 40 km beneath the Jaisalmer region. The deepest Moho at about 42-45 km is found under the Aravalli-Delhi fold belt and Indo-Gangetic plain. A comparable crustal thickness of 41-44 km is observed beneath the Bundelkhand craton. The mean crustal density varies laterally between 2850 and 2890 kg/m-3 in the study area. The lithosphere-asthenosphere boundary (LAB) is 115-145 km deep under Saurashtra and Kachchh regions, which remains comparable under the Barmer-Pokharan region (130-160 km). A moderate LAB depth under the Aravalli-Delhi fold belt (140-180 km) and relatively deeper under the Bundelkhand craton (175 to 185 km) reaches over 200 km below the Indo-Gangetic plain. The thinned lithosphere beneath the Kachchh, Saurashtra, and Barmer-Pokhran regions, crisscrossed by two different rift systems, together with the dominance of the volcanic centres of two major magmatic episodes, and the possible presence of a low-velocity zone in the subcrustal upper mantle region suggest a genetic link between the evolution and the multistage thermal modification of the mantle lithosphere of this region.

A plume origin for hydrous melt at the lithosphere-asthenosphere boundary.

Plate tectonics requires a low-viscosity layer beneath the lithosphere-asthenosphere boundary (LAB), yet the origin of this ductile transition remains debated1,2. Explanations include the weakening effects of increasing temperature3,4, mineral hydration5 or partial melt6. Electrical resistivity is sensitive to all three effects7, including melt volatile content8, but previous LAB constraints from magnetotelluric soundings did not simultaneously consider the thermodynamic stability of the inferred amount of melt and the effect of uncertainty in the estimated resistivity8-14. Here we couple an experimentally constrained parameterization of mantle melting in the presence of volatiles15,16 with Bayesian resistivity inversion17 and apply this to magnetotelluric data sensitive to a LAB channel beneath the Cocos Plate9. Paradoxically, we find that the conductive channel requires either anomalously large melt fractions with moderate volatile contents or moderate melt fractions with anomalously large volatile contents, depending on the assumed mantle temperature. Large melt fractions are unlikely to be mechanically stable and conflict with melt-migration models18. As large volatile contents require a highly enriched mantle source inconsistent with mid-ocean-ridge estimates19, our results indicate that a mantle plume emplaced volatile-rich melts in the LAB channel. This requires the presence of a previously undetected nearby plume or the influence of the distant Galápagos hotspot. Plumes that feed thin, hydrous melt channels9,14,20 may be an unrecognized source of LAB anomalies globally.

Distinct Lithospheric Structure in the Xing'an‐Mongolian Orogenic Belt

AbstractDetailed knowledge of the lithospheric structure is essential for a better understanding of the tectonic evolution of the Xing'an‐Mongolian Orogenic Belt (XMOB). We constrained the lithospheric structure across the XMOB by receiver function imaging and shear wave splitting analysis from a dense seismic array. The lithosphere‐asthenosphere boundary (LAB) is coherently imaged in different tectonic blocks. The mid‐lithospheric discontinuity is also identified at ∼75–100 km depth over the deeper LAB (∼110–130 km) beneath the western side of the North‐South Gravity Lineament (NSGL), roughly at the same depth as the LAB beneath the eastern side of the NSGL. The distinct variations in both the lithospheric structure and fast shear‐wave polarization direction across the NSGL are likely caused by different tectonic processes in the two sides, indicating that the NSGL might represent the western boundary influenced by the subduction of the (Paleo‐)Pacific Plate. Lithosphere thinning in the XMOB is evidenced to have been limited mainly to the east of the NSGL.

Nature of the lithosphere-asthenosphere boundary beneath the Eastern Dharwar Craton of the Indian Shield

• Thickness and nature of lithosphere-asthenosphere boundary has been discussed. • Partial melts may explain the shallow and sharp lithosphere-asthenosphere boundary. • Plume induced lithospheric thinning argued beneath the Eastern Dharwar Craton . • The Eastern Dharwar Craton lithosphere has been substantially modified. The lithosphere-asthenosphere boundary (LAB) is a fundamental element of the plate tectonic hypothesis that accommodates the differential motion of rigid lithosphere over the weaker asthenosphere . In recent times, various usages have been used to define the LAB, depending on the nature of their measurements. Here, we investigate the lithospheric structure beneath the Eastern Dharwar Craton (EDC) of the Indian Shield using geochemical, thermal and seismological data sets. We analysed S-receiver functions from the stations deployed in the EDC along with the surface wave dispersion tomography. We also added thermal measurements from 5 different locations and geochemical data from 34 Kimberlite/Lamproite xenolith samples to constrain the nature of the LAB. The seismological measurements using Rayleigh wave dispersion and receiver function analysis indicate the lithospheric thickness of 98–118 and 94–118 km respectively, with sharp transition across the LAB. The P-T results from xenoliths are interpreted in concurrence with the heat-flow measurements suggesting a thick thermal lithosphere of ∼ 200 km for the normal mantle solidus with cold geotherm. To reconcile our observations, we invoke partial melts or melts enriched in volatiles, which significantly lowers the viscosity of mantle rocks inducing a zone of weakness between the rigid lithosphere (∼125 km) and the convective asthenosphere. Further, we favour the view that the thick lithosphere beneath the Indian plate has been thinned by a plume during the Gondwanaland breakup at ∼ 130 Ma. The presence of younger kimberlites from the Indian shield support that it is further degenerated by the delamination leading to an uneven topography in the LAB.

Three-dimensional crustal and upper-mantle resistivity structure of Alberta, Canada: implications for Precambrian tectonics

SUMMARY This paper describes the first detailed 3-D study of the Alberta basement using the magnetotelluric (MT) method. Long-period MT data were used to generate a 3-D electrical resistivity model of the crust and upper mantle beneath Alberta. The Western Canada Sedimentary Basin was imaged as a low-resistivity layer. A number of crustal and upper-mantle conductors were imaged and include: (1) the Red Deer conductor, (2) the Kiskatinaw conductor, (3) a conductor coincident with the Kimiwan oxygen isotope anomaly, (4) the southern Alberta–British Columbia conductor, (5) the Loverna block conductor and (6) a conductor beneath the Birch Mountain kimberlite field. Conductors (3) and (6) are newly identified in this study. All of these conductors can be related to past tectonic events associated with the assembly of Laurentia and subsequent kimberlite magmatism. The depth of the lithosphere–asthenosphere boundary across Alberta varies from 150 to 300 km as defined from the 100 Ω m contour. Furthermore, the Snowbird tectonic zone in Alberta is characterized by a thick, resistive lithosphere. Similarly, the diamondiferous Buffalo Head Hills kimberlite field is imaged with a thick, resistive lithosphere, both may indicate depletion. In contrast, the Birch Mountain kimberlite field is underlain by a pronounced conductor that extends through the lithosphere and may be the result of metasomatism.

The removal of multiple reflection waves in the P receiver function through parabolic Radon transformation

SUMMARY Receiver function (RF) analysis is an indispensable method that is widely utilized to image the deep structure of the crust and upper mantle. Although the P-wave receiver function (PRF) has the advantages of a higher signal-to-noise ratio and higher resolution than S-wave receiver function, the information about P-to-S converted (Ps) phases from velocity discontinuities of the crust and upper mantle is usually obstructed by multiples from shallower crustal interfaces. Sometimes it is difficult to identify useful information about lithospheric discontinuities. In this paper, we propose a method for suppressing multiples and isolating the Ps phases of RFs. This task is accomplished by exploiting the differences in arrival times resulting from the slowness between the conversion and the multiples that occur at the same interface in the PRFs and by separating them in the Radon domain, which enables removal of crustal and sedimentary multiples. This method can effectively remove multiples and isolate the useful signals of Ps phases contained in the PRFs. We test the method on synthetic PRFs and demonstrate that crustal multiples can be effectively eliminated or suppressed and that the Ps phase of lithospheric discontinuities can be coherently traced. Next, we apply the method to real PRF data collected from the seismic station of the China Array seismic experiment (Phase II and Phase III) in the Ordos block and its adjacent area, and successfully obtain an improved common conversion point stacking image of lithospheric discontinuity structures in the depth domain, particularly for detecting the lithosphere–asthenosphere boundary.

New Insights Into Lithospheric Structure and Melting Beneath the Colorado Plateau

AbstractThe Colorado Plateau and its surroundings serve as an archetypal case to investigate the interaction of mantle melting processes and lithospheric structure. It has been hypothesized that widespread Cenozoic volcanism indicates the encroachment of the convective upwelling of asthenosphere toward the Plateau center. In this study, we generate a Common Conversion Point (CCP) stack of S‐to‐p (Sp) receiver functions to image the locations of lithospheric discontinuities in the southwestern United States. Our results are broadly similar to prior work, showing a strong and continuous Negative Velocity Gradient (NVG) consistent with the Lithosphere‐Asthenosphere Boundary (LAB) over much of the study area. However, with several methodological improvements, we are able to obtain more reliable NVG depth picks below the Colorado Plateau where the LAB becomes weaker, deeper, and broader. We compare the inferred topography of NVGs with the locations of volcanoes, and find that the majority of recent volcanoes are co‐located with lithosphere that is ∼80 km thick. This appears to be the critical depth at which partial melt from upwelling asthenosphere pooling at the base of (or within) the lithosphere may percolate to the surface. We compare our CCP profiles with magma equilibration conditions determined from petrologic analysis and find good agreement between the depth of NVGs and depth of magma equilibration. This analysis provides insight into the progression of magmatism and lithospheric loss toward the center of the Colorado Plateau, and demonstrates how small‐scale processes like melting influence lithosphere‐asthenosphere interactions that persist over large temporal and spatial scales.

Post‐Subduction Tectonics of Sabah, Northern Borneo, Inferred From Surface Wave Tomography

AbstractWe use two‐plane‐wave tomography with a dense network of seismic stations across Sabah, northern Borneo, to image the shear wave velocity structure of the crust and upper mantle. Our model is used to estimate crustal thickness and the depth of the lithosphere‐asthenosphere boundary (LAB) beneath the region. Calculated crustal thickness ranges between 25 and 55 km and suggests extension in a NW‐SE direction, presumably due to back‐arc processes associated with subduction of the Celebes Sea. We estimate the β‐factor to be 1.3–2, well below the initiation of seafloor spreading. The LAB is, on average, at a depth of 100 km, which is inconsistent with models that ascribe Neogene uplift to wholescale removal of the mantle lithosphere. Instead, beneath a region of Plio‐Pleistocene volcanism in the southeast, we image a region 50–100 km across where the lithosphere has thinned to &lt;50 km, supporting recent suggestions of lower lithospheric removal through a Rayleigh‐Taylor instability.

Adjoint Waveform Tomography of South America

AbstractWe used 3D spectral‐element seismic wave simulations and data from 112 earthquakes and 1,311 seismic stations, totalizing 20,884 unique ray paths, to construct an adjoint waveform tomographic model of South America. We performed 23 conjugate‐gradient iterations using exponentiated phase (EP) measurements. Our final model (SAAM23, South American Adjoint Model—iteration 23) shows a 50% decrease in the EP misfit relative to its 3D starting model. We further assessed the phase misfit reduction by using cross‐correlation travel‐time measurements of 53 earthquakes not included in the inversion. We estimated SAAM23 resolution using point‐spread function tests and density coverage analysis. The Nazca Slab is well imaged and is shown to be continuous in the 300–500 km depth range. Beneath northern South America, the slab traverses the mantle transition zone and continues into the lower mantle. In the central and southern part of South America, the slab appears to flatten near the 650 km discontinuity before continuing into the lower mantle. In the stable Precambrian platform, both cratons (Amazonian and São Francisco), as well as covered cratonic blocks beneath the intracratonic Paraná and Parnaíba basins (Paranapanema and Parnaíba, respectively), show high velocities at lithospheric depths. The seismic Lithosphere/Asthenosphere boundary (LAB) agrees well with published values obtained by S‐wave receiver functions. In the Amazonian craton, the positive lithospheric S‐wave velocity anomalies and LAB depth increase with the average age of the geochronological provinces. No lithospheric high‐velocity anomalies were found beneath the Río de la Plata Craton.

Sublithosphere Mantle Crystallization and Immiscible Sulfide Melt Segregation in Continental Basalt Magmatism: Evidence from Clinopyroxene Megacrysts in the Cenozoic Basalts of Eastern China

AbstractThis study explores the effects of high-pressure crystallization and immiscible sulfide melt segregation under mantle conditions on the compositional variation of basaltic magmas, using clinopyroxene megacrysts in the Cenozoic basalts of eastern China. These clinopyroxene megacrysts are large (up to &amp;gt;10 cm in size) and homogeneous at the grain scale. They were crystallized from variably evolved parental magmas and then captured by their host basalts. The large and systematic variations of [Sm/Yb]N, Lu/Hf, Fe/Mn, Sc/La, Ni and Cu with Mg# in the clinopyroxene megacrysts suggest their co-precipitation with garnet and with immiscibility between sulfide and silicate melts. This is consistent with the appearance of garnet megacrysts in the host basalts and abundant sulfide globules in the clinopyroxene megacrysts. The covariation between Ni contents of sulfide globules and Mg# of the clinopyroxene megacrysts suggests a genetic relationship between sulfide globules and clinopyroxene megacrysts. High-pressure crystallization of clinopyroxene and garnet results in decrease of Mg# and concentrations of CaO, MnO and heavy rare earth elements (e.g. Yb) and increase of Fe/Mn and [Sm/Yb]N in the residual melts. Therefore, geochemical characteristics of low Mg#, low CaO and MnO contents and high Fe/Mn and [Sm/Yb]N in basalts do not necessarily indicate a pyroxenite mantle source. In addition, caution is needed when applying the olivine addition method to infer the primary compositions of alkali basalts without considering the effects of high-pressure crystallization of clinopyroxene and garnet. The calculated pressure (P) and temperature (T) conditions of the clinopyroxene megacrysts are close to those of the lithosphere–asthenosphere boundary (LAB) beneath eastern China, and the low primitive [Sm/Yb]N (~4.0) of melts parental to the clinopyroxene megacrysts suggests final equilibration at relatively low pressures most likely beneath the LAB. Hence, a melt-rich layer is expected close beneath the LAB. Melt pools in this melt-rich layer provide a stable and closed environment for the growth of compositionally homogeneous clinopyroxene megacrysts. As a result, melts in these melt pools are compositionally evolved with low and variable Mg#. Subsequent pulses of melt aggregation/supply from depths with primitive compositions and high Mg# will disturb these melt pools, cause magma mixing and trigger the eruption of magmas carrying clinopyroxene and garnet megacrysts.

The thermal structure of subduction zones predicted by plate cooling models with variable thermal properties

SUMMARY Previous modelling studies have investigated the effects of experimentally constrained thermal properties (i.e. thermal conductivity, specific heat and the thermal expansion coefficient) on the thermal structure of subduction zones. However, these studies have not carefully considered whether the assumed thermal structure of the slab before subduction is consistent with geophysical observations. This study investigates the effects of thermal properties on the thermal structure of the Tohoku subduction zone, northeast Japan, by using the slab temperature at the trench determined from plate cooling models. Three types of thermal properties were tested: constant, temperature-dependent and temperature- and lithology-dependent types. For each case, the parameters for the plate cooling models were inferred based on the observed surface heat flow and seafloor depth using Bayes’ theorem. It was found that the predicted temperature and location of phase boundaries in the slab, which are possibly related to intermediate-depth earthquakes, are similar for the three cases. This suggests that, in the Tohoku subduction zone, constant thermal properties can be used in modelling to examine phenomena related to slab dehydration. The depth uncertainties for isotherms in the oceanic plate and slab increase with temperature, and are about ±10 and ±20 km for the 600 and 1200 °C isotherms, respectively. When this uncertainty is considered, the location of the serpentinite-out boundary matches that of the lower plane of double seismic zone, suggesting that dehydration may be important in triggering intermediate-depth seismicity. However, the large uncertainty makes it difficult to discuss in detail the origins of intraplate earthquakes, the lithosphere–asthenosphere boundary, and the lower boundary of the slab in terms of temperature.