new inferences of mantle viscosity from joint inversion of long wavelength mantle convection and post glacial rebound data we perform joint inversions for mantle viscosity of geophysical observables associated with both mantle convection and glacial isostatic adjustment gia our data include coefficients of the non hydrostatic geoid up to degree 8 and decay times associated with post glacial relative sea level rsl variations in hudson bay and fennoscandia we find that both data sets may be reconciled using a single profile of mantle viscosity which is characterized by a significant increase with depth across this region this result weakens previous arguments that the mantle rheology has significant transient effects the viscosity profiles we obtain are found to reconcile the observed free air gravity anomaly over hudson bay

three dimensional instantaneous mantle flow induced by subduction 1 we conduct three dimensional subduction experiments by a finite element approach to study flow around slabs which are prescribed based on a transient stage of upper mantle subduction from a laboratory model instantaneous velocity field solutions are examined focusing on the toroidal vs poloidal components as a function of boundary conditions plate width and viscosity contrast between slab and mantle we show how the slab to mantle viscosity ratio determines the strength of toroidal flow and find that the toroidal flow component peaks for slab mantle viscosity ratios 100 independent of slab width or geometry

buckling instabilities of subducted lithosphere beneath the transition zone a sheet of viscous fluid poured onto a surface buckles periodically to generate a pile of regular folds recent tomographic images beneath subduction zones together with quantitative fluid mechanical scaling laws suggest that a similar instability can occur when slabs of subducted oceanic lithosphere encounter enhanced resistance to penetration at 660 km depth beneath subduction zones such as central america and java the fast anomalies in elastic wavespeed are too broad 400 km to be explained by simple compressional thickening of subducted lithosphere that is initially 50 100 km thick by contrast a periodic buckling mechanism can account for both the characteristic wedge shapes of the anomalies and for their widths which agree closely with those predicted by the scaling laws in addition to providing a plausible physical interpretation of tomographic observations beneath some subduction zones buckling instabilities have important consequences for other aspects of global geodynamics such as tectonic plate reconstructions interpretation of geoid anomalies and modeling of secular changes in earth s moment of inertia

dynamics of retreating slabs 2 insights from three dimensional laboratory experiments 1 a laboratory analogue of a three layer linear viscous slab upper mantle lower mantle system is established in a silicone putty honey and crystallized honey tank experiment the same setup as in the numerical investigation part 1 is used we focus on the interaction of the slab with the induced passive mantle flow by widely varying the mantle volume flux boundary conditions in our numerical experiments the lateral volume flux was set to zero in interpreting the results relative to the real earth the base of the box is taken as the bottom of the mantle convection system while the lateral boundaries may be associated with the presence of other nearby slabs dynamic force equilibrium assessed on the basis of an analytical review of forces is described for four different phases 1 the subduction initiation instability 2 the accelerating dynamic free fall phase of the slab 3 the dynamic interaction with the 660 km discontinuity and 4 a final phase of steady state trench retreat phase 3 is an important feature not observed in the numerical experiments this highly dynamic phase of interrupted trench retreat can therefore be attributed to boundary conditions on mantle volume flux on the basis of integration constants of force equilibrium in phases 2 and 4 we identify two different classes of volume flux one in which the lateral boundary can be considered open and the other class where it is closed closed boundary condition cases are obtained if any of the lateral box boundaries are 600 km away from the slab assuming a one to one relation between trench retreat and back arc spreading enigmatic observations of episodic opening of back arc basins can be explained by our experimental observations

mantle kinematics driving collisional subduction insights from analogue modeling abstract since several decades the processes allowing for the subduction of the continental lithosphere less dense than the mantle in a collision context have been widely explored but models that are based upon the premise that slab pull is the prominent driver of plate tectonics fail the india asia collision where several episodes of continental subduction have been documented constitute a case study for alternative views one of these episodes occurred in the early collision time within the asian plate where continental lithosphere not attached to any oceanic lithosphere subducted southward in front of the indian lithosphere during its northward subduction that followed the oceanic subduction of the tethys ocean this process known as collisional subduction has a counter intuitive behavior since the subduction is not driven by slab pull it has been speculated that the mantle circulation can play an important role in triggering collisional subduction but a detailed qualitative analysis of it is not available yet in this work we explore the southward subduction dynamics of the asian lithosphere below tibet by means of analogue experiments with the aim to highlight how the mantle circulation induces or responds to collisional subduction we found that during the northward oceanic subduction analogue of tethys subduction attached to the indenter indian analogue the main component of slab motion is driven vertically by its negative buoyancy while the trench rolls back in the mantle the convective pattern consists in a pair of wide convective cells on both sides of the slab but when the indenter starts to bend and plunge in the mantle trench motion reverses its advance transmits the far field forces to two upper plates asian analogues the more viscous frontal plate thickens and the less viscous hinterland plate which is attached to the back wall of the box subducts during this transition a pair of sub lithospheric convective cells is observed on both sides of the asian analogue slab driven by the shortening of the frontal plate it favors the initiation of the backwall plate subduction such subduction is maintained during the entire collision by a wide cell with a mostly horizontal mantle flow below tibet passively advecting the asian analogue slab experimental results suggest that once the tectonic far field force related to the forward horizontal motion becomes dominant upon the buoyancy forces trench advancing and the transmission of the tectonic force to the upper and backwall plates are promoted this peculiar condition triggers the subduction of the backwall plate despite it is light and buoyant

numerical modeling of stress and topography coupling during subduction inferences on global vs regional observables interpretation abstract the mechanical strength of the oceanic lithospheres exerts a fundamental control on the dynamics and tectonics of convergent margins here we use numerical models of coupled subducting upper plates and mantle to investigate the impact of the subducting lithosphere properties on the parameters commonly measured at subduction zones the models predict a correlation among parameters such as subducting and convergence velocities slab dip and curvature and stress coupling with the upper plates such signature emerges in a dynamic system through the slab buoyancy similar trends are found when assessing the impact of laboratory constrained non linear rheologies uncertainties on the mechanical strength and mantle flow although some features are more affected when the role of upper plate thickness is tested we find a positive correlation of upper plate subsidence and fore arc topography with the buoyancy of the down going plate however plasticity in the crust limits the development of topography becoming increasingly relevant with thinner upper plates so that no correlation can be established similar control is found on the width of the long wavelength dynamic subsidence due to flexure which increases non linearly with the upper plate thickness the models constrain relations among parameters within the dynamic system driven that allow inferences on the observable parameters the comparisons supports the idea that the subduction dynamics emerges as a global phenomenon in plate motions bending and dynamic subsidence yet andean type margin tectonics is best addressed in a regional context

