controls on back arc basin formation the relationship between subduction and back arc spreading has been well known since the early days of plate tectonics however the reasons why back arc basins are associated with some subduction systems but not all has remained elusive we examine the kinematic controls on subduction and back arc basins for both the present day and cenozoic to differentiate between the major competing hypotheses for back arc basin formation and to explain their temporal and spatial distribution our new data set of subduction and back arc basin parameters uses a new set of paleo oceanic age grids muller et al 2005 associated with a moving atlantic indian ocean hot spot reference frame o neill et al 2005 the plate model includes detailed reconstructed spreading histories of back arc basins based on marine geophysical and satellite gravity data our combined rotation and oceanic paleo age model provides the age distribution of subducting lithosphere through space and time convergence rates and the absolute motion of the downgoing and overriding plates we find that back arc basins develop when the age of subducting normal oceanic lithosphere is greater than 55 million years additionally we establish an age dip relationship showing that the intermediate dip angle of the subducting slab is always greater than 30 with back arc spreading our results suggest that back arc basin formation is always preceded by an absolute motion of the overriding plate away from the subduction hinge thereby creating accommodation space between the overriding and subducting plates once back arc extension is established it continues regardless of overriding plate motion indicating back arc spreading is not a simple consequence of overriding plate behaviour the landward migration of the overriding plate as a precursor to back arc extension may indicate that extension on the overriding plate is influenced by the oceanward lateral flow of the mantle however once back arc extension is established rollback of the subduction hinge appears to be the primary force responsible for the continued creation of accommodation space our analysis indicates the driving mechanism for back arc extension is a combination of surface kinematics properties of the downgoing slab the effect of lateral mantle flow on the slab and mantle wedge dynamics

constraints on the lateral strength of slabs from three dimensional dynamic flow models if the viscosity of subducted lithosphere is described purely by temperature and pressure it should remain considerably more viscous than the surrounding material as it descends through the upper mantle many dynamic models of the long wavelength geoid associated with slabs require an increase in viscosity from the upper to the lower mantle but assume the perturbation to the lateral viscosity structure arising from the slabs themselves can be ignored previous studies indicate that strong and localized viscosity variations should have a dramatic influence on the geoid we present 3d finite element models of the regional geoid of the western pacific subduction zones slab buoyancies and viscosities are defined using the distribution of seismicity the geoid is very sensitive to the lateral strength of the slab very viscous slabs penetrating a low viscosity mantle generate significant 10 50 m long wavelength geoid lows opposite to the geoid high which is observed over slabs to obtain a geoid high comparable to that observed the lower mantle viscosity must be 60 200 times greater than the upper mantle viscosity and the slab must be in contact with the lower mantle these strict requirements suggest that slabs have been weakened and cannot act as stress guides from the deep mantle to the surface

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

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

haskell 1935 revisited the haskell 1935 value of 1021 pa s for mantle viscosity is a classic and enduring constraint on the rheology of the earth s interior we revisit this inference using spherically symmetric self gravitating maxwell viscoelastic earth models our inference is based on both forward and inverse analyses of decay times associated with uplift at two sites considered by haskell angerman river sweden and oslo norway rather than the raw relative sea level rsl data at these sites we demonstrate that predictions of the decay time associated with the angerman river data are insensitive to variations in both the late pleistocene ice load history and the lithospheric thickness of the earth model the predictions at oslo are sensitive to both these inputs and hence decay times at this site provide a remarkably robust constraint on mantle viscosity we derive a constraint on the average viscosity of the mantle of 0 65 1 10 1021 pa s where the average resolved by the data encompasses a region which extends from the base of the lithosphere to a depth of near 1400 km this indicates that many previous analyses which have invoked the haskell value of 1021 pa s as a constraint on the average upper mantle i e above 670 km depth viscosity alone have misinterpreted the resolving power of the inference furthermore our analysis indicates that a number of apparently contradictory inferences of viscosity based on fennoscandian data satisfy the new rigorous interpretation of the haskell constraint finally we demonstrate how the ambiguity in the upper mantle lower mantle viscosity contrast associated with the haskell average may be reduced by invoking decay time constraints estimated from rsl curves in hudson bay a preliminary inversion of decay times at angerman river and richmond gulf in hudson bay suggests a contrast of approximately an order of magnitude between the average viscosities of the upper mantle and the top 1000 km of the lower mantle however a conclusive analysis in this regard must await the determination of consistent decay time estimates for the hudson bay region

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

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

