the analogue shear zone from rheology to associated geometry the geometry of ductile strain localization phenomena is related to the rheology of the deformed rocks both qualitative and quantitative rheological properties of natural rocks have been estimated from finite field structures such as folds and shear zones we apply physical modelling to investigate the relationship between rheology and the temporal evolution of the width and transversal strain distribution in shear zones both of which have been used previously as rheological proxies geologically relevant materials with well characterized rheological properties newtonian strain hardening strain softening mohr coulomb are deformed in a shear box and observed with particle imaging velocimetry piv it is shown that the width and strain distribution histories in model shear zones display characteristic finite responses related to material properties as predicted by previous studies application of the results to natural shear zones in the field is discussed an investigation of the impact of 3d boundary conditions in the experiments demonstrates that quantitative methods for estimating rheology from finite natural structures must take these into account carefully

strain localisation and weakening of the lithosphere during extension we explore the sensitivity of extensional systems to the thermal structure of the lithosphere using numerical simulations that fully couple the energy momentum and continuity equations the rheology of the lithosphere is controlled by weakening processes such as shear heating that localises strain into shear zones and faults numerical models show that during extension of an initially unpatterned lithosphere structures develop spontaneously out of basic thermodynamic energy fluxes and without the imposition of ad hoc rules on strain localisation this contrasts with the classical mohr coulomb theory for brittle localisation which prescribes the angles of faults by a mathematical rule our results show that the mode of extension is sensitive to subtle changes in rheology heat flux and geometry of the system this sensitivity lies at the core of the variety and complexity observed in extensional systems localisation processes make the lithosphere weaker than previously estimated from the brace goetze quasi static approach consequently typical estimates for plate tectonic forces are capable of splitting the lithosphere under extension even without the role of active magmatism

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

complex patterns of fluid and melt transport in the central andean subduction zone revealed by attenuation tomography abstract we present a high resolution 3 d model of p wave attenuation q p 1 for the central andean subduction zone data from 1500 mostly intermediate depth 60 250 km earthquakes recorded at three temporary seismic networks covering the forearc arc and backarc around 23 s were used for tomographic inversion the forearc is characterised by uniformly high q p values indicating low temperature rocks in accordance with low surface heat flow values prominent low q p anomalies are found beneath the magmatic arc and the backarc in the crust and mantle continuous regions of low q p connect earthquake clusters at 100 km and 200 km depth with zones of active volcanism in the arc and backarc fluids fluxed from the subducted oceanic lithosphere into the overlying mantle wedge where they induce melting explain our observations we propose that low q p regions indicate source and ascent pathways of metamorphic fluids and partial melts ascent of fluids and melts as imaged by seismic q p are not vertical as is often implicitely assumed instead sources of fluids are located at different depth levels and ascent paths are complex and exhibit significant variation within the study area the largest quaternary backarc volcano cerro tuzgle is fed by mantle melts which are imaged as a plume of low q p material that reaches to the strong earthquake cluster at 200 km depth

scattering of near inertial internal waves along the japanese coast of the japan sea 1 the scattering of near inertial internal waves niiws induced at the tips of landmasses was investigated by mooring observations using current meters and thermometers in the surface layer along peninsulas and an island located along the japanese coast of the japan sea during the fall of 2006 signals of niiws were found in alongshore currents and in temperatures recorded near the seasonal thermocline the signals appeared after wind stress had reached a maximum and they lasted for more than 2 days relatively high current energy regions due to niiws were localized around the northern parts of the peninsulas and island the phase lag of the temperature fluctuations increased in the southern regions along the eastern coasts of the peninsulas and island the generation and propagation of niiws were investigated by numerical experiments that used a multilayer level model with realistic topography model results showed that niiws induced by the wind struck the tips of landmasses and scattered and that niiws were trapped by the coast and propagated with characteristics of coastal trapped waves ctws and then rapidly attenuated the ctws had a distorted current structure at landmass tips because of periodic beam like current structures from the seafloor these signals show that ctws radiate energy away from the shore as internal waves originating from topographical scattering which is considered to cause the rapid attenuation of the ctws

fore arc deformation at the transition between collision and subduction insights from 3 d thermomechanical laboratory experiments 1 three dimensional thermomechanical laboratory experiments of arc continent collision investigate the deformation of the fore arc at the transition between collision and subduction the deformation of the plates in the collision area propagates into the subduction collision transition zone via along strike coupling of the neighboring segments of the plate boundary in our experiments the largest along strike gradient of trench perpendicular compression does not produce sufficiently localized shear strain in the transition zone to form a strike slip system because of the fast propagation of arc lithosphere failure deformation is continuous along strike but the deformation mechanism is three dimensional progressive along strike structural variations arise because coupling between neighboring segments induces either advanced or delayed failure of the arc lithosphere and passive margin the modeling results suggest that orogenic belts should experience deeper subduction of continental crust and hence higher pressure metamorphism where the two plates first collided than elsewhere along the plate boundary where collision subsequently propagated furthermore during the initial stage of collision the accretionary wedge is partially subducted which leads to lubrication of the interplate zone and a reduction of shear traction therefore a large convergence obliquity angle does not produce a migrating fore arc sliver rather the pressure force generated by subduction of the buoyant continental crust causes fore arc motion it follows that convergence obliquity during collision does not yield trench parallel deformation of the fore arc and its influence on the collision process is limited however convergence obliquity may control the geometry of the active margin during the oceanic subduction stage prior to collision and inherited structures may influence the propagation mechanism