the dynamics of back arc extension an experimental approach to the opening of the tyrrhenian sea the e w opening tyrrhenian sea developed after the cretaceous palaeogene alpine collision nearly perpendicular to the motion of the african plate as a back arc of the adria ionian westward subduction three driving mechanisms have been proposed to explain the dynamic evolution of the tyrrhenian apennine system 1 the northward indentation of the african plate 2 the retreating subduction of the adria ionian lithosphere and 3 the gravitational collapse of the alpine post collisional wedge in order to define the relative contribution of each of these mechanisms in the neogene dynamic of the tyrrhenian apennine system we performed 3 d laboratory experiments in which we simulated a retreating subduction process in a compressional regime oriented perpendicularly to the direction of subduction in this framework we also tested the influence of the gravitational collapse of the overriding plate experiments were constructed using dry sand and silicone putties to simulate brittle upper crust and ductile lower crust upper mantle respectively these layers floated on a high density low viscosity glucose syrup which simulated the asthenosphere the main conclusion of our experiments is that large scale continental extension similar to that observed in the tyrrhenian area could be reproduced perpendicular to the shortening direction induced by the indentation of the african plate in this framework extensional processes are indeed possible if the trench retreat velocity is higher than the rate of shortening induced by the advancing african plate our experimental results indicate that this high trench retreat velocity could be explained by the coexistence of the gravitational collapse of the post alpine wedge with a slab pull process linked to the retreating subduction of the adria ionian plate while the first mechanism is predominant in the northern tyrrhenian area the second one seems to be important in the latest stage of extension and oceanic accretion of the southern tyrrhenian area

predicting plate velocities with mantle circulation models 1 we predict plate motions from a comprehensive inversion of theoretical estimates of tectonic forces in order to evaluate the relative importance of these and the uncertainties of such models plate driving forces from the mantle are calculated using global flow models that are driven by tomography and subduction derived density fields observed and predicted plate velocities agree well for a variety of models leading to varied conclusions about the relative importance of forces the dominance of the subduction related density pattern in the mantle is confirmed it appears that p wave models do not satisfactorily image all of the slab associated anomalies in the upper mantle furthermore lower mantle structure always improves the plate motion fit with respect to models that are based on upper mantle anomalies and lithospheric thickening only we show that the average torques from the lower mantle scale with the radial flow through the 660 km phase transition the amplitude of the lower mantle torques will be significant for a range of models if there is mass flux through 660 km we also evaluate parameterized edge forces and find that the additional inclusion of such torques does not significantly improve the model fit the main reason for the nonuniqueness of the inversions is plate boundary geometry since all plate motions are dominated by the trench ridge system and plates move from ridges to trenches

geodynamic models of continental subduction and obduction of overriding plate forearc oceanic lithosphere on top of continental crust continental subduction takes place in the final stage of subduction when all oceanic lithosphere is consumed and continental passive margin is pulled into the mantle when the overriding plate is oceanic dense forearc oceanic lithosphere might be obducted onto light continental crust forming an ophiolite tethyan style ophiolite obduction four dimensional dynamic analog subduction models have been constructed to evaluate the mechanical feasibility of continental subduction and forearc oceanic lithosphere obduction on top of continental crust the roles of continental crust thickness passive margin length subducting lithosphere thickness and overriding plate thickness were investigated to determine the maximum continental subduction depth maximum forearc obduction distance and forearc deformation during continental subduction our buoyancy driven experiments indicate that deep continental subduction occurs in most circumstances down to 560 km and that obduction of dense oceanic forearc lithosphere on top of light continental crust is mechanically feasible maximum obduction distances are relatively small 26 37 km but are sufficient to explain obduction of short ophiolite sheets such as observed in new caledonia when including the thin 5 10 km thick accretionary wedge of off scraped deep sea sediments oceanic crust and mantle then maximum obduction distances are much larger 60 160 km sufficient to account for the obducted northland allochthon in new zealand results indicate that increasing continental crust thickness decreases continental subduction depth whereas increasing passive margin length and subducting lithosphere thickness increases continental subduction depth notably during continental subduction backarc extension continues while forearc deformation shortening increases moderately compared to the preceding phase of normal oceanic subduction

effects of plate bending and fault strength at subduction zones on plate dynamics for subduction to occur plates must bend and slide past overriding plates along fault zones because the lithosphere is strong significant energy is required for this deformation to occur energy that could otherwise be spent deforming the mantle we have developed a finite element representation of a subduction zone in which we parameterize the bending plate and the fault zone using a viscous rheology by increasing the effective viscosity of either the plate or the fault zone we can increase the rates of energy dissipation within these regions and thus decrease the velocity of a plate driven by a given slab buoyancy we have developed a simple physical theory that predicts this slowing by estimating a convecting cell s total energy balance while taking into account the energy required by inelastic deformation of the bending slab and shearing of the fault zone the energy required to bend the slab is proportional to the slab s viscosity and to the cube of the ratio of its thickness to its radius of curvature the distribution of dissipation among the mantle lithosphere and fault zone causes the speed of a plate to depend on its horizontal length scale using the observation that earth s plate velocities are not correlated to plate size we can constrain the lithosphere viscosity to be between 50 and 200 times the mantle viscosity with higher values required if the fault zone can support shear tractions 50 mpa over 300 km these subduction zone strengths imply that the mantle fault zone and lithosphere dissipate about 30 10 and 60 of a descending slab s potential energy release if the slab is 100 km thick the lithospheric component is highly dependent on slab thickness it is smaller for thin plates but may be large enough to prevent bending in slabs that can grow thicker than 100 km ubduction zone strength should be more stable than mantle viscosity to changes in mantle temperature so the controlling influence of subduction zones could serve to stabilize plate velocities over time as the earth cools because the details of convergent plate boundaries are so important to the dynamics of plate motion numerical models of mantle flow should treat subduction zones in a realistic way

geotechnical modeling at the city scale using statistical and geostatistical tools the pessac case france abstract in the context of the riviera project the building of a 3d geotechnical model at the city scale pessac france has been undertaken from several hundreds of boreholes and geotechnical tests it is first shown how the combination of the lithological information and of geotechnical results can improve thanks to bayesian statistics the knowledge of mechanical characteristics in the various alluvial terraces which can be encountered in this area secondly the upper and lower limits of the 3d model at the city scale are computed by improving an initial digital elevation model for the upper limit and by kriging under inequality constraints for the lower limit these limits border quaternary formations which are of interest for geotechnical applications in a third stage it is focused on the spatial modelling of the pressuremeter modulus the sequential indicator simulation method enables to obtain the spatial probability of occurrence of a given pressiometer modulus class coupled with other information the analysis of these statistical and geostatistical models makes possible to develop decision support tools such as to localise for instance areas more prone to the clay shrinkage swelling hazard

radial profile of mantle viscosity results from the joint inversion of convection and postglacial rebound observables we present new inferences of the radial profile of mantle viscosity that simultaneously fit long wavelength free air gravity harmonics associated with mantle convection and a large set of decay times estimated from the postglacial uplift of sites within previously glaciated regions hudson bay arctic canada and fennoscandia the relative sea level variation at these latter sites is constrained by age height pairs obtained by geological survey rather than the subjective trends which are commonly used in glacial isostatic adjustment gia studies our viscosity inferences are generated using two approaches first we adopt a relative viscosity profile which is known to provide a good fit to the free air gravity harmonics and determine an absolute scaling which yields a best fit to the gia decay time constraints second we perform an iterative nonlinear joint inversion of the two data sets in both cases our inferred profiles are characterized by a significant increase of viscosity 2 orders of magnitude with depth to values of 1022 pa s in the bottom half of the lower mantle the new viscosity profiles are shown to satisfy constraints based on the postglacial uplift of both fennoscandia the classic haskell 1935 number and hudson bay which have commonly been invoked to argue for an isoviscous mantle furthermore the models are used to predict a set of long wavelength signatures of the gia process these include predictions of gia induced variations in 1 the length of day over the late holocene period 2 the earth s precession constant and obliquity over the last 2 6 myr and 3 the present day zonal harmonics of the geopotential jl l 7 the predictions 1 and 3 bound the late holocene and ongoing mass flux between the large polar ice sheets greenland and antarctic and the global oceans to small values 0 4 mm yr equivalent eustatic sea level rise

