Bending-related faulting of the incoming oceanic plate and its effect on lithospheric hydration and seismicity: A passive and active seismological study offshore Maule,Chile

We have studied the dependency between incoming plate structure, bending-related faulting, lithospheric hydration, and outer rise seismic activity offshore Maule, Chile. We derived a 2D Poisson's ratio distribution from P- and S-wave seismic wide angle data collected in the trench-outer rise. High values of Poisson's ratio in the uppermost mantle suggest that the oceanic lithosphere is highly hydrated due to the water infiltration through bending-related normal faults outcropping at the seafloor. This process is presumably facilitated by the presence of a seamount in the area. We conclude that water infiltrates deep into the lithosphere, when it approaches the Chile trench, producing a reduction of crustal and upper mantle velocities, supporting serpentinization of the upper mantle. Further, we observed a mantle Vp anisotropy of 8%, with the fast velocity axis running normal to the abyssal hill fabric and hence in spreading direction, indicating that outer rise processes have yet not affected anisotropy.The first weeks following the megatrust Mw = 8.8 Maule earthquake in 2010 were characterized by a sudden increase of the outer rise seismic activity, located between 34° S and 35°30′ S. We concluded that this phenomenon is a result of an intensification of the water infiltration process in the outer rise, presumably triggered by the main shock, whose epicenter was located some 100 km to the south east of the cluster.     

Ancient Os isotope signatures from the Ontong Java Plateau lithosphere: Tracing lithospheric accretion history

In order to better understand the nature and formation of oceanic lithosphere beneath the Early Cretaceous Ontong Java Plateau, Re–Os isotopes have been analysed in a suite of peridotite xenoliths from Malaita, Solomon Islands. Geological, thermobarometric and petrological evidence from previous studies reveal that the xenoliths represent virtually the entire thickness of the southern part of subplateau lithospheric mantle (< 120 km). This study demonstrates that vertical Os isotopic variations correlate with compositional variations in a stratified lithosphere. The shallowest plateau lithosphere (< 85 km) is dominated by fertile lherzolites showing a restricted range of ¹⁸⁷Os/¹⁸⁸Os (0.1222 to 0.1288), consistent with an origin from ~ 160 Ma Pacific lithosphere. In contrast, the basal section of subplateau lithospheric mantle (~ 95–120 km) is enriched in refractory harzburgites with highly unradiogenic ¹⁸⁷Os/¹⁸⁸Os ratios ranging from 0.1152 to 0.1196, which yield Proterozoic model ages of 0.9–1.7 Ga. Although the whole range of Os isotope compositions of Malaita peridotites is within the variations seen in modern abyssal peridotites, the contrasting isotopic compositions of shallow and deep plateau lithosphere suggest their derivation from different mantle reservoirs. We propose that the subplateau lithosphere forms a genetically unrelated two-layered structure, comprising shallower, typical oceanic lithosphere underpinned by deeper impinged material, which included a component of recycled Proterozoic lithosphere. The impingement of residual but chemically heterogeneous mantle, mechanically coupled to the recently formed, thin lithosphere, may have a bearing on the anomalous initial uplift and late subsidence history of the seismically anomalous plateau root.     

Crustal intrusion beneath the Louisville hotspot track

We report here the first detailed 2D tomographic image of the crust and upper mantle structure of a Cretaceous seamount that formed during the interaction of the Pacific plate and the Louisville hotspot. Results show that at ～ 1.5 km beneath the seamount summit, the core of the volcanic edifice appears to be dominantly intrusive, with velocities faster than 6.5 km/s. The edifice overlies both high lower crustal (> 7.2–7.6 km/s) and upper mantle (> 8.3 km/s) velocities, suggesting that ultramafic rocks have been intruded as sills rather than underplated beneath the crust. The results suggest that the ratio between the volume of intra-crustal magmatic intrusion and extrusive volcanism is as high as ～ 4.5. In addition, the inversion of Moho reflections shows that the Pacific oceanic crust has been flexed downward by up to ～ 2.5 km beneath the seamount. The flexure can be explained by an elastic plate model in which the seamount emplaced upon oceanic lithosphere that was ～ 10 Myr at the time of loading. Intra-crustal magmatic intrusion may be a feature of hotspot volcanism at young, hot, oceanic lithosphere, whereas, magmatic underplating below a pre-existing Moho may be more likely to occur where a hotspot interacts with oceanic lithosphere that is several tens of millions of years old.     

Miocene to Late Quaternary Patagonian basalts (46–47°S): Geochronometric and geochemical evidence for slab tearing due to active spreading ridge subduction