trench parallel shortening in the forearc caused by subduction along a seaward concave plate boundary insights from analogue modelling experiments abstract three dimensional thermo mechanical analogue experiments are employed to test the hypothesis that oceanic subduction along a seaward concave plate boundary can generate trench parallel shortening in the forearc near the axis of curvature the model deformation is analyzed with a particle imaging velocimetry piv system that allows for comparison of forearc deformation along the oblique limbs of the curved plate boundary and near the axis of curvature moreover piv allows for separation of the trench parallel and trench perpendicular components of strain regardless of trench orientation the resulting deformation maps show a remarkable symmetry and indicate drag of the forearc above the interplate coupling area towards the axis of curvature trench perpendicular profiles show that along the oblique limbs of the plate boundary the forearc is submitted to trench normal shortening and trench parallel shearing but not trench parallel shortening or extension this contrasts with the situation near the axis of symmetry where the forearc is submitted to trench parallel and trench perpendicular normal shortening but is not sheared the experimental results confirm that trench normal thrusts observed in the fore arc of the central andes can be a mechanical consequence of subduction along a seaward concave plate boundary if the degree of interplate coupling is large

frontal accretion an internal clock for bivergent wedge deformation and surface uplift 1 we propose a conceptual model of frontal accretion within bivergent wedges which is based on two dimensional sandbox simulations and the analysis of particle displacement fields each frontal accretion cycle consists of a thrust initiation an underthrusting phase and a reactivation phase the location and magnitude of deformation within a bivergent wedge and its associated surface uplift vary systematically with the phase of the frontal accretion cycle and are thus predictable therefore the frontal accretion cycle can be considered as an internal clock for wedge scaled deformation and surface uplift we further demonstrate that the geometry of the deformation front and the spatial distribution of surface uplift can be used to infer the currently active phase within a frontal accretion cycle surface uplift of the axial zone and the retrowedge may reach up to half of the thickness of the incoming sedimentary layer during a frontal accretion cycle we estimate cycle duration to range between 104 and 105 years which is thus similar to the time frame of climatic cycles the proposed conceptual model provides an explanation for the commonly observed transience of deformation and surface uplift and for the discrepancy between geodetic paleoseismologic and geologic estimates of fault slip

rotation of the peruvian block from palaeomagnetic studies of the central andes palaeomagnetic studies were performed on the rocks of the central andes to test the hypothesis of oroclinal bending of the andes around the peru chile border natural remanent magnetization of mesozoic volcanic and sedimentary rocks of peru shows a counterclockwise shift of declination by several tens of degrees relative to the field directions of the same ages in the stable south american craton two mesozoic dyke swarms from arica region northernmost chile also indicate about 20 of declination shift in a counterclockwise sense these results suggest a considerable amount of counterclockwise rotation of large area north of arica that is oroclinal bending of the andes around some hinge near the peru chile border one neogene dyke swarm sampled in central peru also gave declination which deviates counterclockwisely by 15 if this declination shift is a result of the wide area tectonic rotation as suggested by the mesozoic rocks oroclinal bending appears to have continued until relatively recent times

rupture process of the 1964 prince william sound alaska earthquake from the combined inversion of seismic tsunami and geodetic data 1 four great earthquakes 1952 1960 1964 and 2004 have occurred since seismic monitoring began and only two since the installation of a global seismic network a reexamination of the 1964 m 9 2 prince william sound pws alaska earthquake is timely due to the 2004 sumatra earthquake because it adds constraints to the potential range of source parameters for these types of infrequent events and aids in the ability to predict the ground motions for other subduction zone including cascadia we first measured the durations of high frequency energy radiation from teleseismic p wave codas recorded by short period world wide standardized seismographic network wwssn sp instruments to put constraints on the earthquake rupture length the durations ranged between 250 and 600 s and are shorter in time between azimuths of 190 and 260 the fault length is estimated to range from 540 to 740 km rupturing at average speeds of 1 4 2 km s in a 220 238 direction we developed a multiple time window kinematic rupture model for the pws earthquake from the least squares inversion of teleseismic p waves tsunami tide gauge records and geodetic leveling survey observations based on the green s function technique we assume three major fault segments for the subduction zone dipping 6 12 based on geologic data the subfault size on the megathrust was set to 50 50 km the patton bay fault pbf imbricate thrust was assigned a 60 dip with a subfault size of 20 20 km we estimated a seismic moment of 5 52 1022 n m mw9 12 we identified three areas of major moment release on the pws segment near montague and middleton islands extending out to the trench beneath the portlock anticline on the cook segment between 58 n fracture zone and the kodiak bowie seamount chain and finally on the katmai segment off the coast of kodiak island extending to the trench our preferred rupture model has a peak slip of 14 9 m along the megathrust and 17 4 m along the pbf the addition of the pbf resulted in an 8 improvement in residuals and also had a significant effect on the overall slip distribution moment release occurred as deep as 50 km depth slab contour near the 350 c isotherm and mantle fore arc wedge which is about 10 20 km deeper than previous coupling depths

on the strength of interplate coupling and the rate of back arc convergence in the central andes an analysis of the interseismic velocity field 1 we interpret the interseismic crustal velocity field of the central andes using a simple three plate model in which the andean mountain belt is treated as a rigid microplate located between the nazca and south american soam plates we assume that the euler vectors associated with these plates are strictly coaxial and that the surface velocity field can be decomposed into an elastic loading field driven by locking of the main plate boundary and a smaller contribution associated with back arc convergence we obtain our best fit to the geodetic velocities if the main plate boundary is fully 100 locked between depths of 10 50 km and 8 5 of nazca soam plate convergence is achieved in the back arc by underthrusting of the brazilian shield beneath the subandean zone