aeolian dune field geomorphology modulates the stabilization rate imposed by climate 1 the activity of inland aeolian dune fields is typically related to the external forcing imposed by climate active bare dunes are associated with windy and or arid settings and inactive vegetated dunes are associated with humid and or calm environments when a climate shifts the dune field reacts however the behavior rate and potential impact of diverse dune geomorphologies on these transitions are poorly understood here we use a numerical model to systematically investigate the influence of dune field geomorphology dune height organization and collisions on the time a dune field takes to stabilize to generate diverse initial un vegetated dune field geomorphologies under unidirectional winds we varied pre stabilization growth time and initial sediment thickness termed equivalent sediment thickness est following dune field development from a flat bed we introduced vegetation simulating a climate shift and transport vegetation feedbacks slowly stabilized the dune fields qualitatively very young and immature dune fields stabilized quickly whereas older dune fields took longer dune fields with greater est stabilized quicker than those with less est larger dunes stabilized quicker because of low celerity which facilitated higher vegetation growth rates extended stabilization times were associated with the extension of parabolic dunes dune dune collisions resulted in premature stabilization the frequency of collisions was related to dune spacing quantitatively comparing the distribution of deposition rates in a dune field to the deposition tolerance of vegetation provides a promising predictor of relative stabilization time dune fields with deposition rates dominantly above the deposition tolerance of vegetation advanced unimpeded and prolonged stabilization as parabolic dunes paleoenvironmental reconstructions or predictions of dune field activity should not assume that dune activity directly translates to climate considerable lags to stabilizing climate shifts may exist in unidirectional dune forms

dynamic models of downgoing plate buoyancy driven subduction subduction motions and energy dissipation abstract it is much debated whether the forces associated with the downgoing plate the overriding plate passive or active mantle flow are dominant in controlling the paths of plates into the mantle we investigate the dynamics and energetics for a free subduction system driven solely by downgoing plate buoyancy using a finite element model of a viscoelastic plate with a free surface sinking into a passive unbounded mantle represented by drag forces parameters are varied to study effects of an asthenosphere ridge push and a passive overriding plate for a range of subducting plate viscosities and densities such a single free plate achieves subduction mainly through trench retreat most of the energy dissipation occurs in driving the passive mantle response as a result the slab s sinking velocity is its stokes velocity determined by lithospheric buoyancy and mantle viscosity the total subduction velocity and dip adjust to minimize bending dissipation in the lithosphere and are affected by slab rheology as well as buoyancy a low viscosity asthenosphere and ridge push facilitate plate advance increasing plate dips and lowering subduction velocity while suction and buoyancy of a work free passive overriding plate decreases plate dips thus increasing subduction and rollback velocities however the geometrical relation between the different parameters is the same in all model cases because the slabs sink according to their stokes velocity the free subduction models thus provide a reference that can be used to distinguish the signature of downgoing plate buoyancy from that of other driving forces in global compilations of subduction parameters

the driving mechanism of plate tectonics relation to age of the lithosphere at trenches sea floor age depth and age trench depth relations suggest that the oceanic lithosphere continues to thicken and subside with age beyond 80 ma if so slab pull fs and ridge push fr forces can be calculated from the age of the sea floor at trenches t and if the shear stress acting on the base of the plate increases with absolute velocity a the approximate force balance equation is dva2 svat32 rt where d s and r are constants related to drag slab pull and ridge push respectively the relation between velocity and age is va2 s vat12 r we have tested this model by linear regression using rate and age data for 15 trenches va mm a 2 9 20 32 vat12 26 250 the strength of this correlation r 0 98 is strong evidence in favor of the validity of this simple model these results suggest that ridge push is just sufficient to overcome drag at the base of the plate and does not contribute significantly to the motions of oceanic plates though the value of r is consistent with the motions of plates not attached to subducted slabs

the viscosity structure of the mantle during 1991 1994 the four year period covered by this report there has been a shift in emphasis away from one dimensional radial mantle viscosity models and toward an understanding of the three dimensional radial and lateral structure of mantle viscosity our understanding of mantle viscosity comes primarily from two sources studying the creep properties of mantle minerals under appropriate pressures and temperatures and using theoretical flow models to predict surface observables such as plate velocities gravity geoid and heat flow laboratory measurements of deformation indicate that the rheology of upper mantle minerals such as olivine mg fe 2sio4 is a strong function of temperatue grain size and stress cf karato and wu 1993 the deformation of minerals under mantle conditions generally follows a flow law of the form e a nd mexp q qrt 1 where is the deformation rate is the deviatoric stress is the shear modulus d is the grain size of the rock q is the activation energy for the deformation mechanism t is the temperature in kelvin r is the gas constant and a is a constant viscosity is defined as e 2 therefore deformation is directly related to viscosity substituting equation 1 into equation 2 an effective viscosity eff a n 1dmexp qrt 3 can be defined for temperature changes of 100 degrees k the viscosity changes by an order of magnitude at constant stress cf karato and wu 1993 changes of deviatoric stress by a factor of 2 change the viscosity by an order of magnitude cf karato and wu 1993 other factors such as partial pressure of oxygen and water content may also have important effects