Miocene to Quaternary large basaltic plateaus occur in the back-arc domain of the Andean chain in Patagonia. They are thought to result from the ascent of subslab asthenospheric magmas through slab windows generated from subducted segments of the South Chile Ridge (SCR). We have investigated three volcanic centres from the Lago General Carrera–Buenos Aires area (46–47°S) located above the inferred position of the slab window corresponding to a segment subducted 6 Ma ago. (1) The Quaternary Río Murta transitional basalts display major, trace elements, and Sr and Nd isotopic features similar to those of oceanic basalts from the SCR and from the Chile Triple Junction near Taitao Peninsula (e.g., (⁸⁷Sr/⁸⁶Sr)_o = 0.70396–0.70346 and εNd = + 5.5 − + 3.0). We consider them as derived from the melting of a Chile Ridge asthenospheric mantle source containing a weak subduction component. (2) The Plio-Quaternary (< 3.3 Ma) post-plateau basanites from Meseta del Lago Buenos Aires (MLBA), Argentina, likely derive from small degrees of melting of OIB-type mantle sources involving the subslab asthenosphere and the enriched subcontinental lithospheric mantle. (3) The main plateau basaltic volcanism in this region is represented by the 12.4–3.3-Ma-old MLBA basalts and the 8.2–4.4-Ma-old basalts from Meseta Chile Chico (MCC), Chile. Two groups can be distinguished among these main plateau basalts. The first group includes alkali basalts and trachybasalts displaying typical OIB signatures and thought to derive from predominantly asthenospheric mantle sources similar to those of the post-plateau MLBA basalts, but through slightly larger degrees of melting. The second one, although still dominantly alkalic, displays incompatible element signatures intermediate between those of OIB and arc magmas (e.g., La/Nb > 1 and TiO₂ < 2 wt.%). These intermediate basalts differ from their strictly alkalic equivalents by having lower High Field Strength Element (HFSE) and higher εNd (up to + 5.4). These features are consistent with their derivation from an enriched mantle source contaminated by ca. 10% rutile-bearing restite of altered oceanic crust. The petrogenesis of the studied Mio-Pliocene basalts from MLBA and MCC is consistent with contributions of the subslab asthenosphere, the South American subcontinental lithospheric mantle and the subducted Pacific oceanic crust to their sources. However, their chronology of emplacement is not consistent with an ascent through an asthenospheric window opened as a consequence of the subduction of segment SCR-1, which entered the trench at 6 Ma. Indeed, magmatic activity was already important between 12 and 8 Ma in MLBA and MCC as well as in southernmost plateaus, i.e., 6 Ma before the subduction of the SCR-1 segment. We propose a geodynamic model in which OIB and intermediate magmas derived from deep subslab asthenospheric mantle did uprise through a tear-in-the-slab, which formed when the southernmost segments of the SCR collided with the Chile Trench around 15 Ma. During their ascent, they interacted with the Patagonian supraslab mantle and, locally, with slivers of subducted Pacific oceanic crust that contributed to the geochemical signature of the intermediate basalts.     

The thermal state and strength of the lithosphere in the Spanish Central System and Tajo Basin from crustal heat production and thermal isostasy

In this work we have modeled the thermal structure of the lithosphere of the Spanish Central System and the Tajo Basin, and their implications for lithospheric strength. For this, we have used refined heat-producing elements (HPE) values to obtain new estimates of heat production rates in the Spanish Central System and Tajo Basin areas, which have been used joined to the relation between topography and thermal structure of the lithosphere to calculate the best-fit surface heat flows in the study area. Moreover, we have implemented a temperature-dependent thermal conductivity (appropriate for olivine) for the lithospheric mantle to improve the calculations of temperature profiles in the mantle. The geotherms so obtained, together with the implementation of a new rheological law for the upper lithospheric mantle, have been used to calculate refined estimations of the strength and effective elastic thickness of the lithosphere. We have obtained surface heat flow values of 84 mW m⁻² and ∼82 mW m⁻² for the Spanish Central System and the Tajo Basin, respectively. The thermal state of the lithosphere affects mantle temperatures, and hence may be playing an important role in the uplift and maintenance of the Spanish Central System.     

Geochemistry of a new enriched mantle type locality in the northern hemisphere: Implications for the origin of the EM-I source