horizontal lithosphere compression and subduction constraints provided by physical modeling physical modeling of the subduction is performed with a two layered mantle model elastico plastic lithosphere and low viscosity asthenosphere and is governed by the criteria of similarity compression of the lithosphere in the area of a passive continental margin has been shown to produce a buckling instability in the oceanic plate with wavelengths of 200 km on the average later a localization of deformation occurs in sagging at some distance from the margin where a strongly dislocated linear ridge is formed due to the thrusting the plate then experiences a failure along the inclined zone and subduction starts the inner trench slope which forms therefore has a scraped structure and include a block of crushed and dislocated oceanic crust and sediments in the lower section if there is an old inclined fault striking across the compression of an oceanic lithosphere it is on that fault that a subduction zone is initiated the inner trench slope has then a different structure and forms due to normal faulting in the frontal part of the overriding plate the formation of a subduction zone requires a compression that is smaller than that in the preceding case by a factor of 2 or 3 the subducting plate experiences an elastico plastic bending and under specific conditions thrusting along the zone dipping from under the overriding plate oceanward and crossing the entire lithosphere the best agreement between generalized relief in subduction zone in the model and nature is achieved when a shear yield limit s 1 3 108 pa modulus of elasticity e about a few times 1011 pa and a thickness h 60 km are adopted for the real lithosphere

lithospheric domains and controls on kimberlite emplacement slave province canada evidence from elastic thickness and upper mantle composition 1 we have mapped the deep structure of the slave craton by combining analysis the effective elastic thickness te with data on mantle samples from numerous kimberlites three dimensional mapping of the subcontinental lithospheric mantle sclm using mantle derived xenoliths and xenocrysts in kimberlites has shown that much of the craton is underlain by a strongly layered sclm a highly depleted upper layer low in basaltic components ca al fe is separated from a relatively fertile lower layer by a sharp boundary this boundary lies at 140 150 km depth in the lac de gras area and shallows to 100 km in the northern and southern parts of the craton weak lithosphere te 56 km in the younger eastern part of the craton is separated from the older western part by a zone of steep te gradient parallel to the major locus of kimberlite intrusion which may map the deep extension of the boundary between the two domains another strong te gradient across the kilohigok basin accompanies a marked compositional change in the upper layer of the lithospheric mantle the basin probably marks a major translithospheric fault correlations between te and mantle composition suggest that te is strongly influenced by the rheology of the upper mantle

topography and subduction geometry in the central andes clues to the mechanics of a noncollisional orogen the central andean orogen between 12 and 32 s latitude exhibits a high degree of spatial order principally an extraordinary bilateral symmetry that is common to the earth s surface the underlying wadati benioff zone and the nazca south america plate kinematics which has been stable since the mid tertiary this spatial order must reflect the physical mechanisms of mountain building in this noncollisional orogen the shapes of the topography and subduction zone can be reduced to symmetric and antisymmeric components relative to any vertical symmetry plane the particular plane which minimizes the antisymmetry and maximizes the symmetry is well resolved and is essentially coincident with the stable euler equator of nazca south america relative motion since the mid tertiary that the topography subduction geometry and persistent mid tertiary plate kinematics share common spatial and geometric elements suggests that the distribution of topography in this orogen depends strongly on the dynamics of subduction other factors that might affect the topography and underlying tectonics such as climate and inherited structural fabric which have different spatial characteristics must be of less significance at a continental scale furthermore the small components of asymmetry among the various elements of the orogen appear to be mutually related in a simple way it is possible that this coupled asymmetry is associated with a late tertiary change in plate kinematics these observations suggest that there is a close connection between plate tectonics and the form of the earth s surface in this noncollisional setting it follows that the distribution of topography near convergent plate boundaries may provide a powerful constraint for understanding the dynamics of subduction

bending the bolivian orocline in real time global positioning system gps data from the central andes record vertical axis rotations that are consistently counterclockwise in peru and bolivia north of the bend in the mountain belt and clockwise to the south in southern bolivia argentina and chile these geologically instantaneous rotations have the same sense as rotations that have accrued over millions of years and are recorded by paleomagnetic and geologic indicators the change in sign of the rotation at both decadal and million year time scales occurs across the axis of topographic symmetry that defines the bolivian orocline when extrapolated to a common time interval the magnitudes of rotation from geologic features and from gps are surprisingly similar given that a significant part of the instantaneous deformation field is probably elastic and due to interseismic locking of the plate boundary some of the interseismic deformation field must reflect permanent deformation and or some of the current elastic deformation will be converted to upper plate permanent deformation over time rather than be recovered by elastic rebound during interplate earthquakes we suggest that the spatial patterns of the elastic and the permanent modes of bending are similar because they are driven by the same stress field

role of the plate margin curvature in the plateau buildup consequences for the central andes 1 the influence of convergent plate boundary curvature on the stress distribution in an overriding plate is explored using analytical and numerical modeling techniques trench parallel compression can be produced near the symmetry axis of a seaward concave plate boundary if the interplate friction is high and or if the subducting lithosphere has a low flexural rigidity which produces little nonhydrostatic normal stress on the plate boundary this situation favors the formation of trench parallel thrusts with minor trench parallel component of slip trench parallel compression is reduced along the most oblique parts of the plate boundary which favors the formation of strike slip faults with major trench parallel slip both the stress conditions on the interplate zone and the 3 d geometry of this zone control whether the trench parallel stress in the center of a seaward concave curvature is a tension or compression low dip angle and high convergence obliquity angle favor trench parallel compression in the central andes n s minor shortening in the center of the arica bend and strike slip systems north and south of the symmetry axis suggest that the effect of shear traction dominated during cenozoic time when the curvature of the plate boundary was forming this result suggests that the processes responsible for the formation of the plate boundary curvature were assisted by enhanced interplate friction and or reduced compressive nonhydrostatic normal stress for a geometry resembling the present day south american plate margin estimations of normal and shear stresses on the plate boundary suggest that the trench parallel stress in the center of the curvature is compressive

