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

late quaternary benthic foraminifera and deep water paleoceanography in the south china sea abstract recent and late quaternary deep sea benthic foraminifera in the south china sea scs have been quantitatively studied in order to approach some disputable questions about the glacial deep circulation of the pacific from the paleoceanographic records of a marginal sea for example whether or not there existed the north pacific deep water npdw during glacials based on the quantitative data of deep sea benthic foraminiferal faunas from sixty six surface sediment samples of the scs a q mode factor analysis indicates that the globocassidulina subglobosa assemblage the astrononion novozealandicum and bulimina aculeata assemblages and the eggerella bradyi assemblage are associated with the intermediate water mass iwm deep water mass dwm and the deep water below ccd in the scs respectively the factor assemblages resulted from surface sample data were then applied to the quantitative analysis of late quaternary deep sea benthic foraminifera in eight cores from the scs and western pacific the results show a generally increased contribution of the dwm assemblage relative to that of the iwm assemblage during glacials which implies that the influence of the dwm increased and that of the iwm decreased in the scs and western pacific however there were two major deep water paleoceanographic changes during the late quaternary 1 the deep sea carbonate dissolution of the scs and western pacific was intensified during 18 0 stages 4 5 evidenced by the higher arenaceous calcareous ratios of benthic foraminifera and relative abundance of eggerella bradyi 2 the influence of the pacific dwm above 2500 m greatly increased during the last glacial maximum lgm inferred from the increased abundances of the pacific dwm benthic foraminiferal assemblage as compared to those of the pacific iwm and pacific bottom water mass assemblages in the western pacific particularly the water in which favocassidulina favus was distributed could not enter the scs since the lgm meanwhile the great changes had taken place in the dwm of the scs during the lgm reflected in the occurrence of a clearly cooled and ventilated dwm and the replacement of a novozealandicum and b aculeata assemblages by b aculeata melonis barleeanum assemblage at water depths of 1000 2500 m these significant changes during the lgm may be related to the formation of the npdw

deformation of subducted oceanic lithosphere summary the deformation of subducted oceanic lithosphere is revealed by detailed studies of the distribution of seismicity in benioff zones and by the increasingly well resolved images of the upper mantle obtained by seismic p wave tomography by using numerical experiments based on a simplified dynamical model of a subducted lithospheric slab we test the idea that the observed deformation of the subducted slab results from a balance between internal buoyancy forces and the viscous stresses associated with deformation the simplified model assumes that the lithosphere has uniform physical properties it is denser and much more viscous than the upper mantle and its penetration below 670 km depth is resisted by a step like density increase at that level the primary determinant of the style of deformation is the buoyancy number f a ratio of buoyancy stress generated by the mass anomaly in the slab to the stress associated with viscous deformation at a strain rate u l defined in terms of subduction rate u and slab thickness l for a small buoyancy number deformation of the slab is dominated by viscous flexure for large f down dip extension of the slab appears to dominate and a buckling instability of the slab is observed if there is resistance to penetration of the 670 km level the transition value of f at which this change in behaviour is observed is about 0 05 for constant viscosity and about 0 2 for stressdependent viscosity with a stress versus strain rate exponent of n 3 with n 3 a boudinage type instability of the slab is observed for large f penetration of the slab below the 670 km level depends on the density contrast between slab and lower mantle resisted only by the density difference the slab may initially penetrate several hundred kilometres below the 670 km level for plausible density parameters before a buckling instability causes this part of the slab to rotate and ascend comparison of the deforming slab geometry and stress field with seismicity distribution and tomographic images from the tonga subduction zone suggests that the effective buoyancy number f for this slab is approximately equal to the transition value that we determined experimentally we use this constraint to estimate average rheological parameters for the tonga slab which when compared with published creep deformation laws for olivine are consistent with the average temperature of the slab being about 0 4 to 0 45 times the melting temperature based on constitutive laws for olivine at shallow depths 100 km these temperatures correspond to around 590 c for wet olivine or around 700 c for dry olivine

