Mineral grain size in the mantle affects fluid migration by controlling mantle permeability; the smaller the grain size, the less permeable the mantle is. Mantle shear viscosity also affects fluid migration by controlling compaction pressure; high mantle shear viscosity can act as a barrier to fluid flow. Here we investigate for the first time their combined effects on fluid migration in the mantle wedge of subduction zones over ranges of subduction parameters and patterns of fluid influx using a 2‐D numerical fluid migration model. Our results show that fluids introduced into the mantle wedge beneath the forearc are first dragged downdip by the mantle flow due to small grain size (<1 mm) and high mantle shear viscosity that develop along the base of the mantle wedge. Increasing grain size with depth allows upward fluid migration out of the high shear viscosity layer at subarc depths. Fluids introduced into the mantle wedge at postarc depths migrate upward due to relatively large grain size in the deep mantle wedge, forming secondary fluid pathways behind the arc. Fluids that reach the shallow part of the mantle wedge spread trench‐ward due to the combined effect of high mantle shear viscosity and advection by the inflowing mantle and eventually pond at 55–65 km depths. These results show that grain size and mantle shear viscosity together play an important role in focusing fluids beneath the arc.

We present a new approach for the lithosphere-asthenosphere interaction in subduction zones.The lithosphere is modeled as a Maxwell viscoelastic body sinking in the viscous asthenosphere. Both domains are discretized by the finite element method, and we use a staggered coupling method. The interaction is provided by a nonmatching interface method called the fictitious domain method. We describe a simplified formulation of this numerical technique and present 2-D examples and benchmarks. We aim at studying the effect of mantle viscosity on the cyclicity of slab folding at the 660 km depth transition zone. Such cyclicity has previously been shown to occur depending on the kinematics of both the overriding and subducting plates, in analog and numerical models that approximate the 660 km depth transition zone as an impenetrable barrier. Here we applied far-field plate velocities corresponding to those of the South-American and Nazca plates at present. Our models show that the viscosity of the asthenosphere impacts on folding cyclicity and consequently on the slab's dip as well as the stress regime of the overriding plate. Values of the mantle viscosity between 3 and 5 3 10 20 Pa s are found to produce cycles similar to those reported for the Andes, which are of the order of 30-40 Myr (based on magmatism and sedimentological records). Moreover, we discuss the episodic development of horizontal subduction induced by cyclic folding and, hence, propose a new explanation for episodes of flat subduction under the South-American plate.

In subductions where the slab stagnates at the 660-km mantle discontinuity, overriding plate kinematics largely controls slab dip and overriding plate tectonics. Although plates kinematics models suggest frequent velocity changes for most plates, the impact of temporal evolution of overriding plate velocity on subduction dynamics has been relatively little addressed. In the present study, we use 2-d numerical models to assess the effects of changes in overriding plate far-field velocity on subduction geometry and on the horizontal stresses transmitted to the overriding plate. When a change in overriding plate velocity arises during slab stagnation, slab dip evolves during a transient period, called adjustment-time, to reach a state in equilibrium with the new boundary conditions. The models predict a dependency of the adjustment-time on the value of velocity change and on several internal parameters (subducting plate density, thickness, and viscosity, and mantle viscosity). We estimate that the adjustment-times may be ∼ 10 − 35 Myrs in Nature, which suggests that most of present-day subduction zones with stagnating slabs might not be at a steady-state. Further, the models predict that changes in overriding plate velocity generate high temporary variations in the state of stresses of the plate.

No abstract

Slab dip controls the state of stress in an overriding plate and affects mountain building. Analog and numerical models have shown variations in tectonic regime induced by slab folding over the 660 km depth discontinuity zone in orthogonal convergence. Here using a three‐dimensional model of oblique subduction (30°) and accounting for free top surfaces, we show how slab folding generates an along‐strike slab dip segmentation, inducing variations in topography of the overriding plate. When the subducting plate begins to curve forward, the elevation height rises inland and varies along the trench from 5 km to 2 km. The Andes are a suitable natural zone to compare our results with because of its linear margin and well‐constrained plates kinematics. Thus, we provide a new explanation to the general decrease in elevation from the central to southern Andes, which still remains to be combined with other 3‐D mechanisms to explain the actual Andean topography.

No abstract

No abstract

Magmatism and volcanism transfer carbon from the solid Earth into the climate system. This transfer may be modulated by the glacial/interglacial cycling of water between oceans and continental ice sheets, which alters the surface loading of the solid Earth. The consequent volcaniccarbon fluctuations have been proposed as a pacing mechanism for Pleistocene glacial cycles. This mechanism is dependant on the amplitude and lag of the mid-ocean ridge response to sea-level changes. Here we develop and analyse a new model for that response, eliminating some questionable assumptions made in previous work. Our model calculates the carbon flux, accounting for the thermodynamic effect of mantle carbon: reduction of the solidus temperature and a deeper onset of melting. We analyse models forced by idealised, periodic sea level and conclude that fluctuations in melting rate are the prime control on magma and carbon flux. We also discuss a model forced by a reconstruction of eustatic sea level over the past 800 kyr. It indicates that peak-to-trough variations of magma and carbon flux are up to about 20% and 10% of the mean flux, respectively. Peaks in mid-ocean ridge emissions lag peaks in sea-level forcing by less than about 20 kyr and the lag could well be shorter. The amplitude and lag are sensitive to the rate of melt segregation. The lag is much shorter than the time it takes for melt to travel vertically across the melting region.