narrow subducting slabs and the origin of backarc basins abstract both tectonic and seismic data suggest that horizontally wide in the direction parallel to the axis of the trench dipping subducting slabs typically move stably at fixed angles to the horizontal and that their associated trenches migrate very slowly or not at all in contrast where backarc basins babs exist the associated trenches typically migrate at rates of several cm year this behavior may be related to the finite width of the slabs under babs analysis shows that hydrodynamic suction created by the downdip motion of a wide slab is sufficient to balance the pull of gravity on the heavy slab however for a narrow slab sideways asthenospheric flow into the slab overriding plate corner effectively eliminates the suction which results in the rapid sinking of the slab hydrodynamic modeling is used to quantify this hypothesis modeling is used also to show that wide slabs are stable in their inclined position not because of the gravitational wing like sliding but due to the lifting force produced by the negative pressure suction in the slab overriding plate corner finally the model results confirm that for relatively narrow slabs this lifting force is small and cannot counterbalance gravity therefore when bending moments in a narrow slab exceed the plastic yield limit the slab becomes unstable this instability can cause the rapid migration of the slab s hinge and the associated arc at rates of 10 20 cm year and thus form the bab

particle fallout thermal disequilibrium and volcanic plumes we develop a model of a volcanic eruption column that includes the effects of fallout of pyroclasts and thermal disequilibrium we show that clast fallout with no thermal disequilibrium has only a small effect upon the column however disequilibrium changes column behaviour significantly and can even induce collapse our results may explain the lower plume heights and less widely dispersed fallout of cone forming eruptions contrasted with sheet forming eruptions the model also predicts that the transition from the gas thrust to the convective region in a column results in an inflection in dispersal curves that some features of the stratigraphy common to many fall deposits may result from column velocity structure and that there may exist a region near the volcanic vent in which maximum pyroclast size does not decrease significantly with distance

continental subduction and exhumation of high pressure rocks insights from thermo mechanical laboratory modelling abstract thermo mechanical physical modelling of continental subduction is performed to investigate the exhumation of deeply subducted continental crust the model consists of two lithospheric plates made of new temperature sensitive analogue materials the lithosphere is underlain by liquid asthenosphere the continental lithosphere contains three layers the weak sedimentary layer the crust made of a stronger material and of a still stronger lithospheric mantle the whole model is subjected to a constant vertical thermal gradient causing the strength reduction with depth in each lithospheric layer subduction is driven by both push force and pull force during subduction the subducting lithosphere is heating and the strength of its layers reduces the weakening continental crust reaches maximal depth of about 120 km and cannot subduct deeper because its frontal part starts to flow up the subducted crust undergoes complex deformation including indicated upward ductile flow of the most deeply subducted portions and localised failure of the subducted upper crust at about 50 km depth this failure results in the formation of the first crustal slice which rises up between the plates under the buoyancy force this process is accompanied by the delamination of the crustal and mantle layers of the subducting lithosphere the delamination front propagates upwards into the interplate zone resulting in the formation of two other crustal slices that also rise up between the plates average equivalent exhumation rate of the crustal material during delamination is about 1 cm year the crust asthenosphere boundary near the interplate zone is uplifted the subducted mantle layer then breaks off removing the pull force and thereby stopping the delamination and increasing horizontal compression of the lithosphere the latter produces shortening of the formed orogen and the growth of relief the modelling reveals an interesting burial exhumation evolution of the sedimentary cover during initial stages of continental subduction the sediments of the continental margin are dragged to the overriding plate base and are partially accreted at the deep part of the interplate zone at 60 70 km depth these sediments remain there until the beginning of delamination during which the pressure between the subducted crust and the overriding plate increases this results in squeezing the underplated sediments out part of them is extruded upwards along the interplate zone to about 30 km depth at an equivalent rate of 5 10 cm year

the temperature and salinity fine structure of the ocean under the ross ice shelf a series of temperature and salinity profiles was made in the ocean under the ross ice shelf at 82 22 5 s 168 37 5 w where the ice was 420 and the underlying seawater 240 m thick the water structure consisted of a fairly well mixed low salinity layer at the in situ freezing point of the ice water interface about 30 m thick a transition layer characterized by intrusions about 85 m thick a strongly stratified layer with increasing temperature and salinity about 50 m thick another transition layer characterized by intrusions about 45 m thick and an isohaline bottom layer with a nearly abiabatic temperature gradient about 30 m thick the temperature fluctuations in the two transition layers can be attributed mainly to intrusive activity even though internal wave activity was highest in these regions in the central stable layer internal waves probably contribute to the temperature fluctuations equally with intrusions the internal wave dropped displacement spectrum seems to be very similar to those found in the open ocean some profiles showed intrusions in the top transition layer colder than the top boundary layer and some showed intrusions in the bottom transition layer warmer than the bottom boundary layer indicating that horizontal shear must be present internal waves of near inertial frequency excited by the semidiurnal tides combined with shear induced instabilities seem to be the likely causes of the observed fine structure

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

