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

lateral injection of oxygen with the bosporus plume fingers of oxidizing potential in the black sea saline and warm mediterranean water flowing through the bosporus strait maintains a permanent pycnocline with vertical separation of oxic o 2 suboxic absence of o2 and h2s and anoxic h2s zones in the black sea the stable suboxic zone implies restricted vertical mixing of the upper oxic and lower anoxic layers and limited vertical flux of oxygen that cannot balance the upward flux of sulfide we report data that directly confirm massive lateral injections 200 km from the bosporus of oxygen enriched waters of the bosporus plume created by the mixing of shallow cold intermediate layer black sea water with mediterranean water these plume waters are laterally injected into the oxic layer and more importantly into the suboxic and anoxic layers over several small vertical scales fingers of 5 m at water densities st from 15 0 to 16 4 o2 injection oxidizes mn ii to mn iii iv which then oxidizes h2s the onset of h2s detection occurs in deeper waters in the southwest 170 m st 16 4 relative to the west central black sea 110 m st 16 2 and coincides with increased mno2 and s8 formation in the southwest

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

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

effect of three dimensional slab geometry on deformation in the mantle wedge implications for shear wave anisotropy shear wave splitting observations from many subduction zones show complex patterns of seismic anisotropy that commonly have trench parallel fast directions three dimensional flow may give rise to trench parallel stretching and provide an explanation for these patterns of seismic anisotropy along strike variations in slab geometry produce trench parallel pressure gradients and are therefore a possible mechanism for three dimensional flow in this study we quantify the effects of variable slab dip curved slabs oblique subduction and slab edges on flow geometry and finite strain in the mantle wedge of subduction zones temperature dynamic pressure velocity and strain are calculated with high resolution three dimensional finite element models these models include temperature and stress dependent rheology and parameterized slab and trench geometry thick layers 20 60 km with strong trench parallel stretching are observed in the mantle wedge when slab geometry involves a transition to slab dip less than 15 or strong curvature in the slab in these cases strong trench parallel stretching develops when flow lines have an oblique to trench normal orientation this suggests that trench parallel seismically fast directions may not indicate trench parallel flow lines in systems with large along strike variations an oblique component of stretching is confined to a 20 30 km layer above the slab in systems with oblique subduction the effects of slab edges include strong toroidal flow and focusing in the mantle near slab edges and trench parallel flow that extends 50 100 km into the core of the mantle wedge

relative motion of the nazca farallon and south american plates since late cretaceous time by combining reconstructions of the south american and african plates the african and antarctic plates the antarctic and pacific plates and the pacific and nazca plates we calculated the relative positions and history of convergence of the nazca and south american plates despite variations in convergence rates along the andes periods of rapid convergence averaging more than 100 mm a between the times of anomalies 21 49 5 ma and 18 42 ma and since anomaly 7 26 ma coincide with two phases of relatively intense tectonic activity in the peruvian andes known as the late eocene incaic and mio pliocene quechua phases the periods of relatively slow convergence 50 to 55 30 mm a at the latitude of peru and less farther south between the times of anomalies 30 31 68 5 ma and 21 and between those of anomalies 13 36 ma and 7 correlate with periods during which tectonic activity was relatively quiescent thus these reconstructions provide quantitative evidence for a correlation of the intensity of tectonic activity in the overriding plate at subduction zones with variations in the convergence rate

uplift of the western altiplano plateau evidence from the precordillera between 20 and 21 s northern chile 1 we present a quantitative reconstruction of uplift of the western flank of the altiplano plateau central andes one of the largest monoclines on the earth on the basis of an analysis of tectonic structures syntectonic deposits and geophysical data uplift occurred on a west vergent slowly propagating system of high angle reverse faults merging into a joint detachment that ramps down to midcrustal levels below the plateau edge the upper ramps determine local fold geometries while the lower ramp controls large scale surface tilting and uplift at 20 s this fault system was active between 30 ma and 5 10 ma with maximum shortening rates of 0 22 mm yr between 17 and 10 ma it generated some 2600 m of surface uplift with only minor shortening of 3000 m its activity was largely synchronous to eruption of large volume ignimbrites from a midcrustal source geophysical data indicate that the fault system localized deformation at the boundary between a cool strong forearc crust and a presumably fluid rich and or partially molten zone underneath the plateau area the systematic relation between crustal melting and shortening with uplift at the western plateau margin can be traced along most of the plateau flank with a stepwise decrease in age of deformation and magmatism toward the south indicating discontinuous addition of plateau segments crustal thickening to as much as 70 km from westward underthrusting in the back arc parts of the plateau isostatically compensated the tectonic surface uplift and monocline formation with respect to a stable forearc which only reacted with minor tilting

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