Late Cretaceous (66.2 ± 0.5 Ma amphibole and 66.7 ± 0.2 Ma phlogopite ⁴⁰Ar/³⁹Ar ages) nephelinitic volcanic rocks from Godzilla Seamount in the eastern North Atlantic (34°N latitude) have trace element and Sr–Nd–Pb–Hf-isotope compositions similar to the Enriched Mantle I (EM-I) endmember, except for their low ²⁰⁷Pb/²⁰⁴Pb relative to ²⁰⁶Pb/²⁰⁴Pb ratios (²⁰⁶Pb/²⁰⁴Pb_in = 17.7, ²⁰⁷Pb/²⁰⁴Pb_in = 15.34) plotting below the Northern Hemisphere Reference Line on the uranogenic Pb isotope diagram. O isotope data on amphibole separates are mantle-like (δ¹⁸O = 5.6–5.8‰). Age and location of the isolated Godzilla Seamount, however, preclude it from being derived from the Madeira or Canary hotspots, making a lower-mantle origin unlikely. Therefore we propose derivation from a shallow (lithospheric/asthenospheric) melting anomaly. As observed in mid-ocean-ridge and ocean-island basalts, there is a systematic decrease of ²⁰⁷Pb/²⁰⁴Pb ratios (and Δ7/4) in the individual EM-I endmember type localities towards northern latitudes with Godzilla lying on the extension of this trend. This trend is mirrored in ultra-potassic volcanic rocks such as lamproites and kimberlites, which reflect the composition of enriched subcontinental lithospheric mantle. Therefore, a global pattern in ²⁰⁷Pb/²⁰⁴Pb ratios and Δ7/4 is suggested. The geochemical composition of EM-I endmember type localities, including Godzilla lavas, and the enriched (DUPAL) anomaly in the southern hemisphere could reflect derivation from ancient, metasomatized subcontinental lithospheric mantle. We propose a two-stage model to explain the trace element and isotopic composition of the EM-I mantle endmember localities worldwide: 1) during the early history of the Earth, subcontinental lithosphere was metasomatized by melts from subducted slabs along convergent margins generating high μ (²³⁸U/²⁰⁴Pb) sources, and 2) as the Earth cooled, hydrous fluids replaced hydrous melts as the main slab component metasomatizing the subcontinental lithospheric mantle (generating EM-I sources with lower μ). In accordance with this model, the global variations in ²⁰⁷Pb/²⁰⁴Pb ratios and Δ7/4 could reflect geographic differences in μ and/or the age at which the transition from stages 1 to 2 took place in the Archaean lithosphere. The model would require a re-definition of the EM-I endmember to low ²⁰⁶Pb/²⁰⁴Pb, high ²⁰⁸Pb/²⁰⁴Pb (positive Δ8/4) but variable ²⁰⁷Pb/²⁰⁴Pb (positive and negative Δ7/4).     

Subduction erosion of forearc mantle wedge implicated in the genesis of the South Sandwich Island (SSI) arc: Evidence from boron isotope systematics

The South Sandwich volcanic arc is sited on a young oceanic crust, erupts low-K tholeiitic rocks, is characterized by unexotic pelagic and volcanogenic sediments on the down-going slab, and simple tectonic setting, and is ideal for assessing element transport through subduction zones. As a means of quantifying processes attending transfer of subduction-related fluids from the slab to the mantle wedge, boron concentrations and isotopic compositions were determined for representative lavas from along the arc. The samples show variable fluid-mobile/fluid-immobile element ratios and high enrichments of B/Nb (2.7 to 55) and B/Zr (0.12 to 0.57), similar to those observed in western Pacific arcs. δ¹¹B values are among the highest so far reported for mantle-derived lavas; these are highest in the central part of the arc (+ 15 to + 18‰) and decrease toward the southern and northern ends (+ 12 to + 14‰). δ¹¹B is roughly positively correlated with B concentrations and with ⁸⁷Sr/⁸⁶Sr ratios, but poorly coupled with other fluid-mobile elements such as Rb, Ba, Sr and U. Peridotites dredged from the forearc trench also have high δ¹¹B (ca. + 10‰) and elevated B contents (38–140 ppm). Incoming pelagic sediments sampled at ODP Site 701 display a wide range in δ¹¹B (+ 5 to ? 13‰; average = ? 4.1‰), with negative values most common. The unusually high δ¹¹B values inferred for the South Sandwich mantle wedge cannot easily be attributed to direct incorporation of subducting slab materials or fluids derived directly therefrom. Rather, the heavy B isotopic signature of the magma sources is more plausibly explained by ingress of fluids derived from subduction erosion of altered frontal arc mantle wedge materials similar to those in the Marianas forearc. We propose that multi-stage recycling of high-δ¹¹B and high-B serpentinite (possibly embellished by arc crust and volcaniclastic sediments) can produce extremely ¹¹B-rich fluids at slab depths beneath the volcanic arc. Infiltration of such fluids into the mantle wedge likely accounts for the unusual magma sources inferred for this arc.     

Estimates of the transition zone temperature in a mechanically mixed upper mantle