rupture dynamics with energy loss outside the slip zone 1 energy loss in a fault damage zone outside the slip zone contributes to the fracture energy that determines rupture velocity of an earthquake a nonelastic two dimensional dynamic calculation is done in which the slip zone is modeled as a fault plane and material off the fault is subject to a coulomb yield condition in a mode 2 crack like solution in which an abrupt uniform drop of shear traction on the fault spreads from a point coulomb yielding occurs on the extensional side of the fault plastic strain is distributed with uniform magnitude along the fault and it has a thickness normal to the fault proportional to propagation distance energy loss off the fault is also proportional to propagation distance and it can become much larger than energy loss on the fault specified by the fault constitutive relation the slip velocity function could be produced in an equivalent elastic problem by a slip weakening friction law with breakdown slip dc increasing with distance fracture energy g and equivalent dc will be different in ruptures with different initiation points and stress drops so they are not constitutive properties they are determined by the dynamic solution that arrives at a particular point peak slip velocity is however a property of a fault location nonelastic response can be mimicked by imposing a limit on slip velocity on a fault in an elastic medium

india eurasia collision chronology has implications for crustal shortening and driving mechanism of plates the motion of the indian plate is determined in an absolute frame of reference and compared with the position of the southern margin of eurasia deduced from palaeomagnetic data in tibet the 2 600 900 km of continental crust shortening observed is shown to have occurred in three different episodes subduction of continental crust intracontinental thrusting and internal deformation and lateral extrusion the detailed chronology of the collision and plate reorganizations in the indian and pacific oceans supports the hypothesis that slab pull is a dominant driving mechanism of plate tectonics

influence of surrounding plates on 3d subduction dynamics 1 our 3d modelling study shows that the presence of lithospheric plates around a subducting plate has a significant influence on subduction dynamics in particular on trench retreat rate slab dip and lateral shortening of the subducting plate neighbouring plates prevent unrealistic plate behaviour with no need for complex rheologies because at the earth s surface plates form a continuous shell they should not be neglected

the role of the overriding plate thermal state on slab dip variability and on the occurrence of flat subduction 1 slab dip varies significantly both between different and along single subduction zones provided that old subducting plates are colder and denser than young plates variations in the slab dip should correlate with slab age however recent statistical analyses do not show this expected correlation we present the results of non newtonian numerical dynamic models where subduction is driven by means of a kinematic boundary condition we systematically vary the age of both the overriding and subducting plates in order to test these effects on the slab dip at different depth ranges we find that colder overriding plates result in shallower slab dips and episodes of flat slab subduction as a result of the increased suction force in the mantle wedge the influence of the thermal state of the overriding plate on slab dip is shown here to be more important than that of the age of subducting lithosphere modeling results are qualitatively compared to the large dip variability of the cocos slab including a flat slab segment we suggest that this variability is likely related to the change of the thermal state of the overriding plates with flat subduction occurring under cold lithosphere in southwestern mexico and steep subduction under the warmer lithosphere of the northwestern caribbean plate

differential rotation between lithosphere and mantle a consequence of lateral mantle viscosity variations the description of plate motions in the so called hotspot reference frame introduces a global rotation of the lithosphere with respect to the mantle this rotation called toroidal field of degree 1 is roughly westward it reaches an amplitude of about 2 cm yr and has been consistently found in the different generation of plate tectonic models various authors have tried to relate this observation to the deceleration of the earth s rotation to polar wander or to tidal drag however these different physical mechanisms cannot explain the requested amplitude in this paper we compare the values of this rotation vector using different relative plate motion models expressed in the hotspot reference frame in a model earth with lateral viscosity variations a differential rotation is predicted the observed net lithospheric rotation is consistent with the dynamics of a model earth where the asthenospheric viscosity below the oceans is at least one order of magnitude lower than underneath the continents this relative westward drift of the lithosphere may account for the significant structural differences between east or west dipping subduction zones

a laboratory investigation of effects of trench migration on the descent of subducted slabs a laboratory investigation of viscous slabs subducted from a migrating trench reveals a range of possible behaviours the slab dip in a uniform mantle is found to be steady or oscillatory depending on the rates of descent and trench migration in addition density and viscosity interfaces are used to model the increased resistance to sinking of stiff slabs through the seismic discontinuity at a depth of about 660 km in earth s mantle if the slab is denser than the lower layer and its dip in the upper layer is steady it can continue to descend through the lower layer in a tabular form and in either a steady or oscillatory manner or it can be laid horizontally on the interface for a distance before descending in a chain of diapirs resulting from gravitational instability at the interface if the slab dip in the upper layer is unstable the slab sinks into the lower layer in a chain of large diapirs at a spacing determined by the frequency of oscillations set by instability of the slab within the upper layer the style of penetration depends on the trench migration speed and the ratio of sinking velocities in the two layers estimates for mantle slabs indicate that they may range across the major regime transitions the experimental system provides a gross simplification of mantle conditions but the results are a testimony to the possibility of a range of complex behaviour of subducted lithosphere they also indicate that the relationship between the structural features in the lower mantle revealed by seismic imaging and present day tectonic processes at the surface may not be obvious tomographic images of western pacific subduction zones and data for the migration rates of associated trenches suggest a dependence of slab behaviour on migration rate similar to that seen in the experiments

the bending mechanics in a dynamic subduction system constraints from numerical modelling and global compilation analysis abstract lithospheric plates bend at subduction zones where the vertical motions of the slabs are converted to surface plate motions to understand the mechanics of plate bending we derive scaling laws for the deflection at the margin i e radius and dip from numerical models of a subducting viscoelastic plate in such dynamic system we find that the buoyancy and the stiffness of the plates control the radius and the dip as well as the plate motions toward the trench this mechanical model successfully predicts the curvature of published three dimensional laboratory and numerical models for a thorough comparison with the observable we have also implemented forces additional to the slab pull such as the suction force and far field stresses by increasing or resisting the torque applied at the trench by the slab these forces can largely rearrange the dip and the radius of slabs and the inherent plate motions although they do not alter the observed anticorrelation between radius and dip similar inverse correlation relationship and dip radius ranges are shown by most of the subduction zones analysed from a global compilation radii in the range of 100 350 km and dips of 30 70 for slabs that extends to the bottom of the upper mantle are compatible with the models and allow estimating an average lithospheric viscosity contrast of 200 in the bending with respect to the ambient mantle radii and dips outside of this range are in good agreement with the trends and the magnitudes of models that include suction and far field forces in all these subduction zones the correlation between dip radius and plate velocity is found to be compatible with that of the models showing how relevant bending is for the dynamics of earth

