latent heat effects of the major mantle phase transitions on low angle subduction abstract very low to zero shallow dip angles are observed at several moderately young subduction zones with an active trenchward moving overriding plate we have investigated the effects of latent heat for this situation where mantle material is pushed through the major mantle phase transitions during shallow low angle subduction below the overriding plate the significance of the buoyancy forces arising from the latent heat effects on the dynamics of the shallowly subducting slab is examined by numerical modeling when a 32 ma old slab is overridden with 2 5 cm yr by a continent flat subduction occurs with a 4 5 cm yr convergence rate when latent heat is included in the model forced downwellings cause a thermal anomaly and consequently thermal and phase buoyancy forces under these circumstances the flat slab segment subducts horizontally about 350 km further and for about 11 ma longer than in the case without latent heat before it breaks through the 400 km phase transition the style of subduction strongly depends on the mantle rheology increasing the mantle viscosity by one order of magnitude can change the style of subduction from steep to shallow similarly an overriding velocity of less than 1 cm yr leads to steep subduction which gradually changes to flat subduction when increasing the overriding velocity however these model parameters do not change the aforementioned effect of the latent heat provided that low angle subduction occurs in all models latent heat resulted in a substantial increase of the flat slab length by 300 400 km varying the olivine spinel transition clapeyron slope from 1 to 6 mpa k reveals a roughly linear relation between and the horizontal length of the slab based on these results we conclude that buoyancy forces due to latent heat of phase transitions play an important role in low angle subduction below an overriding plate

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

back arc strain in subduction zones statistical observations versus numerical modeling 1 recent statistical analysis by lallemand et al 2008 of subduction zone parameters revealed that the back arc deformation mode depends on the combination between the subducting vsub and upper vup plate velocities no significant strain is recorded in the arc area if plate kinematics verifies vup 0 5 vsub 2 3 cm a in the hs3 reference frame arc spreading shortening occurs if vup is greater lower than the preceding relationship we test this statistical law with numerical models of subduction by applying constant plate velocities far away from the subduction zone the subducting lithosphere is free to deform at all depths we quantify the force applied on the two converging plates to sustain constant surface velocities the simulated rheology combined viscous non newtonian and brittle behaviors and depends on water content the influence of subduction rate vs is first studied for a fixed upper plate after 950 km of convergence steady state slab pull the transition from extensional to compressive stresses in the upper plate occurs for vs 1 4 cm a the effect of upper plate velocity is then tested at constant subduction rate upper plate retreat advance with respect to the trench increases extension compression in the arc lithosphere and increases decreases the subducting plate dip our modeling confirms the statistical kinematic relationship between vsub and vup that describes the transition from extensional to compressive stresses in the arc lithosphere even if the modeled law is shifted toward higher rates of upper plate retreat using our set of physical parameters e g 100 km thick subducting oceanic plate and short term simulations our results make valid the choice of the hs3 reference frame for assessing plate velocity influence on arc tectonic regime the subduction model suggests that friction along the interplate contact and the mantle stokes reaction could be the two main forces competing against slab pull for upper mantle subductions besides our simulations show that the arc deformation mode is strongly time dependent

numerical comparison of different convergent plate contacts subduction channel and subduction fault summary at convergent plate boundaries the properties of the actual plate contact are important for the overall dynamics convergent plate boundaries both mechanically decouple and link tectonic plates and accommodate large amounts of strain we investigate two fundamental physical states of the subduction contact one based on a fault and the other based on a subduction channel using a finite element method we determine the specific signatures of both states of the subduction contact we pay particular attention to the overriding plate in a tectonic setting of converging plates where the subducting plate is freely moving the subduction channel reduces compression relative to the fault model in a land locked basin setting where the relative motion between the far field of the plates is zero the subduction channel model produces tensile stress regime in the overriding plate even though the amount of slab roll back is small the fault model shows a stronger development of slab roll back and a compressive stress regime in the upper plate based on a consistent comparison of fault and channel numerical models we find that the nature of the plate contact is one of the controlling factors in developing or not of backarc extension we conclude that the type of plate contact plays a decisive role in controlling the backarc state of stress to obtain backarc extension roll back is required as an underling geodynamic process but it is not always a sufficient condition

relations among subduction parameters clues to the dynamics of the subduction process are found in the many measurable parameters of modern subduction zones based on a critical appraisal of the geophysical and geological literature 26 parameters are estimated for each of 39 modern subduction zones to isolate causal relationships among these parameters multivariate analysis is applied to this data set this analysis yields empirical quantitative relations that predict strain regime and strike slip faulting in the overriding plate maximum earthquake magnitude benioff zone length slab dip arc trench gap and maximum trench depth excellent correlation is found between length of the benioff zone and the product of convergence rate and age of the downgoing slab this relationship is consistent with the conductive heating model of molnar et al 1979 if the model is modified in one respect the rate of heating of the slab is not constant it is substantially slower during passage of the slab beneath the accretionary prism and overriding plate the structural style in the overriding plate is determined by its stress state though the stress state of overriding plates cannot be quantified one can classify each individual subduction zone into one of seven semiquantitative strain classes that form a continuum from strongly extensional class 1 back arc spreading to strongly compressional class 7 active folding and thrusting this analysis indicates that strain class is probably determined by a linear combination of convergence rate slab age and shallow slab dip interplate coupling controlled by convergence rate and slab age is an important control on strain regime and the primary control on earthquake magnitude arc parallel strike slip faulting is a common feature of convergent margins forming a forearc sliver between the strike slip fault and trench optimum conditions for the development of forearc slivers are oblique convergence a compressional environment and a continental overriding plate the primary factor controlling presence of strike slip faulting is coupling strongly oblique convergence is not required the rate of strike slip faulting is affected by both convergence obliquity and convergence rate maximum trench depth is a response to flexure of the underthrusting plate the amount of flexural deflection at the trench depends on the vertical component of slab pull force which is very sensitive to slab age and shallow slab dip shallow slab dip conforms to the cross sectional shape of the overriding plate which is controlled by width of the accretionary prism and duration of subduction deep slab dip is affected by the mantle trajectory established at shallow depth but may be modified by mantle flow much of the structural complexity of convergent margins is probably attributable to terrane juxtaposition associated with temporal changes in both forearc strike slip faulting and strain regime empirical equations relating subduction parameters can provide both a focus for future theoretical studies and a conceptual and kinematic link between plate tectonics and the geology of subduction zones

a numerical study of the interaction between the mantle wedge subducting slab and overriding plate abstract we have formulated a numerical model with strongly temperature dependent viscosity to calculate thermal structure and flow field in subduction zones one important particularity of the model is that the overriding plate is not fixed over its whole thickness in order to allow material exchange between the wedge and the upper lithosphere numerical problems due to very high viscosity contrasts are avoided by coupling a finite difference method and a finite element method for solving the energy conservation equation and the stockes equation respectively in this model a temperature decrease from 1400 to 1300 c increases the viscosity by an order of magnitude we study the temperature structure and the velocity field of the subducting slab and mantle wedge surface heat flow velocity anomalies and geometry of the partial melting zone are also calculated to study the effect that boundary conditions play on the interaction between the mantle wedge overriding plate and subducting plate we examine models with both fixed and free slip conditions applied to the overriding plate when the overriding plate is allowed to move laterally free slip the subducting slab is thick and both the temperature field and the convective motions in the mantle wedge are similar to those observed when using constant viscosity numerical models or analytical corner flow models if the surface of the overriding plate is fixed the subducting slab is thin and the mantle wedge impinges upon the overriding plate forming a high temperature nose between the overriding plate and subducting lithosphere furthermore viscous decoupling occurs implicitly at shallow depth between the slab and the wedge because hot material from the wedge is entrained close to the trench in that case the subducting slab tectonically erodes the lower lithosphere of the overriding plate leading to high temperatures low seismic velocities high attenuation and high heat flow beneath volcanic arc in agreement with geophysical observations