observational hints for a plume fed suboceanic asthenosphere and its role in mantle convection an asthenosphere layer which is entirely fed from below by plumes and which loses equal mass by accretion to the overlying oceanic lithosphere and at subduction zones may play a critical role in shaping the form of mantle convection in this study we discuss geochemical seismic and geoid depth evidence for lateral flow within this type of asthenosphere in particular we suggest that there are large scale layered horizontal flow structures that connect upward plume input beneath hotspots to near ridge regions of increased asthenosphere accretion into the growing oceanic lithosphere lateral asthenosphere flow is also shaped by oceanic subduction zones with a partial return flow from trenches and by deep continental roots that are migrating barriers to asthenosphere flow this alternative paradigm offers relatively simple explanations for several puzzles about mantle convection for example the low mantle heat flow beneath continents it also offers an explanation for why mid ocean ridges appear to be passive features that migrate with little geochemical or morphological change with respect to the lower mantle and seem to be uncoupled from large scale mantle flow while in contrast trenches appear to be strongly coupled to mantle thick regions of fast colder seismic velocity anomalies we also discuss several implications of this paradigm that should be testable in future studies such as the prediction of cogenetic off axis seamount volcanism that is created between an off axis hotspot and its neighboring ridge axis

geophysics of the pitman fracture zone and pacific antarctic plate motions during the cenozoic multibeam bathymetry and magnetometer data from the pitman fracture zone fz permit construction of a plate motion history for the south pacific over the past 65 million years reconstructions show that motion between the antarctic and bellingshausen plates was smaller than previously hypothesized and ended earlier at chron c27 61 million years ago the fixed hot spot hypothesis and published paleomagnetic data require additional motion elsewhere during the early tertiary either between east antarctica and west antarctica or between the north and south pacific a plate reorganization at chron c27 initiated the pitman fz and may have been responsible for the other right stepping fracture zones along the ridge an abrupt 8 clockwise rotation in the abyssal hill fabric along the pitman flowline near the young end of chron c3a 5 9 million years ago dates the major change in pacific antarctic relative motion in the late neogene

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

strong plate coupling along the nazca south america convergent margin the force balance in plate tectonics is fundamentally important but poorly known here we show that two prominent and seemingly unrelated observations trench parallel gravity anomalies along the nazca south america margin that coincide with the rupture zones of great earthquakes and a rapid slowdown of nazca south america convergence over the past 10 m y provide key insights both result from rapid miocene pliocene uplift of the andes and provide quantitative measures of the magnitude and distribution of plate coupling along the nazca south america margin we compute the plate tectonic force budget using global models of the faulted lithosphere coupled to high resolution mantle circulation models and find that andean related plate margin forces are comparable to plate driving forces from the mantle and they have sufficient magnitude to account for pronounced bathymetry variations along the trench our results suggest that plate coupling gravity anomalies and bathymetry variations along a given trench are all controlled by long term stress variations in the upper portion of plate boundaries and that an explicit budget of driving and resisting forces in plate tectonics can be obtained for the convergent margin considered here spatial variations in the effective coefficient of friction associated with the distribution of lubricating sediments entering the trench are by comparison of minor importance

was the laramide orogeny related to subduction of an oceanic plateau numerous models have been presented to explain the late cretaceous early cenozoic laramide orogeny which affected the foreland region of the western cordillera within the us1 the most attractive models invoke low angle subduction2 4 which can develop for various reasons5 6 evidence from south america indicates that one such reason may be the subduction of buoyant ocean floor we here extend the low angle subduction models by attributing the shallowing of the subduction angle during the laramide orogeny to subduction of a large oceanic plateau of anomalously thick and buoyant oceanic crust pacific plate history indicates that a twin to the hess rise may represent the oceanic plateau that triggered laramide events

