Phase transformations in a harzburgite composition to 26 GPa: implications for dynamical behaviour of the subducting slab

The mineralogy adopted by a depleted harzburgite composition has been studied over the pressure interval 5–26 GPa at temperatures of 1300–1400°C. The pyroxene-garnet component of the harzburgite composition (harzburgite minus 82 wt.% olivine) transforms to majorite garnet by 18–19 GPa, and further disproportionates to the assemblage of garnet + stishovite + Mg₂SiO₄ spinel above 20 GPa. At still higher pressures, first ilmenite (22–24 GPa) and then perovskite MgSiO₃ (24–26 GPa) are found to coexist with garnet. Garnet disappears at 26 GPa and almost complete transition to perovskite is achieved at this pressure. The mineral proportions and density profiles in the subducting oceanic lithosphere, modelled by a combination of 80% harzburgite + 20% primitive MORB compositions are calculated as a function of depth under conditions isothermal with surrounding pyrolite mantle, and also for a temperature distribution in which the slab is substantially cooler than surrounding mantle to below 700 km. Under isothermal conditions, the slab has a density similar to surrounding mantle to a depth of 600 km. However, between 600 and 700 km, the slab is up to 0.08 g/cm³ denser than surrounding mantle. This is caused primarily by the higher alumina content in pyrolite as compared to harzburgite, which causes the transition to perovskite in pyrolite to occur at substantially higher pressures than in harzburgite. The presence of alumina also smears out the garnet-perovskite transition in pyrolite over a depth interval of 50 km, whereas this transformation is much sharper in the harzburgite composition. Calculations based on the observed phase equilibria also show that a subducted cool slab remains much denser (by 0.1–0.3 g/cm³) than surrounding mantle to a depth of 700 km but possesses a density similar to surrounding mantle below this depth. These results have important implications for the dynamical behaviour of slabs possessing different thermal regimes when they encounter the 670 km discontinuity and also for the nature of this discontinuity.     

Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

Formation time of the big mantle wedge beneath eastern China and a new lithospheric thinning mechanism of the North China craton—Geodynamic effects of deep recycled carbon

High-resolution P wave tomography shows that the subducting Pacific slab is stagnant in the mantle transition zone and forms a big mantle wedge beneath eastern China. The Mg isotopic investigation of large numbers of mantle-derived volcanic rocks from eastern China has revealed that carbonates carried by the subducted slab have been recycled into the upper mantle and formed carbonated peridotite overlying the mantle transition zone, which becomes the sources of various basalts. These basalts display light Mg isotopic compositions(δ26 Mg = –0.60‰ to –0.30‰) and relatively low87 Sr/86 Sr ratios(0.70314–0.70564) with ages ranging from 106 Ma to Quaternary, suggesting that their mantle source had been hybridized by recycled magnesite with minor dolomite and their initial melting occurred at 300-360 km in depth. Therefore, the carbonate metasomatism of their mantle source should have occurred at the depth larger than 360 km, which means that the subducted slab should be stagnant in the mantle transition zone forming the big mantle wedge before 106 Ma. This timing supports the rollback model of subducting slab to form the big mantle wedge. Based on high P-T experiment results, when carbonated silicate melts produced by partial melting of carbonated peridotite was raising and reached the bottom(180–120 km in depth) of cratonic lithosphere in North China, the carbonated silicate melts should have 25–18 wt% CO2 contents, with lower Si O2 and Al2 O3 contents, and higher Ca O/Al2 O3 values, similar to those of nephelinites and basanites, and have higher εNdvalues(2 to 6). The carbonatited silicate melts migrated upward and metasomatized the overlying lithospheric mantle, resulting in carbonated peridotite in the bottom of continental lithosphere beneath eastern China. As the craton lithospheric geotherm intersects the solidus of carbonated peridotite at 130 km in depth, the carbonated peridotite in the bottom of cratonic lithosphere should be partially melted, thus its physical characters are similar to the asthenosphere and it could be easily replaced by convective mantle. The newly formed carbonated silicate melts will migrate upward and metasomatize the overlying lithospheric mantle. Similarly, such metasomatism and partial melting processes repeat, and as a result the cratonic lithosphere in North China would be thinning and the carbonated silicate partial melts will be transformed to high-Si O2 alkali basalts with lower εNdvalues(to-2). As the lithospheric thinning goes on,initial melting depth of carbonated peridotite must decrease from 130 km to close 70 km, because the craton geotherm changed to approach oceanic lithosphere geotherm along with lithospheric thinning of the North China craton. Consequently, the interaction between carbonated silicate melt and cratonic lithosphere is a possible mechanism for lithosphere thinning of the North China craton during the late Cretaceous and Cenozoic. Based on the age statistics of low δ26 Mg basalts in eastern China, the lithospheric thinning processes caused by carbonated metasomatism and partial melting in eastern China are limited in a timespan from 106 to25 Ma, but increased quickly after 25 Ma. Therefore, there are two peak times for the lithospheric thinning of the North China craton: the first peak in 135-115 Ma simultaneously with the cratonic destruction, and the second peak caused by interaction between carbonated silicate melt and lithosphere mainly after 25 Ma. The later decreased the lithospheric thickness to about70 km in the eastern part of North China craton.     