the dynamical control of subduction parameters on surface topography the long wavelength surface deflection of earth s outermost rocky shell is mainly controlled by large scale dynamic processes like isostasy or mantle flow the largest topographic amplitudes are therefore observed at plate boundaries due to the presence of large thermal heterogeneities and strong tectonic forces distinct vertical surface deflections are particularly apparent at convergent plate boundaries mostly due to the convergence and asymmetric sinking of the plates having a mantle convection model with a free surface that is able to reproduce both realistic single sided subduction and long wavelength surface topography self consistently we are now able to better investigate this interaction we separate the topographic signal into distinct features and quantify the individual topographic contribution of several controlling subduction parameters results are diagnosed by splitting the topographic signal into isostatic and residual components and by considering various physical aspects like viscous dissipation during plate bending performing several systematic suites of experiments we are then able to quantify the topographic impact of the buoyancy rheology and geometry of the subduction zone system to each and every topographic feature at a subduction zone and to provide corresponding scaling laws we identify slab dip and slightly less importantly slab buoyancy as the major agents controlling surface topography at subduction zones on earth only the island arc high and the back arc depression extent are mainly controlled by plate strength overall his modeling study sets the basis to better constrain deep seated mantle structures and their physical properties via the observed surface topography on present day earth and back through time

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

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

trench migration and subduction zone geometry seismic tomography results show a large variety in upper mantle structure along convergent plate boundaries we numerically investigate the effect of trench migration on the evolution of a slab with temperature dependent viscosity encountering a viscosity interface we find that subduction zone geometry is sensitive to even small rates 1 cm yr of retrograde motion increase in trench migration rate decreases the dip angle of the slab and its ability to penetrate the lower mantle upon including a background mantle flow it turns out that trench migration relative to the upper mantle flow is more decisive than the absolute plate velocities subduction zone geometry appears to be strongly time dependent we conclude that individual tectonic setting and time dependent slab behaviour can account for many different types of observable subduction zone anomalies our model temperature can adequately account for magnitudes and patterns of seismic anomalies as obtained from seismic tomography

slab mantle flow interaction influence on subduction dynamics and duration we investigate the influence of mantle flow relative to the lithosphere on subduction dynamics we use 2d thermo mechanical models assuming incompressible non newtonian fluid rheology different mantle flow velocities consistent with absolute plate motion models are tested as well as both directions of flow either sustaining or opposing slab dip the effects of different inflow outflow velocity profiles slab strengths and upper lower mantle viscosity contrasts are also evaluated slab dip deviations between models with opposite mantle flow directions range from 37 for relatively strong slabs max 1025 pa s to 50 for weaker slabs max 1024 pa s accounting for a significant amount of natural slab dip variability for imposed mantle flow supporting the slab the initial stage of slab steepening is followed by a stage of continuous slab dip decrease this slab shallowing eventually leads to mantle wedge closure subduction cessation and slab break off possibly driving subduction flips

slab dip vs lithosphere age no direct function one paradigm of subduction relates the dip of the slab to the buoyancy of the downgoing lithosphere along subduction zones with the negative buoyancy proportional to the age of the oceanic lithosphere we measured the dip of the slab down to depths of 250 km along 164 sections crossing 13 subduction zones and compared it with the age of the subducting oceanic lithosphere both at the trench and at depth we show here that this relationship is far more irregular than previously suggested and that it is not possible to simply correlate the increase of the slab dip to the increasing age of the downgoing cooler lithosphere younger oceanic lithosphere may show steeper dip than older segments of slabs e g central america vs south america in contrast with predictions of models considering only slab pull the combination of slab age and subduction rate better accounts for slab dip however the correlation is not satisfactory correlation coefficient equal to 0 450 these results suggest that supplemental forces or constraints have to be accounted for such as thickness and shape of the hangingwall plate absolute plate velocity presence of lateral density variations in the hosting upper mantle effects of accretion erosion subduction of oceanic plateaus and slab deformation due to the motion of the mantle relative to the subducting plate d 2005 elsevier b v all rights reserved

tectonic drivers of the wrangell block insights on fore arc sliver processes from 3 d geodynamic models of alaska intracontinental shear zones can play a key role in understanding how plate convergence is manifested in the upper plate in regions of oblique subduction however the relative role of the driving forces from the subducting plate and the resisting force from within intracontinental shear zones is not well understood results from high resolution geographically referenced instantaneous 3d geodynamic models of flat slab subduction at the oblique convergent margin of alaska are presented these models investigate how viscosity and length of the denali fault intracontinental shear zone as well as coupling along the plate boundary interface modulate motion of the wrangell block forearc sliver and slip across the denali fault models with a weak denali fault 1017 pa s and strong plate coupling 1021 pa s were found to produce the fastest motions of the wrangell block 10 mm yr the 3d models predict along strike variation in offs et along the denali fault changing from dextral strike slip motion in the eastern segment to oblique convergence towards the fault apex models further show that the flat slab drives oblique motion of the wrangell block and contributes to 20 models with a short fault and 28 models with a long fault of the observed quaternary slip rates along the denali fault the 3d models provide insight into the general processes of forearc sliver mechanics and also offer a 3d framework for interpreting hazards in regions of flat slab subduction

non steady state subduction and trench parallel flow induced by overriding plate structure abstract the direction of plate tectonic motion and the direction of mantle flow as inferred from observations of seismic anisotropy measurements show a good global correlation far from subduction zones however this correlation is poor near subduction zones where below the slab seismic anisotropy is aligned parallel to the trench and above the slab has a complex pattern which has not been fully explained here we present time dependent three dimensional 3d fully dynamic simulations of subduction to study the effect of overriding plate structure on the evolution of slab geometry and induced mantle flow we find that along strike variation in thermal thickness of the overriding plate causes increased hydrodynamic suction and shallower slab dip beneath the colder portion of the overriding plate the variation in slab geometry drives strong trench parallel flow beneath the slab and a complex flow pattern above the slab this new mechanism for driving trench parallel flow provides a good explanation for seismic anisotropy observations from the middle and south america subduction zones where both slab dip and overriding plate thermal state are strongly variable and correlated and thus may be an important mechanism in other subduction zones the location and strength of trench parallel flow vary with the time dependent evolution of the slab suggesting that the global variability in seismic anisotropy observations in subduction zones is in part due to the non steady state behavior of these systems

dynamics and implications of slab detachment due to ridge trench collision 1 the approach of a buoyant spreading ridge to a subduction zone may lead to detachment of a subducted slab previous work has called upon the detachment process as an explanation for observed ridge abandonment and slab window related magmatism e g in baja california but such a scenario has not been tested using fully dynamic numerical models we use dynamic two dimensional models including a non newtonian rheology to study the approach of a trench parallel spreading ridge to a subduction zone we find that before the ridge approaches within 100 km of the trench detachment of the subducted slab occurs due to the combined effects of increased buoyancy and reduced strength of the lithosphere with ridge proximity in models exploring effects of subducted slab length distance of the ridge from the trench shear zone strength and lithospheric yield strength we find the following 1 a decrease in subduction velocity as the ridge approaches the trench 2 a shrinking surface plate that maintains a uniform plate like subduction velocity 3 ridge abandonment distances 100 275 km from the trench and 4 slab gap distances 195 285 km from the trench these results are consistent with observations in baja california where detachment of the cocos slab may explain abandonment of observed segments of the east pacific rise 50 200 km outboard of the trench and the presence of nonarc related magmatism located 100 250 km inboard of the trench ridge subduction does not occur within any of these two dimensional cases of ridge trench collision indicating that 3 d geometric effects may be responsible for ridge subduction

subduction adjacent to deep continental roots enhanced negative pressure in the mantle wedge mountain building and continental motion abstract subduction near to a deep continental root is modeled with finite elements the root restricts asthenospheric corner flow which increases mantle wedge suction the suction pulls the continental root toward the subduction zone simultaneously promoting trenchward continental motion and compressive tectonics in the lithosphere between the root and the subduction zone in a 2 d model with a root that extends to 250 km depth and lies 300 km from the mantle wedge corner we find two fold increase in mantle wedge suction and nearly a four fold increase in continental compression to 5 10 12 n m compared to a rootless reference model in 3 d modeling we find that distance to the subduction zone is more important than root width in subduction strike direction or reasonable variations in sub lithospheric viscosity structure with the root enhanced suction effects diminishing rapidly as the wedge corner root distance increases beyond 400 km when finite width roots are wider than the distance to the mantle wedge corner and closer than 400 km compression of the continent between the root and the subduction zone is 50 that found in the infinite width root