evolutionary model of the himalaya tibet system geopoem based on new modelling geological and geophysical data a two dimensional thermo mechanical laboratory modelling of continental subduction was performed the subducting continental lithosphere includes a strong brittle upper crust a weak ductile lower crust and a strong upper mantle the lithosphere is underlain by a low viscosity asthenosphere subduction is produced by a piston push force and the pull force from the mantle lithospheric layer which is denser than the asthenosphere the lithospheric layers are composed of material whose strength is sensitive to and inversely proportional to temperature throughout the experiment the model surface was maintained under relatively low temperature and the model base at higher temperature the subduction rate satisfied the peclet criterion modelling confirms that the continental crust can be deeply subducted and shows that slab break off delamination and tectonic underplating are fundamental events with drastic consequences on the subsequent evolution of the convergent system combining these results with previous purely mechanical modelling we elaborate a new evolutionary model for the himalaya tibet convergent system the principal successive stages are 1 subduction of the indian continental lithosphere to 200 250 km depth following subduction of the tethys oceanic lithosphere 2 failure and rapid buoyancy driven uplift of the subducted continental crust from ca 100 km depth to some depth that varies along the mountain belt 20 30 km on average 3 break off of the indian subducted lithospheric mantle with the attached oceanic lithosphere 4 subduction underplating of the indian lithosphere under asia over a few to several hundred kilometers 5 delamination roll back and break off of the indian lithospheric mantle 6 failure of the indian crust in front of the mountain belt formation of the main central thrust and underthrusting of the next portion of indian lithosphere beneath tibet for a few hundred kilometers at the beginning of stage 6 the crustal slice corresponding to the crystalline himalayas undergoes erosion activated uplift and exhumation fl 2000 elsevier science b v all rights reserved

physical modelling of lithosphere subduction in collision zones abstract physical mechanical modeling of continental collision is performed with a 4 layered mantle model compositionally distinct layers of hydrocarbons are used to model a lithosphere which consists of two crustal layers the upper one quasi brittle the lower one plastic and a plastic mantle layer the lower crustal layer the crustal asthenosphere has a low yield limit a low viscosity liquid water represents the asthenosphere under the lithosphere lateral compression of the model plate by a piston is the driving mechanism of collision similarity criteria are satisfied in the modeling the collisional process was investigated with various boundary conditions and values of the model parameters the yield limit and the thickness of the crustal asthenosphere were varied as well as the action of the piston on the plate the force was applied either throughout the plate section or on its mantle layer only in the latter case this layer slipped relative to the upper crust hence different results were obtained a more realistic pattern seems to occur when the yield limit of the crustal asthenosphere is not too small and free sliding of the crustal layers over the mantle lithosphere is not possible the upper crustal layer then experiences failure in the collision zone and thrusting develops there intensive thickening of the lower crust was observed the mantle layer being subducted like oceanic lithosphere in a subduction zone

terrane motion by strike slip faulting of forearc slivers forearc slivers bounded by a trench and an active strike slip fault occur in about 50 of modern subduction zones analysis of earthquake slip vectors indicates that modern sliver terranes typically migrate at rates of 1 2 cm yr tertiary transport of some forearc slivers by 1000 km or more is therefore expected as a consequence of normal subduction active arc parallel strike slip faulting occurs whenever interplate coupling is strong and convergence is somewhat oblique strike slip rate increases with greater convergence obliquity

strain accumulation controls failure of a plate boundary zone linking deformation of the central andes and lithosphere mechanics 1 we make use of observations on orogenic strain accumulation and deformation partitioning in the central andes to explore the backarc strength evolution at the lithospheric scale in plan view the altiplano puna plateaux experienced rapid initial increase of surface area undergoing active deformation during the cenozoic beyond the maximum lateral extent reached around 10 15 ma 40 50 of entire proto andes undergoing deformation at 10 20 total strain rapid localization initiated at the eastern flank of the altiplano inter and subandean thrust belt but not at the puna latitude localization was associated with a significant increase in bulk shortening rate average fault slip rates equally increased by an order of magnitude following a protracted period of stable average rates estimates of strength evolution based on force balance calculations and critical wedge analysis suggest significant backarc weakening driving this change after the middle miocene strain accumulation led to localization and weakening with development of a detachment propagating through crust and upper mantle we find that lithosphere scale failure resulting from strain weakening beyond a critical strain threshold c 20 and fault coalescence with formation of a weak detachment in shales effective coefficient of friction 0 1 plays a key role in the evolution of the andes strain related lithosphere weakening appears to dominate over the impact of external forcing mechanisms such as variations of plate convergence mantle assisted processes or erosion comparison of these orogen scale observations with experimental rock rheology indicates substantial similarity of deformation behavior with similar weakening thresholds across a wide range of scales

experimentally based water budgets for dehydrating slabs and consequences for arc magma generation phase diagrams of hydrous mid ocean ridge mor basalts to 330 km depth and of hydrous peridotites to 250 km depth are compiled for conditions characteristic for subduction zones a synthesis of our experimentally determined phase relations of chlorite lawsonite epidote zoisite amphibole paragonite chloritoid talc and phengite in basalts and of phase relations from the literature of serpentine talc chlorite amphibole and phase a in ultramafics permits calculation of h2o contents in hydrous phase assemblages that occur in natural compositions this yields the information necessary to calculate water budgets for descending slabs starting from low grade blueschist conditions 10 20 km depth with h2o contents between 5 and 6 wt for hydrated oceanic crust complete dehydration is achieved between 70 and 300 km depth as a function of individual slab geotherms hydrous phases which decompose at depth below volcanic arcs are lawsonite zoisite chloritoid and talc phengite in mafic compositions and chlorite and serpentine in peridotite approximately 15 35 of the initially subducted h2o are released below volcanic arcs the contribution of amphibole dehydration to the water budget is small 5 20 and occurs at relatively shallow depth 65 90 km in any predicted thermal structure dehydration is a combination of a stepwise and a continuous process through many different reactions which occur simultaneously in the different portions of the descending slab such a dehydration characteristic is incompatible with single phase dehydration models which focus fluid flow through a unique major dehydration event in order to explain volcanic fronts as a consequence of continuously progressing dehydration water ascending from the slab will be generally available to depth of ca 150 200 km the fluid rising from the subducting lithosphere will cause partial melting in the hot portion of the mantle wedge we propose that the volcanic front simply forms above the mantle wedge isotherm where the extent of melting is sufficient to allow for the mechanical extraction of parental arc magmas thermal models show that such an isotherm ca 1300oc locates below volcanic fronts slab surface depths below such an isotherm are compatible with the observed depths of the slab surface below volcanic fronts 1998 elsevier science b v all rights reserved