a new self organizing mechanism for deep focus earthquakes the mechanism of deep focus earthquakes has been a puzzle since their discovery almost 70 years ago1 2 because brittle fracture and frictional sliding at depths in excess of 100 200 km would require unrealistic rock strengths3 4 rock strength does increase with pressure but a few hundred mpa equivalent to 10 20 km depth suffices to inhibit most fracture and elevated temperature activates ductile mechanisms that operate at stresses less than the fracture stength a range of mechanisms has been proposed for deep earthquakes including plastic instabilities5 7 shear induced melting8 11 and instabilities accompanying recrystallization12 13 or polymorphic phase transformation14 23 each of these proposed mechanisms has exhibited certain inherent weaknesses see kirby22 for review and discussion here we report experimental observations of high pressure faulting of metastable mg2geo4 olivine as it undergoes incipient transformation to a spinel 24 25 structure we present a model in which this faulting arises from the generation propagation and linking up of spinel filled anticracks26 when applied to the olivine spinel transformation in the earth s mantle the anticrack model27 satisfactorily accounts for the similarities and differences between shallow and deep earthquakes and removes the problems associated with frictional sliding

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

coupled effects of phase transitions and rheology in 2 d dynamical models of subduction author s arredondo km billen mi abstract 2017 american geophysical union all rights reserved observations of seismicity and seismic tomography provide a present day constraint on the geometry of slabs within the mantle however it is challenging to directly relate the observed shape of slabs to their time evolution phase transitions modify the buoyancy forces within slabs affecting how slabs deform as they sink into the lower mantle and how forces are transmitted to the surface plates here we show that if the crustal shear zone is sufficiently weak 1020 pa s then the coupled effects of a compositionally dependent phase transition model and a stress dependent rheology lead to oscillatory plate speeds and trench motion in response to buckling and folding of the slab on average subducting plate velocity is 4 cm yr and trench motion is 0 7 cm yr however during folding subducting plate speed and trench motion increases by a factor of 3 and trench motion switches from advance to retreat in models with a higher viscosity shear zone 1021 pa s average plate speed is lower 2 cm yr and there is only trench advance stress dependent weakening due to increasing driving stresses can also cause unexpected changes in plate velocity i e decrease or no change in addition due to the added negative buoyancy from the phase transitions slab breakoff occurs in models with a yield stress 500 mpa because of the oscillatory motion of the trench long flat slabs do not form in the transition zone this suggests that other factors such as interaction of the slab with deep mantle flow may be required to create long flat slabs

folding of viscous plumes impinging on a density or viscosity interface summary the possibility that slabs of subducted oceanic lithosphere may buckle and fold under longitudinal compressive stresses when they encounter the major seismic discontinuity near a depth of 670 km in the mantle has led us to explore experimentally the behaviour of very viscous plumes falling onto a density interface or a viscosity step we begin with the case of axisymmetric and planar streams falling through air onto the free surface of a viscous fluid which has the same or smaller density and observe the conditions under which they coil or fold in order to remove surface tension and to introduce effects of a viscous environment we study plumes falling onto a density interface between two generally very viscous liquid layers stability to buckling is found to depend largely on a geometrical criterion based on the length to thickness ratio of plumes but also the tendency to buckle is greater for large density differences across the interface and larger viscosity contrasts between plume and upper layer when plumes are unstable coils or folds are laid down with a horizontal scale in a fixed proportion to the thickness of the plume and the frequency of folding is determined by the velocity and thickness of the slab as it enters the region of compressive stress entrainment of upper layer fluid can have a marked effect on the fate of the folded plume material which may spread at the interface or continue to descend into the lower layer similar behaviour is observed at a viscosity interface application of results to the mantle depends on the nature of the 670 km discontinuity and is therefore inconclusive however cool slabs are likely to be unstable if mantle circulation involves two layers separated by a density interface or if whole mantle circulation involves a large viscosity increase with depth near 700 km