horizontal subduction and truncation of the cocos plate beneath central mexico 1 based on analysis of data from a trans mexico temporary broadband seismic network centered on mexico city we report that the subducting cocos plate beneath central mexico is horizontal and tectonically underplates the base of the crust for a distance of 250 km from the trench it is decoupled from the crust by a very thin low viscosity zone the plate plunges into the mantle near mexico city but is truncated at a depth of 500 km probably due to an e w propagating tear in the cocos slab unlike the shallow slab subduction in peru and chile there is active volcanism along the trans mexican volcanic belt tmvb that lies much further inland than regions to either side where subduction dip is not horizontal geodynamical modeling indicates that a thin weak layer such as imaged by the seismic experiment can explain the flat 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

trench migration net rotation and slab mantle coupling article i nfo laboratory models have been conducted to improve our understanding of the role that the resistance of the slab to bending and its coupling to the ambient mantle play in subduction dynamics over geological time scales our models are set up with a viscous plate of silicone lithosphere subducting under negative buoyancy in a viscous layer of glucose syrup mantle for our study the lithosphere upper mantle viscosity contrast has been systematically varied from 10 to 10 5 in order to explore the parameter space between weak and strong slab dynamics we found that subduction is characterized by a retreating mode for viscosity ratios n10 4 by the coexistence of a retreating mode and an advancing mode for viscosity ratios between 10 4

the influence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity results of high resolution seismic tomography showing subducting slabs deflected in the transition zone and thickened in the lower mantle seem to call for slab material weaker than inferred from mineral physics deformation mechanisms a possible mechanism suggested by several authors could be the weakening due to grain size reduction which should occur in the cold portion of fast slabs after an exothermic phase transition at a depth of 400 km since the amount of weakening as well as the rate of subsequent strengthening due to the grain growth are not precisely known we present here a parametric study of slab behavior in the transition zone and upper part of the lower mantle we simulate a subducting slab in a two dimensional 2 d cartesian box in the numerical model with composite rheology including diffusion creep dislocation creep and a general stress limiting rheology approximating peierl s creep we concentrate on two rheologic effects the dynamic effect of slab weakening due to grain size reduction at the phase boundary and the effect of yield stress of stress limiting rheology the effect of trench migration on slab deformation is also included in our study results show that the slab ability to penetrate into the lower mantle is not significantly affected by a trench retreat in the absence of grain size weakening however in case of a 4 cm yr trench retreat grain size weakening provides a viable mechanism to deflect the slab in the transition zone provided that stress limiting deformation mechanism would limit the effective viscosity outside the areas of grain size weakening to about 10 24 pa s 2002 elsevier science b v all rights reserved

the mesoproterozoic rayner province in the lambert glacier area its age origin isotopic structure and implications for australia antarctica correlations within the rayner province in the lambert glacier area two proterozoic lithotectonic zones are distinguished which differ in age and lithology in the fisher zone a significant juvenile component 1nd t 1 4 2 to 4 is represented by mafic intermediate and felsic plutonic and volcanic rocks the tectonomagmatic processes were concentrated at c 1400 1200 ma in the beaver zone orthogneisses record multiple emplacement metamorphism and ductile deformation at c 1150 930 ma new u pb shrimp ii zircon ages c 1140 ma and 1095 ma protolith emplacement c 950 850 and c 540 ma metamorphism chemical and sm nd isotopic data indicate that both zones have their continuations in the eastern amery ice shelf coast the rayner province has some geological features in common with the albany fraser orogen in western australia suggesting that these regions evolved in similar tectonic environments at c 1400 1300 ma followed by steady closure from east to west of the ocean that separated the mawson continent from western australia supplementary material chemical compositions of rocks a summary of trace element ratios a summary of published u pb dates field photo views new u pb analyses by sims shrimp ii and new sm nd data are available at www geolsoc org uk sup18621 the crystalline basement of east antarctica comprises a few relatively small archaean cratons enclosed in a vast proterozoic circum east antarctic mobile belt ceamb yoshida kizaki 1983 ivanov kamenev 1990 tingey 1991 yoshida 1992 black et al 1992a development of the mobile belt occurred largely between c 1150 ma and c 950 ma fitzsimons 2000a 2003 and references therein the ceamb was correlated with the grenville province of north america and thought to have been conjugated within the rodinia supercontinent hoffman 1991 dalziel 1991 condie 2001 in some areas the ceamb was affected by major tectono thermal reworking at c 550 500 ma and this suggested a model of continental collision resulting in the amalgamation of gondwana at that time shiraishi et al 1992 fitzsimons 2000a boger et al 2001 typical rock assemblages of the ceamb were first described in western enderby land by kamenev et al 1965 and kamenev 1972 these authors distinguished the rayner metamorphic complex considering it to be a proto geosynclinal formation other investigators considered the rayner complex to be strongly and probably repeatedly reworked earliest precambrian crust with the rayner structural episode c 1000 950 ma being the most prominent event sandiford wilson 1984 sheraton et al 1987 u pb zircon age data obtained from various parts of the ceamb showed considerable spread e g fitzsimons 2000b and three separate tectonic provinces have now been from harley s l fitzsimons i c w zhao y eds 2013 antarctica and supercontinent evolution geological society london special publications 383 35 57 first published online may 20 2013 http dx doi org 10 1144 sp383 1 the geological society of london 2013 publishing disclaimer www geolsoc org uk pub ethics by guest on december 12 2013 http sp lyellcollection org downloaded from

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

newtonian versus non newtonian upper mantle viscosity implications for subduction initiation 1 the effect of rheology on the evolution of the slab tip during subduction initiation is analyzed using 2 d numerical flow models experimentally determined flow laws have both strong temperature and stress dependence which leads to large local variations in viscosity with direct consequences for subduction initiation we find that models with newtonian viscosity lead to flat or coupled subduction due to hydrodynamic stresses that pull the slab tip up towards the overriding plate non newtonian rheology reduces these hydrodynamic stresses by decreasing the wedge viscosity and the slab coupling to wedge corner flow rendering the small negative slab buoyancy of the slab tip sufficient to maintain its dip during the early stages of subduction

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

predicting trench and plate motion from the dynamics of a strong slab the motion of oceanic plates is commonly related to the subduction of cold and dense oceanic material into the mantle these models predict plate velocities from subduction velocities but the trench motion is not directly included in the computation here using a recent compilation of a global data set we found that the motion of trenches either advancing or retreating with respect to upper plates scales with their corresponding subducting plates motion based on simple experimental tests we found that subduction of strong slabs inside the upper mantle correctly predicts these kinematic relationships we deduce that the motion of the trenches represent the surface manifestation of the resistance encountered by the subducting lithosphere to bend and penetrate within the upper mantle

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

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

geometry of subduction and depth of the seismogenic zone in the guerrero gap mexico the western coast of guerrero in southern mexico has been identified as a seismic gap on the middle american trench in which no large earthquakes have occurred at least since 1908 it has been suggested that the seismic energy accumulated since 1908 will eventually be released by a large earthquake a permanent seismic network was installed to monitor the seismicity of this mature seismic gap and to understand the geometry of the subducted slab beneath this region the seismicity defines an unusual distribution along two bands of seismic activity parallel to the coast the resulting geometry of the subduction zone shows that the cocos plate dips at a shallow angle beneath the north american plate to a depth of 40 km from there the subducted slab is bent upward following a subhorizontal trajectory extending inland at a depth of 50 km this plate geometry is reminiscent of that found in peru and central argentina two other regions where a young oceanic plate is being subducted in mexico however the slab underplates an overriding plate which is only half as thick as that observed in south america a possible oceanic origin of the allochtonous terranes comprising southern mexico may explain the presence of the anomalously thin lithosphere in this region