Partial melting of mantle peridotite generates a physically and chemically layered oceanic lithosphere that is cycled back into the mantle in subduction zones. Stirring times of the mantle are too long to allow for complete re-homogenization of subducted basalt and harzburgite, given the low chemical diffusivity of the solid mantle. This suggests that the Earth's mantle is a mechanical mixture of basaltic and harzburgitic components. Using a recently developed thermodynamic formulism we determine the phase equilibria and the seismic properties of a mantle comprised of a mechanical mixture of basalt and harzburgite (MM) and a homogeneous mantle (EA) with identical pyrolitic bulk chemistry. We use the theoretical shear velocity profiles as a new thermometer of the mantle below the magma-genetic zone by modeling the difference ΔT_410-660 between traveltimes of shear wave reflections off the 410-km and 660-km with the potential temperature T_P. ΔT_410-660 are measured from waveform stacks. They indicate that, over 1000+ km wave lengths, the temperature varies by about 200 K. Lowest and highest temperatures are resolved for the western Pacific subduction zones and the central Pacific, respectively. This variation is similar for the EA and MM and is in excellent agreement with estimates of transition zone thickness and shear velocity variations. The median value of T_P for the EA is 1720 K. It is about 1625 K for the MM, a value that is in better agreement with the Normal-MORB values of 1610 ± 40 K inferred from olivine-liquid equilibria given that our sampling region encompasses the Western Pacific subduction zones and the oldest parts of the Pacific Plate. We argue therefore that a mechanical mixed mantle, with generally higher velocities and steeper velocities gradients, represents a better physical reference model than a model based on a fully equilibrated assemblage.     

High temperatures and inverted metamorphism in the schist of Sierra de Salinas,California

New field and thermobarometric work in the Californian Salinian block clarifies current and pre-Tertiary relationships between the schist of Sierra de Salinas and Cretaceous arc-related granitic rocks. The contact is variably preserved as a brittle fault and high-temperature mylonite zone, the Salinas shear zone, which represents the contact between North America and sediments accreted above the Farallon slab between ∼ 76 Ma and ∼ 70 Ma. Near granulite facies, prograde replacement of hornblende with clinopyroxene is associated with deformation of plutonic rocks at the base of the upper plate. In the lower plate, the schist of Sierra de Salinas, garnet–biotite thermometry indicates decreasing temperatures down-section from at least 714 °C to ∼ 575 °C over an exposed thickness of ∼ 2.5 km, consistent with petrologic evidence of an inverted metamorphic gradient. The measured temperatures are significantly higher than observed at shallow levels above subducting slabs or predicted by 2D computational models assuming low shear stresses. Previous workers have called upon shear heating to explain similar observations in the correlative Pelona schist, an unlikely scenario given the results of recent rock deformation experiments which predict that feldspar–quartz–mica aggregates are far too weak to withstand stresses of ∼ 70 MPa required by the shear heating hypothesis. As an alternative, we propose that high temperatures resulted from conductive heating while the leading edge of the schist traveled ∼ 150 km beneath the recently active Salinian continental arc during the initiation of shallow subduction. Weakening of the schist due to high temperatures helped facilitate the collapse of the Salinian arc as the schist was emplaced. Schist emplacement coincided with loss of lower, mafic portions of the arc, and therefore evolution of the Southern California crust towards a more felsic composition.     

A new conceptual model for whole mantle convection and the origin of hotspot plumes

A new conceptual model of mantle convection is constructed for consideration of the origin of hotspot plumes, using recent evidence from seismology, high-pressure experiments, geodynamic modeling, geoid inversion studies, and post-glacial rebound analyses. This conceptual model delivers several key points. Firstly, some of the small-scale mantle upwellings observed as hotspots on the Earth's surface originate at the base of the mantle transition zone (MTZ), in which the Archean granitic continental material crust (TTG; tonalite-trondhjemite-granodiorite) with abundant radiogenic elements is accumulated. Secondly, the TTG crust and the subducted oceanic crust that have accumulated at the base of MTZ could act as thermal or mechanical insulators, leading to the formation of a hot and less viscous layer just beneath the MTZ; which may enhance the instability of plume generation at the base of the MTZ. Thirdly, the origin of some hotspot plumes is isolated from the large low shear-wave velocity provinces (LLSVPs) under Africa and the South Pacific. I consider that the conceptual model explains why almost all the hotspots around Africa are located above the margins of the African LLSVP. Because a planetary-scale trench system surrounding a “Pangean cell” has been spatially stable throughout the Phanerozoic, a large amount of the oceanic crustal layer is likely to be trapped in the MTZ under the Pangean cell. Therefore, under Africa, almost all of the hotspot plumes originate from the base of the MTZ, where a large amount of TTG and/or oceanic crusts has accumulated. This conceptual model may explain the fact that almost all the hotspots around Africa are located on margins above the African LLSVP. It is also considered that some of the hotspot plumes under the South Pacific thread through the TTG/oceanic crusts accumulated around the bottom of the MTZ, and some have their roots in the South Pacific LLSVP while others originate from the MTZ. The numerical simulations of mantle convection also speculate that the Earth's mantle convection is not thermally double-layered at the ringwoodite to perovskite + magnesiowüstite (Rw → Pv + Mw) phase boundary, because of its gentle negative Clapeyron slope. This is in contrast with some traditional images of mantle convection that have independent convection cells between the upper and lower mantle. These numerical studies speculate that the generation of stagnant slab at the base of the MTZ (as seismically observed globally) may not be due to the negative Clapeyron slope, and may instead be related to a viscosity increase (i.e., a viscosity jump) at the Rw → Pv + Mw phase boundary, or to a chemically stratified boundary between the upper and the lower mantle, as suggested by a recent high-pressure experiment.     