Subduction zone seismicity and the thermo-mechanical evolution of downgoing lithosphere

In this paper we discuss characteristic features of subduction zone seismicity at depths between about 100 km and 700 km, with emphasis on the role of temperature and rheology in controlling the deformation of, and the seismic energy release in downgoing lithosphere. This is done in two steps. After a brief review of earlier developments, we first show that the depth distribution of hypocentres at depths between 100 km and 700 km in subducted lithosphere can be explained by a model in which seismic activity is confined to those parts of the slab which have temperatures below a depth-dependent critical valueT _cr.Second, the variation of seismic energy release (frequency of events, magnitude) with depth is addressed by inferring a rheological evolution from the slab's thermal evolution and by combining this with models for the system of forces acting on the subducting lithosphere. It is found that considerable stress concentration occurs in a reheating slab in the depth range of 400 to 650–700 km: the slab weakens, but the stress level strongly increases. On the basis of this stress concentration a model is formulated for earthquake generation within subducting slabs. The model predicts a maximum depth of seismic activity in the depth range of 635 to 760 km and, for deep earthquake zones, a relative maximum in seismic energy release near the maximum depth of earthquakes. From our modelling it follows that, whereas such a maximum is indeed likely to develop in deep earthquake zones, zones with a maximum depth around 300 km (such as the Aleutians) are expected to exhibit a smooth decay in seismic energy release with depth. This is in excellent agreement with observational data. In conclusion, the incoroporation of both depth-dependent forces and depth-dependent rheology provides new insight into the generation of intermediate and deep earthquakes and into the variation of seismic activity with depth.Our results imply that no barrier to slab penetration at a depth of 650–700 km is required to explain the maximum depth of seismic activity and the pattern of seismic energy release in deep earthquake zones.     

Time dependence of negative buoyancy and the subduction of continental lithosphere

The time evolution of negative buoyancy of a subducting slab is modelled from the beginning of subduction under various kinematic conditions (dip angle and subduction velocity). The calculations take into account the thermal and density effects of the variations of the thermophysical parameters with temperature and pressure, and of phase transitions. The magnitude of the negative buoyancy increases during subduction of oceanic lithosphere, up to values in the (2–4) × 10¹³ N m⁻¹ range when the tip of the slab reaches a depth of 600–700 km. If continental material arrives at the trench and is subducted, the downward buoyancy decreases by an amount proportional to the volume of the subducted continental crust. Assuming that subduction stops when the buoyancy becomes zero, and that delamination of the continental crust or slab breakoff do not occur, the maximum downdip length of the subductable continental crust is estimated as a function of the dip angle, subduction velocity and geometry of the margin. In most cases, subduction of continental material down to depths of 100–250 km is possible, and continental subduction can continue for times up to 10–15 Ma if the velocity is low. These estimates are not significantly affected by the hypothetical occurrence of a metastable olivine wedge within the slab, and could be lower bounds if the lower continental crust is mafic and transforms to eclogite.     