subduction zones observations and geodynamic models abstract this review of subduction and geodynamic models is organized around three central questions 1 why is subduction asymmetric 2 are subducted slabs strong or weak 3 how do subducted slabs interact with phase transformations changes in mantle rheology and possibly chemical boundaries in the mantle based on laboratory measurements of the temperature dependence of olivine one would conclude that the core of a subducting slab is at least 10 000 times more viscous than ambient mantle however there are a number of complementary but independent observations that suggest that slabs are much weaker than this slabs undergo significant deformation in the upper mantle and may thicken to twice their original width by the time they reach the base of the transition zone the lack of a clear correlation between the observed dip angle of deep slabs and plate velocity rate of trench migration and slab age in modern subduction zones is consistent with hypothesis that subduction is a time dependent phenomenon both tank and numerical convection experiments with plates conclude that subduction is not a steady phenomenon but that slabs bend thicken stretch and change dip through time this is at odds with the assumptions used in steady state slab thermal models where slab deformation is not considered

epeirogenic uplift above a detached slab in northern central america p wave tomographic images reveal that the northern central america highlands east of the modern volcanic arc overlie a detached slab hypsometric analysis of the highlands in honduras demonstrates that the region is a dissected plateau that is disrupted by normal faults near the north american caribbean plate margin the dissected central american plateau contains a network of superimposed rivers with meanders cut into bedrock such a geomorphic character indicates that the regional uplift occurred in the absence of tilting we propose that the epeirogenic uplift of northern central america is the buoyant upper plate response to the influx of mantle asthenosphere following the break off and sinking of the slab

two stage subduction history under north america inferred from multiple frequency tomography eastward subduction of oceanic tectonic plates has shaped the geologic history of western north america over the past 150 million years the mountain building and volcanism that brought forth the spectacular landscapes of the west are credited to the vast ancient farallon plate which interacted mechanically and chemically with the overlying continent as it plunged back into the mantle here we use finite frequency travel time and amplitude measurements of teleseismic p waves in seven frequency bands to obtain a high resolution tomographic image to approx1 800 km depth we discover several large previously unknown pieces of the plate which show that two distinct stages of whole mantle subduction are present under north america the currently active one descends from the pacific northwest coast to 1 500 km depth beneath the great plains whereas its stalled predecessor occupies the transition zone and lower mantle beneath the eastern half of the continent we argue that the separation between them is linked to the laramide era 70 50 myr ago a time of unusual volcanism and mountain building far inland generally explained by an episode of extremely flat subduction

the influence of olivine metastability on the dynamics of subduction abstract both seismological observations and mineral physics experiments suggest that old and cold subducting slabs might contain a wedge of metastable olivine mo dynamically consistent models of subducting slabs in which the kinetics of the olivine to spinel wadsleyite or ringwoodite transition is treated in a simplified way are carried out to evaluate the buoyancy effect of mo on the subduction velocity assuming slab thicknesses according to the cooling half space model we find that mo forms in slabs with ages greater than approximately 70 myr and that a significant slowing down of subduction velocity occurs for ages greater than 100 myr because a decrease in the subduction velocity for old lithospheric ages is not observed our models suggest that the amount of mo in slabs is less than approximately 5000 km2 in cross section and that for old lithosphere the cooling half space model is not applicable

bending mechanics and mode selection in free subduction a thin sheet analysis summary to elucidate the dynamics of free buoyancy driven subduction of oceanic lithosphere i study a model in which a 2 d sheet of viscous fluid with thickness h and viscosity 1 subducts in an infinitely deep ambient fluid with viscosity 1 numerical solutions for the sheet s evolution are obtained using the boundary element method bem starting from an initial configuration comprising a short protoslab attached to a longer horizontal plate that is free to move laterally beneath an impermeable traction free surface interpretation of the solutions using thin viscous sheet theory shows that the fundamental length scale controlling the subduction is the bending length lb defined at each instant as the length of the portion of the sheet s midsurface where the rate of change of curvature is significant geophysically speaking lb is the sum of the lengths of the slab and of the region seaward of the trench where flexural bulging occurs the bending length in turn enters into the definition of the sheet s dimensionless stiffness s h lb 3 which controls whether the sinking speed of the slab is determined by the viscosity of the sheet itself s 1 or by that of the ambient fluid s 1 motivated by laboratory observations of different modes of subduction retreating versus advancing trench folding versus no folding etc in fluid layers with finite depth i calculate numerically the dip d of the slab s leading end as a function of and the normalized depth d h to which it has penetrated the contours of the function d d h strongly resemble the intermode boundaries in the laboratory based regime diagram of schellart supporting the hypothesis that the mode of subduction observed at long times in experiments is controlled by the dip of the slab s leading end when it reaches the bottom of the layer in particular the bem solutions explain why trenches advance in the laboratory only when lies in an intermediate range and why they retreat when is either smaller or larger than this application of the bem model to wu et al s compilation of the minimum curvature radii of subducted slabs suggests 140 510 for the earth this is too small to permit the laboratory type trench advancing mode in agreement with the lack of tomographic evidence for slabs that are bent over backwards

the dynamics of subduction and trench migration for viscosity stratification summary although subduction plays a dominant role in plate tectonics related processes such as trench migration are not well understood we investigate the influence of viscosity stratification and different flow boundary conditions on the trench migration dynamics of a subducting lithosphere in a 2 d model using a finite difference code subduction and mantle flow are driven by a compositionally prescribed density contrast between the high viscosity slab and the ambient mantle no overriding plate is assumed the model rheology is viscoplastic at shallow and viscous at greater depth the mantle is stratified with a higher viscosity in the lower part of the computational domain the plastic rheology allows the plate to decouple from the top of the box and is described by a simplified friction law to mimic the brittle deformation of the shallow lithosphere the flow patterns and slab morphologies differ as a function of the velocity boundary conditions bcs for reflective bcs we model either a plate which is fixed on the box side or a plate that can freely move a slab with periodic bcs is always laterally free we find that when the slabs interact with the highly viscous lower mantle the slab with reflective bcs is more folded than the periodic one due to the confinement of mantle flow around the slab in general the trench is observed to roll back toward the oceanic plate for free slabs with reflective bcs the trench retreat velocity deceases after the interaction with the lower mantle in contrast for the free slab with periodic bcs it is nearly constant or even increases with time if described in a reference frame fixed to the lower mantle the free slab with periodic bcs exerts a force on the lower mantle which causes a net differential flow between the upper and lower mantle we also carried out additional experiments with different mechanical slab properties and thicknesses stiffer free slabs subduct more steeply and are able to penetrate straight into the lower mantle after penetration into the lower mantle the trench moves toward the continent thinner slabs deform more easily and are folded for all boundary conditions we can generalize the rollback systematics for free slabs by normalizing trench velocities by the plate speed free slabs being more representative of the west pacific subduction zones show large temporal variations of normalized trench velocities ranging between a retreat value of 1 0 to an advancing value of 0 3

the impact of the south american plate motion and the nazca ridge subduction on the flat subduction below south peru 1 flat subduction near peru occurs only where the thickened crust of the nazca ridge subducts furthermore the south america continent shows a westward absolute plate motion both the overriding motion of south america and the subduction of the nazca ridge have been proposed to explain the flat slab segment below south peru we have conducted a series of numerical model experiments to investigate the relative importance of both mechanisms results suggest that the average upper mantle viscosity should be about 3 5 1020 pa s or less and basaltic crust should be able to survive 600 to 800 c ambient temperature before transforming into eclogite to explain the slab geometry below peru the effect of the overriding plate is estimated to be as large or twice as large as that of the plateau subduction

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