the development of slabs in the upper mantle insights from numerical and laboratory experiments we have performed numerical and laboratory experiments to model subduction of oceanic lithosphere in the upper mantle from its beginnings as a gravitational instability to the fully developed slab a two dimensional finite element code is applied to model newtonian creep in the numerical experiments scaled analog media are used in the laboratory a sand mixture models the brittle crust silicone putty simulates creep in the lower crust and mantle lithosphere and glucose syrup is the asthenosphere analog both model approaches show similar results and reproduce first order observations of the subduction process in nature based on density and viscosity heterogeneities in a stokes flow model subduction nucleates slowly and a pronounced slab forms only when the viscosity contrast between oceanic plate and mantle is below a threshold we find that the subduction velocity and angle are time dependent and increase roughly exponentially over tens of millions of years before the slab reaches the 670 km discontinuity the style of subduction is controlled by the prescribed velocity of convergence the density contrast between the plates and the viscosity contrast between the oceanic plate and the mantle these factors can be combined in the buoyancy number f which expresses the ratio between driving slab pull and resisting viscous dissipation in the oceanic plate variations in f control the stress in the plates the speed and the dip of subduction and the rate of trench retreat reproducing the contrasting styles of subduction observed in nature the subduction rate is strongly influenced by the work of bending the lithosphere as it subducts

control of seafloor aging on the migration of the izu bonin mariana trench article i nfo recent global kinematic studies reveal that most of the trenches roll back but a significant number of them advance toward the upper plate those advancing trenches are mostly located in the western pacific and correspond to the subduction of very old mesozoic oceanic lithosphere while retreating trenches are commonly explained by the slab pull action of the descending lithosphere the origin of advancing trenches is still debated since this relationship is dependent upon the adopted reference frame we select region where geological studies show the variability of trench migration style with time the izu bonin mariana ibm region represents a key example the detailed reconstruction of the trench migration of the ibm subducting system reveals that after a long episode of asymmetric rollback the ibm trench recently started advancing we propose that this change from retreating to advancing trench mode results from the subduction of progressively older and stiffer lithospheric material we test this hypothesis by means of two dimensional 2 d numerical models reproducing the effects of the lithospheric aging during subduction the result of our numerical tests shows that the entrance of old and stiff lithosphere forces the trench to advance because the increasing stiffness of the slab prevents the slab to unbend once it has subducted we adapt this physical result to the ibm evolution showing that the switch from trench retreat to trench advance can reconcile the shape of ibm slab as previously suggested finally we discuss the possibility that the subduction of the cretaceous lithosphere might have triggered the change of the motion of the pacific plate around 5 ma

mantle rheology and the scaling of bending dissipation in plate tectonics 1 plate tectonics on earth involves the bending deformation of plates at subduction zones and because plates are generally considered to be stiff owning to the rheology of mantle minerals the role of energy dissipation by plate bending in the global energy balance has been frequently debated in the recent literature here we consider how bending dissipation should scale with slab parameters such as dip angle plate age the radius of curvature and plate velocity by systematically exploring the parameter space with instantaneous stokes flow calculations we derive the scaling of bending dissipation for a range of mantle viscosity functions including pseudoplastic rheology with olivine flow laws our results indicate that as we move away from the isoviscous case the scaling gradually deviates from what has commonly been assumed in previous studies most notably for the radius exponent which exhibits more than threefold reduction and even a sign reversal in some cases these modifications in scaling exponents originate in the complication of the deformation field caused by viscosity variations within the bending plate approximating the lithospheric rheology by a single effective viscosity in the dynamical models of subduction has been a common practice but we suggest that such approximation may limit the geological relevance of modeling studies in particular when estimating the significance of bending dissipation

plate kinematics slab shape and back arc stress a comparison between laboratory models and current subduction zones a combination of statistical studies on present day subduction zones and three dimensional 3d laboratory models is performed with the aim to clarify the way that plate kinematics control the geometry of the slab and the overriding plate deformation in subduction zones in 3d laboratory models the analogue of a two layer linearly viscous lithosphere upper mantle system is achieved by means of silicon putty glucose syrup tank experiment the subducting and overriding plate velocities are systematically changed by exploring the variability field of natural plate kinematics both statistical and modelling approaches recognize the importance of overriding plate motion on subduction process behavior 1 trenches migrate at a rate close to the overriding plate motion but always move slower than the overriding plates the mechanism at work is a direct consequence of slab anchoring opposed by both lithosphere and mantle viscous resistance and is responsible for overriding plate deformation and slab geometry variability 2 an overriding plate shortens when the overriding plate moves toward the trench and conditions that are favourable for overriding plate extension are created when the overriding plate moves away from the trench 3 shallow and steep dips are found if the overriding plate moves toward and away from the trench respectively

impact of regional mantle flow on subducting plate geometry and interplate stress insights from physical modelling summary physical models of subduction investigate the impact of regional mantle flow on the structure of the subducted slab and deformation of the downgoing and overriding plates the initial mantle flow direction beneath the overriding plate can be horizontal or vertical depending on its location with respect to the asthenospheric flow field imposed mantle flow produces either over or underpressure on the lower surface of the slab depending on the initial mantle flow pattern horizontal or vertical respectively overpressure promotes shallow dip subduction while underpressure tends to steepen the slab horizontal mantle flow with rates of 1 10 cm yr 1 provides sufficient overpressure on a dense subducting lithosphere to obtain a subduction angle of 60 while the same lithospheric slab sinks vertically when no flow is imposed vertical drag force due to downward mantle flow exerted on a slab can result in steep subduction if the slab is neutrally buoyant but fails to produce steep subduction of buoyant oceanic lithosphere the strain regime in the overriding plate due to the asthenospheric drag force depends largely on slab geometry when the slab dip is steeper than the interplate zone the drag force produces negative additional normal stress on the interplate zone and tensile horizontal stress in the overriding plate when the slab dip is shallower than the interplate zone an additional positive normal stress is produced on the interplate zone and the overriding plate experiences additional horizontal compressive stress however the impact of the mantle drag force on interplate pressure is small compared to the influence of the slab pull force since these stress variations can only be observed when the slab is dense and interplate pressure is low

laboratory experiments on the structure of subducted lithosphere abstract mantle seismic tomography images indicate a variety of shapes for slabs of oceanic lithosphere subducted into the transition zone in the vicinity of the 660 km discontinuity both horizontally deflected slabs and nearly vertical slabs have been inferred at subduction zones on the pacific and indian plates here we determine experimentally the behavior of a cold negatively buoyant and highly viscous slab incident on a fluid interface with a viscosity and a density increase the experiments are scaled to mantle conditions with realistic subduction peclet number ratio of trench migration to subduction velocities and physical properties of the interface a wide range of slab deformation styles is observed when slight changes in the parameters are imposed sinking slab stagnant slab spreading slab sinking pile and stagnant pile represent the five different modes identified we show that these deformation modes are governed by two velocity ratios which characterize the horizontal and the vertical components of the slab velocity near the interface slab penetration into the lower layer takes the form of folded slab piles for large values of the velocity ratios and stagnant behavior is promoted by rapid retrograde trench motion the sensitivity of slab behavior to transition zone density and viscosity changes and to trench migration can account for the coexistence of different modes of deformation in subducted lithosphere