grain size evolution in subducted oceanic lithosphere associated with the olivine spinel transformation and its effects on rheology we investigate the role of grain size reduction during the olivine spinel transformation on rheological properties of subducting slabs on the basis of a scaling model for microstructural development during nucleation and growth in this model the size of spinel grains nucleating at olivine grain boundaries is controlled by the relative rates of nucleation and growth taking into account the impingement through the collision of grains due to growth when the volume fraction of spinel reaches a certain threshold value critical volume fraction 1 10 depending on the p t conditions in the slab the new phase will form a continuous film and will significantly reduce the strength of the two phase aggregate if spinel grain size is small the size of spinel grains d at this stage is calculated and is shown to be highly sensitive to 0 size and 2 a subducting slab has a complicated rheological structure containing a weak region below the tip of a metastable olivine bearing wedge in a cold slab possible implications of these anomalous rheological structures of slabs on the dynamics of subduction are discussed including the mechanisms of deep earthquakes

finite amplitude convection through a phase boundary summary a finite amplitude numerical model of convective motions through a phase boundary is studied the phase change is assumed to occur over a finite vertical transition zone in the fluid and thus a continuously stratified model is considered the motions are driven by horizontal temperature gradients on the boundaries or super critical vertical stratification in the case of rayleigh instability the numerical results reproduce the neutral curve previously obtained in the analytic solution of the linearized problem when the parameters of the problem are chosen to simulate the olivine spinel transition in the upper mantle the following results are obtained 1 the vertical scale of motion is the entire depth of the fluid 2 the horizontal scale is not significantly changed from the case of a single phase fluid 3 the amplitude of the motion is not significantly changed due to buoyancy changes at the phase boundary being offset by corresponding sources and sinks of latent heat and 4 the phase boundary varies in depth by as much as 30 km when the vertical velocities are of the order of 10 1 cm yr 1

crustal shear wave velocity and radial anisotropy in the xiaojiang fault zone system se tibet revealed by ambient noise interferometry abstract the xiaojiang fault zone system xjfs is located in the southeastern tibet with high seismicity in this study we invert rayleigh and love wave dispersion curves obtained from three dense seismic arrays jointly for high resolution 3 d crustal average shear wave velocity and radial anisotropy models simultaneously in xjfs our model reveals that the upper crust and mid lower crust generally exhibit negative and positive radial anisotropy respectively implying that the deformation pattern is depth dependent to the east of the lvzhijiang fault most of the low velocity zones in the mid crust correspond to the positive radial anisotropy vsh vsv the channelized weak zone seems to be continuous across the red river fault at depth of 20 km probably within a thin layer west of the lvzhijiang fault where it is inferred to be the inner zone of the permian emeishan large igneous province the high velocity zone and positive radial anisotropy in the mid lower crust can be attributed to the underplating and intruded magmatic rocks in the upper crust the high and low velocity zones are mostly belt shape and follow the north south direction parallel to the regional major faults the lateral variation of the radial anisotropy spatially corresponds to the activity of different segments of the xiaojiang fault zone our models depict detailed geometry of the channelized weak zone in the mid lower crust and provide new insight into the role of major faults in regional tectonics

controls of the structure of subducted slabs numerical simulations of subducting slabs are formulated in which the shape and dip of the slab are determined by the dynamics of the flow rather than imposed a priori the dip of slabs is a function of the time since the initiation of subduction slabs fold develop a kink in dip and thicken on entry into a high viscosity lower mantle comparison of the simulations with seismic observations suggest that the lower mantle is at least 10 30 times more viscous than the upper mantle