modeling the dynamics of subducting slabs cold dense subducting lithosphere provides the primary force driving tectonic plates at earth s surface the force available to drive the plates depends on a balance between the buoyancy forces driving subduction and the mechanical and buoyancy forces resisting subduction because both the buoyancy and rheology of the slab and mantle depend on temperature composition grain size water content and melt fraction unraveling which of these processes exert a first order control on slab dynamics and under what circumstances other processes become first order effects can be challenging laboratory and numerical models of slab dynamics provide a powerful method for testing the combined effects of buoyancy and strength changes that accompany the slab evolution in the upper mantle transition zone and lower mantle recent studies have focused on understanding how rheologic variations newtonian versus non newtonian viscosity or water content geometry 2d versus 3d and plate motions trench roll back or advance influence the evolution of slabs in the upper mantle and how they sink into the lower mantle these models suggest that spatial and temporal variations in slab strength and the history of subduction determine whether slabs sink directly into the lower mantle or are trapped in the transition zone

the yakutat terrane dramatic change in crustal thickness across the transition fault alaska we present new constraints on the crustal structure of the yakutat terrane and evidence of the role of the transition fault in southern alaska the yakutat terrane south of yakutat bay includes crystalline crust that is 24 27 km thick overlain by sedimentary units that are 4 5 7 5 km thick the yakutat terrane crustal thickness and velocity structure are consistent with an oceanic plateau origin the southern edge of the yakutat terrane is bounded by the transition fault which is imaged as a near vertical fault zone 1 km wide the transition fault is coincident with a dramatic change in moho depth from 32 km for yakutat oceanic plateau crust to 11 5 km for pacific ocean crust occurring over a horizontal distance of 0 5 km there is no evidence for underthrusting of the pacific ocean crust beneath the yakutat terrane at the transition fault we argue that the yakutat terrane formed on the kula or farallon plate and was later juxtaposed next to younger pacific ocean crust by the transition fault

orogens and slabs vs their direction of subduction abstract subduction zones appear primarily controlled by the polarity of their direction i e w directed or e to nne directed probably due to the westward drift of the lithosphere relative to the asthenosphere the decollement planes behave differently in the two end members in the w directed subduction zone the decollement of the plate to the east is warped and subducted whereas in the e to nne directed it is ramping upward at the surface there are w directed subduction zones that work also in absence of active convergence like the carpathians or the apennines w directed subduction zones have shorter life 30 40 ma than e or ne directed subduction zones even longer than 100 ma the different decollements in the two end members of subduction should control different ptt paths and therefore generate variable metamorphic assemblages in the associated accretionary wedges and orogens these asymmetries also determine different topographic and structural evolutions that are marked by low topography and a fast eastward migrating structural wave along w directed subduction zones whereas the topography and the structure are rapidly growing upward and expanding laterally along the opposite subduction zones the magmatic pair calc alkaline and alkaline tholeiitic volcanic products of the island arc and the back arc basin characterise the w directed subduction zones magmatic rocks associated with e or ne directed subduction zones have higher abundances of incompatible elements and mainly consist of calc alkaline shoshonitic suites with large volumes of batholithic intrusions and porphyry copper ore deposits the subduction zones surrounding the adriatic plate in the central mediterranean confirm the differences among subduction zones as primarily controlled by the geographic polarity of the main direction of the slab the western margin of the adriatic plate contemporaneously overridden and underthrust europe toward the west to generate respectively the alps and the apennines while the eastern margin subducted under the dinarides hellenides these belts confirm the characters of the end members of subduction zones as a function of their geographic polarity similarly to the pacific subduction zones

the potential influence of subduction zone polarity on overriding plate deformation trench migration and slab dip angle abstract a geodynamic model exists the westward lithospheric drift model in which the variety of overriding plate deformation trench migration and slab dip angles is explained by the polarity of subduction zones the model predicts overriding plate extension a fixed trench and a steep slab dip for westward dipping subduction zones e g mariana and predicts overriding plate shortening oceanward trench retreat and a gentle slab dip for east to northeastward dipping subduction zones e g chile this paper investigates these predictions quantitatively with a global subduction zone analysis the results show overriding plate extension for all dip directions azimuth 180 to 180 and overriding plate shortening for dip directions with 90 to 110 the wide scatter in data negate any obvious trend and only local mean values in overriding plate deformation rate indicate that overriding plate extension is somewhat more prevalent for west dipping slabs west dipping subduction zones are never fixed irrespective of the choice of reference frame while east to northeast dipping subduction zones are both retreating and advancing in five out of seven global reference frames in addition westward dipping subduction zones have a range in trench migration velocities that is twice the magnitude of that for east to northeastward dipping slabs finally there is no recognizable correlation between slab dip direction and slab dip angle east to northeast dipping slabs 30 to 120 have shallow 0 125 km slab dip angles in the range 10 60 and deep 125 670 km slab dip angles in the range 40 82 while west dipping slabs 60 to 120 have shallow slab dip angles in the range 19 50 and deep slab dip angles in the range 25 86 local mean deep slab dip angles are nearly identical for east and west dipping slabs while local mean shallow slab dip angles are lower by only 4 7 8 1 for east to northeast dipping slabs it is thus concluded that overall there is no observational basis to support the three predictions made by the westward drift model and for some sub predictions the observational basis is very weak at most alternative models which incorporate and underline the importance of slab buoyancy driven trench migration slab width and overriding plate motion are better candidates to explain the complexity of subduction zones including the variety in trench migration velocities overriding plate deformation and slab dip angles

a numerical investigation of lutocline dynamics and saturation of fine sediment in the oscillatory boundary layer 1 understanding the state of the muddy seabed is critical to sediment transport hydrodynamic dissipation and seabed properties however this endeavor is challenging because the availability and settling velocity of sediment in muddy environments are highly variable for a given reynolds number typical of fine sediment settling in a moderately energetic muddy shelf recent 3d numerical simulations have revealed four distinct regimes of wave induced fine sediment transport these regimes depend on the availability or concentration of sediment and range from well mixed condition to the formation of lutocline and eventually a complete flow laminarization by keeping the sediment availability unchanged this study further demonstrates the existence of these flow regimes for a range of sediment settling velocities simulation results suggest that when settling velocity is larger the location of the lutocline becomes lower closer to the bed and the flow eventually laminarizes when there is further increase in the settling velocity hence the dynamics of lutocline is clearly related to the transition between these flow regimes the vertical flow structure in the presence of lutocline is revealed through the budgets of sediment flux and turbulent kinetic energy the suppressed mixing in the lutocline layer is further illustrated from a new perspective i e the budgets of turbulent suspension and concentration fluctuation variance the concept of saturation commonly used for tidal boundary layer is extended here for wave boundary layer

influence of trench width on subduction hinge retreat rates in 3 d models of slab rollback subduction of tectonic plates limited in lateral extent and with a free trailing tail i e free subduction is modeled in a three dimensional 3 d geometry the models use a nonlinear viscoplastic rheology for the subducting plate and exhibit a wide range of behaviors depending on such plate characteristics as strength width and thickness we investigate the time evolution of this progressive rollback subduction measure the accompanying return flow in the upper mantle and quantify the plate kinematics due to the 3 d geometry flow is allowed to accompany slab rollback around the lateral edges of the slab the toroidal component as opposed to 2 d geometry where material is forced to flow underneath the slab tip the poloidal component a simple force balance is provided which relates the speed of backward trench migration to the resistive forces of generating flow and weakening the plate our results indicate most of the gravitational energy of the system i e the negative buoyancy of the subducting slab is converted into a toroidal flow 69 a much smaller amount goes into weakening the plate 18 and the remaining amount goes into driving flow parallel to displacement of the slab 13 for the trench widths w we investigate 1500 km a maximum trench retreat rate occurs for trenches 600 km wide which is attributed to the interaction between a plate of finite width and the induced flow which has a lengthscale in the horizontal direction these numerical results quantitatively agree with comparable 3 d laboratory experiments using analogue models with a purely viscous plate material schellart 2004a 2004b including correlations between increasing retreat rate with increasing plate thickness trench width for maximum retreat rate 500 km and estimated amount of slab buoyancy used to drive rollback induced flow 70 several implications for plate tectonics on earth result from these models such as rollback subduction providing a physical mechanism for ephemeral slab graveyards situated above the more viscous lower mantle and endothermic phase transition prior to a flushing event into the lower mantle mantle avalanche