Subduction zone earthquakes and stress in slabs

Summary The pattern of seismicity as a function of depth in the world, and the orientation of stress axes of deep and intermediate earthquakes, are explained using viscous fluid models of subducting slabs, with a barrier in the mantle at 670 km. 670 km is the depth of a seismic discontinuity, and also the depth below which earthquakes do not occur. The barrier in the models can be a viscosity increase of an order of magnitude or more, or a chemical discontinuity where vertical velocity is zero. LongN versus depth, whereN is the number of earthquakes, shows (1) a linear decrease to about 250–300 km depth, (2) a minimum near that depth, and (3) an increase thereafter. Stress magnitude in a subducting slab versus depth, for a wide variety of models, shows the same pattern. Since there is some experimental evidence thatN is proportional toe , where is a constant and is the stress magnitude, the agreement is encouraging. In addition, the models predict down-dip compression in the slab at depths below 400 km. This has been observed in earlier studies of earthquake stress axes, and we have confirmed it via a survey of events occurring since 1977 which have been analysed by moment tensor inversion. At intermediate depths, the models predict an approximate but not precise state of down-dip tension when the slab is dipping. The observations do not show an unambiguous state of down-dip tension at intermediate depths, but in the majority of regions the state of stress is decidedly closer to down-dip tension than it is to down-dip compression. Chemical discontinuities above 670 km, or phase transitions with an elevation of the boundary in the slab, predict, when incorporated into the models, stress peaks which are not mirrored in the profile of seismicity versus depth. Models with an asthenosphere and mesosphere of appropriate viscosity can not only explain the state of stress observed in double Benioff zones, but also yield stress magnitude profiles consistent with observed seismicity. Models where a nonlinear rheology is used are qualitatively consistent with the linear models.     

Intraplate seismicity in the western Bohemian Massif (central Europe): A possible correlation with a paleoplate junction

Locations of the Eger Rift, Cheb Basin, Quaternary volcanoes, crustal earthquake swarms and exhalation centers of CO₂ and ³He of mantle origin correlate with the tectonic fabric of the mantle lithosphere modelled from seismic anisotropy. We suggest that positions of the seismic and volcanic phenomena, as well as of the Cenozoic sedimentary basins, correlate with a “triple junction” of three mantle lithospheres distinguished by different orientations of their tectonic fabric consistent within each unit. The three mantle domains most probably belong to the originally separated microcontinents – the Saxothuringian, Teplá-Barrandian and Moldanubian – assembled during the Variscan orogeny. Cenozoic extension reactivated the junction and locally thinned the crust and mantle lithosphere. The rigid part of the crust, characterized by the presence of earthquake foci, decoupled near the junction from the mantle probably during the Variscan. The boundaries (transitions) of three mantle domains provided open pathways for Quaternary volcanism and the ascent of ³He- and CO₂-rich fluids released from the asthenosphere. The deepest earthquakes, interpreted as an upper limit of the brittle–ductile transition in the crust, are shallower above the junction of the mantle blocks (at about 12 km) than above the more stable Saxothuringian mantle lithosphere (at about 20 km), probably due to a higher heat flow and presence of fluids.     

Geodynamic implications of deep mantle upwelling in the source of Tertiary volcanics from the Veneto region (South-Eastern Alps)

Major and trace element and Sr–Nd–Hf–Pb isotopic data for the most primitive Tertiary lavas from the Veneto region (South-Eastern Alps, Italy) show the typical features of HIMU hotspot volcanism, variably diluted by a depleted asthenospheric mantle component (⁸⁷Sr/⁸⁶Sr_i=0.70306–0.70378; _Ndi=+3.9 to +6.8; _Hfi=+6.4 to +8.1, ²⁰⁶Pb/²⁰⁴Pb_i=18.786–19.574). P-wave seismic tomography of the mantle below the Veneto region shows the presence of low-velocity anomalies at depth, which is consistent with possible upwellings of plume material. Between the depths of 100–250 km the velocity anomalies are approximately 2–2.5% slower than average, implying a temperature excess of about 220–280 K, in agreement with estimates for other mantle plumes in the world. In this context, the Veneto volcanics may represent the shallow expression of a mantle upflow. The presence of a HIMU-DM component in a collision environment has significant geodynamic implications. Slab detachment and ensuing rise of deep mantle material into the lithospheric gap is proposed to be a viable mechanism of hotspot magmatism in a subduction zone setting.     