kinetics of the olivine spinel transformation in subducting lithosphere experimental constraints and implications for deep slab processes abstract the persistence of metastable olivine to depths greater than 400 km in subducting slabs has implications for the generation of deep focus earthquakes the magnitude of buoyancy forces driving plate motion and the state of stress in the slab the depth to which metastable olivine can survive in a subduction zone and the depth interval over which transformation to or mg fe 2sio4 occurs have been evaluated from experimental kinetic data for the mg2geo4 and ni2sio4 transformations and the mg2sio4 transformation the data were extrapolated using a kinetic model for grain boundary nucleation and interface controlled growth under non isothermal and non isobaric conditions the results predict that metastable mg1 8fe0 2sio4 olivine survives to depths greater than 500 km in the cold interior of rapidly subducting slabs of old lithosphere the onset of transformation to mg fe 2sio4 spinel depends only on growth kinetics and coincides with the 550 50 c isotherm including the effects of latent heat production causes the transformation to occur by a runaway process over a very narrow depth interval at the onset of transformation high nucleation rates and low growth rates are consistent with the formation of very fine grained reaction products which are required for the transformational faulting mechanism of deep focus earthquakes when olivine crosses the equilibrium boundaries at higher temperatures e g higher than 700 c at 400 km depth transformation to or occurs much closer to equilibrium at a depth controlled by nucleation kinetics in this case the effect of latent heat production on the transformation kinetics is small and microstructural evolution is unlikely to result in transformational faulting below the depth of olivine breakdown cold slabs are likely to have a complex rheological structure owing to temperature dependent microstructural evolution during the phase transformation

effects of relative plate motion on the deep structure and penetration depth of slabs below the izu bonin and mariana island arcs an increasing number of seismological studies indicate that slabs of subducted lithosphere penetrate the earth s lower mantle below some island arcs but are deflected or rather laid down in the transition zone below others recent numerical simulations of mantle flow also advocate a hybrid form of mantle convection with intermittent layering we present a multi disciplinary analysis of slab morphology and mantle dynamics in which we account explicitly for the history of subduction below specific island arcs in an attempt to understand what controls lateral variations in slab morphology and penetration depth central in our discussion are the izu bonin and mariana subduction zones we argue that the differences in the tectonic evolution of these subduction zones in particular the amount and rate of trench migration can explain why the slab of subducted oceanic lithosphere seems to be at least temporarily stagnant in the earth s transition zone below the izu bonin arc but penetrates into the lower mantle below the mariana arc we briefly speculate on the applicability of our model of the temporal and spatial evolution of slab morphology to other subduction zones although further investigation is necessary our tentative model shows the potential for interpreting seismic images of slab structure by accounting for the plate tectonic history of the subduction zones involved we therefore hope that the ideas outlined here will stimulate and direct new research initiatives

effects of slab mineralogy on subduction rates although velocities of subducting slabs should be controlled primarily by their negative thermal buoyan cies their mineralogy can also have signicant eects we explore this by using thermo kinetic modeling to predict mineralogy and compare the resultant buoyant driving force to the opposing viscous drag phase transitions of mg fe 2sio4 in subducting slabs depend on thermal struc ture in two ways first equilibrium phase boundaries should be deflected causing local buoyancy anomalies whose sign depends upon that of the clapeyron slope as slabs rst enter the transition zone negative anomalies should acceler ate them but positive anomalies that arise when they fully penetrate the transition zone should slow them such ef fects may induce geologically abrupt changes in plate mo tions second olivine that persists metastably in slabs will form regions of positive buoyancy which should reduce slab velocities the coldest and fastest slabs should be slowed more greatly thus narrowing the range of feasible subduc tion rates decreased descent rates however allow slabs to warm and metastable wedges to thermally erode such neg ative feedback mechanisms may serve to regulate subduction rates