the role of boundary conditions in numerical models of subduction zone dynamics abstract numerical studies of subduction zone dynamics on a regional scale usually use a limited spatial extent for their models and therefore need to define boundary conditions on all model edges these boundary conditions imply a choice for the mechanical and thermal state of the surrounding regions which may influence the evolution of the model system we investigate the role of the surface and lateral boundary conditions for simple mechanical subduction models using a standard finite element method we find that subduction is only possible if the slab can decouple from the surface this decoupling can be achieved by a true free surface a free slip surface with a weak crust for the subducting plate or a free slip surface with a layer of low viscosity low density material sticky air between the model surface and the crust models of slab dynamics that employ a free slip surface reproduce trench migration slab sinking velocities and slab geometry of models with a free surface as long as they use either a weak crust which can be viscous viscoelastic and or brittle or a sticky air layer the vertical topography will however not be reproduced for free slip models without a sticky air layer for ocean ocean convergent models we find that the application of inflow boundary conditions at the edges of the subducting or overriding lithosphere controls trench motion and the geometry of the subducting slab inflow on the overriding side causes trench retreat and a slab resting on the lower mantle whereas inflow restricted to the subducting side can cause trench advance and a slab which folds on the lower mantle

subduction zone evolution and low viscosity wedges and channels dehydration of subducting lithosphere likely transports fluid into the mantle wedge where the viscosity is decreased such a decrease in viscosity could form a low viscosity wedge lvw or a low viscosity channel lvc on top of the subducting slab using numerical models we investigate the influence of low viscosity wedges and channels on subduction zone structure slab dip changes substantially with the viscosity reduction within the lvws and lvcs for models with or without trench rollback overthickening of slabs is greatly reduced bylvws or lvcs two divergent evolutionary pathways have been found depending on the maximum depth extentofthelvwandwedgeviscosity assumingaviscositycontrastof0 1withbackgroundasthenosphere modelswithalvwthat extendsdown to 400km depthshow asteeply dipping slab while models with anlvw that extends to muchshallower depth such as 200 km can produce slabs that are flat lying beneath the overriding plate there is a narrow range of mantle viscosities that produces flatslabs 5to10 10 19 pas andthe slab flatteningprocessisenhanced bytrenchrollback slabcanbedecoupled fromthe overriding platewithalvcifthe thicknessisatleastafew10s ofkm the viscosityreduction isatleastafactor oftwoandthedepthextentofthe lvc is several hundred km these models have important implications for the geochemical and spatial evolution of volcanic arcs and the state of stress within the overriding plate the models explain the poor correlation between traditional geodynamic controls subductingplateageandconvergencerates onslabdip wepredictthatwhenvolcanicarcschangetheirdistancefromthetrench they could be preceded by changes in arc chemistry we predict that there could be a larger volatile input into the wedge when arcs migrate toward the trench and visa versa the transition of a subduction zone into the flat lying regime could be preceded by changes in the volatile budget such that the dehydration front moves to shallower depths our flat slab models shed some light on puzzling flat subductionsystems likeincentralmexico wherethereisnodeformationwithintheoverridingplateabovetheflatsegment thelack

slab stiffness control of trench motion insights from numerical models 1 subduction zones are not static features but trenches retreat roll back or advance here we investigate the dominant dynamic controls on trench migration by means of two and three dimensional numerical modeling of subduction this investigation has been carried out by systematically varying the geometrical and rheological model parameters our viscoplastic models illustrate that advancing style subduction is promoted by a thick plate a large viscosity ratio between plate and mantle and a small density contrast between plate and mantle or an intermediate width w 1300 km advancing slabs dissipate 45 to 50 of the energy in the system thin plates with relatively low viscosity or relatively high density or wide slabs w 2300 km on the other hand promote subduction in the retreating style i e slab roll back the energy dissipated by a retreating slab is 35 to 40 of the total dissipated energy most of the energy dissipation occurs in the mantle to accommodate the slab motion whereas the lithosphere dissipates the remaining part to bend and unbend with a simple scaling law we illustrate that this complex combination of model parameters influencing trench migration can be reduced to a single one plate stiffness stiffer slabs cause the trench to advance whereas more flexible slabs lead to trench retreat the reason for this is that all slabs will bend into the subduction zone because of their low plastic strength near the surface but stiff slabs have more difficulty unbending at depth when arriving at the 660 km discontinuity those bent slabs tend to cause the trench to advance in a similar way variation of the viscoplasticity parameters in the plate may change the style of subduction a low value of friction coefficient weakens the plate and results in a retreating style while higher values strengthen the plate and promote the advancing subduction style given the fact that also on earth the oldest and therefore probably stiffest plates have the fastest advancing trenches we hypothesize that the ability of slabs to unbend after subduction forms the dominant control on trench migration

influence of the zone of weakness on dip angle and shear heating of subducted slabs 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 2004 elsevier b v all rights reserved

a regime diagram for subduction styles from 3 d numerical models of free subduction abstract previous models of subduction both analogue and numerical have observed a number of distinct styles of subduction each with particular subduction motions partitioned between slab rollback and forward plate advance and associated slab morphologies we use 3 d numerical models to investigate subduction dynamics by varying the strength of slabs as well as the buoyancy and propose a new classification based on these parameters the slab strength is specified both through the ratio of viscosities between the subducting plate and upper mantle plate um as well as the plate thickness h plate only a very restrictive range of plates strong plates with smaller buoyancy tend to favor modes of subduction which are exclusively advancing plates which have greater negative buoyancy will eventually transition into a retreating style we find that the flexural strength and the buoyancy determine the subduction style as distinguished by a characteristic slab morphology and control several subduction characteristics including the partitioning between slab rollback and plate advance the trench curvature and the slab s radius of curvature plates that are 80 100 km thick with plate um 100 300 are classified here as weak and are the only plates to exhibit slab geometry with several recumbent folds atop the more viscous lower mantle this regime of weak plates with their associated slab morphologies predominant folding is argued to be most similar to slabs on earth based on the presence of folded slab piles in earth s upper mantle as interpreted from seismic tomography

shape of the subducted rivera and cocos plates in southern mexico seismic and tectonic implications the geometry of the subducted rivera and cocos plates beneath the north american plate in southern mexico was determined based on the accurately located hypocenters of local and teleseismic earthquakes the hypocenters of the teleseisms were relocated and the focal depths of 21 events were constrained using a body wave inversion scheme the suduction in southern mexico may be approximated as a subhorizontal slab bounded at the edges by the steep subduction geometry of the cocos plate beneath the caribbean plate to the east and of the rivera plate beneath north america to the west the dip of the interplate contact geometry is constant to a depth of 30 km and lateral changes in the dip of the subducted plate are only observed once it is decoupled from the overriding plate on the basis of the seismicity the focal mechanisms and the geometry of the downgoing slab southern mexico may be segmented into four regions 1 the jalisco region to the west where the rivera plate subducts at a steep angle that resembles the geometry of the cocos plate beneath the caribbean plate in central america 2 the michoacan region where the dip angle of the cocos plate decreases gradually toward the southeast 3 the guerrero oaxaca region bounded approximately by the onshore projection of the orozco and o gorman fracture zones where the subducted slab is almost subhorizontal and underplates the upper continental plate for about 250 km and 4 the southern oaxaca and chiapas region in southeastern mexico where the dip of the subduction gradually increases to a steeper subduction in central america these drastic changes in dip do not appear to take place on tear faults suggesting that smooth contortions accommodate these changes in geometry the inferred 80 and 100 km depth contours of the subducted slab lie beneath the southern front of the trans mexican volcanic belt suggesting it is directly related to the subduction thus the observed nonparallelism with the middle american trench is apparently due to the changing geometry of the rivera and cocos plates beneath the north american plate in southern mexico and not to zones of weakness in the crust of the north american plate as some authors have suggested

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

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