episodic back arc extension during restricted mantle convection in the central mediterranean abstract in the central mediterranean during the last 80 myr africa has been slowly converging toward eurasia with an average relative velocity of less than 1 cm yr 1 for 50 myr this slow convergence led to the build up of the alpine orogenic belt and to the initiation of consumption of the former tethys ocean in the last 30 myr however the whole central western mediterranean has been swept by rapid episodes of trench migration up to 6 cm yr 1 and back arc opening which consumed the whole former tethys ocean we combine plate tectonic history geological timing and laboratory modelling to reconstruct the subduction history and the evolution of the central mediterranean over the last 80 myr we find that the dynamic evolution of the subducting oceanic lithosphere can reconcile the rapid episodic back arc migration with the slow african convergence but only if the subduction process is restricted to the upper mantle

quantifying the energy dissipation of overriding plate deformation in three dimensional subduction models in a subduction system the force and the energy required to deform the overriding plate are generally thought to come from the negative buoyancy of the subducted slab and its potential energy respectively such deformation might involve extension and back arc basin formation or shortening and mountain building how much of the slab s potential energy is consumed during overriding plate deformation remains unknown in this work we present dynamic three dimensional laboratory experiments of progressive subduction with an overriding plate to quantify the force fopd that drives overriding plate deformation and the associated energy dissipation rate opd and we compare them with the negative buoyancy fbu of the subducted slab and its total potential energy release rate bu respectively we varied the viscosity ratio between the plates and the sublithospheric upper mantle with sp um 157 560 and the thickness of the overriding plate with top 0 5 2 5 cm scaling to 25 125 km in nature the results show that fopd fbu has average values of 0 5 2 0 with a maximum of 5 3 and opd bu has average values of 0 05 0 30 with a maximum of 0 41 the results indicate that only a small portion of the negative buoyancy of the slab and its potential energy are used to deform the overriding plate our models also suggest that the force required to deform the overriding plate is of comparable magnitude as the ridge push force furthermore we show that in subduction models with an overriding plate bending dissipation at the subduction zone hinge remains low 3 15 during steady state subduction

plate motions slab dynamics and back arc deformation back arc extension or compression is often related to a particular behavior of the slab dip change and or forward backward migration with respect to the upper plate slabs act either passively when anchored in the deep mantle whereas back arc deformation accommodates the backward or forward motion of the upper plate with respect to the trench or actively when slab pull increasing with slab age forces the trench to rollback in addition to these two mechanisms i e slab anchor and slab rollback numerous observations support the existence of dynamic mantle flow that can exert an overpressure on one side of the slab causing its forward or rearward migration with respect to the arc based on a compilation of upper plate absolute motion trench absolute motion back arc deformation rate upper plate strain regime and slab age for all oceanic subduction zones excluding any kind of collision with continents arcs or plateaus we have examined how the combined effects of these parameters can account for the observed back arc deformations our main results are 1 a global correlation exists between upper plate absolute motion and back arc deformation i e back arc extension when upper plate retreats and vice versa 2 there are as many advancing trenches as retreating ones with trench motion globally limited to 50 mm y 1 furthermore there is no positive correlation between trench retreat and slab age and 3 upper plate absolute motion often fails to explain the back arc deformation rates and the trench motions observed at several subduction zones as tonga new hebrides sandwich or ryukyu we propose that the trench migration of these subduction zones are on the influence of mantle flows we conclude that spontaneous trench rollback related to slab pull is negligible with respect to upper plate motion back arc deformation regime is mostly controlled by the upper plate absolute motion relative to a partly anchored slab trench motion is limited to 50 mm y 1 even if locally mantle flows force slab to move

mantle flow and dynamic topography associated with slab window opening insights from laboratory models abstract we present dynamically self consistent mantle scale laboratory models that have been conducted to improve our understanding of the influence of slab window opening on subduction dynamics mantle flow and associated dynamic topography over geological time scales the adopted setup consists of a two layer linearly viscous system simulating the subduction of a fixed plate of silicone lithosphere under negative buoyancy in a viscous layer of glucose syrup mantle our experimental setting is also characterized by a constant width rectangular window located at the center of a laterally confined slab modeling the case of the interaction of a trench parallel spreading ridge with a wide subduction zone we found that the opening of a slab window does not produce consistent changes of the geometry and the kinematics of the slab on the contrary slab induced mantle circulation quantified both in the vertical and horizontal sections using the feature tracking image analysis technique is strongly modified in particular rollback subduction and the opening of the slab window generate a complex mantle circulation pattern characterized by the presence of poloidal and toroidal components with the importance of each evolving according to kinematic stages mantle coming from the oceanic domain floods through the slab window indenting the supra slab mantle zone and producing its deformation without any mixing between mantle portions the opening of the slab window and the upwelling of sub slab mantle produce a regional scale non isostatic topographic uplift of the overriding plate that would correspond to values ranging between ca 1 and 5 km in nature assuming that our modeling results can be representative of the natural behavior of subduction zones we compared them to the tectonics and volcanism of the patagonian subduction zone we found that the anomalous backarc volcanism that has been developing since the middle miocene could result from the lateral flow of sub slab mantle and that part of the patagonian uplift could be dynamically supported

subducted slabs and the geoid constraints on mantle rheology and flow geoid anomalies are primarily the result of the density contrasts driving mantle convection and plate motion the total geoid anomaly resulting from a given density contrast in a convecting earth is affected by the mass anomalies associated with the flow induced deformation of the upper surface and internal compositional boundaries as well as by the density contrast itself these boundary deformations and hence the total gravity field depend on the radial distribution of effective viscosity if the internal density contrasts can be estimated as is the case for subducted slabs useful constraints can be placed on the depth and on the variation of viscosity with depth of the convecting system the degree 4 9 components of the observed long wavelength geoid are highly correlated with those predicted by a density model for seismically active subducted slabs the positive sign of the correlation requires that the effective viscosity increases with depth by a factor of 30 or more the amplitude of the correlation cannot be explained by the density contrasts associated with just the seismically active parts of subducted slabs however the amplitude can be explained if the density contrasts associated with subduction extend into the lower mantle or if old lithosphere is piled up at the base of the upper mantle beneath subduction zones to a thickness in excess of 350 km over horizontal distances of thousands of kilometers mantlewide convection in a mantle that has a viscosity increasing with depth provides a simple explanation of the long wavelength geoid anomalies over subduction zones