slab stagnation and detachment under northeast china abstract results of tomography models around the japanese islands show the existence of a gap between the horizontally lying stagnant slab extending under northeastern china and the fast seismic velocity anomaly in the lower mantle a simple conversion from the fast velocity anomaly to the low temperature anomaly shows a similar feature this feature appears to be inconsistent with the results of numerical simulations on the interaction between the slab and phase transitions with temperature dependent viscosity such numerical models predict a continuous slab throughout the mantle i extend previous analyses of the tomography model and model calculations to infer the origins of the gap beneath northeastern china results of numerical simulations that take the geologic history of the subduction zone into account suggest two possible origins for the gap 1 the opening of the japan sea led to a breaking off of the otherwise continuous subducting slab or 2 the western edge of the stagnant slab is the previous subducted ridge which was the plate boundary between the extinct izanagi and the pacific plates origin 2 suggesting the present horizontally lying slab has accumulated since the ridge subduction is preferable for explaining the present length of the horizontally lying slab in the upper mantle numerical models of origin 1 predict a stagnant slab in the upper mantle that is too short and a narrow or non existent gap preferred models require rather stronger flow resistance of the 660 km phase change than expected from current estimates of the phase transition property future detailed estimates of the amount of the subducted izanagi plate and the present stagnant slab would be useful to constrain models a systematic along arc variation of the slab morphology from the northeast japan to kurile arcs is also recognized and its understanding may constrain the 3d mantle flow there

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

layered convection induced by phase transitions we report a systematic study on the conditions under which an endothermic phase transition can enforce layered convection two dimensional numerical calculations of convection in a domain containing a divariant phase change were performed in the framework of the extended boussinesq approximation i e considering the effects of adiabatic gradient latent heat and frictional heating in the energy equation we find that the critical value of the negative clapeyron slope which must be surpassed in order to induce layered convection decreases in magnitude with increasing rayleigh number ra in the range 104 ra 2 106 near the critical clapeyron slope vacillations between double and single layer convection or strongly leaking double layer convection are possible the breakdown into layers is influenced very little by the latent heat release but depends solely on the phase boundary deflection caused by lateral temperature differences the value of the critical clapeyron slope also seems little affected by the width of the transition zone or by its depth a possible superplastic rheology within the transition zone would tend to favor layered convection scaling the model results to the 670 km discontinuity in the earth s mantle as a possible endothermic phase boundary we estimate the critical clapeyron slope to be in the range of 4 to 8 mpa k 40 to 80 bar k the possibility that the spinel perovskite periclase transition is within this range appears to be remote but certainly cannot be neglected

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

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

structure and evolution of lithospheric slab beneath the sunda arc indonesia tomographic imaging reveals seismic anomalies beneath the sunda island arc indonesia that suggest that the lithospheric slab penetrates to a depth of at least 1500 kilometers the sunda slab forms the eastern end of a deep anomaly associated with the past subduction of the plate underlying the mesozoic tethys ocean in accord with previous studies the lithospheric slab was imaged as a continuous feature from the surface to the lower mantle below java with a local deflection where the slab continues into the lower mantle the deep slab seems to be detached from the upper mantle slab beneath sumatra this complex slab structure is related to the tertiary evolution of southeastern asia and the indian ocean region

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

geodynamic models of deep subduction abstract numerical and laboratory models that highlight the mechanisms leading to a complex morphology of subducted lithospheric slabs in the mantle transition zone are reviewed an increase of intrinsic density with depth an increase of viscosity or phase transitions with negative clapeyron slope have an inhibiting influence on deep subduction the impingement of slabs on a viscosity and density interface has been studied in laboratory tanks using corn syrup slab interaction with equilibrium and non equilibrium phase transitions has been modelled numerically in two dimensions both the laboratory and the numerical experiments can reproduce the variety of slab behaviour that is found in tomographic images of subduction zones including cases of straight penetration into the lower mantle flattening at the 660 km discontinuity folding and thickening of slabs and sinking of slabs into the lower mantle at the endpoint of a flat lying segment aside from the material and phase transition properties the tectonic conditions play an important role in particular the retrograde motion of the point of subduction trench rollback has an influence on slab penetration into the lower mantle a question that still needs to be clarified is the mutual interaction between plate kinematics and the subduction process through the transition zone