time dependent interaction between subduction dynamics and phase transition kinetics summary numerical subduction models are presented that account for the kinetics of the olivine spinel transformation and the equilibrium transformation of spinel into perovskite and magnesiowustite temperature depth and stress dependent rheology is used buoyancy due to temperature and the olivine spinel and magnesiowustite perovskite phases leads to free subduction into the transition zone and retarded penetration into the lower mantle a dynamic feed back mechanism is found between the different buoyancy forces the kinetics of the phase transition and the release of latent heat this mechanism modulates the parachute effect of the metastable olivine mo in a subducting slab and leads to characteristic temporal variations of the depth of mo and with a time delay variations of the subduction velocity time periods are of the order of 3 4 myr this mechanism allows to explain variations of the depths of the deepest earthquakes in otherwise similar subduction zones further conclusions are as follows the maximum depth of mo increases roughly linearly with lithospheric thickness and the maximum depth may reach as deep as 720 750 km different subduction zones may represent different stages of velocity depth phase loops of mo

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

how valid are dynamic models of subduction and convection when plate motions are prescribed a detailed comparison between fully dynamic and kinematic plate formulations has been made in models of mantle convection plate velocity is computed self consistently from fully dynamic plate models with temperature and stress dependent viscosity and preexisting mobile faults in fully dynamic models the flow is driven solely by internal buoyancy while in kinematic models the flow is driven by a combination of the prescribed surface velocity and internal buoyancy only a temperature dependent viscosity close to the effective viscosity determined from the fully dynamic models is used in the kinematic models the two types of models give very similar temperature structures and slab evolutionary histories when the effective viscosity and surface velocity are nearly identical in kinematic plate models the additional work introduced by the prescribed velocity boundary condition is apparently dissipated within the lithosphere and has little influence on the convection under the lithosphere in models with periodic lateral boundary conditions slabs sink into the lower mantle at an oblique angle and this contrasts with the vertical sinking which occurs with reflecting boundary conditions models show that we can simulate fully dynamic models with kinematic models under either periodic boundary conditions or reflecting boundary conditions

the influence of temperature and depth dependent viscosity on geoid and topography profiles from models of mantle convection abstract we consider the geoid and topography profiles from convection models with depth and temperature dependent viscosity models in a 2 d cartesian geometry we use a very low viscosity region in the top center of the geometry to generate plate like viscosity models we focus on subduction regions which have narrow 100 200 km trenches the narrow trench is only a minor feature in the regional geoid in depth dependent viscosity models we change the geometry of the weak zone around the ocean trench to study the impact of the weak zone on trench topography and the geoid we also vary the vertical viscosity structure of the models the weak zone geometry has very little effect on the surface velocity of the ocean plate or the large scale pattern of convection the width of the weak zone within a reasonable range has little effect on topography and geoid profiles there is a narrow trough in the geoid and topography profiles this trough becomes wider and deeper with a high viscosity upper mantle this may suggest that the high viscosity upper mantle is not suitable for our models in temperature dependent viscosity models we vary the viscosity values for different layers temperature dependent viscosity is also important to producing realistic subduction models however the activation energy must be weaker than laboratory estimates of olivine under mantle pressures and temperatures from both depth and temperature dependent viscosity models the viscosity in the lower mantle should be at least ten times greater than in the upper mantle consistent with previous studies of the geoid in regions of subduction the topography and geoid profiles can match the observations when the viscosity increases with depth

kinematics of subduction and subduction induced flow in the upper mantle 1 results of fluid dynamical experiments are presented to model the kinematics of lithospheric subduction in the upper mantle the experiments model a dense high viscosity plate subducting lithosphere overlying a less dense low viscosity layer upper mantle the overriding lithosphere is not incorporated several important features of slab behavior were investigated including the temporal variability of hinge line migration the kinematic behavior of the slab and the subduction induced upper mantle flow both fixed and free trailing edge boundary conditions of the subducting plate were investigated results show that hinge line retreat is a natural consequence of subduction of a negatively buoyant slab the migration rate increases until the slab approaches the upper lower mantle discontinuity resulting in a decrease in migration rate followed by a renewed increase and finally approaching a steady state slab retreat results in mantle flow with material initially located underneath the slab flowing around the lateral slab edges toward the mantle wedge experimental results indicate that all rollback induced flow occurs around the lateral slab edges forcing the hinge line to attain a convex shape toward the direction of retreat no signs for poloidal flow underneath the slab tip have been detected only a small component of toroidal type flow was observed underneath slanting slab tips for a fixed trailing edge the slab does not sink vertically downward but sinks at an angle in a regressive manner for a free trailing edge slab sinking is oriented more vertically while the surface part of the subducting plate is pulled into the subduction zone

upper plate controls on deep subduction trench migrations and deformations at convergent margins abstract thus far relatively simplistic models of free subduction in which the trench and plate motions are emergent features completely driven by the negative buoyancy of the slab have investigated the dynamics of a single isolated subducting plate here we extend such models to incorporate an overriding plate and present the results of how such an overriding plate feedbacks into the dynamics of free subduction in this study we address three fundamental aspects of these dynamics 1 how does the presence of an overriding plate change the force balance at the convergent margins 2 how are the forces from deep subduction propagated to the surface and 3 what controls the stress regime in a system of coupled upper and subducting plates and how is it expressed in the deformations and plate motions in general we find that the evolution of subduction zones is strongly controlled by both the interactions between the slab and the upper lower mantle discontinuity as well as the strength of the upper plate when either the subducting or upper plates are unable to move subduction motions are steady state and partitioned entirely into either slab rollback or plate advance respectively when conditions favour a quasi stationary trench the subducted lithosphere can form into a pile with multiple recumbent folds of slab material atop the lower mantle alternating between forwards and backwards draping slab the corresponding horizontal trench motions at the surface are frontward and rearward respectively resulting in either a compressive or extensional regime in the back arc time dependent forcing arising from the slab piling behaviour can have a feedback with upper plate and produce strongly non steady state intermittent phases of upper plate deformation as those commonly observed on earth two types of discontinuous back arc strain evolution are identified 1 periodic when recurrent phases of strain over finite durations are accommodated by viscous stretching thickening of the plate and 2 episodic when upper plate deformation localizes plastic strain and allows for punctuated episodes these phases can include extension quiescence and compression giving rise to a large variety of possible tectonic evolutions the models presented here provide insight into the dynamics behind the non steady state evolution of subduction which can help unravel seemingly erratic motions of major convergent margins and back arc deformations around the pacific and indian oceans during the cenozoic

effective elastic thickness variations along the andean margin and their relationship to subduction geometry 1 we present a new map of the spatial variations in effective elastic thickness te along the andes estimated using bouguer coherence the te variations reflect interactions between subducting slab and preexisting terrane structure in the forearc conductive cooling of the continent by the subducting slab exerts primary control on rigidity resulting in te that is highest 40 km where the oceanic lithosphere is oldest and coldest 20 s in the central andes te is relatively low 20 km along the volcanic chain and the altiplano and puna plateaus we interpret this weakening to reflect a high geothermal gradient maintained by advective magmatic processes a shallow and hot asthenosphere and a very weak lower crust throughout this region east of the plateaus high te delineates underthrusting of the brazilian shield finally north and south of the plateaus flat subduction areas are characterized by high te high shear wave velocity thick thermal lithosphere and low heat flow indicating that continental lithosphere there is thicker colder and stronger on the basis of these relationships we suggest that variations in slab dip along the margin relate to variations in structure of the continental lithosphere in particular we propose that upper plate structure influences the width and viscosity of the asthenospheric wedge which control the suction moment responsible for the subduction angle at depths 70 100 km for example when oceanic lithosphere subducts beneath thin continental lithosphere the low viscosity asthenosphere allows the slab to detach from the continent and sink into the mantle at normal angles however when oceanic lithosphere subducts close or beneath thick and strong continental lithosphere the asthenospheric wedge narrows and corner flow drags high viscosity mantle from the base of the thick 150 km cold continent into the wedge suction forces increase with both narrowing of the wedge and its increasing viscosity we estimate the asthenospheric viscosity underneath thick cold continents to be 1020 pa s sufficient to induce flat subduction later after prolonged hydration and weakening of the continent s base asthenospheric flow into the wedge may resume allowing the slab to sink again into the mantle at normal angles