a low viscosity wedge in subduction zones geochemical petrologic and seismological observations indicate that there may be high concentrations of water in the region above a subducting slab the mantle wedge which could decrease the viscosity of the mantle locally by several orders of magnitude using numerical models we demonstrate that a low viscosity wedge has a dramatic influence on the force balance in a subduction zone and leads to an observable signal in the topography gravity and geoid a regional dynamic model of the tonga kermadec subduction zone shows that the viscosity of the wedge is at least a factor of 10 smaller than surrounding mantle lithosphere and asthenosphere consistent with estimates from seismic dissipation and deformation experiments

coupling of south american and african plate motion and plate deformation although the african plate s northeastward absolute motion slowed abruptly 30 million years ago the south atlantic s spreading velocity has remained roughly constant over the past 80 million years thus requiring a simultaneous westward acceleration of the south american plate this plate velocity correlation occurs because the two plates are coupled to general mantle circulation the deceleration of the african plate due to its collision with the eurasian plate diverts mantle flow westward increasing the net basal driving torque and westward velocity of the south american plate one result of south america s higher plate velocity is the increased cordilleran activity along its western edge beginning at about 30 million years ago

age spreading rates and spreading asymmetry of the world s ocean crust we present four companion digital models of the age age uncertainty spreading rates and spreading asymmetries of the world s ocean basins as geographic and mercator grids with 2 arc min resolution the grids include data from all the major ocean basins as well as detailed reconstructions of back arc basins the age spreading rate and asymmetry at each grid node are determined by linear interpolation between adjacent seafloor isochrons in the direction of spreading ages for ocean floor between the oldest identified magnetic anomalies and continental crust are interpolated by geological estimates of the ages of passive continental margin segments the age uncertainties for grid cells coinciding with marine magnetic anomaly identifications observed or rotated to their conjugate ridge flanks are based on the difference between gridded age and observed age the uncertainties are also a function of the distance of a given grid cell to the nearest age observation and the proximity to fracture zones or other age discontinuities asymmetries in crustal accretion appear to be frequently related to asthenospheric flow from mantle plumes to spreading ridges resulting in ridge jumps toward hot spots we also use the new age grid to compute global residual basement depth grids from the difference between observed oceanic basement depth and predicted depth using three alternative age depth relationships the new set of grids helps to investigate prominent negative depth anomalies which may be alternatively related to subducted slab material descending in the mantle or to asthenospheric flow a combination of our digital grids and the associated relative and absolute plate motion model with seismic tomography and mantle convection model outputs represents a valuable set of tools to investigate geodynamic problems

bison bonasus skull from the bihor mountains romania isotopic and morphological investigations in this study we investigated morphologically and geochemically a skull from an open pit situated on the c lineasa esu g rzii plateau bihor mountains the study presents the first 14c dating and stable isotope composition of bison skeletal remains from the romanian carpathians the result of 14c determination yielded two relevant maxima one is around ad 1550 and the second around ad 1645 the data indicate that the investigated bison died during the little ice age lia approximately 200 years before the last ones were hunted to extinction stable isotope composition of carbon from bison tooth enamel suggests that the bison diet consisted mainly of c3 grasses compatible with a high altitude habitat vegetation distribution and low mobility oxygen isotope composition of both phosphate and carbonate groups of tooth enamel indicates that the bison drank from stagnant water sources such as lakes or puddles formed from rain water this further suggests low mobility of the bison as well as the presen

what drives orogeny in the andes the andes the world9s second highest orogenic belt were generated by the cenozoic tectonic shortening of the south american plate margin overriding the subducting nazca plate we use a coupled thermomechanical numerical modeling technique to identify factors controlling the intensity of the tectonic shortening from the modeling we infer that the most important factor was accelerated westward drift of the south american plate changes in the subduction rate were less important other important factors are crustal structure of the overriding plate and shear coupling at the plates9 interface the model with a thick 40 45 km at 30 ma south american crust and relatively high friction coefficient 0 05 at the nazca south american interface generates 300 km of tectonic shortening during 30 35 m y and replicates the crustal structure and evolution of the high central andes the model with an initially thinner