late cenozoic right lateral strike slip faulting in southern tibet field observations collected during the sino french himalaya tibet cooperation project 1980 1982 show that right lateral faulting between southern and northern tibet roughly along the chord joining the eastern and western syntaxes of the himalayan arc is a major component of the active tectonics north of the himalayas just north of the well developed rifts of southern tibet we studied three right lateral 100 200 km long strike slip faults forming an en echelon array along that chord evidence for earthquake surface breaks was found on each of these faults for the m 8 beng co or damxung strike slip earthquake of november 18 1951 rupture parameters could be constrained to be 8 2 m and l 90 km offset glacial and postglacial morphological markers 105 years and 104 years b p imply displacement rates of the order of 1 2 cm yr on the three faults the en echelon fault zone appears to extend westward 2000 km across the plateau to the karakorum fault and is thus referred to as the karakorum jiali fault zone kjfz fast rates of slip on it and on faults along the northern rim of tibet appear to contrast with the slower rates of deformations in the interior of the plateau our estimates of slip rates on the strike slip and normal faults of southeastern tibet suggest that eastward extrusion of north tibet between the right lateral kjfz and the left lateral altyn tagh and kun lun faults absorbs at least 30 of the present convergence between india and asia we propose an instantaneous velocity model of the active tectonics of asia compatible with such faulting rates and geometries evidence for tertiary right lateral strike slip movements in southern tibet and around the eastern syntaxis of the himalayas is also found in the structural record offsets of various large scale geological features imply that these movements have been large

great earthquakes and slab pull interaction between seismic coupling and plate slab coupling great earthquakes the few largest earthquakes that account for most of the earth s seismic energy release have occurred at only a few subduction zones around the world strong locking or seismic coupling of the interface between plates at certain subduction zones is often invoked to explain these great earthquakes although past studies have correlated strong seismic coupling with a compressional stress environment that is characterized by back arc compression and caused by trenchward motion of the overriding plate the consequences of this compressional environment for the tectonic forces that drive global plate motions are not yet clear to examine these consequences we compared subduction zone earthquake magnitudes to tectonically constrained estimates of the degree to which each slab transmits its excess weight as a direct pull force on a subducting plate at seismically uncoupled subduction zones that generate only moderate sized earthquakes we find that slabs must transmit nearly their entire upper mantle weight as a pull force on the subducting plate at seismically coupled subduction zones that produce great earthquakes however we find that slabs must be nearly completely detached from their subducting plates this suggests that slabs subducting in a compressional environment experience stress induced weakening that prevents the effective transmission of the slab pull force convergent mantle flow above a descending slab that becomes decoupled from its surface plate may induce additional surface compression that further locks the subduction zone and leads to additional slab decoupling and detachment the resulting redistribution of plate driving forces may be responsible for rapid changes in plate motion c 2003 elsevier b v all rights reserved

parameters controlling dynamically self consistent plate tectonics and single sided subduction in global models of mantle convection recent advances in numerical modeling allow global models of mantle convection to more realistically reproduce the behavior at convergent plate boundaries in particular the inclusion of a free surface at the outer boundary has been shown to facilitate self consistent development of single sided subduction this allows for a more extensive study of subduction in the context of global mantle convection as opposed to commonly used regional models our first study already indicated important differences between mantle convection with single sided subduction and mantle convection with double sided subduction here we further investigate the effect of various physical parameters and complexities on inducing earth like plate tectonics and its evolution in time results reinforce the previous finding that using a free surface instead of a free slip outer boundary dramatically changes subduction style with free surface cases displaying many episodes of single sided subduction which leads to more realistic slab dip stress state trench retreat rate and slab induced mantle flow longevity of single sided subduction is promoted by a layer of hydrated crust with a low yield strength to lubricate the subduction channel a low viscosity asthenosphere and a high strength of the slab determined by a combination of high diffusion creep viscosity and intermediate friction coefficient although its effective viscosity is in the observationally constrained range in the bending region the time evolution displays interesting events including subduction polarity reversals subduction shut off and slab break off

the interaction of a subducting lithospheric slab with a chemical or phase boundary we have used a finite element model of time dependent convection to determine the conditions for penetration of the subducted plate into the lower mantle a temperature dependent and non newtonian rheology is applied to achieve platelike behavior of the upper and sinking thermal boundary layer of convection the 650 km discontinuity is taken as either a chemical or phase boundary or as a combination of both it is represented by a marker chain which effects additional buoyancy when distorted out of its equilibrium position when the compositional density contrast is greater than about 5 the descending slab is deflected sidewards at the boundary and two layer convection prevails a resulting depression of the boundary in the range of 50 200 km should be detectable with seismic methods below 5 density difference the slab plunges several hundred kilometers into the lower mantle and below 2 it will probably not stop before reaching the core mantle boundary and extensive mixing would be expected with a pure phase change a negative clapeyron slope of about 6 mpa k 60 bar k is required to establish a type of leaky double layer convection a more moderate slope can aid a small compositional density difference to prevent slab penetration into the lower mantle with the present uncertainties about the physical nature of the 650 km discontinuity a variety of convective styles appears possible on dynamical grounds

how mantle slabs drive plate tectonics the gravitational pull of subducted slabs is thought to drive the motions of earth s tectonic plates but the coupling between slabs and plates is not well established if a slab is mechanically attached to a subducting plate it can exert a direct pull on the plate alternatively a detached slab may drive a plate by exciting flow in the mantle that exerts a shear traction on the base of the plate from the geologic history of subduction we estimated the relative importance of pull versus suction for the present day plates observed plate motions are best predicted if slabs in the upper mantle are attached to plates and generate slab pull forces that account for about half of the total driving force on plates slabs in the lower mantle are supported by viscous mantle forces and drive plates through slab suction

quantifying the net slab pull force as a driving mechanism for plate tectonics 1 it has remained unclear how much of the negative buoyancy force of the slab fb is used to pull the trailing plate at the surface into the mantle here i present three dimensional laboratory experiments to quantify the net slab pull force fnsp with respect to fb during subduction results show that fnsp increases with increasing slab length and dip up to 8 12 of fb making fnsp up to twice as large as the ridge push force the remainder of fb is primarily used to drive rollback induced mantle flow 70 to bend the subducting plate at the trench 15 30 and to overcome shear resistance between slab and mantle 0 8

the influence of trench migration on slab penetration into the lower mantle abstract a two dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth the transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface the viscosity depends on temperature and depth a variety of styles of slab behaviour is found depending predominantly on the trench velocity when trench retreat is faster than 2 4 cm a the descending slab flattens above the phase boundary at slower rates it penetrates straight into the lower mantle although flattening in the transition zone may occur later leading to a complex slab morphology the slab can buckle independent of whether it penetrates or not especially when there is a localised increase in viscosity at the phase boundary flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle the results offer a framework for understanding the variety in slab geometry revealed by seismic tomography

an experimental study of subduction and slab migration results of fluid dynamical experiments on the behavior of subducted slabs are presented that address two important characteristics of subduction slab penetration through the transition zone and horizontal slab migration cold negatively buoyant molded slabs of concentrated sucrose solution with viscosities of 3 5 106 p were introduced into a more dilute two layered sucrose solution representing the upper and lower mantle the transition zone was modeled by a step increase in both density and viscosity ratios of slab thickness to upper layer depth the relative viscosities in each layer and the effective rayleigh number were close to mantle values the initial configuration consisted of two plates separated by a trench gap with one plate attached to a dipping slab simulating a developing subduction zone with an overriding plate two different boundary conditions were used at the trailing ridge end of the subducting plate 1 fixed end corresponding to zero ridge motion and 2 free end corresponding to zero ridge resistance based on results from 15 experiments with various combinations of upper layer lower layer and slab densities u l s and viscosities we find that slab penetration depth is sensitive to both the normalized slab density anomaly r s l s u and the dip angle three penetration regimes are recognized r 0 5 in which the slab sinks into the lower layer without distortion 0 2 r 0 5 which permits limited penetration the amount dependent on dip angle in this regime the slab develops a thick root just beneath the transition zone and r 0 2 in which the slab is deflected by the transition and little or no penetration occurs retrograde slab motion and trench migration occurred in nearly every case with the overriding plate moving parallel to but faster than the trench causing the gap between subducting and overriding plates to close with time prograde trench motion was never observed and the trench remained stationary only in the extreme case free ridge condition with equal upper and lower layer densities these results suggest that the transient extension found in many back arc basins may be a direct consequence of the tendency for retrograde subduction and the amount of extension may be governed in part by the degree of slab penetration through the transition zone