On the heat flux related to stretching in the NW-Mediterranean continental margins

Summary The surface thermal flux of the continental margins of the northwestern Mediterranean Sea is interpreted on the basis of a 1-D instantaneous pure shear stretching model of the lithosphere in terms of three components: the background heat flowing out from the asthenosphere (38 mW m^–2), the transient contribution depending on the rift age and extension amount (35 mW m^–2 at the most), and the contribution due to the radiogenic elements of the lithosphere. The radiogenic component is estimated at the continental margins of the Ligurian-Provençal basin and Valencia trough, and in the surrounding mainland areas by means of available data of surface heat generation from Variscan Corsica, Maures-Estérel and the Central Massif along with a geophysical-petrological relationship between heat production and seismic velocity. The lithosphere radiogenic heat contribution q_l decreases with the thinning factor according to the exponential law: q_l() = a exp(-b), in which factor b is greater for that part of the lithosphere below the uppermost 10 km. Considering also the heat generated by radioactive isotopes in sediments, the stable Variscan lithosphere produces an average thermal flux of 30 mW m^–2 which decreases by about one half where the lithosphere is thinned by one third. Although the surface heat generation is 2·1 – 3·3 µW m^–3 in the Maures-Estérel massif — excepting small outcrops of dioritic rocks with lower heat production — and 1·8 µW m^–3 for most of Corsica, the radiogenic heating within the lithosphere for such areas is nearly the same and does not explain the higher heat flux of the Corsica margin. This asymmetric thermal pattern with surface heat flux which is 10 – 15 mW m^–2 higher than predictions is probably of upper mantle origin, or can be ascribed to penetrative magmatism.     

Depth distribution of stresses in the Hokkaido Wadati-Benioff zone as deduced by inversion of earthquake focal mechanisms

This study is based on the detailed geometry of the Hokkaido Wadati-Benioff zone and the paleosubduction zone as delineated by Hanus and Vanek (1984). The used data includes 217 CMT Harvard solutions for earthquakes, which belong to the Wadati-Benioff zone and 13 for the paleosubduction zone. The inverse technique by Gephart and Forsyth (1984) was incorporated for determining the best fit principal stress directions σ₁, σ₂, σ₃ and the ratio (R=σ₂−σ₁/σ₃−σ₁) for 20 km depth intervals in the Wadati-Benioff zone and for the paleosubduction zone considered as a single body. In almost all the considered depth layers, the maximum compressive stress σ₁ is normal to the strike of the slab and dips less than 25°, indicating the NW-SE convergence between the Pacific and Eurasian lithospheric plates. Exceptions are in the depth layer 81–120 km, the paleosubduction zone with steeply dipping along-strike σ₁, and the lower part of the subduction zone (161–220 km) where σ₁ is almost horizontal and of E trend. The minimum compressive stress σ₃ is mostly along-strike and of a different dip with the exception of the 21–60 km layer wher they are down-dipping. The results obtained for the depth ranges 0–20 km, 81–100 km, 121–160 km, and the paleosubduction zone indicate heterogeneous stress fields. These results show that the slab pull and the mantle resistance, acting on the slab edge, are not the main forces which control the contemporary plate tectonics in the Hokkaido region. Along-strike compression at depths 81–120 km and along-strike extension at 0–20 and 61–220 km are involved in the slab dynamics. These can be related to horizontal bending of the subducting Pacific plate.     

Mesozoic mafic magmatism in North China: Implications for thinning and destruction of cratonic lithosphere