closing the gap between regional and global travel time tomography recent global travel time tomography studies by zhou 1996 and van der hilst et al 1997 have been performed with cell parameterizations of the order of those frequently used in regional tomography studies i e with cell sizes of 1 2 these new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations 5 cells the inferred structures are however of larger scale than is usually obtained in regional models and it is not clear where and if individual cells are actually resolved this study aims at resolving lateral heterogeneity on scales as small as 0 6 in the upper mantle and 1 2 3 in the lower mantle this allows for the adequate mapping of expected small scale structures induced by for example lithosphere subduction deep mantle upwellings and mid ocean ridges there are three major contributions that allow for this advancement first we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the earth s mantle by seismic waves spakman and bijwaard 1998 second we exploit the global data set of engdahl et al 1998 which is a reprocessed version of the global data set of the international seismological centre their reprocessing included hypocenter redetermination and phase reidentification finally we combine all data used p pp and pwp phases into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal to noise ratio is improved in the approximate solution of the huge inverse problem we obtain a variance reduction of 57 1 synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells 0 6 or 1 2 in the shallow parts of subduction zones decreasing to approximately 2 3 resolution in well sampled regions in the lower mantle vertical resolution can be worse up to several hundreds of kilometers in subduction zones with rays predominantly pointing along dip important features of the solution are as follows 100 200 km thick high velocity slabs beneath all major subduction zones sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle large high velocity anomalies in the lower mantle that have been attributed to subduction of the tethys ocean and the farallon plate and low velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under iceland east africa the canary islands yellowstone and the society islands our findings corroborate that the 660 km boundary may resist but not prevent present day large scale mass transfer from upper to lower mantle or vice versa this observation confirms the results of previous global mantle studies that employed coarser parameterizations

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

elasticity composition and temperature of the earth s lower mantle a reappraisal the interpretation of seismological models for the earth s lower mantle in terms of chemical composition and temperature is sensitive both to the details of the methodology employed and also to uncertainties in key thermoelastic parameters especially for the dominant mg fe sio3 perovskite phase here the alternative approaches adiabatic decompression of the lower mantle for comparison at zero pressure with laboratory data and the projection of laboratory data to lower mantle p t conditions for direct comparison with seismological observations are assessed along with the equations of state on which they are based it is argued that adiabatic decompression of the lower mantle is best effected by a strategy in which consistent third order finite strain expressions are required simultaneously to fit the strain dependence of both the seismic parameter and the density this procedure accords due weight to the most robust seismological observations namely the wave speeds and reduces to an acceptable level the otherwise very strong covariance among the fitted coefficients it is demonstrated that this approach can readily be adapted to include the effects of relaxation from the hill average to the reuss lower bound of the aggregate bulk modulus which governs the radial variation of density for the projection of laboratory thermoelastic data to lower mantle p t conditions the preferred equation of state is of the mie gruneisen type involving the addition at constant volume of the pressure along a finite strain principal 300 k isotherm and the thermal pressure calculated from the debye approximation to the lattice vibrational energy a high degree of consistency is demonstrated between these alternative equations of state possible departures of the lower mantle from conditions of adiabaticity and large scale homogeneity are assessed these are negligible for the prem model bullen parameter b 0 99 0 01 but significant for ak135 b 0 94 0 02 reviving the possibility of a substantially superadiabatic temperature gradient experimentally determined thermoelastic properties of the major mg fe sio3 perovskite and mg fe o magnesiowustite phases are used to constrain the equation of state parameters employed in the analysis of seismological information although the scarcity of information concerning the pressure and temperature dependence of the elastic moduli for the perovskite phase remains a serious impediment nevertheless it is demonstrated that the combination of a pyrolite composition simplified to the three component system sio2 mgo feo with molar xpv 0 67 and xmg 0 89 and a lower mantle adiabat with a potential temperature of 1600 k consistent with the preferred geotherm for the upper mantle and transition zone is compatible with the prem seismological model within the residual uncertainties of the thermoelastic parameters for the perovskite phase in particular such consistency requires values for the perovskite phase of ks ks p t and q ln ln v t of approximately 3 8 and 2 and g t near 0 022 gpa k 1 for the lower mantle assemblage more silicic models which have sometimes been advocated can also be reconciled with the seismological data but require lower mantle temperatures which are much higher by about 700 k for pyroxene stoichiometry however the absence of seismological and rheological evidence for the pair of thermal boundary layers separating convection above and below which is implied by such models remains a formidable difficulty under these circumstances the simplest possible model that of grossly uniform chemical composition throughout the mantle is preferred

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

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

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

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