Depth constraints on azimuthal anisotropy in the Great Basin from Rayleigh-wave phase velocity maps

We present fundamental-mode Rayleigh-wave azimuthally anisotropic phase velocity maps obtained for the Great Basin region at periods between 16 s and 102 s. These maps offer the first depth constraints on the origin of the semi-circular shear-wave splitting pattern observed in central Nevada, around a weak azimuthal anisotropy zone. A variety of explanations have been proposed to explain this signal, including an upwelling, toroidal mantle flow around a slab, lithospheric drip, and a megadetachment, but no consensus has been reached. Our phase velocity study helps constrain the three-dimensional anisotropic structure of the upper mantle in this region and contributes to a better understanding of the deformation mechanisms taking place beneath the western United States. The dispersion measurements were made using data from the USArray Transportable Array. At periods of 16 s and 18 s, which mostly sample the crust, we find a region of low anisotropy in central Nevada coinciding with locally reduced phase velocities, and surrounded by a semi-circular pattern of fast seismic directions. Away from central Nevada the fast directions are ～ N–S in the eastern Great Basin, NW–SE in the Walker Lane region, and they transition from E–W to N–S in the northwestern Great Basin. Our short-period phase velocity maps, combined with recent crustal receiver function results, are consistent with the presence of a semi-circular anisotropy signal in the lithosphere in the vicinity of a locally thick crust. At longer periods (28–102 s), which sample the uppermost mantle, isotropic phase velocities are significantly reduced across the study region, and fast directions are more uniform with an ～ E–W fast axis. The transition in phase velocities and anisotropy can be attributed to the lithosphere–asthenosphere boundary at depths of ～ 60 km. We interpret the fast seismic directions observed at longer periods in terms of present-day asthenospheric flow-driven deformation, possibly related to a combination of Juan de Fuca slab rollback and eastward-driven mantle flow from the Pacific asthenosphere. Our results also provide context to regional SKS splitting observations. We find that our short-period phase velocity anisotropy can only explain ～ 30% of the SKS splitting times, despite similar patterns in fast directions. This implies that the origin of the regional shear-wave splitting signal is complex and must also have a significant sublithospheric component.     

Mantle transition zone beneath the Caribbean-South American plate boundary and its tectonic implications

We analyzed receiver-function data recorded by a temporary broadband array deployed as part of the BOLIVAR project and the permanent seismic network of Venezuela to study the mantle transition zone structure beneath the Caribbean-South American plate boundary and Venezuela. Significant topography on both the 410-km and the 660-km discontinuities was clearly imaged in the CCP (common-conversion-point) stacked images. Beneath the southeastern Caribbean, the 410-km is featured by a narrow (～ 200 km EW) ～ 25-km uplift extending in the NS direction around 63° west, while the 660-km is depressed by ～ 20 km in a narrow region slightly west to the uplift, a scenario that is more consistent with westward descent of the oceanic South American plate rather than a break-off of NNW dipping proto-Caribbean oceanic lithosphere along the El Pilar Fault. We also found a thick transition zone beneath the Falcon region in northwestern Venezuela, possibly associated with the subducted Nazca plate. A flat 410-km was observed beneath the Guayana shield, suggesting that the shield has a stable and moderately deep keel, which has little effect on the underlying transition zone structure.     

Thermal conditions of the central Sinai Microplate inferred from new surface heat-flow values and continuous borehole temperature logging in central and southern Israel