asymmetric mantle dynamics in the melt region of the east pacific rise abstract the mantle electromagnetic and tomography melt experiment found a surprising degree of asymmetry in the mantle beneath the fast spreading southern east pacific rise melt seismic team science 280 1998 1215 1218 forsyth et al science 280 1998 1235 1238 toomey et al science 280 1998 1224 1227 wolfe and solomon science 280 1998 1230 1232 scheirer et al science 280 1998 1221 1224 evans et al science 286 1999 752 756 pressure release melting of the upwelling mantle produces magma that migrates to the surface to form a layer of new crust at the spreading center about 6 km thick canales et al science 280 1998 1218 1221 seismic and electromagnetic measurements demonstrated that the distribution of this melt in the mantle is asymmetric forsyth et al science 280 1998 1235 1238 toomey et al science 280 1998 1224 1227 evans et al science 286 1999 752 756 at depths of several tens of kilometers melt is more abundant beneath the pacific plate to the west of the axis than beneath the nazca plate to the east melt investigators attributed the asymmetry in melt and geophysical properties to several possible factors asymmetric flow passively driven by coupling to the faster moving pacific plate interactions between the spreading center and hotspots of the south pacific an off axis center of dynamic upwelling and or anomalous melting of an embedded compositional heterogeneity melt seismic team science 280 1998 1215 1218 forsyth et al science 280 1998 1235 1238 toomey et al science 280 1998 1224 1227 wolfe and solomon science 280 1998 1230 1232 evans et al science 286 1999 752 756 here we demonstrate that passive flow driven by asymmetric plate motion alone is not a sufficient explanation of the anomalies asthenospheric flow from hotspots in the pacific superswell region back to the migrating ridge axis in conjunction with the asymmetric plate motion can create many of the observed anomalies

intraslab seismicity and thermal stress in the subducted cocos plate beneath central mexico we present a model of the subducting cocos slab beneath central mexico that provides an explanation for stresses causing the occurrence of the majority of the intraslab earthquakes which are concentrated in a long flat segment based on the recently developed thermal models for the central mexico subduction zone the thermal stresses due to non uniform temperature contrast in the subducting slab are calculated using a finite element approach the slab is considered purely elastic but due to high temperature at its bottom the behavior is considered as ductile creep the calculation results show a 20 km slab core characterized by a tensional state of stress with stresses up to 70 mpa on the other hand the top of the slab experiences high compressive thermal stresses up to 110 mpa depending on the elastic constants used and location along the flat part of the subducting plate these compressive stresses at the top of the slab are not consistent with the exclusive normal fault intraslab earthquakes and two different sources of stress are proposed the trenchward migration of the mexican volcanic arc for the last 7 ma indicates an increase of the slab dip through time this observation suggests that the gravity torque might exceed the suction torque considering the flat slab as an embedded plate subject to an applied clockwise net torque of 0 5 10 16 1 5 10 16 n m the upper half would exhibit tensional stresses of 40 110 mpa that can actually balance the compressive thermally induced stresses an alternative stress source might come from the slab pull force caused by the slab positive density anomaly based on our density anomaly estimations 75 20 kg m 3 a 350 km slab length dipping at 20 into the asthenosphere induces a slab pull force of 1 7 10 12 4 6 10 12 n m this force produces a tensional stress of 41 114 mpa sufficient to balance the compressive thermal stresses at the top of the flat slab the linear superposition of the thermally and torque or slab pull induced stresses shows tensile stresses up to 60 180 mpa inside the flat slab core also our results suggest that the majority of the intraslab earthquakes inside the flat slab are situated where the resultant stresses are larger than 40 80 mpa this study provides a reasonable explanation for the existence of exclusively normal fault intraslab earthquakes in the flat slab beneath central mexico and also it shows that thermal stresses due to non uniform reheating of subducting slabs play a considerable role in the total stress field

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