trench motion slab geometry and viscous stresses in subduction systems summary a semi analytic 3 d model for subduction within a newtonian viscous upper mantle provides a dynamically consistent means of computing viscous stress trench motion and slab geometry in subduction systems although negative slab buoyancy provides the basic driving force for subduction slabs that extend from the surface to the base of the upper mantle are oversupported by viscous stresses in the shallow 100 150 km mantle and undersupported by viscous stresses at greater depth in the upper mantle these deeper parts of the subduction system act as an engine for subduction while shallower parts act as a brake on trench motion trench migration rates and slab geometry reflect a competition between these two effects during steady state subduction trench migration rates vary approximately linearly with slab buoyancy and model rates of trench motion are in good agreement with the range of observed rates for a two layer upper mantle viscosity of 2 10 20 pa s above 300 km and 5 10 20 pa s below steady state slab dip increases as slab density decreases especially for very low density slabs which dip significantly more steeply than high density slabs the horizontal velocity at the top of the lower mantle measured relative to the foreland has a very large effect on trench migration rates rivalling or even exceeding that of slab buoyancy slab width parallel to the trench also has a significant effect on trench migration rates due to the viscous pressure of toroidal flow around the slab the stiffness of the subducting lithosphere does not exert a significant effect on trench migration rates or slab geometries for rigidities compatible with oceanic lithosphere very stiff slabs with elastic plate thicknesses more that 40 km or viscosities in the range of 10 25 10 26 pa s subduct significantly more slowly than weak slabs with trench migration rates in the range of half to a third that of weak slabs large unexpected effects on trench migration rates and slab geometry are exerted by the structure and density of the frontal prism and overriding plate indicating that local geology can exert important constraints on slab dynamics during non steady state subduction rates of trench migration respond rapidly as variably buoyant lithosphere penetrates into the asthenosphere in the absence of other driving forces for convergence trench migration rates can change by a factor of two or more in as little as 2 3 myr for example from 35 to 70 mm yr 1 when an oceanic piece of slab follows a continental one into the subduction system subduction of variable buoyancy lithosphere is accompanied by changes in slab dip with depth and through time

curvature of subducted lithosphere from earthquake locations in the wadati benioff zone estimates of slab curvature from 133 profiles across 12 major subduction zones are obtained by fitting splines through earthquakes locations several geometric quantities are computed from the spline fits including the minimum radius of curvature we also introduce a quantity called the equivalent radius of curvature to precisely define the bending dissipation in either a viscous or a plastic slab the minimum and equivalent radius of curvature are approximately equal in slabs with simple geometry i e a single bend but deviate substantially when the geometry is more complicated a compilation of estimates shows that the equivalent radius of curvature is insensitive to the age of the plate at the trench mean values are 175 km or 191 km depending on the rheology of the slab when bending occurs by viscous deformation the resulting dissipation can exceed the power available from subduction with reasonable choices for viscosity a more consistent description of the energetics is obtained by assuming a plastic rheology

global paleo lithospheric models for geodynamical analysis of plate reconstructions we present a general framework to generate time dependent global subduction history models from kinematic plate reconstructions and explore their associated coupled plate mantle dynamic behaviour slabs are constructed by advecting material into the mantle by prescribing its radial velocity and following the absolute tangential motion of the subducting plate a simple geodynamic scenario where plates and slabs define isopycnic and isoviscous regions in an homogeneous or layered mantle was explored using the boundary element method based software bemearth the resulting dynamic behaviour was used to predict the absolute plate motion directions for the present day and a particular mid cretaceous 125 ma kinematic model we show how the methodology can be used to compare and revise kinematic reconstructions based on their effect on the balance of plate driving forces and the resulting euler poles of subducting plates as an example we compare the farallon plate dynamics at 125 ma in a global model with two reconstructions in the context of the evolution of the western north american cordillera our results suggest a method to identify episodes of absolute plate motions that are inconsistent with the expected plate dynamics

capture of the canary mantle plume material by the gibraltar arc mantle wedge during slab rollback summary recent evidence suggests that a portion of the canary plume travelled northeastwards below the lithosphere of the atlas mountains in north africa towards the alboran domain and was captured 10ma ago by the gibraltar subduction system in the western mediterranean the capture would have been associated with the mantle return flow induced by the westwardretreating slab that would have dragged and trapped a portion of the plume material in the mantle wedge of the gibraltar subduction zone such material eventually contaminated the subduction related volcanism in the alboran region in this work we use scaled analogue models of slab plume interaction to investigate the plausibility of the plume capture an upper mantle scaled model combines a narrow 400 km edge fixed subduction plate with a laterally offset compositional plume the subduction dominated by slab rollback and toroidal mantle flow is seen to increasingly impact on the plume dynamics as the area of influence of the toroidal flow cells at the surface is up to 500 1350 km 2 while the plume head initially spreads axisymmetrically it starts being distorted parallel to the plate in the direction of the trench as the slab trench approaches the plume edge at a separation distance of about 500 km before getting dragged towards mantle wedge when applied to the canary plume gibraltar subductionsystem ourmodelsupportstheobservationallybasedconceptualmodelthatmantle plume material may have been dragged towards the mantle wedge by slab rollback induced toroidal mantle flow using a scaling argument for the spreading of a gravity current within a channel we also show that more than 1500 km of plume propagation in the sublithospheric atlas corridor is dynamically plausible

seismicity and rheology of subducted slabs since the incorporation of the inclined zones of intermediate and deep seismic activity into the theory of plate tectonics1 2 it has been generally accepted that earthquakes with focal depths 100 km occur in subducted slabs of predominantly oceanic lithosphere yet our understanding of the processes which generate intermediate and in particular deep earthquakes is still incomplete this is largely due to our limited knowledge of the rheological behaviour of subducted oceanic lithosphere at depths between 100 and 700 km however recently considerable insight has been gained into the rheology of oceanic lithosphere in near surface conditions and i investigate here whether the seismic activity at depths 100 km can be understood on the basis of this information i first present results of temperature calculations performed for all subduction zones for which adequate data on seismicity relative plate motion and age of the descending oceanic lithosphere exist sub sequently i adapt the temperature conditions governing the rheology near the surface by including depth dependence this leads to a depth dependent critical temperature tcr z above which the subducted lithosphere cannot sustain the stresses necessary to generate seismic events seismicity data and calculated temperature versus depth relationships are in good agreement with the rheological model predictions this implies that the absence of seismic activity at depths 700 km should not be interpreted as direct evidence for the hypothesis that the 650 km discontinuity in the mantle acts as a barrier to vertical motion of subducted slabs