The North China Craton (NCC) has been thinned from >200 km to <100 km in its eastern part. The ancient subcontinental lithospheric mantle (SCLM) has been replaced by the juvenile SCLM in the Meoszoic. During this period, the NCC was destructed as indicated by extensive magmatism in the Early Cretaceous. While there is a consensus on the thinning and destruction of cratonic lithosphere in North China, it has been hotly debated about the mechanism of cartonic destruction. This study attempts to provide a resolution to current debates in the view of Mesozoic mafic magmatism in North China. We made a compilation of geochemical data available for Mesozoic mafic igneous rocks in the NCC. The results indicate that these mafic igneous rocks can be categorized into two series, manifesting a dramatic change in the nature of mantle sources at ~121 Ma. Mafic igneous rocks emplaced at this age start to show both oceanic island basalts (OIB)-like trace element distribution patterns and depleted to weakly enriched Sr-Nd isotope compositions. In contrast, mafic igneous rocks emplaced before and after this age exhibit both island arc basalts (IAB)-like trace element distribution patterns and enriched Sr-Nd isotope compositions. This difference indicates a geochemical mutation in the SCLM of North China at ~121 Ma. Although mafic magmatism also took place in the Late Triassic, it was related to exhumation of the deeply subducted South China continental crust because the subduction of Paleo-Pacific slab was not operated at that time. Paleo-Pacific slab started to subduct beneath the eastern margin of Eruasian continent since the Jurrasic. The subducting slab and its overlying SCLM wedge were coupled in the Jurassic, and slab dehydration resulted in hydration and weakening of the cratonic mantle. The mantle sources of ancient IAB-like mafic igneous rocks are a kind of ultramafic metasomatites that were generated by reaction of the cratonic mantle wedge peridotite not only with aqueous solutions derived from dehydration of the subducting Paleo-Pacific oceanic crust in the Jurassic but also with hydrous melts derived from partial melting of the subducting South China continental crust in the Triassic. On the other hand, the mantle sources of juvenile OIB-like mafic igneous rocks are also a kind of ultramafic metasomatites that were generated by reaction of the asthenospheric mantle underneath the North China lithosphere with hydrous felsic melts derived from partial melting of the subducting Paleo-Pacific oceanic crust. The subducting Paleo-Pacific slab became rollback at ~144 Ma. Afterwards the SCLM base was heated by laterally filled asthenospheric mantle, leading to thinning of the hydrated and weakened cratonic mantle. There was extensive bimodal magmatism at 130 to 120 Ma, marking intensive destruction of the cratonic lithosphere. Not only the ultramafic metasomatites in the lower part of the cratonic mantle wedge underwent partial melting to produce mafic igneous rocks showing negative ε_Nd(t) values, depletion in Nb and Ta but enrichment in Pb, but also the lower continent crust overlying the cratonic mantle wedge was heated for extensive felsic magmatism. At the same time, the rollback slab surface was heated by the laterally filled asthenospheric mantle, resulting in partial melting of the previously dehydrated rocks beyond rutile stability on the slab surface. This produce still hydrous felsic melts, which metasomatized the overlying asthenospheric mantle peridotite to generate the ultramafic metasomatites that show positive ε_Nd(t) values, no depletion or even enrichment in Nb and Ta but depletion in Pb. Partial melting of such metasomatites started at ~121 Ma, giving rise to the mafic igneous rocks with juvenile OIB-like geochemical signatures. In this context, the age of ~121 Ma may terminate replacement of the ancient SCLM by the juvenile SCLM in North China. Paleo-Pacific slab was not subducted to the mantle transition zone in the Mesozoic as revealed by modern seismic tomography, and it was subducted at a low angle since the Jurassic, like the subduction of Nazca Plate beneath American continent. This flat subduction would not only chemically metasomatize the cratonic mantle but also physically erode the cratonic mantle. Therefore, the interaction between Paleo-Pacific slab and the cratonic mantle is the first-order geodynamic mechanism for the thinning and destruction of cratonic lithosphere in North China.     

