We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated. 相似文献
We use lithosphere-scale gravity models to calculate gravity anomalies resulting from oceanic subduction, continental collision, slab steepening, delamination, and break-off. Local isostasy was assumed for determining vertical movements caused by mass changes related to these tectonic processes. Our results show that subduction is accompanied by basin subsidence on the upper plate caused by the heavy lithospheric root of the subducting slab. The basin evolution goes parallel with the slab evolution and shows considerable modifications when the processes at depth change (slab steepening, delamination, break-off). Characteristic gravity anomaly curves were acquired for the different tectonic scenarios. These curves together with other data (e.g. basin evolution on the upper and the lower plate) were used for the reconstruction of the tectonic evolution of the SE-Carpathians which includes Tertiary subduction and collision followed by slab steepening and delamination. 相似文献
Fluid-undersaturated experiments were conducted to determine the phase relations in the simplified peridotite system MgO-SiO2-H2O (MSH) at 11.0-14.5 GPa and 800-1400 °C. Stability relations of dense hydrous magnesium silicates (DHMSs) under fluid-undersaturated conditions were experimentally examined. From the fluid-absent experimental results, we retrieved thermodynamic data for clinohumite, phase A, phase E, and hydrous wadsleyite, consistent with the published data set for dry mantle minerals. With this new data set, we have calculated phase equilibria in the MSH system including dehydration reactions. The dehydration reactions calculated with lower water activities of 0.68-0.60 match the fluid-present experiments of this study above 11.0 GPa and 1000 °C, indicating that considerable amounts of silicate component were dissolved into the fluid phase. The calculated phase equilibria illustrate the stability of the post-antigorite phase A-bearing assemblages. In the cold subducting slab peridotite, phase A + enstatite assemblage survives into the transition zone, whereas phase A + forsterite + enstatite assemblage forms hydrous wadsleyite at a much shallower depth of about 360-km. The slab is subducted with no dehydration reactions occurring when entering the transition zone. The phase equilibria also show the high temperature stability of phase E. Phase E is stable up to 1200 °C at 13.5 GPa, a plausible condition in the mantle of relatively low temperature, i.e., beneath subduction zones. Phase E is a possible water reservoir in the mantle as well as wadsleyite and ringwoodite. 相似文献
In-situ Hf isotope analyses and U–Pb dates were obtained by laser ablation-MC-ICP-MS for a zircon-bearing mantle eclogite xenolith from the diamondiferous Jericho kimberlite located within the Archean Slave Province (Nunavut), Canada. The U–Pb zircon results yield a wide range of ages (2.0 to 0.8 Ga) indicating a complex geological history. Of importance, one zircon yields a U–Pb upper intercept date of 1989 ± 67 Ma, providing a new minimum age constraint for zircon crystallization and eclogite formation. In contrast, Hf isotope systematics for the same zircons display an intriguing uniformity, and corresponding Hf depleted mantle model ages range between 2.1 ± 0.1 and 2.3 ± 0.1 Ga; the youngest Hf model age is within error to the oldest U–Pb date.
The Jericho eclogites have previously been interpreted as representing remnants of metamorphosed oceanic crust, and their formation related to Paleoproterozoic subduction regimes along the western margin of the Archean Slave craton during the Wopmay orogeny. Hf isotope compositions and U–Pb results for the Jericho zircons reported here are in good agreement with a Paleoproterozoic subduction model, suggesting that generation of oceanic crust and eclogite formation occurred between 2.0 and 2.1 Ga. The slightly older Hf depleted mantle model ages (2.1 to 2.3 Ga) may be reconciled with this model by invoking mixing between ‘crustal’-derived Hf from sediments and more radiogenic Hf associated with the oceanic crust during the 2 Ga subduction event. This results in intermediate Hf isotope compositions for the Jericho zircons that yield ‘fictitiously’ older Hf model ages. 相似文献
Seismic tomography studies in the northeastern Japan arc have revealed the existence of an inclined sheet-like seismic low-velocity and high-attenuation zone in the mantle wedge at depths shallower than about 150 km. This sheet-like low-velocity, high-attenuation zone is oriented sub-parallel to the subducted slab, and is considered to correspond to the upwelling flow portion of the subduction-induced convection. The low-velocity, high-attenuation zone reaches the Moho immediately beneath the volcanic front (or the Ou Backbone Range) running through the middle of the arc nearly parallel to the trench axis, which suggests that the volcanic front is formed by this hot upwelling flow. Aqueous fluids supplied by the subducted slab are probably transported upward through this upwelling flow to reach shallow levels beneath the Backbone Range where they are expelled from solidified magma and migrate further upward. The existence of aqueous fluids may weaken the surrounding crustal rocks, resulting in local contractive deformation and uplift along the Backbone Range under the compressional stress field of the volcanic arc. A strain-rate distribution map generated from GPS data reveals a notable concentration of east–west contraction along the Backbone Range, consistent with this interpretation. Shallow inland earthquakes are also concentrated in the upper crust of this locally large contraction deformation zone. Based on these observations, a simple model is proposed to explain the deformation pattern of the crust and the characteristic shallow seismic activity beneath the northeastern Japan arc. 相似文献
We present a new three-dimensional SV-wave velocity model for the upper mantle beneath South America and the surrounding oceans, built from the waveform inversion of 5850 Rayleigh wave seismograms. The dense path coverage and the use of higher modes to supplement the fundamental mode of surface waves allow us to constrain seismic heterogeneities with horizontal wavelengths of a few hundred kilometres in the uppermost 400 km of the mantle.The large scale features of our tomographic model confirm previous results from global and regional tomographic studies (e.g. the depth extent of the high velocity cratonic roots down to about 200–250 km).Several new features are highlighted in our model. Down to 100 km depth, the high velocity lid beneath the Amazonian craton is separated in two parts associated with the Guyana and Guapore shields, suggesting that the rifting episode responsible for the formation of the Amazon basin has involved a significant part of the lithosphere. Along the Andean subduction belt, the structure of the high velocity anomaly associated with the sudbduction of the Nazca plate beneath the South American plate reflects the along-strike variation in dip of the subducting plate. Slow velocities are observed down to about 100 km and 150 km at the intersection of the Carnegie and Chile ridges with the continent and are likely to represent the thermal anomalies associated with the subducted ridges. These lowered velocities might correspond to zones of weakness in the subducted plate and may have led to the formation of “slab windows” developed through unzipping of the subducted ridges; these windows might accommodate a transfer of asthenospheric mantle from the Pacific to the Atlantic ocean. From 150 to 250 km depth, the subducting Nazca plate is associated with high seismic velocities between 5°S and 37°S. We find high seismic velocities beneath the Paraná basin down to about 200 km depth, underlain by a low velocity anomaly in the depth range 200–400 km located beneath the Ponta Grossa arc at the southern tip of the basin. This high velocity anomaly is located southward of a narrow S-wave low velocity structure observed between 200 and 500–600 km depth in body wave studies, but irresolvable with our long period datasets. Both anomalies point to a model in which several, possibly diachronous, plumes have risen to the surface to generate the Paraná large igneous province (LIP). 相似文献
Understanding the dominant force responsible for supercontinent breakup is crucial for establishing Earth's geodynamic evolution that includes supercontinent cycles and plate tectonics. Conventionally,two forces have been considered: the push by mantle plumes from the sub-continental mantle which is called the active force for breakup, and the dragging force from oceanic subduction retreat which is called the passive force for breakup. However, the relative importance of these two forces is unclear. Here we model the supercontinent breakup coupled with global mantle convection in order to address this question. Our global model features a spherical harmonic degree-2 structure, which includes a major subduction girdle and two large upwelling(superplume) systems. Based on this global mantle structure,we examine the distribution of extensional stress applied to the supercontinent by both subsupercontinent mantle upwellings and subduction retreat at the supercontinent peripheral. Our results show that:(1) at the center half of the supercontinent, plume push stress is ~3 times larger than the stress induced by subduction retreat;(2) an average hot anomaly of no higher than 50 K beneath the supercontinent can produce a push force strong enough to cause the initialization of supercontinent breakup;(3) the extensional stress induced by subduction retreat concentrates on a ~600 km wide zone on the boundary of the supercontinent, but has far less impact to the interior of the supercontinent. We therefore conclude that although circum-supercontinent subduction retreat assists supercontinent breakup, sub-supercontinent mantle upwelling is the essential force. 相似文献