compositional vs thermal buoyancy and the evolution of subducted lithosphere we formulate 2 d cartesian finite element models that explore the fate of compositionally defined lithosphere as it encounters a viscosity increase at the boundary between the upper and lower mantle subducted lithosphere is represented as a cold stiff layered composite of denser eclogite underlain by more buoyant harzburgite slabs impinging on a lower mantle 30 and 100 times more viscous than the upper mantle thicken and fold strongly as they penetrate the lower mantle approximately a factor of two thickening occurs via pure shear just above the discontinuity with additional enhancement due to folding by over a factor of two no separation of the individual slab components occurs at the discontinuity and direct comparison with models in which compositional buoyancy is explicitly ignored indicates that slab evolution is largely controlled by the thermal buoyancy these results are at odds with hypotheses about slab evolution in which the compositional buoyancy contributions lead to component separation and the formation of slab megaliths or a compositionally layered upper mantle

cenozoic plate driving forces past studies of plate driving forces have concluded that the forces due to subducted slabs in the upper mantle and those due to the thickening of the oceanic lithosphere are the principal driving forces we reexamine the balance of driving forces for the present day and extend our analysis through the cenozoic using an analytical torque balance method which accounts for interactions between plates via viscous coupling to the induced mantle flow we use an evolving mantle density heterogeneity field based on the last 200 myr of subduction to drive plate motions an approach which has proven successful in predicting the present day mantle heterogeneity field we find that for plausible upper mantle viscosities the forces due to subducted slabs in the cenozoic and mesozoic account for in excess of 90 of plate driving forces and those due to lithospheric thickening for less than 10

effects of mantle and subduction interface rheologies on slab stagnation and trench rollback abstract trench rollback has been a widely discussed phenomenon in recent years and multiple studies have concentrated on various parameters that may influence trench migration and related aspects of slab deformation in the upper mantle here we concentrate on the effects of rheological description yield stress lower mantle viscosity viscosity of crust in controlling the rollback and associated stagnation of slabs in the transition zone 410 660 km depth we perform numerical simulations of slab evolution in a 2d cartesian model with strongly nonlinear rheology combining diffusion creep dislocation creep and a power law stress limiter we demonstrate that trench retreat develops in most models considered regardless of the subducting plate age or prescribed strength rollback then mostly produces slabs that are horizontally deflected at the 660 km phase boundary and remain subhorizontal at the bottom of the transition zone slab morphologies are in agreement with stagnant horizontally deflected structures reported in the transition zone by seismic tomography furthermore if the strength of the slab is limited to less than 0 5 gpa the slab experiences a significant amount of horizontal buckling the amplitude of the rollback velocity is sensitive to several model parameters as one might expect it increases with the age of the subducting plate thus reflecting its increasingly negative buoyancy on the other hand rollback velocity decreases if we increase the viscosity of the crust and thereby strengthen the coupling between the subducting and overriding plates high friction on the contact between the subducting and overriding plates may even result in slabs penetrating into the lower mantle after a period of temporary stagnation also reducing the additional negative buoyancy imparted by the 410 km exothermic phase transition suppresses trench rollback interpretation of the controls on slab rollback and stagnation may be rather complex in strongly nonlinear rheological models where for example buoyancy effects may be counteracted by associated yield stress weakening