Space Variations of Stress Along the Tyrrhenian Wadati-Benioff Zone

—Gephart and Forsyth’s (1984) algorithm for stress inversion of earthquake fault-plane solutions has been applied to a set of ninety intermediate and deep events occurring in the southern Tyrrhenian region between 1976 and 1995. P- and S-wave data from local seismic networks in southern Italy, the Italian National Network and international bulletins, have been used for hypocenter and focal mechanism computations. Stress inversion runs performed after accurate selection and weighting of fault-plane solutions have allowed us to identify stress space variations at a higher level of detail than available from all previous investigations carried out in the study area. The maximum compressive stress has been shown to follow the depth-decreasing dip of the Wadati-Benioff zone, along the entire zone from a depth of 90 km, to the depth of the deepest events (about 500 km). Variations to such a stress pattern have been found, possibly related to mantle dynamics and the complex composition of the subducting structure. The diffused state of down-dip compression suggests that the Tyrrhenian subduction has already evolved to the point where the lower end of the slab has reached high-strength mantle materials, the load of the excess mass is entirely supported from below and most of the subducted slab is under compression. In agreement with the lack of large, shallow thrusting events in the immersion zone, the findings of the present study appear to agree well with geodynamic models assuming a passive subduction process with eastward roll-back of the Ionian lithosphere in the study area. In this context, the depth-decrease of the slab dip may also find a reasonable explanation.     

Composition,temperature, and thickness of the lithosphere of the Archean Kaapvaal craton

A new model is proposed for the structure of the Kaapvaal craton lithosphere. Based on chemical thermodynamics methods, profiles of the chemical composition, temperature, density, and S wave velocities are constructed for depths of 100–300 km. A solid-state zone of lower velocities is discovered on the S velocity profile in the depth interval 150–260 km. The temperature profiles are obtained from absolute values of P and S velocities, taking into account phase transformations, anharmonicity, and anelastic effects. The examination of the sensitivity of seismic models to the chemical composition showed that relatively small variations in the composition of South African xenoliths result in lateral temperature variations of ～200°C. Inversion of some seismic profiles (including IASP91) with a fixed bulk composition of garnet peridotites (the primitive mantle material) leads to a temperature inversion at depths of 200–250 km, which is physically meaningless. It is supposed that the temperature inversion can be removed by gradual fertilization of the mantle with depth. In this case, the craton lithosphere should be stratified in chemical composition. The depleted lithosphere composed by garnet peridotites exists to depths of 175–200 km. The lithospheric material at depths of 200–250 km is enriched in basaltoid components (FeO, Al₂O₃, and CaO) as compared with the material of garnet peridotites but is depleted in the same components as compared with the fertile substance of the underlying primitive mantle. The material composing the craton root at a depth of ～275 km does not differ in its physical and chemical characteristics from the composition of the normal mantle, and this allows one to estimate the thickness of the lithosphere at 275 km. The results of this work are compared with data of seismology, thermal investigations, and thermobarometry.     

2-D Crustal Velocity Structure Along Hirapur-Mandla Profile in Central India: An Update

The 2-D crustal velocity model along the Hirapur-Mandla DSS profile across the Narmada-Son lineament in central India (Murty et al., 1998) has been updated based on the analysis of some short and discontinuous seismic wide-angle reflection phases. Three layers, with seismic velocities of 6.5–6.7, 6.35–6.40 and 6.8 km s^–1, and upper boundaries located approximately at 8, 17 and 22 km depth respectively, have been identified between the basement (velocity 5.9 km s^–1) and the uppermost mantle (velocity 7.8 km s^–1). The layer with 6.5–6.7 km s^–1 velocity is thin and is less than 2-km deep between the Narmada north (at Katangi) and south (at Jabalpur) faults. The upper crust shows a horst feature between these faults, which indicates that the Narmada zone acts as a ridge between two pockets of mafic intrusion in the upper crust. The Moho boundary, at 40–44 km depth and the intra-crustal layers exhibit an upwarp suggesting that the Narmada faults have deep origins, involving deep-seated tectonics. A smaller intrusive thickness between the Narmada faults, as compared to those beyond these faults, suggests that the intrusive activities on the two sides are independent. This further suggests that the two Narmada faults may have been active at different geological times. The seismic model is constrained by 2-D gravity modeling. The gravity highs on either side of the Narmada zone are due to the effect of the high velocity/high density mafic intrusion at upper crustal level.     