paleogeographic maps of the continents bordering the indian ocean since the early jurassic paleomagnetic data from the african plate from the early jurassic onward are critically re evaluated and updated resulting in a set of 30 paleopoles from which an apparent polar wander path is derived using 20 ma intervals this path is then used in conjunction with kinematic models of the indian and atlantic oceans to place other continents bordering the indian ocean eurasia india islands and plateaus madagascar seychelles mascarene and blocks involved in the collision of india with eurasia south tibet indochina in an absolute paleogeographic frame fixed with respect to africa paleogeographic maps are presented for the early jurassic pre breakup the jurassic cretaceous boundary anomaly ml5 the early cretaceous m2 the late cretaceous a 34 the cretaceous tertiary boundary a 29 the early eocene a 24 the middle eocene a 21 and the oligocene miocene boundary a 6 observed and predicted paleomagnetic directions for landmasses other than africa i e for western europe china korea india and australia which are not used in deriving the maps are found to be in very satisfactory agreement the quality of reconstructions as measured by rms discrepancies between observed and predicted paleolatitudes is of the order of 5 the main features of the paleogeographic maps are discussed with the aim of relating oceanic and on land geodynamic events the constraints provided by tibetan data on the india asia collision and the 20 ma reorganization of geology and plate boundaries in a band extending from the south china sea through tibet to the red sea system are discussed in more detail we emphasize the episodic nature of the breakup of gondwanaland and formation of the indian ocean and subsequent interaction with eurasia in which six major discrete events are identified separated by 30 ma phases of rather uniform plate motion

role of lateral mantle flow in the evolution of subduction systems insights from laboratory experiments summary we present 3 d laboratory experiments constructed to investigate the pattern of mantle flow around a subducting slab under different boundary conditions in particular we present a set of experiments characterized by different conditions imposed at the trailing edge of the subducting plate that is plate fixed in the far field plate detached in the far field imposed plate motion experiments have been performed using a silicone slab floating inside a honey tank to simulate a thin viscous lithosphere subducting in a viscous mantle for each set we show differences between models that do or do not include the possibility of out of plane lateral flow in the mantle by varying the lateral boundary conditions our results illustrate how a subducting slab vertically confined over a 660 km equivalent depth can be influenced in its geometry and in its kinematics by the presence or absence of possible lateral pathways on the basis of these results we show implications for natural subduction systems and we highlight the importance of suitable simulations of lateral viscosity variations to obtain a realistic simulation of the history of subduction

variations of slab dip and overriding plate tectonics during subduction insights from analogue modelling we present small scale laboratory models of oceanic subduction in which plates motion is imposed by lateral boundary conditions the oceanic plate moves trenchward at constant speed and subducts below a fixed overriding plate in this configuration the long term process of subduction is not steady state slab interaction with the upper mantle lower mantle boundary results in periods of slab flattening during which the dip of the slab diminishes followed by periods of slab steepening the overriding plate tectonic regime is influenced by the dynamics of subduction slab anchoring favouring trench perpendicular shortening when the slab is anchored slab flattening further favours shortening while slab steepening favours extension or smaller shortening rates non steady state long term subduction may explain part of the variability of slab geometries evidenced by statistical analyses of present day subduction zones experiments suggest that despite boundary conditions applied on the converging plates do not change tectonics pulses within the overriding plate may be caused by this non steady state behaviour

dynamics of subduction and plate motion in laboratory experiments insights into the plate tectonics behavior of the earth 1 three dimensional laboratory experiments have been designed to investigate the way slab bearing plates move during subduction inside the mantle in our experiments a viscous plate of silicone lithosphere subducts under its negative buoyancy in a viscous layer of pure honey mantle varying thickness width viscosity and density of the plate and mantle three characteristic modes of subduction are observed a retreating trench mode mode i a retreating trench mode following a transient period of advancing trench mode ii and an advancing trench mode mode iii these modes are characterized by different partitioning of the amount of subduction into plate and trench motion our experiments show that the velocity of subduction can be modeled by the dynamic interaction between acting and resisting forces where lithospheric bending represents 75 95 of the total resisting forces however our experimental results also show the impossibility to predict a priori the plate velocity only from the velocity of subduction without considering trench migrations we find experimentally that the lithospheric radius of curvature which depends upon plate characteristics stiffness and thickness and the mantle thickness exerts a primary control on the trench behavior our results suggest that the complexity of the style of subduction could be controlled by geometrical rules of a plate bending inside a stratified mantle the earth system is in the crucial range for the interplay between the rigidity of the plate and the mantle stratification this setting may be the responsible for the complexity of the past and present tectonic styles

the influence of olivine metastability on deep subduction of oceanic lithosphere abstract the dynamics of subduction of oceanic lithosphere is mainly controlled by thermal buoyancy forces and mineralogy in self consistent models the influence of metastable olivine mo on the dynamics and evolution of subducting oceanic lithosphere is investigated we consider both phase boundaries olivine ol to spinel sp 410 km and sp to perovskite pv magnesiowustite mw 660 km treating the kinetics of the phase transition in a simplified way while models do not allow the formation of mo penetrate straight into the lower mantle slabs with a sufficient amount of mo flatten out at the 660 km boundary before penetrating at a later stage furthermore our models show that the effect of mo to slow down subduction velocities within the transition zone is rather moderate

mantle flow mechanisms for the large scale subsidence of continental interiors evidence in the geological record shows that continental interiors periodically undergo enigmatic episodes of large scale subsidence we propose that mantle flow associated with the descent of cold plumes and slabs and their interaction with the endothermic phase change at 660 km depth may provide a plausible mechanism for these epeirogenic events simulations of mantle convection that incorporate the thermodynamic effects of the phase change and both depth and temperature dependent viscosity are used to model the descent of plumes and slabs through the 660 km boundary we find that the plume slab flow scenarios are capable of supporting topographic deflections of amplitudes of 1 km and horizontal wavelengths of 1000 km that persist over time scales of 100 150 m y

the relationship between plate velocity and trench viscosity in newtonian and power law subduction calculations convection with a newtonian temperature dependent rheology leads to little or no surface velocity unless zones of weakness are introduced plate like features are observed in calculations both with newtonian rheology employing imposed weak zones and with power law non newtonian rheology where high stresses at the trench reduce the effective viscosity since deformation at subduction zones involves faulting both of these parameterizations should be treated with some skepticism it is important to understand how the parameterizations affect the model results we study the relationship between trench viscosity and plate velocity using a newtonian rheology by varying the viscosity at the trench the plate velocity is a function of the trench viscosity for fixed rayleigh number and plate slab viscosity slab velocities for non newtonian rheology calculations are significantly different from slab velocities from newtonian rheology calculations at the same effective rayleigh number both models give reasonable strain rates for the slab when compared with estimates of seismic strain rate non newtonian rheology eliminates the need for imposed weak zones and provides a self consistent fluid dynamical mechanism for subduction in numerical convection models