influence of continental roots and asthenosphere on plate mantle coupling the shear tractions that mantle flow exerts on the base of earth s lithosphere contribute to plate driving forces and lithospheric stresses we investigate the sensitivity of these tractions to sub lithospheric viscosity variations by comparing shear tractions computed from a mantle flow model featuring laterally varying lithosphere and asthenosphere viscosity with those from a model with layered viscosity lateral viscosity variations generally do not change the direction of shear tractions but deeply penetrating continental roots increase traction magnitudes by a factor of 2 5 compared to 100 km thick lithosphere a low viscosity asthenosphere decreases traction magnitudes by a smaller amount and is important only if 100 km thick increased plate mantle coupling beneath thick continental lithosphere may increase plate driving forces surface deformation and mantle derived lithospheric stresses in these regions by contrast a low viscosity asthenosphere does not decouple the lithosphere from mantle flow highlighting the geological importance of mantle tractions on the lithosphere

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

uncertainties and implications of the late cretaceous and tertiary position of north america relative to the farallon kula and pacific plates we present updated global plate reconstructions and calculated uncertainties of the pacific kula and farallon vancouver plates relative to north america for selected times since 68 ma improved magnetic data from the indian ocean decrease the uncertainties in the global plate circuit approach these uncertainties are now considerably smaller than those inherent in equivalent reconstructions based on the assumption of fixed hotspots major differences between these results and those of others are due to our use of more detailed africa north america reconstructions separate vancouver and farallon plate reconstructions and the assumption of a rigid antarctica plate during cenozoic time the uncertainties in the relative positions of the pacific and north america plates since the time of anomaly 7 26 ma range up to 100 km in position or from 1 to 3 m y in time if the mendocino triple junction initiated at about 28 5 ma its position would have been at 31 3 n 130 km relative to fixed north america unacceptable overlap of oceanic crust of the pacific plate with continental crust of western north america in the anomaly 10 30 ma reconstruction is a minimum of 340 200 km along an azimuth of n60 e and may be accounted for by basin and range extension pacific north america displacement in the past 20 ma is found to be considerably less than that calculated by fixed hotspot reconstructions farallon vancouver north america convergence velocity decreased greatly between the times of anomalies 24 and 21 56 to 50 ma prior to the 43 ma age of the hawaiian emperor bend and the often quoted 40 ma end of the laramide orogeny a change in direction of farallon north america convergence occurred sometime between 50 and 42 ma and also may not correlate with the time of the hawaiian emperor bend the lack of data from subducted parts of the farallon and kula plates permits many possibilities regarding the position of the kula farallon ridge the age of subducted crust or the position of oceanic plateaus during the laramide orogeny leaving open the question of the relationship between plate tectonic scenarios and tectonic style during laramide time displacements of points on the various oceanic plates along the west coast of an arbitrarily fixed north america during the interval between anomalies 30 31 and 18 68 to 42 ma are found to be pacific plate 1700 200 km northward farallon plate 3200 400 km northeastward vancouver plate 3000 400 km northeastward kula plate if attached to the pacific plate after a24 time 2500 400 km northward

cenozoic global plate motions cenozoic global plate motions relative to the hot spots are investigated and compared to plate motions in a mean lithosphere reference frame plate motions were analyzed over six time intervals divided by ages 10 25 43 48 and 56 ma chosen as much as possible to coincide with key plate reorganizations alternative motion circuits and rotational parameters were considered and evaluated with paleomagnetic data from the pacific and north american plates the circuit found to be in best agreement with the paleomagnetic data is one in which the hot spots in the atlantic region are assumed to be fixed relative to the hot spots in the pacific region throughout the cenozoic the hot spot and mean lithosphere reference frames have been in continual slow relative motion the rate of motion is nonuniform however most of the motion having occurred during the middle cenozoic the net cenozoic rotation of the lithosphere relative to the hot spots is described by a right handed rotation of 7 about a euler pole at 46 s 87 e which yields a 5 displacement of the north poles of the two reference frames this motion is small enough that inferences drawn about plate speeds in one reference frame should be valid in the other analysis of the global motions resulting from our preferred model showed that many characteristics of current plate motions have persisted throughout the cenozoic plate speeds correlate with latitude plates moving faster near the equator than near the poles throughout the cenozoic as at present continental plates except for the indian plate moved slower than oceanic plates throughout the cenozoic even the structure of the velocity fields as revealed in a contour of root mean square velocities in equatorial bands persists throughout the cenozoic the migration of the paleomagnetic axis over time is also compared to the hot spot and mean lithosphere reference frames the paleomagnetic axis has shifted 5 10 relative to the hot spot frame and a lesser amount relative to the mean lithosphere frame

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

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