The thermal structure of subduction zones constrained by seismic imaging: Implications for slab dehydration and wedge flow

Thermal models of subduction zones often base their slab–wedge geometry from seismicity at mantle depths and, consequently, cannot be used to evaluate the relationship between seismicity and structure. Here, high-resolution seismic observations from the recent Broadband Experiment Across the Alaska Range (BEAAR) constrain, in a rare instance, the subducting slab geometry and mantle wedge temperature independent of seismicity. Receiver functions reveal that the subducting crust descends less steeply than the Wadati-Benioff Zone. Attenuation tomography of the mantle wedge reveals a high Q and presumably cold region where the slab is less than 80 km deep. To understand these two observations, we generate thermal models that use the improved wedge geometry from receiver functions and that incorporate temperature- and strain-rate-dependent olivine rheology. These calculations show that seismicity within the subducting crust falls in a narrow belt of pressure–temperature conditions, illuminating an effective Clapeyron slope of 0.1 K/MPa at temperatures of 450–750 °C. These conditions typify the breakdown of high-pressure hydrous minerals such as lawsonite and suggest that a single set of dehydration reactions may trigger intermediate-depth seismicity. The models also require that the upper, cold nose of the mantle wedge be isolated from the main flow in the mantle wedge in order to sustain the cold temperatures inferred from the Q tomography. Possibly, sufficient mechanical decoupling occurs at the top of the downgoing slab along a localized shear zone to 80 km depth, considerably deeper than inferred from thrust zone seismicity.     

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.     

亚稳态橄榄石对俯冲带负浮力的影响及其动力学意义

利用有限元方法,计算了不同俯冲速度及热传导系数俯冲岩石层的负浮力.在h-lt;400km和h-gt;740km的深度范围,低温高密俯冲带的负浮力随深度单调增加.因为尖晶石相到后尖晶石相的相变有负的克兰帕龙斜率,俯冲带冷的物质在660km间断面以下不到100km的深度范围内仍以低密度的低压物质相存在.所以在该深度范围,俯冲带受到了周围地幔阻止其插入下地幔的浮力作用.在400km-lt;h-lt;660km深度范围,由于受橄榄石相变的影响,不同计算模型负浮力随深度的变化有明显的不同.对于可能的相变动力学模型,亚稳态橄榄石的存在使负浮力随深度的增加值减小.其作用是不利于俯冲带直接穿透660km间断面,并引起该深度范围俯冲带沿俯冲方向压应力分布的变化.     

The transport of water in subduction zones

The transport of water from subducting crust into the mantle is mainly dictated by the stability of hydrous minerals in subduction zones. The thermal structure of subduction zones is a key to dehydration of the subducting crust at different depths. Oceanic subduction zones show a large variation in the geotherm, but seismicity and arc volcanism are only prominent in cold subduction zones where geothermal gradients are low. In contrast, continental subduction zones have low geothermal gradients, resulting in metamorphism in cold subduction zones and the absence of arc volcanism during subduction. In very cold subduction zone where the geothermal gradient is very low(?5?C/km), lawsonite may carry water into great depths of ?300 km. In the hot subduction zone where the geothermal gradient is high(25?C/km), the subducting crust dehydrates significantly at shallow depths and may partially melt at depths of 80 km to form felsic melts, into which water is highly dissolved. In this case, only a minor amount of water can be transported into great depths. A number of intermediate modes are present between these two end-member dehydration modes, making subduction-zone dehydration various. Low-T/low-P hydrous minerals are not stable in warm subduction zones with increasing subduction depths and thus break down at forearc depths of ?60–80 km to release large amounts of water. In contrast, the low-T/low-P hydrous minerals are replaced by low-T/high-P hydrous minerals in cold subduction zones with increasing subduction depths, allowing the water to be transported to subarc depths of 80–160 km. In either case, dehydration reactions not only trigger seismicity in the subducting crust but also cause hydration of the mantle wedge. Nevertheless, there are still minor amounts of water to be transported by ultrahigh-pressure hydrous minerals and nominally anhydrous minerals into the deeper mantle. The mantle wedge overlying the subducting slab does not partially melt upon water influx for volcanic arc magmatism, but it is hydrated at first with the lowest temperature at the slab-mantle interface, several hundreds of degree lower than the wet solidus of hydrated peridotites. The hydrated peridotites may undergo partial melting upon heating at a later time. Therefore, the water flux from the subducting crust into the overlying mantle wedge does not trigger the volcanic arc magmatism immediately.     