subduction of the lithosphere and back arc dynamics insights from physical modeling physical modeling of the subduction of an elastico plastic lithosphere underlained by a low viscosity asthenosphere and driven by both horizontal compressional force and gravitational sinking of the subducting plate has shown the subduction regime to be strongly affected by 1 a where a and 1 are the asthenosphere and lithosphere densities at 0 a steady state process develops which requires a horizontal lithosphere compression of h 0 3 0 4 s 108 pa where s is the lithosphere yield limit the central part of the leading wedge of the overriding plate is nonisostatically elevated in this case especially in the region of the frontal arc which corresponds to a free air gravity anomaly maximum of 200 400 mgal a decrease in 0 compression of the overriding plate is considerably less from the onset of subduction and then decreases during subduction at some point compression changes into tension tensional nonhydrostatic stresses grow during the process to a value near h 0 1 0 2 s 0 4 108 pa both subsidence of the overriding wedge below the isostatic equilibrium level and an increase in the benioif dip angle accompany this process the gravity anomaly maximum is displaced into the volcanic arc in this subduction regime where it is normally less than 150 mgal if some weaknesses are present in the overriding plate tensional stresses result in a failure of the lithosphere at these locations and the broken lithosphere segment the arc plate starts to move counter to the subduction plate the arc plate slides along the underthrusting plate s surface without loss of contact and causes back arc opening hydrostatic suction prevents the plates from being separated this effect is caused by the peculiarity of a mechanical system whereby a low viscosity asthenosphere and a plastic lithosphere which effective yield limit is considerably lower than the hydrostatic stresses in the plate if the overriding lithosphere is not weakened enough to fail under the existing tension then the hanging slab breaks down under its own weight when it reaches a specific length both interacting plates undergo isostatic recovery in this case and the whole cycle is repeated again the hydrostatic suction effect also provides back arc tension under other conditions such as when asthenospheric flow forces the subducted plate to move away from the island arc

on the relationship between the formation of shear zones and the form of the flow law for rocks undergoing dynamic recrystallization abstract the formation and maintenance of shear zones in rocks at high temperatures when viscous flow takes place at constant volume requires strain weakening several processes can lead to a strain weakening characteristic in the flow of such rocks to high strains but conventional experimental rock mechanics techniques cannot normally provide high quality mechanical data to high strains recent experimental data for calcite marble have demonstrated strain weakening accompanied by dynamic recrystallization a power law with strain exponent r can be used to describe the dependence of the flow stress on strain together with the conventional description of the dependence of strain rate on stress described by a stress exponent n the stable localization of strain into shear zones is favoured by more negative values of r e g r 0 8 and more positive values of n e g n 5 to 7 a further important factor is the critical strain e0 for the onset of recrystallization and hence strain weakening shear zone formation is favoured by smaller values for e0 shear zone formation weakens the crustal section in which they occur but the amount of weakening may not be great perhaps no more than 50 of the previous steady state flow stress certain rock types such as mafic igneous rocks may through their mechanical properties be more susceptible to shear zone formation than others such as mica schists it is emphasized that the analysis presented is based upon a strain dependent flow law that describes single phase dynamic recrystallization by grain boundary migration that leads to enhanced recovery without any change in deformation mechanism no consideration is given to possibly important effects arising from the activation of grain size sensitive flow or to potential weakening effects arising from syntectonic metamorphic reactions

shocks in a box an analogue model of subduction earthquake cycles with application to seismotectonic forearc evolution 1 we introduce and test an experimental approach to simulate elastoplastic megathrust earthquake cycles using an analogue model and apply it to study the seismotectonic evolution of subduction zones the quasi two dimensional analogue model features rate and state dependent elastic frictional plastic and viscoelastic material properties and is scaled for gravity inertia elasticity friction and viscosity the experiments are monitored with a high resolution strain analysis tool based on digital image correlation particle imaging velocimetry piv providing deformation time series comparable to seismologic geodetic and geologic observations in order to separate elastic and nonelastic effects inherent the experimental deformation patterns we integrate elastic dislocation modeling edm into a hybrid approach we use the analogue earthquake slip and interseismic locking distribution as edm dislocation input and forward model the coseismic and interseismic elastic response the residual which remains when the edm prediction is subtracted from the experimental deformation pattern highlights the accumulation of permanent deformation in the model the setup generates analogue earthquake sequences with realistic source mechanisms and elastic forearc response and recurrence patterns and reproduces principal earthquake scaling relations by applying the model to an accretionary type plate margin we demonstrate how strain localization at the rupture peripheries may lead to a seismotectonically segmented forearc including a tectonically stable shelf and coastal high 20 plate convergence accommodated by internal shortening overlying the area of large megathrust earthquake slip fifty to 75 of plate convergence is accommodated by internal shortening in the slope region where earthquake slip tapers out toward the trench the inner forearc region remains undeformed and represents a basin