intermittent layered convection in a model mantle with an endothermic phase change at 670 km it is not yet known whether convection in the earth s mantle occurs in cells that extend throughout the entire depth of the mantle or if the upper and lower mantle convect in independent layers both possibilities seem to be consistent with geochemical and seismologi cal observations1 4 it is generally accepted however that the seismologically observed boundary between the upper and lower mantle at 670 km depth arises from the endothermic phase change spinel perovskite mgo here we investigate the results of including this phase change in a model of mantle convection for realism and to eliminate scaling problems our calculations use commonly accepted values of geophysical parameters for a clapeyron slope y 2 106 a value close to experimental results and theoretical expectations we observe local intermittent mixing between the upper and lower mantle such behaviour may offer a way to reconcile the existence of geophysical evidence for both whole mantle and layered convection

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

metastable mantle phase transformations and deep earthquakes in subducting oceanic lithosphere earth s deepest earthquakes occur as a population in subducting or previously subducted lithosphere at depths ranging from about 325 to 690 km this depth interval closely brackets the mantle transition zone characterized by rapid seismic velocity increases resulting from the transformation of upper mantle minerals to higher pressure phases deep earthquakes thus provide the primary direct evidence for subduction of the lithosphere to these depths and allow us to investigate the deep thermal thermodynamic and mechanical ferment inside slabs numerical simulations of reaction rates show that the olivine spinel transformation should be kinetically hindered in old cold slabs descending into the transition zone thus wedge shaped zones of metastable peridotite probably persist to depths of more than 600 km laboratory deformation experiments on some metastable minerals display a shear instability called transformational faulting this instability involves sudden failure by localized superplasticity in thin shear zones where the metastable host mineral transforms to a denser finer grained phase hence in cold slabs such faulting is expected for the polymorphic reactions in which olivine transforms to the spinel structure and clinoenstatite transforms to ilmenite it is thus natural to hypothesize that deep earthquakes result from transformational faulting in metastable peridotite wedges within cold slabs this consideration of the mineralogical states of slabs augments the traditional largely thermal view of slab processes and explains some previously enigmatic slab features it explains why deep seismicity occurs only in the approximate depth range of the mantle transition zone where minerals in downgoing slabs should transform to spinel and ilmenite structures the onset of deep shocks at about 325 km is consistent with the onset of metastability near the equilibrium phase boundary in the slab even if a slab penetrates into the lower mantle earthquakes should cease at depths near 700 km because the seismogenic phase transformations in the slab are completed or can no longer occur substantial metastability is expected only in old cold slabs consistent with the observed restriction of deep earthquakes to those settings earthquakes should be restricted to the cold cores of slabs as in any model in which the seismicity is temperature controlled via the distribution of metastability however the geometries of recent large deep earthquakes pose a challenge for any such models transformational faulting may give insight into why deep shocks lack appreciable aftershocks and why their source characteristics including focal mechanisms indicating localized shear failure rather than implosive deformation are so similar to those of shallow earthquakes finally metastable phase changes in slabs would produce an internal source of stress in addition to those due to the weight of the sinking slab such internal stresses may explain the occurrence of earthquakes in portions of lithosphere which have foundered to the bottom of the transition zone and or are detached from subducting slabs metastability in downgoing slabs could have considerable geodynamic significance metastable wedges would reduce the negative buoyancy of slabs decrease the driving force for subduction and influence the state of stress in slabs heat released by metastable phase changes would raise temperatures within slabs and facilitate the transformation of spinel to the lower mantle mineral assemblage causing slabs to equilibrate more rapidly with the ambient mantle and thus contribute to the cessation of deep seismicity because wedge formation should occur only for fast subducting slabs it may act as a parachute and contribute to regulating plate speeds wedge formation would also have consequences for mantle evolution because the density of a slab stagnated near the bottom of the transition zone would increase as it heats up and the wedge transforms to denser spinel favoring the subsequent sinking of the slab into the lower mantle