Geodynamic evolution of a forearc rift in the southernmost Mariana Arc

The southernmost Mariana forearc stretched to accommodate opening of the Mariana Trough backarc basin in late Neogene time, erupting basalts at 3.7–2.7 Ma that are now exposed in the Southeast Mariana Forearc Rift (SEMFR). Today, SEMFR is a broad zone of extension that formed on hydrated, forearc lithosphere and overlies the shallow subducting slab (slab depth ≤ 30–50 km). It comprises NW–SE trending subparallel deeps, 3–16 km wide, that can be traced ≥ ∼30 km from the trench almost to the backarc spreading center, the Malaguana‐Gadao Ridge (MGR). While forearcs are usually underlain by serpentinized harzburgites too cold to melt, SEMFR crust is mostly composed of Pliocene, low‐K basaltic to basaltic andesite lavas that are compositionally similar to arc lavas and backarc basin (BAB) lavas, and thus defines a forearc region that recently witnessed abundant igneous activity in the form of seafloor spreading. SEMFR igneous rocks have low Na₈, Ti₈, and Fe₈, consistent with extensive melting, at ∼23 ± 6.6 km depth and 1239 ± 40°C, by adiabatic decompression of depleted asthenospheric mantle metasomatized by slab‐derived fluids. Stretching of pre‐existing forearc lithosphere allowed BAB‐like mantle to flow along the SEMFR and melt, forming new oceanic crust. Melts interacted with pre‐existing forearc lithosphere during ascent. The SEMFR is no longer magmatically active and post‐magmatic tectonic activity dominates the rift.     

On the origin of El Chichón volcano and subduction of Tehuantepec Ridge: A geodynamical perspective

The origin of El Chichón volcano is poorly understood, and we attempt in this study to demonstrate that the Tehuantepec Ridge (TR), a major tectonic discontinuity on the Cocos plate, plays a key role in determining the location of the volcano by enhancing the slab dehydration budget beneath it. Using marine magnetic anomalies we show that the upper mantle beneath TR undergoes strong serpentinization, carrying significant amounts of water into subduction. Another key aspect of the magnetic anomaly over southern Mexico is a long-wavelength (∼ 150 km) high amplitude (∼ 500 nT) magnetic anomaly located between the trench and the coast. Using a 2D joint magnetic-gravity forward model, constrained by the subduction P–T structure, slab geometry and seismicity, we find a highly magnetic and low-density source located at 40–80 km depth that we interpret as a partially serpentinized mantle wedge formed by fluids expelled from the subducting Cocos plate. Using phase diagrams for sediments, basalt and peridotite, and the thermal structure of the subduction zone beneath El Chichón we find that ∼ 40% of sediments and basalt dehydrate at depths corresponding with the location of the serpentinized mantle wedge, whereas the serpentinized root beneath TR strongly dehydrates (∼90%) at depths of 180-200 km comparable with the slab depths beneath El Chichón (200-220 km). We conclude that this strong deserpentinization pulse of mantle lithosphere beneath TR at great depths is responsible for the unusual location, singularity and, probably, the geochemically distinct signature (adakitic-like) of El Chichón volcano.