development of the gulf of guayaquil ecuador during the quaternary as an effect of the north andean block tectonic escape 1 interpretation of industrial multichannel seismic profiles and well data are used to identify the main tectonic features of the gulf of guayaquil gg area these include two e w trending major detachments the posorja and the jambeli detachment systems which represent half grabens with oppositely dipping detachments to the south and to the north respectively the ne sw trending puna santa clara fault system developed as a transfer fault system between the posorja and the jambeli detachments the esperanza and the jambeli basins exhibit 3 4 km of sediment that accumulated during the past 1 6 1 8 myr the puna santa clara fault system bounds the esperanza and the jambeli basins evidencing that the evolution of these basins is tightly controlled by the two detachments at depth to the west the n s trending domito fault system bounding the posorja detachment system and the esperanza basin to the west acted as a transfer zone between the shelf area and the continental slope the pliocene series show no significant variations in thickness throughout the gulf of guayaquil area suggesting that no important tectonic deformation occurred from 5 2 to 1 8 1 6 ma the major period of tectonic deformation in the gulf of guayaquil area occurred during the pleistocene times three main tectonic steps are identified from early pleistocene to 180 ka major subsidence occurred along the esperanza and jambeli basins from 180 to 140 ka most of the gulf of guayaquil area was above sea level during the isotope substage 6 low stand from 140 ka to present tectonic activity is restricted along the normal faults bounding the esperanza basin the tenguel fault and the puna santa clara and domito fault systems a n s trending tensional stress regime characterizes the pleistocene times throughout the northward drifting of the north andean block is proposed to control the tectonic evolution and associated subsidence of the gulf of guayaquil area it is also accepted that the collision of the carnegie ridge with the trench axis has to play a major role in controlling the north andean block northward drift because the carnegie ridge subduction is possibly an ongoing process which began prior to the pliocene we postulate the along strike morphology of the ridge at the origin of interplate coupling variations the subduction of an along strike positive relief of the ridge is proposed at the origin of the major tectonic reorganization of the gg area occurring at 1 8 1 6 ma

oceanic plateaus the fragmentation of continents and mountain building many anomalous rises in today s oceans may be submerged continental fragments detached from previous continents ancient island arcs or basaltic piles formed by hot spots and spreading centers these rises are embedded in their respective moving oceanic plates and are fated to be consumed at active margins where such rises are being consumed at present e g the nazca ridge they cause cessation of volcanism disruption of the downgoing slab and possible shifts in plate boundary configuration many past rises including numerous continental fragments have been recognized within mountain belts as allochthonous terranes they constitute a large portion of the orogenic belts in the north pacific from mexico through western north america alaska east siberia japan and in new zealand the orogenic deformation in these belts is possibly the result of the accretion of the allochthonous terranes many terranes have been accreted with substantial deformation also in the alpine chain well before major continent continent collisions it is suggested therefore that the accretion of fragments may be the common process of the deformation phase of mountain building subduction of normal oceanic crust may be insufficient for deformation whereas full continent continent collision may not be necessary the general validity of this conclusion depends critically on whether allochthonous terranes caused orogenic deformation in the andes or not most of the accreted fragments with continental affinities in the mesozoic cenozoic orogenic belts of the world can be traced back to the breakup of gondwana beginning with a pacifica domain in the permian through a larger india domain in the early mesozoic and continuing through the separation of the somalia plate in the near future the reasons for this 250 million year breakup process are not known but some kind of thermal process possible of mantle wide scale is implied

lithosphere tearing at step faults response to edges of subduction zones abstract slab edges are a relatively common feature in plate tectonics two prominent examples are the northern end of the tonga subduction zone and the southern end of the new hebrides subduction zone near such horizontal terminations of subduction trenches ongoing tearing of oceanic lithosphere is a geometric consequence we refer to such kinks in the plate boundary as a subduction transform edge propagator or step other steps are the north and south ends of the lesser antilles trench the north end of the south sandwich trench the south end of the vrancea trench and both ends of the calabria trench volcanism near steps is distinct from typical arc volcanism in some cases slab edges appear to coincide with mantle plumes using 3d mechanical models we establish that step faults are stable plate tectonic features in most circumstances in the probably rare cases that the resistance to fault propagation is high slab break off will occur relative motion along the transform segment of the plate boundary often is non uniform and the step is not a transform plate boundary in the rigid plate tectonics sense of the phrase step propagation may result in substantial deformation rotation topography and sedimentary basins with a very specific time space evolution surface velocities are substantially affected by nearby steps

magmatic processes under arcs and formation of the volcanic front temperature structure and stress under arcs are simulated in a two dimensional cross section taking into account the flow induced by the subducting slab in the mantle wedge results of the calculations show three important features with respect to magmatic processes under arcs 1 temperature structure in the crust and the mantle wedge under arcs is insensitive to the angle and velocity of slab subduction the temperature structure of the slab and that of the back arc region this indicates that physical conditions such as temperature and pressure are similar under various arcs it is thus inferred that primary magmas generated under various arcs should have similar chemical compositions if chemical composition and the flux rate of fluid from the slab are similar and the chemical compositions of mantle wedge materials are the same 2 calculated deviatoric stress magnitude is relatively large more than a few tens of megapascals in the partially molten mantle cracks may open under high differential stress and magma can easily segregate and accumulate through interconnected cracks while the buoyancy driven compaction of partially molten mantle proceeds 3 the deviatoric stress values in the region over the partially molten mantle are relatively large and the direction of the principal stress changes horizontally the direction of the maximum compressional stress is nearly vertical under the volcanic zone and is nearly horizontal under the fore arc region it is considered that magma segregated in partially molten mantle migrates upward through the brittle mantle and crust by the magma fracturing mechanism the propagation direction of magma filled fissures is controlled by the stress field in the crust and mantle and is parallel to the maximum principal stress the calculated stress is highly compressive horizontally on the trench side while it becomes tensile on the back arc side the location of this stress transition coincides with that of the volcanic front the location of this transition indicates that the volcanic front marks a change in the ease of upward migration of the magma filled cracks under relatively high differential stress field

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