This paper reports ten new surface heat-flow density (qs) values for central and southern Israel (central Sinai Microplate), whose crystalline crust and lithosphere formed as part of the Neoproterozoic Arabian-Nubian Shield. Heat flow was calculated in Mesozoic sediments using the classical approach of heat-flow determination by implementing in the analysis high-precision continuous temperature logs obtained in air- and/or water-filled boreholes. Thermal conductivity (TC) measured for a large suite of rock samples of lithotypes making up the sequence was assigned to temperature gradients in intervals for which the lithology was known. The heat-flow values obtained for different depth intervals in a borehole as well as the average values for the individual borehole locations cover a narrow range, attesting heat-conduction conditions. A steady-state thermal model along an E–W crustal cross section through the area shows that the observed systematic spatial distribution of the qs values, which range between 50 and 62 mW m⁻², can primarily be explained by variations in the thickness of the upper crust and in the ratio between sedimentary and crystalline rocks therein. Given the time lapse of thermal heat transfer through the lithosphere, the qs data monitor the crustal thermal conditions prior to rift- and plume-related lithospheric thermal perturbations that have started in the larger area ca. 30 Ma ago. Observed and modeled qs display the best fit for a pre-Oligocene lithosphere–asthenosphere boundary (LAB) at ∼150 km, which would be at the upper end of LAB depths determined from stable areas of the Arabian Shield (150–120 km) not affected by the young, deep-seated thermal processes that have caused a further uprise of the LAB. Our data imply or predict that the surface heat flow of the Sinai Microplate generally tends to increase along N–S and W–E traverses, from ∼45–50 mW m⁻² to ∼55–60 mW m⁻². Surface heat flows on the order of 55–60 mW m⁻² may be common in the northern Arabian Shield, where it exhibits typical lithosphere structure and composition and is unaffected by young heating processes, compared to values of ≤45 mW m⁻² recently determined in the southern Arabian Plate for the Arabian Platform.     

Geological,tomographic, kinematic and geodynamic constraints on the dynamics of sinking slabs

We use geodynamic models with imposed plate velocities to test the forward-modeled history of subduction based on a particular plate motion model against alternative seismic tomography models. We utilize three alternative published reference frames: a hybrid moving hotspot-palaeomagnetic, a hybrid moving hotspot-true polar wander corrected-palaeomagnetic, and a Subduction Reference Frame, a plate model including longitudinal shifts of subduction zones by matching subduction volumes imaged by P-wave tomography, to assess which model best predicts present day mantle structure compared with seismic tomography and volumetrically derived subduction history. Geodynamic modeling suggests paleo-longitudinal corrections applied to the Subduction Reference Frame result in lower mantle slab material beneath North America and East Asia accumulating up to 10–15° westward of that imaged by tomography, whereas the hybrid models develop material offset by 2–9°. However, the Subduction Reference Frame geodynamic model produces slab material beneath the Tethyan Domain coinciding with slab volumes imaged by tomography, whereas the hybrid reference frame models do not, suggesting regional paleo-longitudinal corrections are required to constrain slab locations. We use our models to test inferred slab sinking rates in the mantle focusing on well-constrained regions. We derive a globally averaged slab-sinking rate of 13 ± 3 mm/yr by combining the ages of onset and cessation of subduction from geological data and kinematic reconstructions with images of subducted slabs in the mantle. Our global average slab-sinking rate overlaps with the 15–20 mm/yr rate implied by mantle convection models using a lower mantle viscosity 100 times higher than the upper mantle.     

Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units,Lyon Turin transect,Western Alps)

In the central Western Alps, a combined structural, petrological and ⁴⁰Ar–³⁹Ar geochronological study of the Modane-Aussois and Southern Vanoise units yields important constraints on the timing of deformation and exhumation of the Briançonnais zone. These data help to decipher the respective roles of oceanic subduction, continental subduction and collision in the burial and exhumation of the main units through time. In the Modane-Aussois unit top to the NW thrusting (D1) was followed by top to the east shearing (D2) interpreted by some as normal faulting and by others as backthrusting. Pseudosection calculations imply that D1 deformation occurred at 1.0 ± 0.1 GPa and 350 ± 30 °C. Analysis of chlorite–phengite pairs yield P–T estimates between 0.15 and 0.65 GPa and between 220 and 350 °C for the D2 event. Phengites along the D1 schistosity (sample M80) yields an ⁴⁰Ar–³⁹Ar age of 37.12 ± 0.39 Ma, while D2 phengites yield ages of 35.42 ± 0.38 (sample M173) and 31.60 ± 0.33 Ma (sample M196). It was not possible to test whether these ages are altered by excess argon or not. Our interpretation is that the D1/D2 transition occurred at ∼37 Ma at the beginning of decompression, and that D2 lasted until at least ∼32 Ma. Pseudosection calculation suggests that the Southern Vanoise unit was buried at 1.6 ± 0.2 GPa and 500–540 °C. D1 deformation occurred during exhumation until 0.7–10.5 GPa and 370 ± 30 °C. Published ages suggest that D1 deformation possibly started at ∼50 Ma and lasted until ∼37 Ma. D2 deformations started at P–T conditions close to that recorded in Modane-Aussois unit and lasted until 0.2 ± 0.1 GPa and 280 ± 30 °C at ∼28 Ma. The gap of 0.6 ± 0.3 GPa and 150 ± 130 °C between peak metamorphic conditions in the two units was concealed by thrusting of the South Vanoise unit on top of the Modane-Aussois unit during D1 Deformation. Top to the east deformation (D2) affects both units and is interpreted as backthrusting.Based on these data, we propose a geodynamic reconstruction where the oceanic subduction of the Piedmont unit until ∼50 Ma, is followed by its exhumation at the time of continental subduction of the continental Southern Vanoise unit until ∼45 Ma. The Southern Vanoise is in turn underthrusted by the Modane-Aussois unit until ∼37 Ma (D1). Between 37 and 31 Ma the Modane-Aussois and Southern Vanoise units exhume together during backthrusting to the east (D2). This corresponds to the collision stage and to the activation of the Penninic Thrust. In the ∼50 Ma to ∼31 Ma time period the main thrusts propagated westward as the tectonic context switched from oceanic to continental subduction and finally to collision. During each stage, external units are buried while internal ones are exhumed.     