rheologic controls on slab dynamics several models have been proposed to relate slab geometry to parameters such as plate velocity or plate age however studies on the observed relationships between slab geometry and a wide range of subduction parameters show that there is not a simple global relationship between slab geometry and any one of these other subduction parameters for all subduction zones numerical and laboratory models of subduction provide a method to explore the relative importance of different physical processes in determining subduction dynamics employing 2 d numerical models with a viscosity structure constrained by laboratory experiments for the deformation of olivine we show that the observed range in slab dip and the observed trends between slab dip and convergence velocity subducting plate age and subduction duration can be reproduced without trench motion i e slab roll back for locations away from slab edges successful models include a stiff slab that is 100 1000 times more viscous than previous estimates from models of plate bending the geoid and global plate motions we find that slab dip in the upper mantle depends primarily on slab strength and plate boundary coupling with a small dependence on subducting plate age once the slab sinks into the lower mantle the primary processes controlling slab evolution are 1 the ability of the stiff slab to transmit stresses up dip 2 resistance to slab descent into the higher viscosity lower mantle and 3 subduction induced flow in the mantle wedge corner

on the relationships between slab dip back arc stress upper plate absolute motion and crustal nature in subduction zones 1 statistical analysis of modern oceanic subduction zone parameters such as the age of a downgoing plate or the absolute plate motions is performed in order to investigate which parameter controls the dip of a slab and conversely what the influence of slab geometry is on upper plate behavior for that purpose parameters have been determined from global databases along 159 transects from all subduction zones that are not perturbed by nearby collision or ridge plateau seamount subduction on the basis of tomographic images slabs that penetrate through or lie on the 670 km discontinuity are also identified the results of the statistical analysis are as follows 1 back arc stress correlates with slab dip i e back arc spreading is observed for deep dips deeper than 125 km larger than 50 whereas back arc shortening occurs only for deep dips less than 30 2 slab dip correlates with absolute motion of the overriding plate the correlation is even better when the slab lies on or even more penetrates through the 670 km discontinuity 3 slabs dip more steeply by about 20 on average beneath oceanic overriding plates than beneath continental ones 4 slabs dip more steeply on average by about 10 near edges 5 slab dip does not correlate with the magnitude of slab pull the age of subducting lithosphere at the trench the thermal regime of the subducting lithosphere the convergence rate or the subduction polarity east versus west the present study provides evidence that the upper plate absolute motion plays an important role on slab dip as well as on upper plate strain retreating overriding plates are often oceanic ones and thus may partially explain the steeper slab dips beneath oceanic upper plates one can infer that low slab dips correlate well with compression in continental advancing upper plates whereas steep dips are often associated with extension in oceanic retreating upper plates excess weight of old slabs is often counterbalanced by other forces probably asthenospheric in origin such as lateral mantle flow near slab edges or anchor forces to determine slab dip

can earth s rotation and tidal despinning drive plate tectonics abstract we re evaluate the possibility that earth s rotation contributes to plate tectonics on the basis of the following observations 1 plates move along a westerly polarized flow that forms an angle relative to the equator close to the revolution plane of the moon 2 plate boundaries are asymmetric being their geographic polarity the first order controlling parameter unlike recent analysis the slab dip is confirmed to be steeper along w directed subduction zones 3 the global seismicity depends on latitude and correlates with the decadal oscillations of the excess length of day lod 4 the earth s deceleration supplies energy to plate tectonics comparable to the computed budget dissipated by the deformation processes 5 the gutenberg richter law supports that the whole lithosphere is a self organized system in critical state i e a force is acting contemporaneously on all the plates and distributes the energy over the whole lithospheric shell a condition that can be satisfied by a force acting at the astronomical scale assuming an ultra low viscosity layer in the upper asthenosphere the horizontal component of the tidal oscillation and torque would be able to slowly shift the lithosphere relative to the mantle

formation of the rocky mountains western united states a continuum computer model one hypothesis for the information of the rocky mountain structures in late cretaceous through eocene time is that plate of oceanic lithosphere was underthrust horizontally along the base of the north american lithosphere the horizontal components of the motion of this plate are known from paleomagnetism and the edge of the region of flat slab can estimated from reconstructed patterns of volcanism new techniques of finite element modeling allow prediction of the thermal and mechanical effects of horizontal subduction on the north american plate a model that has a realistic temperature dependent rheology and a simple plane layered initial condition is used to compute the consequences of horizontal underthrusting in the time interval 75 million to 30 million years before present successful prediction of this model include i the location amount and direction of horizontal shortening that has been inferred from laramide structures ii massive transport of lower crust from southwest to northeast iii the location and timing of the subsequent extension in metamorphic core complexes and the rio grande rift and iv the total area eventually involved in basin and range style extension in a broad sense this model has predicted the belt of laramide structures the transport of crust from the coastal region to the continental interior the subsequent extension in metamorphic core complexes and the rio grande rift and the geographic region of late tertiary basin and range extension its principal defects are that i many events are predicted about 5 million to 10 million years too late and ii the wave of crustal thickening does not travel far enough to the east reasonable modifications to the oceanic plate kinematics and rheologies that were assumed may correct these defects the correspondence of model predictions to actual geology is already sufficiently close to show that the hypothesis that horizontal subduction caused the laramide orogeny is probably correct the rocky mountain thrust and reverse faults formed in an environment of east west to northeast southwest compressive stress that was caused by the viscous coupling between the oceanic plate and the base of the north american crust nonuniform crustal thickening by simple shear transport also caused relative uplifts therefore this model is consistent with both of the range forming mechanisms that have been inferred 1 a new proposal that arises from this simulation is that horizontal subduction also caused the subsequent extensional basin and range taphrogeny by stripping away the mantle lithosphere so that the crust was exposed to hot asthenosphere after the oceanic slab dropped away

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

plate bending at subduction zones consequences for the direction of plate motions abstract bending of lithospheric plates at subduction zones is thought to be an important source of dissipation for convection in the earth s mantle however the influence of bending on plate motion is uncertain here we use a variational description of mantle convection to show that bending strongly affects the direction of plate motion subduction of slabs and subsidence of oceanic lithosphere with age provide the primary driving forces dissipation is partitioned between plate bending and various sources of friction at plate boundaries and in the interior of the mantle due to viscous flow we determine the poles of rotation for the pacific and nazca plates by requiring the net work to be stationary with respect to small changes in the direction of motion the best fit to the observed rotation poles is obtained with an effective lithospheric viscosity of 6 1022 pa s bending of the pacific plate dissipates roughly 40 of the energy released by subduction through the upper mantle

role of the overriding plate in the subduction process insights from numerical models active convergent margins are primarily shaped by the interplay among the subducting plate overriding plate and mantle the effect of important forces like far field mantle flow overriding plate motion and inter plate coupling however remains partially ambiguous in a preliminary attempt to clarify their role a self consistent viscoelastic plane strain mechanical finite element model in which subducting plate overriding plate and mantle interact dynamically is developed in this quasi static framework with a freely moving slab trench and inter plate fault the role of a compressive overriding plate on subduction zone kinematics morphology and stress state is characterized a slab interacting solely with a semi analytical three dimensional mantle flow formulation shows that local non induced mantle flow influences slab geometry and kinematics adding an important dynamic term to the system the impact of an overriding plate on this system is determined completely by overriding plate trench ward motions and is only pertinent if the overriding plate actively advances the trench a trench ward moving overriding plate indents the slab and thereby enforces trench retreat and decreases slab dip it also stimulates over thrusting of the overriding plate onto the slab and thereby permits mountain building within the overriding plate frictional resistance is observed to have a dominant local effect within the overriding plate as it is increasingly dragged down thereby inhibiting the growth of overriding plate topography a distinguishable effect on large scale trench motions and deep slab dip is however absent for re normalized friction coefficients ranging up to about 0 2 minor additional effects include a decrease in plate motions of about 15 and slab bending stresses of about 10