a unified description of localization for application to large scale tectonics 1 localized regions of deformation such as faults and shear zones are ubiquitous in the earth s lithosphere however we lack a simple unified framework of localization that is independent of the mechanism or scale of localization we address this issue by introducing the effective stress exponent ne a parameter that describes how a material responds to a local perturbation of an internal variable being tested for localization the value of ne is based on micromechanics a localizing regime has a negative ne indicating a weakening behavior and localization is stronger for more negative 1 ne we present expressions for the effective stress exponent associated with several mechanisms that trigger localization at large scale brittle failure with loss of cohesion elastoplasticity rate and state dependent friction shear heating and grain size feedback in ductile rocks in most cases localization does not arise solely from the relation between stress and deformation but instead requires a positive feedback between the rheology and internal variables brittle mechanisms failure and friction are generally described by ne of the order of 100 shear heating requires an already localized forcing which could be provided by a brittle fault at shallower levels of the lithosphere grain size reduction combined with a transition from dislocation to diffusion creep leads to localization only if the grain size departs significantly from its equilibrium value because either large scale flow moves rocks through different thermodynamic environments or new grains are nucleated when shear heating or grain size feedback produce localization 1 ne can be extremely negative and can control lithospheric scale localization

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

active tectonics of the andes nearly 90 mm a 1 of relative plate convergence is absorbed in the andean plate boundary zone the pattern of active tectonics shows remarkable variations in the way in which the plate slip vector is partitioned into displacement and strain and the ways in which compatibility between different segments is solved along any traverse across the plate boundary zone the sum of relative velocities between points must equal the relative plate motion we have developed a kinematic synthesis of displacement and strain partitioning in the andes from 47 s to 5 n relevant for the last 5 ma based upon 1 relative plate motion deduced from oceanic circuits giving a roughly constant azimuth between 075 and 080 2 moment tensor solutions for over 120 crustal earthquakes since 1960 3 structural studies of deformed plio pleistocene rocks 4 topographic geomorphic studies 5 palaeomagnetic data and 6 geodetic data we recognize four neotectonic zones with subzones and boundary transfer zones that are partitioned in different ways these zones are not coincident with the classic zones defined by the presence or absence of a volcanic chain or differences in finite displacements and strains and tectonic form the long term segmentation and finite evolution of the andes may not occur in constantly defined segments in space and time in segment 1 47 39 s the slip vector is partitioned into roughly orthogonal benioff zone slip with large magnitude large slip surface earthquakes and both distributed dextral shear giving clockwise rotations of up to 50 and dextral slip in the curved liquine ofqui fault system giving 5 10 of anticlockwise fore arc rotation in segment 2 39 20 s the slip vector is partitioned into benioff zone slip roughly parallel with the slip vector andean crustal shortening and a very small component of dextral slip including that on the atacama fault system between 39 and 34 s a cross strike dextral transfer which deflects the chile trench and the volcanic arc absorbs the shortening contrast between segments 1 and 2 in segment 3 20 6 s the slip vector is partitioned into roughly orthogonal benioff zone slip crustal shortening sinistral trench parallel faulting and northeast southwest extension compatibility between segments 2 and 3 is maintained by the sinistral east southeast trending cochabamba shear zone and north trending dextral faults in segment 4 6 s to 5 n the slip vector is partitioned into roughly orthogonal benioff zone slip and dextral strike slip faulting in the fore arc and volcanic chain

paleogene building of the bolivian orocline tectonic restoration of the central andes in 2 d map view using available information on the magnitude and age of tectonic shortening as well as paleomagnetically determined tectonic rotations we have run a series of 2 d map view restorations of the central andes neogene shortening in the foreland belt induced only slight orogenic curvature of the central andes the constraints on the ages of the large observed fore arc rotations average of 37 counterclockwise in southern peru and 29 clockwise in northern chile indicate that the bolivian orocline formed during the eocene oligocene as a consequence of differential shortening focused in the eastern cordillera to minimize local block rotations along the fore arc the restoration that best fits the central andean rotation pattern requires about 400 km of total paleogene plus neogene shortening near the arica bend this value corresponds to the upper bound of estimates of maximum horizontal shortening for the central andes along strike variations in horizontal shortening in the back arc induced bending of the continental margin block rotations with fore arc along strike extension and or orogen parallel transport of upper crustal material toward the symmetry axis of the orocline

strain partitioning and interplate friction in oblique subduction zones constraints provided by experimental modeling physical modeling of oblique subduction is performed to study the mechanism of strain partitioning the model is two layer and includes the elasto plastic lithosphere the overriding and subducting plates and the low viscosity liquid asthenosphere the subduction is driven by a push force from a piston and a pull force when the density contrast ap between the subducting plate and the asthenosphere is positive we vary both ap and the interplate friction frictional stresses slip partitioning is obtained only in the models with high interplate friction and only when the overrid ing plate contains a weak zone this zone in the models corresponds either to locally thinned litho sphere or to cut fault the horizontal trench normal component of the interplate friction force ffh can be comparable with the absolute value of the horizontal component of the nonhydrostatic interplate pressure force fph in the subduction zone ffh is always negative compression while fph can be either negative compressional subduction regime or positive extensional regime high friction which promotes partitioning can completely cancel the extensional suction force fph back arc tension and strike slip faulting appear thus as conflicting processes although they can coexist in the same subduction zone depending on the relative values of relevant forces it ap pears that high friction can exist only in compressional subduction zones where partitioning should develop more easily this conclusion is supported by the comparison of two oblique subduction zones having similar geometry the compressional southern kurile zone strong partitioning and extensional southern ryukyu zone no lithospheric scale partitioning other factors controlling the strain partitioning are the length of the oblique subductlon zone the boundary conditions at the transverse limits of the forearc sliver and of course the obliquity of subduction