dynamics of retreating slabs 1 insights from two dimensional numerical experiments from the effects of mantle flux part 2 therefore in this paper we apply forces to the slab using simple analytical functions related to buoyancy and viscous forces in order to isolate the role of rheology on slab dynamics we analyze parameters for simplified elastic viscous and nonlinear viscoelastoplastic single layer models of slabs and compare them with a stratified thermomechanical viscoelastoplastic slab embedded in a thermal solution the near surface behavior of slabs is summarized by assessing the amplitude and wavelength of forebulge uplift for each rheology in the complete thermomechanical solutions vastly contrasting styles of slab dynamics and force balance are observed at top and bottom bends however we find that slab subduction can be modeled using simplified rheologies characterized by a narrow range of selected benchmark parameters the best fit linear viscosity ranges between 5 10 22 pa s and 5 10 23 pa s the closeness of the numerical solution to nature can be characterized by a deborah number 0 5 indicating that elasticity is an important ingredient in subduction index terms 8120 tectonophysics dynamics of lithosphere and mantle general 8159 tectonophysics rheology crust and lithosphere 8160 tectonophysics rheology general keywords subduction numerical models lithospheric rheology

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

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

on the mechanism of seismic decoupling and back arc spreading at subduction zones to address the problem of the great variability of the mechanical state of subduction zones we investigate the mechanics of back arc spreading and seismic decoupling back arc spreading is assumed to be due to rifting of the upper plate and hence occurs when trench normal tension reaches a critical value seismic decoupling is assumed to occur when the normal stress at the frictional interface is decreased by an amount sufficient to cross the friction stability transition two forces are important in this problem the first is a small component of the slab pull force which remains unbalanced by the subduction resistance and exerts a vertical suction force at the trench the second is a sea anchor force exerted on the slab that resists its lateral motion assumed to occur at the upper plate velocity both forces contribute to the coupling problem only the sea anchor force is responsible for back arc spreading the unbalanced slab pull force is determined from a force balance for subduction the sea anchor force is computed as the hydrodynamic resistance to the facewise translation of an elliptical disc through a viscous fluid the model predicts three regimes seismically coupled compressional arcs with advancing upper plates seismically decoupled extensional arcs with retreating upper plates and strongly extensional arcs which also have back arc spreading this model is applied in detail to the izu bonin marianas system it shows that back arc spreading occurs when the integrated tension in the upper plate exceeds a value of about 1 1013 n m 1 and requires a residual tension of about a third that to drive the back arc spreading once rifting is completed it shows why the plate interface near guam is seismically coupled while the plate boundary everywhere farther north is decoupled when applied globally it successfully predicts the state of seismic coupling and back arc spreading in more than 80 of the world s subduction zones of the remaining half can be seen to contain additional complications not included in the model about 10 of cases remain unexplained but some of these may have incorrectly determined seismic coupling coefficients

retracted signatures of downgoing plate buoyancy driven subduction in motions and seismic coupling at major subduction zones this article has been retracted at the request of the editor in chief and authors please see elsevier policy on article withdrawal http www elsevier com locate withdrawalpolicy reason after publication errors were discovered in the plate motion database that it was based on this dataset was an updated version of the dataset presented in sdrolias and muller 2006 provided to us by the first author the errors in this version were in the away from trench towards trench assignment for subduction zones with back arcs and also due to the fact that the next generation plate model had only partially been completed these errors affect the conclusions about seismic coupling they also change some of the points in most of the other plots and although this does not invalidate the other conclusions the discussion to reach them would be altered

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

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

towards understanding the early gondwanan margin in southeastern australia abstract this review synthesizes the proterozoic and early paleozoic geology of tasmania bass strait and western and central victoria we examine the many different conflicting hypotheses that have been proposed to solve the paradoxical relationships between tasmanian geology and that of mainland australia most notably the prevalence of proterozoic basement of western and central tasmania while immediately across bass strait evidence of proterozoic rocks is much more cryptic we conclude that the selwyn block model is the most satisfactory hypothesis to date since it fits best with the obvious patterns in the magnetic and gravity data this model proposes that the central victorian melbourne zone is underlain by the northern extension of thin tasmanian proterozoic and cambrian crust under bass strait and that the silurian to middle devonian melbourne zone was shortened along a decollement during the tabberabberan orogeny the ordovician rocks of eastern tasmania correlate more closely with the tabberabbera zone than the melbourne zone in victoria however the silurian and devonian correlations are less certain major unresolved issues are the origins of the proterozoic and early cambrian lithostratigraphic packages tectonic models for their assembly during the tyennan orogeny and how these models fit with those for mainland australia

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

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

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

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