Relation between mantle shear zone deformation and metasomatism in spinel peridotite xenoliths of Jeju Island (South Korea): Evidence from olivine CPO and trace elements

This study presents geochemical and fabric analysis of coarse-grained, porphyroclastic and mylonitic spinel peridotite xenoliths derived from the shallow subcontinental lithospheric mantle of Jeju Island (South Korea). Fabric analysis of olivines in the studied peridotites shows activation of the high temperature (0 1 0)[1 0 0] slip system; however, in the porphyroclastic and mylonitic peridotites, activation of (0 k l)[1 0 0] slip system results in a weaker fabric. Formation of porphyroclastic and mylonitic fabrics are thought to form in a shear-dominated environment. The results of the trace element analysis reveal that the smaller the grain size and weaker the fabric, the more enriched in LREE and HFSE are the peridotites, which indicates a strong relationship between metasomatic agents and mantle shear zones.     

Critical phenomena in thermal conductivity: Implications for lower mantle dynamics

Microscopic mechanisms for heat transport in dense minerals (phonon scattering and photon attenuation) exhibit aspects of threshold behavior, discussed qualitatively here. For all minerals examined so-far using laser-flash analysis, the lattice component of the thermal conductivity of the mantle asymptotes to a constant above a critical temperature of ∼1500 K. Radiative transfer calculated from absorption spectra has thresholds in both grain-size and Fe content, and a rather complex dependence on temperature. These critical phenomena impact convection of the lower mantle, because the lattice contribution tends to destabilize the cold boundary layers, whereas radiative transfer mostly promotes stability in the lower mantle, unless the grains are large and Fe-rich, which makes convection chaotic and time-dependent. The specific behavior suggests that flow in the lower mantle is sluggish, whereas flow in the upper mantle-transition zone is time-dependent. The decrease in k_rad as Fe/(Fe + Mg) increases beyond ∼0.1 may be connected with formation of lower mantle, thermo-chemical plumes through positive feedback.     

Early crust on top of the Earth's core

In an effort to resolve the current conflict between geochemical requirements for an apparently isolated mantle reservoir and recent geophysical evidence for whole-mantle convection, we investigate the possibility that the region above the core-mantle boundary, termed D″, serves as an early-isolated, rare-gas- and incompatible-element-bearing reservoir, and we propose a mechanism for its formation that is a likely outcome of Earth accretion models. In these models, the most cataclysmic event in Earth history, the moon-forming giant impact on the proto-Earth of a Mars-size planet (perhaps as early as 4540 Ma ago) was followed by accretion of smaller bodies long afterwards (until ∼3900 Ma). Some collisions probably triggered melting, metal segregation and degassing. However, the small bodies, fragments, particles, dust, including those of chondrite-like composition, existed on near-earth orbits, were irradiated by intense solar wind, and finally fell on an early-formed, incompatible-element-enriched basaltic crust without causing extensive melting. The respective regions of the crust, loaded with chondrite-like debris, were therefore significantly enriched in iron. When this mixed material was subducted, the bulk density of the subducted lithosphere exceeded that of the bulk silicate mantle, which had already lost its metallic iron to the core. Segregation of this denser material at the base of the mantle was facilitated by the high temperatures at the core-mantle boundary, which greatly reduce the viscosity, as was quantitatively modelled by Christensen and Hofmann (Christensen, U.R., Hofmann, A.W., 1994. Segregation of subducted oceanic-crust in the convecting mantle. J. Geophys. Res.-Solid Earth 99 (B10), 19867–19884). Assuming a basalt/chondrite mass ratio of about 4/1, we obtain a density contrast of ∼7%, which would stabilize the subducted material between the metal core and silicate mantle.Mass balance considerations and preliminary results of geochemical modelling of the above scenario (similar to that performed by Tolstikhin and Marty [Tolstikhin, I.N., Marty, B., 1998. The evolution of terrestrial volatiles, A view from helium, neon, argon and nitrogen isotope modeling. Chem. Geol. 147, 27–52]) show the potential geochemical importance of D″. (1) Modelling of Pu–U–I–Xe isotope systematics predicts formation of this reservoir early in Earth history, ∼100 Ma after formation of the Solar system. (2) The total amount of heat-generating U, Th, K (and other highly incompatible elements) in D″ exceeds 20% of the Earth inventory, and a similar portion of terrestrial heat is being transferred from the core + D″ into the base of the overlying convecting mantle. (3) D″ is enriched in solar implanted rare gases because the small (re)-accreting debris with high surface/mass ratios will have been subjected to intense radiation by the early sun. (4) Rare gases diffuse from D″ into the overlying mantle and are then transferred into upwelling plumes, which originate above D″. In addition, small amounts of D″ material may be entrained by the mantle convective flow as was recently discussed by Schott et al. [Schott, B., Yuen, D.A., Braun, A., 2002. The influences of composition and temperature-dependent rheology in thermal-chemical convection on entrainment of the D″ layer. Physics Earth Planet. Inter. 129, 43–65]. From the rare-gas modelling it follows that initially (∼4500 Ma ago) D″ could have been more massive by a factor of ∼1.2 than at present (about 2 × 10²⁶ g). The present-day mass flux from D″ into the convecting mantle is estimated to be ≤0.05 × 10¹⁶ g year⁻¹, a factor of ∼100 less than the rate of ridge magmatism. This small contribution of D″ material makes it difficult to trace fingerprints of D″ even using such sensitive tracers as Pb isotope ratios. (5) The density contrast that stabilizes D″ is maintained by its higher intrinsic density due to the iron-rich chondrite-like component. Subduction of this material, its entrainment by convective mantle flow and mixing could also account for the preservation of the chondritic relative abundances of siderophile elements in the mantle. If D″ is partially molten, the density contrast may be caused by a high-density melt fraction.     

Scattering coefficients in the mantle revealed from the seismogram envelope analysis based on the multiple isotropic scattering model

At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10^− 4∼1.32 × 10^− 3 km^− 1 and 2.08 × 10^− 4∼6.23 × 10^− 4 km^− 1 at 4 s, and 4.51 × 10^− 4∼7.37 × 10^− 4 km^− 1 and 2.80 × 10^− 5∼2.71 × 10^− 4 km^− 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.     

Mantle dynamics of the Paleoproterozoic North China Craton: A perspective based on seismic tomography

We investigate the mantle dynamics beneath the North China Craton (NCC) and surrounding regions based on a synthesis of recent P-wave mantle tomographic data down to depths of 600–800 km and their correlation with the surface geological features, with particular reference to the Paleoproterozoic tectonic events associated with the incorporation of the NCC within the Columbia supercontinent amalgam. From the tomographic images, we identify a hot corridor in the mantle transition zone beneath the central region of the Western Block of the NCC sandwiched between two cold corridors. This scenario is similar to the donut-shaped high-velocity anomaly surrounding a region of low-velocity anomaly in the lowermost mantle under the Pacific and suggests that the cold regions might represent slab graveyards which provide the fuel for the plumes rising from the center. A tomographic transect along the collisional suture of the NCC with the Columbia supercontinent, covering the Yinshan-Ordos Blocks in the Western Block through the Central Orogenic Belt and into the Eastern Block of the NCC reveals a ca. 250 km thick lithospheric keel below the Ordos Block defined by a prominent high-velocity anomaly. We identify slab break-off and asthenospheric upwelling in this region and suggest that this process probably initiated the thermal and material erosion of the tectosphere beneath the Eastern Block from the Paleoproterozoic, which was further intensified during the Mesozoic when a substantial part of the sub-continental mantle lithosphere was lost. We visualize heat input from asthenosphere and interaction between asthenosphere and overlying carbonated tectosphere releasing CO₂-rich fluids for the preservation of ultra-high temperature (ca. 1000 °C) metamorphic rocks enriched in CO₂ as well as high-pressure mafic granulites as a paired suite in this region. We also identify a hot swell of the asthenosphere rooted to more than 200 km depth and reaching up to the shallow mantle in the tomographic section along 35°N latitude at a depth of 800 km. This zone represents a cross-section through the southern part of the NCC. The surface distribution of Paleoproterozoic Xiong’er lavas and mafic dykes in this region would indicate that this region might have evidenced similar upwellings in the past. Our study has important implications in understanding the evolution of the NCC and suggests that the extensive modification of the mantle architecture and lithospheric structure beneath one of the fundamental Precambrian nuclei of Asia had a prolonged history probably dating from the Paleoproterozoic suturing of the NCC within the Columbia supercontinent amalgam.