no flat wadati benioff zone in the central and southern central andes abstract the wadati benioff zone wbz is an approximate plane defined by earthquakes hypocentres observed in convergent plate boundaries and that usually dips at angles greater than 30 in some areas of the andes where there are gaps in volcanic activity and where heat flow is abnormally low this plane in most studies has nearly horizontal dip at a depth of about 75 100 km and it has been associated to flat subduction of the oceanic lithosphere this situation has been taken as the present day analogue of the laramide orogeny of western north america for which a flat slab episode has been proposed in the past years in this work the observed low heat flow in areas of the andes is assumed to be due to low radiogenic heat generation in geologically old and allochthonous terranes constituting large regions of western south america on the basis of geotherms obtained for areas of ecuador peru chile and argentina and of rheological results describing the partition between brittle and ductile regimes the seismic activity observed both in the lower crust and at depths of about 75 100 km is thoroughly explained at these depths earthquakes occur within the subcontinental upper mantle and then there is no flat wbz associated to subduction of the oceanic lithosphere there is evidence from recent seismological observations that the real wbz lies not horizontally and deeper in the tectonosphere

timing of the acadian orogeny in northern new hampshire abstract new u pb geochronology constrains the timing of the acadian orogeny in the central maine terrane of northern new hampshire sixteen fractions of one to six grains each of zircon or monazite have been analyzed from six samples 1 an early syntectonic diorite that records the onset of the acadian 2 a schist a migmatite and two granites that together record the peak of the acadian and 3 a postkinematic pluton that records the end of the acadian zircon from the syntectonic wamsutta diorite gives a 207pb 206pb age of circa 408 ma the time at which the boundary between the deforming orogenic wedge and the foreland basin was in the vicinity of the presidential range this age agrees well with the emsian position of the northwest migrating acadian orogenic front and records the beginning of the acadian in this part of the central maine terrane we propose a possible acadian tectonic model that incorporates the geochronologic structural and stratigraphic data monazite from the schist migma

evolving geographic patterns of cenozoic magmatism in the north american cordillera the temporal and spatial association of magmatism and metamorphic core complexes four maps are presented here that show the location and extent of magmatic fields between eastern alaska and northern mexico during the successive time intervals of 55 40 40 25 25 10 and 10 0 ma and four others show the distribution of metamorphic core complexes during the same cenozoic time intervals the maps are based on u s geological survey and canadian cordilleran data bases contining about 6000 isotopic dates and extensive literature review for nearly 60 ma the development of metamorphic core complexes has coincided with the locus of a really extensive and voluminous intermediate composition magmatic fields the association is suggestive of a close link between magmatism and core complex formation namely that magma directly and indirectly lowers the strength of the crust magmatism thus controls the location and timing of core complex formation the stresses responsible may be inherited from mesozoic crustal thickening locally created by uplift and magmatic thickening of the crust and imposed by the global pattern of plate motions and driving forces since the miocene rates of magmatism extension and core complex formation have declined the modern basin and range province is not a suitable model for the situation that existed during major magmatic culminations the singular event of early miocene time the merging of two large magmatic fields extinguishing the laramide magmatic gap explains several disconnected observations the hyperextension episode of the colorado river corridor rapid reorientation of stress patterns across much of western north america and subsequent rapid tectonic movements in california magma triggered breakup of western north america lithosphere coincided with development of the san andreas transform system thermal destruction of the laramide magmatic gap created a california microplate about 22 ma ago that moved rapidly away from north america thus two plate tectonic processes thermal destruction of the lithosphere bridge and northward growth of a transform system interacted to produce miocene and later tectonic patterns and events

the late cretaceous to recent tectonic history of the pacific ocean basin abstract a vast ocean basin has spanned the region between the americas asia and australasia for well over 100 myr represented today by the pacific ocean its evolution includes a number of plate fragmentation and plate capture events such as the formation of the vancouver nazca and cocos plates from the break up of the farallon plate and the incorporation of the bellingshausen kula and aluk phoenix plates which have been studied individually but never been synthesised into one coherent model of ocean basin evolution previous regional tectonic models of the pacific typically restrict their scope to either the north or south pacific and global kinematic models fail to incorporate some of the complexities in the pacific plate evolution e g the independent motion of the bellingshausen and aluk plates thereby limiting their usefulness for understanding tectonic events and processes occurring in the pacific ocean perimeter we derive relative plate motions with 95 uncertainties for the pacific farallon vancouver kula pacific bellingshausen pacific and early pacific west antarctic spreading systems based on recent data including marine gravity anomalies well constrained fracture zone traces and a large compilation of magnetic anomaly identifications we find our well constrained relative plate motions result in a good match to the fracture zone traces and magnetic anomaly identifications in both the north and south pacific in conjunction with recently published and well constrained relative plate motions for other pacific spreading systems e g aluk west antarctic pacific cocos recent pacific west antarctic spreading we explore variations in the age of the oceanic crust seafloor spreading rates and crustal accretion and find considerable refinements have been made in the central and southern pacific asymmetries in crustal accretion within the overall pacific basin where both flanks of the spreading system are preserved have typically deviated less than 5 from symmetry and large variations in crustal accretion along the southern east pacific rise i e pacific nazca farallon spreading appear to be unique to this spreading corridor through a relative plate motion circuit we explore the implied convergence history along the north and south america s where we find that the inclusion of small tectonic plate fragments such as the aluk plate are critical for reconciling the history of convergence with onshore geological evidence

mesozoic aseismic ridges on the farallon plate and southward migration of shallow subduction during the laramide orogeny we propose that an aseismic ridge of late jurassic and early cretaceous age on the farallon plate was subducted beneath the unites states cordillera during the laramide orogeny the relative buoyancy of the aseismic ridge caused shallow subduction which resulted in a magmatic lull 70 40 ma in the existing near trench magmatic arc and laramide uplift and faulting of crystalline basement 1000 1500 km inland from the trench the timing and latitudinal limits of the subduction of the aseismic ridge are similar to the timing and latitudinal limits of both the laramide orogeny and the magmatic lull in the western cordillera because the north northeast trending aseismic ridge was moving northeast relative to the north american plate the point of collision of the aseismic ridge with north america and therefore the locus of shallow subduction migrated southward with time causing the well documented southward migration of the magmatic lull

miocene quaternary structural evolution of the uyuni atacama region andes of chile and bolivia abstract we describe the miocene quaternary geological structural evolution of the region between the salar de uyuni and de atacama andes of chile and bolivia we recognized four main tectonic events based on fold geometry fault kinematics and stratigraphic data the oldest event of miocene age is characterized by folding and reverse faulting of the sedimentary successions with an e w direction of shortening in the northern part of the studied area and a wnw ese shortening in the southern part the following two events of pliocene age are characterized by lower shortening amounts they occurred first by reverse faulting with a nw se trending greatest principal stress 1 computed with striated fault planes and a vertical least principal stress 3 followed by pervasive strike slip faulting with the same nw se trending 1 and a horizontal ne sw 3 the fourth event dating to the late pliocene quaternary is characterized by normal faulting the 3 still trends ne sw whereas the intermediate principal stress 2 exchanged with 1 volcanism accompanied both the contractional transcurrent and extensional tectonic phases the mio pliocene compression appears directly linked to a rapid convergence and an apparently important coupling between the continental and oceanic plates the e w to wnw ese direction of shortening of the miocene structures and the nw se 1 of the pliocene structures seem to be more linked to an intra andean re orientation of structures following the wnw directed absolute motion of the south american plate the extensional deformations can be interpreted as related to gravity forces affecting the highest parts of the volcanic belt in a sort of asymmetrical sw ward collapse of the belt

crust and upper mantle structure beneath the pacific northwest from joint inversions of ambient noise and earthquake data 1 we perform a joint inversion of phase velocities from both earthquake and ambient noise induced rayleigh waves to determine shear wave velocity structure in the crust and upper mantle beneath the pacific northwest we focus particularly on the areas affected by mid miocene to present volcanic activity the joint inversion combined with the high density seismic network of the high lava plains seismic experiment and data from the earthscope transportable array provides outstanding resolution for this area in oregon we find that the pattern of low velocities in the crust and uppermost mantle varies between the high lava plains physiographic province and the adjacent northwestern basin and range these patterns may be due to the presence of the brothers fault zone which separates the clockwise rotating northwest basin and range from the relatively undeformed areas further north further to the east the owyhee plateau snake river plain srp and northeastern basin and range are characterized by high crustal velocities though the depth extent of these fast wave speeds varies by province of particular interest is the mid crustal high velocity sill previously only identified within the srp we show this anomaly extends significantly further south into utah and nevada we suggest that one possible explanation is lateral crustal extrusion due to the emplacement of the high density mafic mid crustal sill structures within the srp

permo triassic tectonism in volcanic arc sequences of the western u s cordillera and implications for the sonoma orogeny the permo triassic sonoma orogeny western u s cordillera involved tight folding and imbricate thrust faulting of basinal strata of the golconda allochthon this event is generally interpreted in terms of accretion of upper paleozoic volcanic arc rocks of the western u s cordillera mccloud arc to the north american continent either via back arc basin closure or arc continent collision the structural history of the mccloud arc during the sonoma orogeny is poorly understood however as are the implications of this history for tectonic models this paper presents a synthesis of the permian to triassic stratigraphic and structural evolution of the mccloud arc incorporating new data from the pine forest range northwest nevada the synthesis indicates that on a regional scale the arc was affected by widespread uplift but only local tilting folding and faulting in the late permian to early triassic the timing of permo triassic tectonic activity in the mccloud arc is thus similar to the middle permian to early triassic age of deformation in the golconda allochthon but the structures produced in the two groups of rocks differ markedly the significance of this contrasting structural response to sonoma age tectonism can be evaluated by comparison with other better understood examples of arc continent accretion back arc basin closure in the southern andes compressional arc continent collision in new guinea and extensional arc continent collision in the mediterranean this comparison indicates that the contrasting structural history of the mccloud arc and golconda allochthon is entirely compatible with an arc continent accretion model for the sonoma orogeny and cannot be used as a basis for concluding that the mccloud arc was not involved in the orogenic event the comparison further indicates that either back arc basin closure or compressional arc continent collision are the most likely settings for the sonoma orogeny in contrast an extensional mediterranean style arc continent collision as has been recently proposed for the earlier middle paleozoic antler orogeny is probably not an appropriate analogue for the sonoma orogeny the weak structural response of the mccloud arc during sonoman compression also provides new constraints on the specific processes involved in arc accretion either via arc continent collision or back arc basin closure these constraints suggest that there may be few fundamental differences between the two models during the final stages of arc accretion regardless of tectonic setting permo triassic uplift and associated deformation in the mccloud arc should be considered as much a response to the sonoma orogeny as is tight folding and imbricate faulting in the golconda allochthon

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

hypocentral trend surface analysis probing the geometry of benioff zones a hypocentral trend surface is a continuous function of latitude and longtitude fitted by least squares to a set of hypocenters so that it predicts depth to the middle of a benioff zone in this paper we take a relatively simple approach to hypocentral trend surface analysis the hypocentral trend surface is constructed from a spherical surface harmonic expansion whose coefficients are selected so as to minimize the standard vertical deviation between hypocenters and trend surface the vertical deviation of a hypocenter from the hypocentral trend surface is called its residual consideration of a hypocentral trend surface cannot be divorced from consideration of the associated residuals and their spatial distribution any analysis of benioff zone geometry trades off presumed thickness of the zone with its presumed shape this trade off can be investigated by examining suites of candidate hypocentral trend surfaces together with their associated residual distributions generated by varying the number of degrees of freedom available to the trend surface hypocentral trend surfaces are generated for three high quality sets of hypocenters obtained by local seismic networks in honshu japan cook inlet alaska and south peru hypocenters beneath honshu generate a bimodal distribution of residuals about their trend surface the vertical separation of the upper and lower sheets of this double benioff zone averages about 33 km no double benioff zone configuration is discovered below cook inlet or south peru regional trend surfaces are established for intermediate depth teleseismic data from south and middle america in south america 0 40 s teleseismic data suggest that along strike transitions between flat and moderately steeply dipping sections of the benioff zone are achieved by flexure of a coherent slab rather than by fragmentation of the slab into tear bound flaps or separate tongues of lithosphere hypocentral trend surface analysis of teleseismic data from south peru suggests a geometry similar to that indicated by hypocentral trend surface analysis of the local data of hasegawa and sacks an eastward shallowing benioff zone beneath southwest mexico is closely associated with the geometry of the trans mexican volcanic arc the benioff zone beneath central america is more steeply dipping the flat slab steep slab transition near 96 w may relate to the difference in age of oceanic lithosphere on either side of the tehuantepec ridge while the strong flexure of the subducted nazca plate beneath ecuador is closely associated with the carnegie ridge the nazca ridge extropolates slightly north of the flexure in subducted lithosphere beneath south peru the flexure appears instead more closely associated with the pronounced inflection in the map view shape of the leading edge of the overriding south american plate

geodynamic evolution of southern costa rica related to low angle subduction of the cocos ridge constraints from thermochronology abstract the late tertiary shallow subduction of the cocos ridge under the caribbean plate controlled the evolution of the cordillera de talamanca in southeast costa rica which is a mountain range that consists mainly of granitoids formed in a volcanic arc setting fission track thermochronology using zircon and apatite as well as 40ar 39ar and rb sr age data of amphibole and biotite in granitoid rocks constrain the thermal history of the cordillera de talamanca and the age of onset of subduction of the cocos ridge shallow intrusion of granitoid melts resulted in fast and isobaric cooling a weighted mean zircon fission track age 13 analyses and rb sr biotite ages of about 10 ma suggest rapid cooling and give minimum ages for granitoid emplacement in some cases 40ar 39ar and rb sr apparent ages of amphibole and biotite are younger than the zircon fission track ages which can be attributed to partial resetting by hydrothermal alteration apatite fission track ages range from 4 8 to 1 7 ma but show no correlation with the 3090 m elevation span over which they were sampled the apatite ages seem to indicate rapid exhumation caused by tectonic and isostatic processes the combination of the apatite fission track ages with subduction parameters of the cocos plate such as subduction angle plate convergence rate and distance of the cordillera de talamanca to the trench implies that the cocos ridge entered the middle america trench between 5 5 and 3 5 ma

model for late mesozoic early tertiary tectonics of coastal california and western mexico and speculations on the origin of the san andreas fault cretaceous and early tertiary arc and forearc rocks found along the coast in southern california and baja california have been shown by paleomagnetic measurements to have originated many hundreds or even several thousands of kilometers south of their current locations northward transport also is found in cretaceous batholithic rocks near the edge of the continent in washington and british columbia the consistency of this pattern suggests that slices of arc and forearc rock originating on the western edge of north america have been translated along the coast by strike slip faulting strike slip faults are shown to be an expected response to north oblique subduction of the farallon and or kula plates faulting and subduction probably were concurrent activities i e the arc and forearc grew and were displaced simultaneously because displacement was dominantly parallel to subduction related lithic belts arc and forearc rocks may have traveled large distances without producing either significant hiatuses in the rock record or conspicuous disruptions in the distribution of lithic belts subduction related strike slip faulting is shown to be favored by high angles of oblique convergence shallow dip of the subducting slab relatively easy slip on potential faults within the arc or forearc and the existence of a place for the moving sliver to go it is suggested that these conditions are not difficult to satisfy and that subduction related strike slip faulting may thus be a common feature of convergent orogens the location of the san andreas fault inboard from the continental margin may be attributable to localization of the pacific north american transform by a system of precollision subduction related faults

three dimensional numerical simulation of the subduction dynamics in the sunda arc region southeast asia abstract the sunda arc has been modelled as a test subduction zone with geometrical complexities and varying dynamics we simulated spatially the structural heterogeneities of the subduction zone and carried out our computations by a three dimensional finite element approach the aim is to estimate the lateral variation in the effective horizontal compressive stress in the interplate region which can produce the present day observed seismogenic stress field in our simplified approach the stress field has been modelled by varying the traction ridge push while keeping the age dependent gravitational body force slab pull constant our result suggests that there is a considerable amount of stress dissipation in the decoupled segment of java at sumatra the required traction is much greater than that expected from the age dependency relations in our modelling the transition from the dominance of the horizontal compressive stress to that of the gravitational pull in the regional stress field occurs around a depth range of 150 200 km lateral variation in the subducted slab volume emerges crucial to the seismogenic stress field and the derived pattern of displacement subduction of the lighter oceanic rise material should enhance the stress level the effect of abrupt lateral variation in structure e g sudden changes in the slab dip or slab depth may also be significant on the regional stress field viscoelastic modelling attempted for estimating the temporal variation in subduction dynamics suggests that the thrust level will increase southward of the present day thrust zone in java while strike slip fault type stress field will be developed south of sumatra

plio quaternary vertical motions and the subduction of the nazca ridge central coast of peru abstract the coastal region of south central peru displays raised marine terraces indicating that unlike the central peruvian coast it underwent a strong uplift since the late pliocene these recent vertical movements which were probably the most rapid motion experienced along the pacific coast of south america are commonly considered to be closely related to the subduction of the aseismic nazca ridge however several questions regarding the geometry and the precise mechanism s responsible for this deformation have not yet been fully addressed the kinematic constraints of the andean convergence suggest that the deformation of the coastal area between pisco and lomas is characterized by a quite continuous sequence of uplift events progressing southwards in spite of some difficulties encountered in precisely dating the uplifted quaternary terraces and other remnants of pliocene marine platforms it is proposed that uplift rates during the late pleistocene reached maximum values of the order of 700 mm 10 3 y as against a recently proposed estimate of 470 mm 10 3 y the reconstructed finite deformation as recorded by the study area for the last 3 ma presents a nw se oriented asymmetric dome shaped pattern it is emphasized that the strongest uplift has been occurring immediately south of the inland projection of the nazca ridge and not along the ridge axis as envisioned by the previously constructed physical models

the effect of oceanic crustal structure on phase changes and subduction abstract recent numerical models for the thermal regime of subduction zones are in disagreement with observations of depth to formation of eclogite which require higher temperatures than those found in the models this discrepancy may be reconciled if more accurate geometric depictions of the subducted slab are incorporated in the modelling studies the transition from 15 dip along the thrust plane to 30 along the benioff zone starts at depths of about 40 km dehydration reactions and the gabbro eclogite phase change are presumed to increase the average density of the lithosphere and to cause the steepening of the subducted lithosphere the benioff zone seismicity which starts at about 40 km depth may itself be a result of stresses caused by the phase changes since the gabbro eclogite transformation occurs in gabbroic sections of the crust thickened crustal sections beneath oceanic plateaus will affect subduction first the width of the thrust zone is reduced due to the increased dip angle because of early transformation of lower and warmer crustal sections second incomplete conversion of the entire crustal thickness results in an average density intermediate between unsubducted oceanic lithosphere and fully transformed normal lithosphere the first effect leads to decreased thrust mechanism seismicity because the coupling between the plates is reduced the second effect leads to flat lying subcontinental subduction or arcward advance of the trench in suboceanic subduction

ablative subduction a two sided alternative to the conventional subduction model subduction is conventionally modelled as a convergence and overlapping of two semirigid plates such a convergence is one sided in that one lithospheric plate descends while the other remains on the surface we suggest an alternative model one based on the dynamics of fluids in this model the viscous lower lithosphere flows downward and the brittle upper lithosphere deforms in passive response this process is potentially double sided for we find that even a buoyant plate can be dragged downward by a dense descending neighbor thus an apparent overriding plate may be worn away by a process of viscous ablation with the rate of ablation a function of plate buoyancy we call this process ablative subduction ablative subduction allows us to simply interpret observations concerning slab profiles interplate seismicity back arc tectonics and complex processes such as double subduction and subduction polarity reversal in performing experiments modelling the evolution of simple fluid slabs we find that slab profile is strongly influenced by ablation in the overriding plate when ablation is weak as when a buoyant continent borders the trench deformable slabs adopt shallow andean style profiles these profiles develop over time from an initially steep shape more vigorous ablation rates as might occur in an ocean ocean convergence yield steeper marianas style profiles thus differences between marianas style and andean style slabs may result from differences in ablation rate and subduction duration the occurrence of ablation might not be easily detected at depth because material from subducting and ablating plates adhere closely as a single slab these slabs show downdip deformations that are consistent with seismic focal mechanisms plate deformations associated with ablative subduction are also consistent with observed patterns of seismicity in particular the relative aseismicity of marianas style plate boundaries is consistent with the ablative model we may thus explain how plates can converge aseismically without requiring that these plates be decoupled ablation should be weakest in andean style subduction although pulses of rapid ablation might lead to episodes of crustal shortening in oceanic regions the manifestation of ablation may depend upon tectonic regime in an extensional area like the marianas ablation may be temporally self limiting and ablative cycles may explain conflicting observations of tectonic erosion and accretion in a compressional regime vigorous ablation might lead to observations of double subduction and subduction reversal in each case the behavior of the lower lithosphere is relatively simple what varies is the rate of ablation and the response of the upper lithosphere

an evaluation of proposed mechanisms of slab flattening in central mexico central mexico is the site of an enigmatic zone of flat subduction the general geometry of the subducting slab has been known for some time and is characterized by a horizontal zone bounded on either side by two moderately dipping sections we systematically evaluate proposed hypotheses for shallow subduction in mexico based on the spatial and temporal evidence and we find no simple or obvious explanation for the shallow subduction in mexico we are unable to locate an oceanic lithosphere impactor or the conjugate of an impactor that is most often called upon to explain shallow subduction zones as in south america japan and laramide deformation in the us the only bathymetric feature that is of the right age and in the correct position on the conjugate plate is a set of unnamed seamounts that are too small to have a significant effect on the buoyancy of the slab the only candidate that we cannot dismiss is a change in the dynamics of subduction through a change in wedge viscosity possibly caused by water brought in by the slab

an early miocene arc tholeiitic magmatic dike event from the alboran sea evidence for precollisional subduction and back arc crustal extension in the westernmost mediterranean a discussion is made on the general significance of an early miocene 22 23 ma k ar dating basaltic dike swarm that outcrops within large areas of the central and western betic zone of the betic cordilleras southern spain this dike swarm is thought to represent the remains of one of the earliest magmatic episodes within the the neogene volcanic province that is associated with the late orogenic evolution of the alboran region and undoubtedly offers the best preserved igneous material related to this early stage that is known at present the dike rocks hypabyssal equivalents of andesitic basalts and basaltic andesites have both major and trace element abundance patterns that are akin to those of arc tholeiitic suites the field relationships and areal distribution of the dike swarm suggest on the other hand that a major e w trending rifting belt could have been developing within the alboran block in oligo miocene times previous to its collision with adjacent continental borderlands in the early middle miocene

constraints on subduction geodynamics from seismic anisotropy 1 much progress has been made over the past several decades in delineating the structure of subducting slabs but several key aspects of their dynamics remain poorly constrained major unsolved problems in subduction geodynamics include those related to mantle wedge viscosity and rheology slab hydration and dehydration mechanical coupling between slabs and the ambient mantle the geometry of mantle flow above and beneath slabs and the interactions between slabs and deep discontinuities such as the core mantle boundary observations of seismic anisotropy can provide relatively direct constraints on mantle dynamics because of the link between deformation and the resulting anisotropy when mantle rocks are deformed a preferred orientation of individual mineral crystals or materials such as partial melt often develops resulting in the directional dependence of seismic wave speeds measurements of seismic anisotropy thus represent a powerful tool for probing mantle dynamics in subduction systems here i review the observational constraints on seismic anisotropy in subduction zones and discuss how seismic data can place constraints on wedge slab and sub slab anisotropy i also discuss constraints from mineral physics investigations and geodynamical modeling studies and how they inform our interpretation of observations i evaluate different models in light of constraints from seismology geodynamics and mineral physics finally i discuss some of the major unsolved problems related to the dynamics of subduction systems and how ongoing and future work on the characterization and interpretation of seismic anisotropy can lead to progress particularly in frontier areas such as understanding slab dynamics in the deep mantle

flat slab thermal structure and evolution beneath central mexico recent seismic and magnetotelluric experiments aimed at better characterizing the shape and state of the subducting slab and continental crust beneath central mexico exposed significant differences with conclusions of previous studies a new slab geometry is revealed in which the subducting cocos slab is perfectly flat between 120 to 290 km from the trench after which it plunges into the asthenosphere at a dip angle of 65 in sharp contrast with the previously proposed 20 dip angle seismic tomography studies show negative p wave velocity anomalies 2 to 4 in the mantle wedge beneath the mexican volcanic belt and positive anomalies 2 to 3 for the subducted cocos slab magnetotelluric experiments exposed a very low resistivity area 1 10 m located within the continental crust just below the mexican volcanic arc finally several spots of non volcanic tremors nvts have been recorded inside the continental crust above the flat slab segment while all these experiments provide a better picture of the subduction system beneath central mexico several key processes need further investigation in this study we take advantage of these new observations to better constrain the thermal structure beneath central mexico two different thermal models are computed for a mantle potential temperature t p of 1 350 and 1 450 c respectively the new thermal structures are then converted into p wave velocity anomalies and compared with the observed v p anomalies we found that a t p of 1 450 c produced larger v p anomalies that do not fit the observations however using a t p of only 1 350 c our predicted v p anomalies are positive 2 to 3 for the cold slab and negative 2 to 4 in the mantle wedge these v p estimates are consistent with the observed seismic tomography from p wave arrivals and therefore we conclude that a t p of 1 350 c is a better estimate for the mantle potential temperature beneath central mexico the new thermal model in conjunction with phase diagrams for sediments hydrated basalt and lithospheric mantle have been used to estimate the amount and location of fluids released from the subducting cocos slab several dehydration pulses have been identified along the slab interface where most of the fluids stored in sediments and oceanic crust are released into the overlying continental crust above the flat slab we found a good correlation between the pattern of these dehydration pulses and the location of nvts suggesting that slab dehydration is responsible for triggering the tremors we suggest that nvt bursts localized above the flat slab segment represent the manifestation of ongoing continental crust hydration and weakening a process that has been going on since 15 ma ago when the cocos slab entered into a flat slab regime such continuous weakening would have reduced the suction forces that kept the slab in a flat regime in the last 15 ma allowing the slab to easily roll back the continuous low resistivity region recorded beneath the volcanic front in central mexico might represent the evidence of slab dehydration and crust weakening over time

k ar age and geochemistry of the sw japan paleogene cauldron cluster implications for eocene oligocene thermo tectonic reactivation abstract systematic k ar dating and geochemical analyses of paleogene cauldrons in the sanin belt of sw japan have been made to explore the relationship between the timing of their formation and the paleogene subduction history of sw japan documented in the shimanto accretionary complex we also examine the magma sources and tectonics beneath the backarc region of sw japan at the eastern plate boundary of eurasia fifty eight new k ar ages and 19 previously reported radiometric age data show that the cauldrons formed during middle eocene to early oligocene time 43 30 ma following a period of magmatic hiatus from 52 to 43 ma the hiatus coincides with absence of an accretionary prism in the shimanto belt resumption of the magmatism that formed the cauldron cluster in the backarc was concurrent with voluminous influx of terrigenous detritus to the trench as a common tectono thermal event within a subduction system the cauldrons are composed of medium k calc alkaline basalts to rhyolites and their plutonic equivalents these rocks are characterized by lower concentrations of large ion lithophile elements lile including k2o ba rb th u and li lower la yb n ratios lower initial sr isotopic ratios 0 7037 0 7052 and higher end t values 0 5 to 3 5 relative to late cretaceous to early paleogene equivalents there are clear trends from enriched to depleted signatures with decreasing age from the late cretaceous to the paleogene the same isotopic shift is also confirmed in lower crust derived xenoliths and is interpreted as mobilization of pre existing enriched lithospheric mantle by upwelling depleted asthenosphere relatively elevated geothermal gradients are presumed to have prevailed over wide areas of the backarc and forearc of the sw japan arc trench system during the eocene to oligocene newly identified late eocene low silica adakites and high mg andesites in the sanin belt and early eocene a type granites in the sw korea peninsula probably formed due to upwelling of hot asthenosphere and subduction of a young plate the backarc region was an extensional tectonic setting and some paleogene rift basins and sanin belt cauldrons occur in linear arrays the eocene oligocene sanin se korea continental arc lies on the ne extension of the east china sea basin the initial stage of which probably formed by continental arc rifting this rifting may have been triggered by upwelling of hot asthenosphere into the wedge space created by rollback of the subducted slab in response to decreased convergence rate between the pacific and eurasian plates

seismicity and shape of the subducted nazca plate an updated compilation of earthquake locations and focal mechanism solutions from the international seismological centre and preliminary determination of earthquakes is the basis of a comprehensive study of the geometry of the wadati benioff zone beneath western south america the new data support previous mapping of a sharp flexure rather than a tear in the subducted nazca plate beneath southern peru and provide evidence for a similar flexure in the southward transition from nearly horizontal subduction to a slab with 30 dip at latitude 33 s in contrast the transition from 30 slab dip beneath bolivia to a nearly horizontal dip in the region between 28 s 32 s is more gradual occurring over several hundred kilometers of along strike distance between 20 s and 32 s this southward flattening corresponds to a broadening of a horizontal benchlike part of the subducted plate formed between 100 and 125 km depth the transition in continental tectonic style near 27 s 28 s from a wide volcanically active plateau to a narrow nonvolcanic cordillera appears not to be associated with the main slab flattening which begins to the north of these latitudes but with a more abrupt change in curvature of the subducted slab from convex upward to concave upward immediately below the plate boundary interface the concept of gaussian curvature is applied to slab bending to explain how subduction geometry is affected by the shape of the south american plate we hypothesize that the polarity of vertical curvature in the subducting slab is governed by the orientation of lateral curvature of the plate margin focal mechanism solutions for intermediate and deep earthquakes are grouped by geographic region and inverted for the orientation and relative magnitudes of the principal stresses results of the inversion indicate that downdip extension dominates in the slab above 350 km while downdip compression dominates at greater depths

structure and seismogenic properties of the mentawai segment of the sumatra subduction zone revealed by local earthquake traveltime tomography 1 on 12 september 2007 an mw8 4 earthquake occurred within the southern section of the mentawai segment of the sumatra subduction zone where the subduction thrust had previously ruptured in 1833 and 1797 traveltime data obtained from a temporary local seismic network deployed between december 2007 and october 2008 to record the aftershocks of the 2007 event was used to determine two dimensional 2 d and three dimensional 3 d velocity models of the mentawai segment the seismicity distribution reveals significant activity along the subduction interface and within two clusters in the overriding plate either side of the forearc basin the downgoing slab is clearly distinguished by a dipping region of highvp 8 0 km s which can be a traced to 50 km depth with an increased vp vs ratio 1 75 to 1 90 beneath the islands and the western side of the forearc basin suggesting hydrated oceanic crust above the slab a shallow continental moho of less than 30 km depth can be inferred suggesting that the intersection of the continental mantle with the subducting slab is much shallower than the downdip limit of the seismogenic zone despite localized serpentinization being present at the toe of the mantle wedge the outer arc islands are characterized by low vp 4 5 5 8 km s and high vp vs greater than 2 0 suggesting that they consist of fluid saturated sediments the very low rigidity of the outer forearc contributed to the slow rupture of the mw 7 7 mentawai tsunami earthquake on 25 october 2010

effect of aseismic ridge subduction on slab geometry and overriding plate deformation insights from analogue modeling abstract we present analogue models simulating the subduction of a buoyant ridge oriented perpendicularly or obliquely with respect to the trench beneath an advancing overriding plate the convergence velocity is imposed by lateral boundary conditions in this experimental set we analyze the three dimensional geometry of the slab the deformation and topography of the overriding plate experiments suggest that ridge subduction diminishes the dip of the slab eventually leading to the appearance of a horizontal slab segment in case boundary conditions impose a rapid convergence this result contrasts with that obtained in free subduction experiments in which ridge subduction diminishes the convergence velocity which in turn increases the dip of the slab beneath the ridge the slab dip decrease is accompanied by the indentation of the overriding plate by the ridge resulting in arc curvature experiments suggest that indentation is larger for small convergence velocity and large slab dip ridge subduction also uplifts the overriding plate uplift first occurs close to the trench fore arc area and is accompanied by the flexural subsidence of the overriding plate behind the uplifted area back arc subsidence the uplifted area migrates within the overriding plate interiors following the appearance of a horizontal slab segment these results are compared with natural examples of ridge subduction in the circum pacific area they explain why ridge subduction may have contrasted effects on the overriding plate dynamics depending on the global conditions that constrain the converging system

orientation and nature of active crustal stresses determined by electromagnetic measurements in the patagonian segment of the south america plate based on electromagnetic measurements we determined the current stress directions in the uppermost continental crust of patagonia between the active plate margin of the chilean pacific coast and the argentinean passive atlantic margin regional variations of the observed stress pattern are giving details onto the acting tectonic processes we distinguish five regional stress domains with different prevailing horizontal stress directions sh 1 southern coastal cordillera and longitudinal valley sh ssw nne 2 chiloe island sh sw ne 3 northern patagonian andes sh wsw ene 4 argentinean pampa and atlantic margin sh wnw ese and 5 southern patagonian andes sh wnw ese these stress regimes can be related to the geometry of the subducting nazca and antarctic plates to the transform fault between the south america and scotia plates and to passive margin processes along the atlantic coast absolute plate motion and rapid relative plate convergence control the subduction geometry and therefore the stress directions along the convergent margin of the south america plate and the structural style within and landward of the magmatic arc the knowledge of current local stress directions permits the characterisation of potential fault kinematics by in situ measuring of electromagnetic emissions from rocks we determined the maximum horizontal stress orientation in the uppermost crust using a new geophysical tool our investigations on the orientation of the stress regimes also allow conclusions about the causative forces of either tectonic or gravitational origin in this part of the south america plate

seismotectonics of the central segment of the indonesian arc abstract in this paper a revision of seismicity affecting the central segment of the indonesian island arc ranging between 110 and 126 e is presented using the areal and in depth distribution of seismic activity strain release maps and focal mechanisms lateral changes in the wadati benioff zone have been analyzed to determine possible boundaries between portions of lithosphere with different subduction geometries the seismicity pattern indicates that the sumbawa flores wetar sector shows different forms of behaviour with respect to the adjacent sectors these include driving mechanism inclination and continuity of the subducting slab and subduction features this area therefore seems to be isolated from the sunda and banda arcs by two principal boundaries these having a nearly n s trend in the bali region and with a nearly e w trend in the region ranging between wetar northern timor and tanibar the first boundary characterized by an absolute minimum of seismic activity at all ranges of depth has been interpreted in terms of subduction of the roo rise aseismic bathymetric ridge for the second boundary dividing a northwardly steeply dipping slab from an e w subducting slab dipping with an angle of 30 40 a tear in the upper part depth less than 300 km and a hinge fault system in the deepest part of the lithosphere have been proposed from the analysis of focal mechanisms of shallow earthquakes it was inferred that the central part of the indonesian arc is subject to a vortex shaped stress field centred on the savu basin in this model the compressive axes appear to rotate counterclockwise from sw to nne in the sumba sumbawa western flores region and clockwise from w to nnw in the timor eastern flores zone to interpret these features on the basis of seismological evidence a lateral discontinuity in the arc trench system close to sumba a collision between sumba and sumbawa and a rotation towards the north northeast of sumba have been suggested the proposed structural discontinuity trending nw se may represent a major transcurrent fault zone through which the australian continental lithosphere comes into contact with the indian oceanic lithosphere

intermediate and deep seismicity and lateral structure of subducted lithosphere in the circum pacific region in this paper we present a region by region review of the wadati benioff zone structure of most of the world s seismically active subduction zones focusing primarily on the intermediate and deep seismicity lateral changes in wadati benioff zone structure are common in every major subduction zone in this study we use these changes to define possible boundaries between portions or segments of lithosphere with differing subduction geometries although earthquake data seldom have the resolution to show conclusively whether these boundaries separate independent blocks of lithosphere the available data indicate that the active process at most of these segment boundaries is ductile deformation of the subducting plate rather than tearing we found the strongest evidence for the existence of tears where tears are geometrically necessary such as where a transform boundary terminates a trench as at the new hebrides tonga and south sandwich trenches weak evidence suggesting other tears does exist in some regions such as taiwan the japan izu bonin corner and the philippines the causes of these changes in structure are in most cases unclear only about 34 of the possible segment boundaries coincide with subducting bathymetric features some boundaries occur where there is apparent lateral strain caused by anomalous trench geometry we have designed a simple modeling procedure which incorporates published plate motion and the observed geometry of trenches and wadati benioff zones to estimate the lateral strain in subduction zones throughout the circum pacific region although no observed subduction zone has a perfectly strain free geometry there is a broad range of geometries for which the lateral strain is small indeed the observed geometry of most subduction zones involve relatively little lateral strain comparison with centroid moment tensor focal mechanisms indicates that in zones where the modeling predicted little lateral strain the mechanisms of intermediate and deep earthquakes show no effects of lateral stress and downdip stresses are clearly dominant in regions such as the mariana arc where the model predicts very large lateral extension lateral tension is very evident in the focal mechanisms in regions such as the hokkaido corner where the modeling predicts large compressional strains the plate appears to buckle and bending stresses parallel to the trench are evident in general this study finds that subducted lithosphere is remarkably cohesive and rigid and only rarely deforms by breaking or stretching

cenozoic deformation and tectonic style of the andes between 33 and 34 south latitude 1 the andes of argentina and chile between latitudes 33 and 34 s are composed from west to east of an oligocene to miocene volcanic arcs and the neogene east vergent aconcagua fold and thrust belt of the cordillera principal and the basement block faulted cordillera frontal a regional cross section suggests that shortening across the andes was achieved by thrusting along detachments at several levels in the crust while thin skinned deformation along newly formed thrusts occurred in mesozoic sequences of the eastern cordillera principal reactivation of preexisting jurassic and oligocene normal faults has resulted in additional hybrid thick and thin skinned structures in the western cordillera principal five major thrusting events are recognized in this part of the andes 1 early to middle miocene tectonic inversion of the extensional faults in the western cordillera principal 2 middle to late miocene development of the aconcagua fold and thrust belt 3 late miocene uplift of cordillera frontal 4 late miocene early pliocene out of sequence thrusting in the cordillera principal and 5 pliocene to present deformation of the foreland

constraints on pacific plate kinematics and dynamics with global positioning system measurements global positioning system gps receivers are capable of generating precise geodetic data that can yield important constraints on both plate kinematics and dynamics geodetic measurements can determine plate motion rates over time scales for which little data are currently available high quality geodetic data may also allow the investigation of plate driving forces because changes in the intraplate stress field can potentially be inferred from measurement of the resulting crustal strain a measurement program designed to investigate kinematic and dynamic aspects of plate tectonics should be concentrated in the pacific region this area contains the largest and fastest moving plate up to 17 cm year furthermore subduction zones which play a quantitatively important role in plate driving forces are largely restricted to the pacific region we summarize accuracy studies showing that for short 100 km baselines centimeter level accuracy can be expected using only mobile stations for longer baselines uncertainty in the orbit ephemeris of the gps satellites is a major error source performing simultaneous observations at widely separated 3000 km fiducial stations near the region of interest however should allow centimeter level accuracy for baselines up to several thousand kilometers in length this performance level is predicated upon the assumption that fiducial baselines are known a priori to the centimeter level for example from very long baseline interferometry and that corrections for tropospheric path delay are accurate to the centimeter level the fiducial network location is flexible limited mainly by the requirement for mutual satellite visibility

andesite and dacite genesis via contrasting processes the geology and geochemistry of el valle volcano panama the easternmost stratovolcano along the central american arc is el valle volcano panama several andesitic and dacitic lava flows which range in age 5 10 ma are termed the old group after a long period of quiescence approximately 3 4 ma volcanic activity resumed approximately 1 55 ma with the emplacement of dacitic domes and the deposition of dacitic pyroclastic flows 0 9 0 2 ma these are referred to as the young group all of the samples analyzed are calc alkaline andesites and dacites the mineralogy of the two groups is distinct two pyroxenes occur in the old group rocks but are commonly absent in the young group in contrast amphibole has been found only in the young group samples several disequilibrium features have been observed in the minerals e g oscillatory zoning within clinopyroxenes these disequilibrium textures appear to be more prevalent among the old as compared with the young group samples and are most likely the result of magma mixing assimilation and or polybaric crystallization mass balance fractionation models for major and trace elements were successful in relating samples from the old group but failed to show a relationship among the young group rocks or between the old and young group volcanics we believe that the old group volcanics were derived through differentiation processes from basaltic magmas generated within the mantlewedge the young group however does not appear to be related to more primitive magmas by differentiation the young group samples cannot be related by fractionation including realistic amounts of amphibole distinctive geochemical features of the young group including la yb ratios 15 yb 1 sr y 150 and y 6 suggest that these rocks were derived from the partial melting of the subducted lithosphere these characteristics can be explained by the partial melting of a source with residual garnet and amphibole dacitic material with the geochemical characteristics of subducted lithosphere melting is generated apparently only where relatively hot crust is subducted based on recent work the young dacite genesis at el valle volcano is related to the subduction of relatively hot lithosphere

the role of transcurrent intra arc tectonics in the configuration of a volcanic arc summary field investigations show that plio quaternary mexican and quaternary ecuadorian volcanic arcs are coeval with intra arc transcurrent or transtensional tectonics with motions along faults parallel or subparallel to the arc this relationship can also be observed in the quaternary kamchatka arc these arcs have an anomalously wide distribution of volcanic centres which cannot be explained simply by a low dip in the subducting lithosphere comparison with the intra arc tectonics and geometry of other quaternary arcs reveals that strike slip fault zones parallel to the arc favour magma rising in a broad surficial zone and may channel magma feeding mainly huge stratovolcanoes

reconstruction of subduction zone paleogeometries and quantification of upper plate material losses caused by tectonic erosion we propose in this paper an original method for quantifying the material losses at subduction zones mainly based on vertical movements of the active margins we first test the method using isostatic and elastic models applying it to 15 different well constrained convergent margins the main conclusion is that most of them are presently close to the isostatic equilibrium consequently we consider as a first hypothesis that such was also the case in the past and that erosion or accretion are natural mechanisms restoring the equilibrium in response to a perturbation then we apply this model to compute the geometry of the palcoslab knowing the palcotopography of the margin the comparison between present and past ocean margin lithospheric wedge geometries provides estimates of the removed material that constituted the initial wedge two cases were studied for palcoreconstructions japan and peru the results show that a considerable amount of material was removed during the neogene along both margins a detailed interpretation of the neogene history of the japan margin is proposed on the basis of our results

internal tectonic structure of the central american wadati benioff zone based on analysis of aftershock sequences 1 relocated engdahl et al 1998 global seismological data for 10 aftershock sequences were used to analyze the internal tectonic structure of the central american subduction zone the main shocks of several of these were the most destructive and often referenced earthquakes in the region e g the 1970 chiapas 1983 osa 1992 nicaragua 1999 quepos 2001 el salvador earthquakes the spatial analysis of aftershock foci distribution was performed in a rotated cartesian coordinate system x y z related to the wadati benioff zone and not in a standard coordinate system h are latitude longitude focal depth respectively available fault plane solutions were also transformed into the plane approximating the wadati benioff zone the spatial distribution of earthquakes in each aftershock sequence was modeled as either a plane fit using a least squares approximation or a volume fit with a minimum thickness rectangular box the analysis points to a quasi planar distribution of earthquake foci in all aftershock sequences manifesting the appurtenance of aftershocks to fracture zones geometrical parameters of fracture zones strike dip and dimensions hosting individual sequences were calculated and compared with the seafloor morphology of the cocos plate the smooth character of the seafloor correlates with the aftershock fracture zones oriented parallel to the trench and commonly subparallel to the subducting slab whereas subduction of the cocos ridge and seamounts around the quepos plateau coincides with steeply dipping fracture zones transformed focal mechanisms are almost exclusively 90 of normal character

mesozoic cenozoic denudation history of the patagonian andes southern chile and its correlation to different subduction processes fission track ft analysis is applied to assess the mesozoic and cenozoic thermal and denudational history of the patagonian andes between 44 and 51 s and the geologic and geomorphic response of late cenozoic subduction of the active chile rise mid oceanic spreading center on the overriding plate seventy two ft ages from 43 samples are presented zircon ft ages indicate fast post intrusion cooling of cretaceous parts of the patagonian batholith and previously unreported miocene magmatic rocks south of 48 s metamorphic basement rocks to the east of the batholith are constrained as having been deposited and metamorphosed in the early carboniferous and late permian apatite ft data reveal initiation of accelerated cooling and denudation at ca 30 ma at the western margin of southern continental south america followed by an up to 200 km eastward migration of the locus of maximum denudation that ceased at ca 12 8 ma at the position of the present day main topographic divide this migration is proposed to be related to either coeval eastward migration of the retroarc deformation the effects of subduction erosion in the overriding plate at the peru chile trench or less likely shallowing of the angle of subduction east of the divide 3 km of denudation has occurred since the late cretaceous enhanced denudation is interpreted to be the result of increased tectonic uplift driven by a large increase in convergence rates at ca 28 26 ma that triggered orographically enhanced precipitation on the west side of the patagonian andes allowing increased erosion by fluvial incision and mass transport processes the actual process of spreading center subduction had remarkably little influence on denudation in the upper plate and indeed coincides with a slowdown in denudation

cretaceous to early paleogene tectonic evolution of the northern central andes 0 18 s and its relations to geodynamics the tectonic phases subsidence patterns and location of the magmatic arc along the peruvian margin during late cretaceous and paleogene times are compared with the main geodynamic parameters absolute plate motion convergence velocity and direction age of the subducted slab and accretional events for the same periods during the late cretaceous long termed compression seems to be controlled by the absolute trenchward motion of the overriding plate and to a minor extent by the young age of the subducted lithosphere short lived contractional or extensional phases are mainly linked to the acceleration or deceleration of the convergence between the oceanic and continental plates and probably to changes in the convergence direction periods of high and low convergence rates coincide with increased and decreased subsidence rates along the margin respectively and seem to be independent of compressional tensional regimes during the paleogene in addition to the mentioned processes the eastward shift of the magmatic belt and subsidence in the forearc zone are interpreted as the result of the subduction erosion of the margin this process seems to be related to the decrease of the dip of the benioff zone expressed by the drastic widening of the magmatic arc

cenozoic plate kinematics subduction and magmatism south american andes revised plate reconstructions and enlarged isotopic age data sets strongly support three mechanisms for controlling the configuration and associated volcanism and tectonism of the andean subduction zone these include age of subducted plate convergence rate and interaction with hotspot traces from 50 to 25 ma in the andes the age of oceanic plate subducted progressively decreased apparently resulting in a9gradual decrease in subduction angle and consequent eastward shift in the loci of magmatic activity a major reorganization of plate motions at 25 ma resulted in an increased convergence rate normal to the chilean andes and was accompanied by subsequent broadening of the locus of magmatic activity into adjacent portions of bolivia and argentina convergence rates changed only slightly across the peruvian margin and no clear effect on subsequent magmatism is indicated interaction of the andean subduction zone with the juan fernandez and easter nazca hotspot traces apparently resulted in lower angle subduction and consequent gaps in magmatic activity as predicted continuations of the traces correspond with regions lacking dated volcanic rocks the current gaps in igneous activity correspond with contemporary subduction of the two traces and this association apparently occurred in the past as well

the tectonic expression slab pull at continental convergent boundaries examination of five thrust belt systems developed at continental subduction boundaries suggests that they comprise two distinct groups that display pronounced and systematic differences in structural style topographic elevation denudation metamorphism postcollisional convergence and foredeep basin geometry and facies the distinctive geological features developed within each thrust belt group appear to be causally linked to the relative rates of subduction and convergence via the magnitude of horizontal compressional stress transmitted across the subduction boundary at subduction boundaries where the rate of overall plate convergence is less than the rate of subduction termed here retreating subduction boundaries the transmission of horizontal compressive stress across the plate boundary is small and regional deformation of the overriding plate is by horizontal extension the tectonic expression of these retreating subduction boundaries includes topographically low mountains little erosion or denudation low grade to no metamorphism little to no involvement of crystalline basement in shortening little to no postcollisional convergence anomalously deep foredeep basins and a protracted history of flysch deposition within the adjacent foredeep basin analysis of deflection and gravity data across three retreating subduction boundaries apennine carpathian and hellenic systems shows that subduction is driven by gravitational forces acting on dense subducted slabs at depths between about 40 and 80 km carpathians 50 and 150 km apennines and 50 and 250 km hellenides the total mass anomalies represented by the slabs are approximately 3 1012 6 1012 and 12 1012 n m respectively the slabs are partially supported by flexural stresses transmitted through the subducted lithosphere to the foreland and partially supported by dynamic viscous stresses in the asthenosphere at subduction boundaries where the rate of overall plate convergence is greater than the rate of subduction termed here advancing subduction boundaries the transmission of horizontal compressive stress across the plate boundary is large and regional deformation of the overriding plate is by horizontal shortening the tectonic expression of these advancing subduction boundaries includes topographically high mountains antithetic thrust belts large amounts of erosion and denudation exposure of high grade metamorphic rocks at the surface extensive deformation of crystalline basement to midcrustal depths protracted postcollisional convergence tens of millions of years and a protracted history of molasse deposition within the adjacent foredeep basins analysis of gravity and deflection data across two advancing subduction boundaries developed within the continental lithosphere western to eastern and southern alps and himalayas shows that the thrust sheets have been translated for great distances over the foreland lithosphere relative to the point at which the subduction forces are applied thus obscuring any flexural and gravity signals from the subducted slab however it appears that far field stresses presumably related to global plate motions drive most of the convergent motion across these subduction boundaries the concept that orogenic belts formed above retreating subduction boundaries have recognizable tectonic signatures that differ from those of orogenic belts formed above advancing subduction boundaries suggests that it may be possible to interpret the plate boundary settings in which ancient orogenic belts evolved appendix b is available with entire article on microfiche order from the american geophysical union 2000 floridaavenue n w washington d c 20009 document t92 004 2 50 payment must accompany order

controls of subduction geometry location of magmatic arcs and tectonics of arc and back arc regions most variation in geometry and angle of inclination of subducted oceanic lithosphere is caused by four interdependent factors combinations of 1 rapid absolute upper plate motion toward the trench and active overriding of the subducted plate 2 rapid relative plate convergence and 3 subduction of intraplate island seamount chains aseismic ridges and oceanic plateaus anomalously low density oceanic lithosphere cause low angle subduction under conditions of low angle subduction the upper surface of the subducted plate is in contact with the base of the overlying plate the wedge of low density asthenosphere is replaced by subducted lithosphere and the width of the arc trench gap either is significantly increased or a magmatic arc is not developed within the overlying plate the fourth factor is age of the subducting lithosphere subduction of young lithosphere produces two opposing tendencies 1 low angle subduction and increased arc trench distance owing to its low density and 2 decreased arc trench distance owing to its higher temperature two factors of secondary importance contribute to variation in subduction zone geometry and arc trench distance accretion of sediment in trenches depresses the upper portion of the subducting oceanic plate and causes the trench axis to migrate seaward prolonged subduction thickens the upper plate depresses the isotherms in the subducted plate and may create a broader arc both factors increase the arc trench gap the four primary factors also control development of other tectonic elements such as regional subsidence for example the amazon basin and a portion of the cretaceous interior seaway of western united states intra arc extension for example the basin and range province foreland fold and thrust belts and laramide style tectonics

the development of slabs in the upper mantle insights from numerical and laboratory experiments we have performed numerical and laboratory experiments to model subduction of oceanic lithosphere in the upper mantle from its beginnings as a gravitational instability to the fully developed slab a two dimensional finite element code is applied to model newtonian creep in the numerical experiments scaled analog media are used in the laboratory a sand mixture models the brittle crust silicone putty simulates creep in the lower crust and mantle lithosphere and glucose syrup is the asthenosphere analog both model approaches show similar results and reproduce first order observations of the subduction process in nature based on density and viscosity heterogeneities in a stokes flow model subduction nucleates slowly and a pronounced slab forms only when the viscosity contrast between oceanic plate and mantle is below a threshold we find that the subduction velocity and angle are time dependent and increase roughly exponentially over tens of millions of years before the slab reaches the 670 km discontinuity the style of subduction is controlled by the prescribed velocity of convergence the density contrast between the plates and the viscosity contrast between the oceanic plate and the mantle these factors can be combined in the buoyancy number f which expresses the ratio between driving slab pull and resisting viscous dissipation in the oceanic plate variations in f control the stress in the plates the speed and the dip of subduction and the rate of trench retreat reproducing the contrasting styles of subduction observed in nature the subduction rate is strongly influenced by the work of bending the lithosphere as it subducts

paleoproterozoic positive δ13ccarb excursion in the northeastern sino korean craton evidence of the lomagundi event abstract the 2 33 2 06 ga positive 13 c carb excursion associated with environmental change and the breakup of the kenorland or superia supercontinent is called the lomagundi or jatulian event or great oxidation event and has been reported in many early precambrian cratons but not yet in the sino korean craton the guanmenshan formation of the liaohe group occurs in the northeastern part of the sino korean craton 13 c carb and 18 o values in 42 samples from this formation range from 3 5 5 9 v pdb and 15 4 24 8 v smow respectively showing a clear positive 13 c carb excursion that characterizes the lomagundi event thirty five of the 42 samples with less hydrothermal alteration have higher 13 c carb and 18 o carb values than the other 7 samples obviously affected by fluid flow confirming that it was fluid flow that reduced the 13 c carb and 18 o carb values this positive 13 c carb excursion places deposition of the guanmenshan formation within the age range of 2 33 2 06 ga

palaeozoic arc growth deformation and migration across the lachlan fold belt southeastern australia abstract the palaeozoic lachlan fold belt lfb is a low pressure metamorphic belt characterized by greenschist facies rocks and infolded volcanic sequences which indicate that the belt has not undergone substantial uplift at any stage porphyroblast matrix relationships in one of the highest grade metamorphic zones of the lfb the omeo complex indicate that the thermal peak was reached before or was synchronous with the earliest deformation similar to that in high grade metamorphic zones of the northern arunta inlier central australia in both areas the metamorphism is of low p high t type centred on major granite intrusions as regional aureoles and characterized by anti clockwise p t t paths however in the outer lower grade parts of the aureoles peak metamorphic conditions post date the earliest deformation this diachronous relationship indicates that 1 outward migration of a thermal front was centred on major granitoid intrusions batholiths 2 deformation propagated outward more rapidly than migration of the thermal front and 3 plastic and ductile deformation began only after heating around the batholiths thus orogeny was initiated after and localized by heat focusing in the mid crust associated with batholith emplacement therefore the early deformation is a result of thermal softening it is typically subhorizontal ductile shear at mid crustal levels but is characterized by upright to inclined folds at upper crustal levels in the latter environment the typically discordant contact aureole plutons are often interpreted as post tectonic granitoids but they were part of the ongoing mid crustal thermal perturbation that induced regional greenschist facies metamorphism which overprinted the early formed upright folds as it migrated to upper crustal levels orogeny migrated generally eastward through the lfb it began with the early silurian benambran orogeny centred on the wagga metamorphic belt wmb and terminated with the carboniferous kanimblan orogeny which was most intense in the hill end trough eastward migration of orogeny across central victoria part of the west lfb began in the ballarat bendigo zone during the early devonian and terminated in the melbourne zone in the middle devonian orogeny was centred on the meridional batholiths of the lfb each of which are considered to represent the transitory axis of an ancient magmatic arc stepwise but generally eastward arc migration caused the eastward migration of orogeny a tectonic model based on the modern southwest pacific arc system can be applied to the lfb following cambrian intra oceanic arc growth associated with west dipping subduction an ordovician marginal sea developed partly on stretched cambrian crust and was flooded by an extensive turbiditic wedge closure of the eastern part the east lfb began in the early silurian after a magmatic arc developed over the wmb possibly associated with an east dipping subduction zone after re establishment of west dipping subduction in the middle silurian absolute motion retreat of the upper australian plate caused by oblique plate convergence resulted in dextral transtension and the generation of siluro devonian inter arc and back arc basins transient compression resulted in deformation along the magmatic arc and induced formation of a new outboard eastward arc system periodic compression within the evolving arc and back arc system resulted in a series of east younging east verging linear siluro devonian fold and thrust belts which were localized adjacent to elongate batholiths that represent the relict arc systems local westward jumps in deformation and plutonism occurred the most significant is the 600 km westward jump into the west lfb caused by closure of the remainder of the passive ordovician back arc basin melbourne and bendigo ballarat zones west dipping subduction beneath the basin produced a low pressure fold thrust belt similar to the east lfb though the granitoids have a more primitive isotopic character indicating closure of thin dominantly oceanic crust in the early to middle devonian localized deformation induced by thermal softening of the crust within magmatic arcs that have developed on earlier passive margins is likely to produce anticlockwise p 7 t paths and may be a typical feature of southwest pacific style tectonism

geochronology petrogenesis and tectonic significance of peraluminous granites from the chinese altai nw china abstract the paleozoic granites in the chinese altai are important for the study of tectonic evolution and crustal growth in the central asian orogenic belt caob four representative peraluminous granitic intrusions were selected for systematic studies of zircon u pb and hf isotopic compositions and whole rock geochemical and nd sr isotopic analyses these rocks have high asi alumina saturation index al 2 o 3 cao na 2 o k 2 o 1 01 1 46 molecular ratios with 0 6 5 6 wt of normative corundum and are characterized by moderately negative eu anomalies eu eu 0 38 0 98 and strong depletion in ba nb and sr elements our data suggest that these intrusions were emplaced from 419 to 393 ma consistent with a period of intensive magmatic activities and high temperature metamorphism in the chinese altai while in situ zircon hf isotopic analyses for these granites give predominantly positive e hf t values 0 8 to 12 8 a few inherited zircons yield negative e hf t values from 12 5 to 1 53 the u pb age and hf isotopic data of these inherited zircons are similar to that of the widespread metasediments in addition the peraluminous granitic rocks have near zero or negative e nd t values 3 3 to 0 5 and relatively high initial 87 sr 86 sr ratios 0 7079 0 7266 distinct from those of the i type granites in the study region but similar to the early paleozoic habahe sediments these isotopic compositions suggest that the newly accreted metasediments of habahe group may be the major source rock of the peraluminous granites the geochemical compositions indicate that their precursor magmas were derived from a relatively shallow crustal level p 5 kbar and zircon saturation temperatures suggest that these granitic intrusions were emplaced at 672 861 c the peraluminous granitic magmas were generated by dehydration melting of newly accreted materials which were possibly brought to at least middle crustal depth by subduction related processes in an active margin and were subsequently molten by strikingly high ambient temperature probably caused by upwelling of the hot asthenosphere associated with ridge subduction in the paleozoic

dynamical and thermal structure of a subduction zone influence of slab geometry on the convective state of the earth s upper mantle preliminary results received 29 march 1996 accepted 3 july 1996 abstract the aim of this study is to investigate the influences of plate kinematics and of the downgoing slab s geometry on internal dynamics and thermal structure in the vicinity of a subduction zone a systematic study depending on these parameters has been carried out based on a two dimensional finite differences code in a convective domain representing the upper mantle plate displacements are prescribed by kinematic boundary conditions in particular the subducting plate moves and dips at a constant velocity under the slower overriding one typical subduction dips 30 45 and 70 have been studied for a realistic rayleigh number of about l05 and for plate velocities around a coupling velocity of 2 3 cm year 1 this limit value was inferred from the balance between surface forces and buoyancy forces for any subduction dip and plate velocities lower than this coupling velocity plate kinematics variations induce significant effects on thermal and mineralogical patterns at greater plate velocities these patterns stabilize for plate velocities of the order of the coupling velocity a subduction dip decrease of 40 involves a heating of the downgoing plate of 100degrees and a depth increase of 20 km of the olivine spinel phase transition in this plate in the particular case of a non moving overriding plate the balance between this plate and underlying internal dynamics requires a moving plate velocity lower than 2 3 cm year j for a lower dip moreover a shallow subduction enhances a lateral large scale convection under the upper plate even if this plate is motionless in contrast to the case of a steep subduction this is consistent with observed back arc extension at the surface which is often associated with steep subductions our model suggests a spreading zone of between 500 and 1000km in the upper plate for a subduction velocity lower than 2 3 cm year i and a rayleigh number of about l05 1 introduction the earth s mantle undergoes a rayleigh benard thermal convection coupled to oceanic lithospheric plate movements rigid and moving lithospheric corresponding author plates form the thermal and mechanical limit bound aries of the earth s mantle convection cells one of the concerns of numerical modelling is to account for a real dynamic convection model coupled to plate tectonics to force the movement of rigid plates one may either model a dynamic plate by prescribing a lower viscosity at plate margins or impose the horizontal surface velocity davies 1989 0031 9201 97 17 00 copyright 1997 elsevier science b v all rights reserved ph s0031 9201 96 03209 8

deformed continental arc sequences in the south tianshan new constraints on the early paleozoic accretionary tectonics of the central asian orogenic belt abstract the tectonic affinity and structural evolution of the south tianshan are key issues to understand the paleozoic accretionary or collisional orogenesis in the southwestern central asian orogenic belt but remain controversial geological and structural investigations on the early paleozoic low grade metamorphic volcanic sedimentary sequences around the bayinbuluk basin demonstrate that they were prevalently reworked by north vergent ductile deformation and overlain unconformably by devonian to carboniferous carbonates zircon u pb dating indicates that the volcanic rocks erupted during the middle ordovician to early devonian time 460 410 ma and contain neoarchean to neoproterozoic zircon xenocrysts 2 75 ga 2 46 ga 1 60 ga 997 963 ma and 827 ma their whole rock elemental sr nd isotopic end t 5 0 to 3 6 and zircon lu hf isotopic ehf t 13 5 to 11 6 features resemble those of andean type continental arc volcanic rocks the coexisting tuffaceous sedimentary rocks indicate a depositional environment amidst a continental volcanic arc at 421 to 404 ma the newly obtained and available geochronological geochemical and structural data suggest that 1 a wide active continental margin covering the south and central tianshan and northern tarim was formed due to early mid paleozoic southward subduction of the paleo tianshan ocean beneath the tarim craton and 2 a mid paleozoic transition from advancing to retreating accretionary orogeny was followed by opening of the south tianshan back arc basins and break up of the northern tarim continental margin the south tianshan was subsequently subjected to regional deformation during the successive diachronous closure of the major and back arc oceanic basins and locally reworked by the ductile strike slip faulting 307 255 ma

geodynamic evolution of the northern and central andes during early to middle mesozoic times a tethyan model the early to middle mesozoic sedimentary and magmatic history of the northern and central andes i e the nne trending colombian ecuadorian segment and the nw trending peruvian segment exhibits a spatially contrasted evolution involving several successive tectonic and geodynamic settings the late triassic to late liassic period began with a widespread marine transgression from the latest triassic the colombian marine shelf was progressively destroyed by southward propagating extensional tectonic activity but marine sedimentation continued in peru during this time no significant magmatic activity is recorded except in the emerging colombian area this period is interpreted as the result of rifting of a tethyan oceanic arm which separated the colombian and palaeo mexican margins from latest liassic an important calc alkaline magmatic arc developed along the emergent colombian segment further south the north peruvian shelf probably emerged but marine sedimentation continued in southern peru during middle and early late jurassic times the colombian segment was characterized by important magmatic activity and by coarse clastic continental sedimentation along the peruvian segment a turbiditic trough emergent areas and continental basins were created and the scarcity of calc alkaline magmatism suggests that only very local subduction took place this period is regarded as one of the southeastward subduction beneath the colombian segment and of tethyan oceanic crust originating in the newly formed colombian oceanic arm this pattern would have induced a chiefly left lateral transform motion along the peruvian segment by kimmeridgian tithonian times the palaeogeographic framework had drastically changed along the colombian segment magmatic activity ceased and continental accretions occurred along dextral strike slip sutures in peru tectonic activity was recorded by the creation of a new turbiditic trough and by the resumption of detrital sedimentation in the coastal area arc related volcanism indicated that subduction took place beneath this segment this geodynamic change is interpreted as the result of a sharp decrease in the spreading activity of the tethyan ridges and replacement by pacific spreading centres inducing a roughly northeastward convergence direction

cenozoic evolution of the savu basin indonesia forearc basin response to arc continent collision abstract a seismo stratigraphic analysis of the savu basin relates the evolution of a complexly evolved forearc basin to the late miocene recent collision with australia initial early miocene subsidence of the outer forearc basin and the subsequent development of a middle late miocene volcanic proto arc located in the south of the basin both resulted from the late oligocene early miocene eastward propagation of the java timor trench block faulting of the relatively thick continental forearc basement probably determined the rather irregular and scattered pattern of middle late miocene volcanic centres just north of the accretionary prism the forearc basement is flexed down and reshaped into a trenchward dipping backstop that facilitates backthrusting of the accretionary prism the southern forearc basement probably acted as a barrier against compression the thickness of this continental basement is suggested to be a critical factor in the response of the forearc to continent collision owing to its strength the savu basin did not react to collision by the development of northward verging inter arc thrusts in contrast it responded to the underthrusting of continental crust by reactivation of the basement ridges this resulted in the differentiation of the forearc basin into a presently extinct and uplifted miocene south savu basin and a pliocene recent active north savu basin the simultaneous late miocene recent uplift of large segments of the outer forearc and subsidence of the north savu basin is interpreted to result from tectonic underplating of continental material and uplift caused by the buoyancy of the subducted continental crust the relatively recent decoupling of this crust from the oceanic slab and the subsequent isostatic rebound may superimpose the initial effects of post collision extension on arc continent collision

spherical shell tectonics effects of sphericity and inextensibility on the geometry of the descending lithosphere the shape of a deep seismic zone is thought to represent that of the descending slab of lithosphere the lithosphere before subduction is a spherical shell and the shape of the descending slab is the result of the deformation of the spherical lithosphere at the subduction zone upon bending a spherical shell often deforms in a very different way from a simple plate we examine whether the actual shape of the descending slab can be explained by a simple bending of an inextensible spherical shell which shows little surface deformation under moderate stress this examination is made region by region for most of the subduction zones in the world by means of an analogue method the lithosphere is simulated by an inextensible spherical shell made of polyvinyl chloride resin the wadati benioff zone is shaped by plaster by referring to the reliable hypocentral data selected from the international seismological centre isc bulletins the spherical shell is forced to fit the miniature of the wadati benioff zone fitting is first attempted only by bending if a good fit is not attainable and if a discontinuity or gap in seismic activity is observed in the relevant region the spherical shell is torn along this discontinuity or gap and the goodness of fit is reexamined the results of the analysis are summarized as follows 1 the shape of the wadati benioff zone can be simulated largely by a simple bending of a spherical shell without surface deformation 2 in almost all of the regions of poor fit with bending a good fit can be achieved by tearing the spherical shell along the trace of low seismicity the sphericity of the lithosphere and the inextensibility upon deformation are the two essential factors in controlling the slab shape this means that the lateral constraint is most important for understanding the geometry of the downgoing slab of lithosphere and the stress state within it further several problems related to the deformational characteristics of the spherical lithosphere are also reviewed and discussed in connection with subduction tectonics

mesozoic episodic magmatism in south korea and its tectonic implication 1 ne trending granitoid batholiths in south korea are part of the huge phanerozoic granitoid belt in east asia the crystallization ages of granitoid plutons show three episodes of mesozoic magmatism in south korea triassic 248 210 ma jurassic 197 158 ma and late cretaceous early tertiary 110 50 ma with a late jurassic early cretaceous 50 m y hiatus the granitoid ages and kinematic analyses of the mesozoic structures in south korea suggest that the middle jurassic 180 170 ma dextral strike slip ductile shear zones formed by the northwest orthogonal subduction of the izanagi plate whereas the cretaceous sinistral faults formed by the northward oblique subduction of the same plate which generally agree with proposed reconstruction of plate motions albeit uncertain around east asia among the models proposed to explain magmatic lull in a continental arc shallow angle subduction associated with high convergence rate and subduction of oceanic plateaus might be the probable cause for the magmatic hiatus in south korea and japan the huge width 500 km of the east asian granitoid belt may be explained by oceanward retreat of subduction zone in the cretaceous in association with change in the movement direction of the subducting plate

tibet uplift and erosion abstract the 5 km high tibetan plateau is an outstanding topographic feature on the earth today its horziontal extent elevation and location cause significant effects on modern atmospheric circulation and climate so the history of uplift of the surface of tibet is linked to cenozoic climate changes at local regional and perhaps global scales geological and geophysical studies of the plateau are contributing data on the present and past deformation of the tibetan lithosphere that has formed the plateau primarily during the cenozoic the principal of isostasy then can be used to estimate the elevation history of the surface for a given deformation history different parts of tibet probably had different uplift histories but presently available data are not sufficient for distinguishing many contrasts in one scenario cretaceous and early cenozoic north south distributed shortening of tibetan crust and mantle lithosphere probably caused significant uplift of the surface relative to sea level to perhaps half of the present elevation by the early miocene thinning of the high density mantle portion of the lithosphere during the miocene may then have allowed the thick tibetan crust to rise close to its present elevation perhaps higher before 8 ma since then slow east west extension of tibet probably reduced the crustal thickness slightly and may have caused the elevation of the plateau to decrease during the late cenozoic erosion of tibet unlike narrow mountain belts has been unable to match the uplift of a broad plateau orographic precipitation and efficient river networks concentrate erosion on the edges while the interior is protected from significant erosion despite its lofty elevation the southern edge of the plateau the himalaya has suffered a minimum of 25 km of denudation since the miocene while central tibet shows little or no sign of major erosion since that time the gangdese arc in southern tibet was rapidly eroded during the mid miocene when 4 km of rock apparently were removed from the surface as shown by mineral cooling ages this pulse of erosion was probably caused by a combination of local thrust system movement and changes in base level and precipitation due to relative elevation changes between the gangdese and the himalaya to the south the modern long wavelength flatness of tibet is unlikely to have been caused by erosion and indicates viscous flow at some level of the lithosphere has been acting to level the plateau surface

geological constraints on the kinematics and fault plane solution of the quiches fault zone reactivated during the 10 november 1946 ancash earthquake northern peru the ancash earthquake of 10 november 1946 mw 6 8 produced surface faulting along the quiches normal slip fault field studies reported discontinuous nw se to nnw sse striking seismic scarps that displace quaternary glacial moraines and outwash as well as mesozoic sedimentary rocks structural analyses show that the 1946 normal displacement was accompanied by a small left lateral slip component morphological analysis of the fault scarp and observations in a trench excavated across the scarp show evidence of a single normal fault displacement before the 1946 seismic event that offset late pleistocene glacial moraines this indicates an approximate recurrence interval of more than 13 000 1 000 years for the major events suggesting a rate of net vertical displacement on the fault of less than 0 25 mm yr kinematic analysis along the quiches fault in comparison with other regional studies suggests that the state of stress in the high andes is regionally and locally homogeneous with only minor changes in the stress directions adjacent to major structures such as the maranon thrust and fold belt the minimum horizontal stress 3 axis is oriented n10 25 e the intermediate horizontal stress 2 axis is oriented n100 115 e and the maximum stress 1 axis is vertical we suggest that the 60 sw dipping quiches surface fault plane is connected at depth to a low angle dipping fault plane this seismic fault is a reactivated thrust fault of the tertiary maranon thrust and fold belt

seismicity and tectonics of the subducted cocos plate we have examined teleseismic earthquake locations reported by the international seismological centre isc for the middle america region and selected 220 as the most reliable these hypocenters and other data are used to delineate the deep structure of the subducted cocos plate the results indicate that the subducted plate consists of three major segments segment i extends from the panama fracture zone to the nicoya peninsula the structure of this segment is poorly defined segment ii is the largest and best defined segment this segment consists of two parts iia and iib part iia extends from the nicoya peninsula to western guatemala and is very well defined and continuous in structure its strike follows the curvature of the trench and dips at about 60 part iib extends from western guatemala to orizaba mexico the dip of this part of the segment decreases slightly toward the northwest and its strike is more northward than that of the trench segment iii extends from orizaba to the rivera fracture zone and is not well defined due to a lack of earthquake activity beneath about 100 km its orientation differs markedly from segment ii and strikes somewhat more westward than the trench between parts iia and iib of segment ii the subducted plate seems to be continuous bending smoothly to accommodate the change in geometry local network data from costa rica suggest there may be a tear between segments i and ii between segments ii and iii there is a gap in the hypocenters which makes it difficult to define the boundary the change in geometry between these two segments indicates that there may be a tear and two strike slip focal mechanisms in the region support this conclusion we find no convincing evidence supporting the existence of segments smaller than the three described above if there is smaller scale segmentation in the shallow part of the subducting plate the plate must still maintain enough continuity to appear continuous at greater depths there is no evidence for any major tear in the subducted plate associated directly with either the tehuantepec ridge or the orozco fracture zone the shallow subduction at the northwestern end of segment ii may be related to the bouyancy of the tehuantepec ridge the cocos ridge is probably directly responsible for the change in geometry between segments i and ii and may even be slowing or stopping subduction in segment i the structure of the subducted plate in segment ii and the changes in the character of volcanism along the arc can be related to the relative motion of the north american and caribbean plates the present geometry of part iib of segment ii is more consistent with the probable configuration of the trench about 7 ma ago than with the present configuration indicating that the north american plate is overriding the subduction zone appendices 2 3 and 4 are available with entire article on microfiche order from american geophysical union 2000 florida avenue n w washington dc 20009 document b84 009 2 50

evidence for the subduction and underplating of an oceanic plateau beneath the south peruvian margin during the late cretaceous structural implications abstract uring late cretaceous times santonian to maastrichtian the southern peruvian and northernmost chilean andes show a series of tectonic magmatic and sedimentological features which appear to be specific to this area when compared with northern and southern parts of the andes a gap in the magmatic activity in the coastal range and the western cordillera between 84 and 70 ma huge overthrusting faults involving the precambrian basement in the arequipa region and the syntectonic filling of the cuzco and sicuani foreland continental basins these particular features may be regarded as consequences of the subduction and underplating of a paleo oceanic plateau the mollendo ridge which would have occurred in this area between 85 and 70 ma b p with an almost n s direction of convergence between the farallon and south american plates this inferred mollendo ridge would have been an eastern equivalent of the present day manihiki plateau of the western pacific ocean the proposed model provides new elements for the interpretation of the subsequent structural and magmatic evolution of the central andes it permits the explanation of the particular location and the petrochemical and metallogenetic features of the andahuaylas yauri batholith the segment of abnormally thick crust created by the underplating of the mollendo ridge and by overthrusting tectonics will act as an accumulation point for further cenozoic deformation along the andean margin then the subduction of the mollendo ridge would mark the outset of the development of the oroclinal bending of the central andes the arica elbow and of the associated altiplano and its anomalous thick crust

tectonic and magmatic responses to the subduction of high bathymetric relief subduction of high bathymetric relief such as aseismic ridges and magmatic plateaus is considered to be responsible for dramatic changes in the dynamics and kinematics of the subduction zone for example the buoyancy of high bathymetric relief is thought to flatten the dip of the subducting slab modifying the structural and magmatic evolution of the overriding plate and terminating arc volcanism in addition the effect of ridge subduction in retreating plate boundaries can inhibit subduction rollback a process that could locally pin the subduction hinge and lead to the development of cusps and slab tearing here we discuss the tectonic response to subduction of high bathymetric relief using examples from the circum pacific subduction systems we demonstrate that flattening of the subduction dip angle is only significant in the eastern pacific where the average slab dip angle is relatively shallow in the western pacific in contrast the average subduction dip angle is steeper and there is no significant flattening of the dip angle in areas of ridge subduction subduction of high bathymetric relief in the circum pacific is commonly associated with reduced arc volcanism and in many cases the area of ridge subduction coincides with a volcanic gap in the overriding plate ridge subduction is associated with pronounced changes in the style of deformation involving uplift reactivation of basement thrusts development of orogen perpendicular tear faults and block rotations leading to oroclinal bending the discussed characteristic patterns associated with ridge subduction provide important guidelines for reconstructing past plate tectonic processes and could help constraining the geodynamics of ancient subduction systems

neocene and quaternary extent and geometry of the subducted pacific plate beneath north island new zealand implications for kaikoura tectonics abstract the extent and geometry of the obliquely subduced oceanic pacific plate beneath north island new zealand for five million year intervals through the mid miocene to quaternary are presented in a series of maps and cross sections these show that the subducted plate progressively increased its extent from ne to sw beneath the north island and in the more northern regions where it was first emplaced concomitantly increased its dip from 10 to 50 the changing extent and geometry of the subducted slab has been established from the age pattern of orogenic andesites and from the geochemical k 2 o h parameter of depth of magma generation the radiometric dates show a migration of the volcanic front back towards the trench at an average rate of 20 km my the trenchward migration is explained by a model of increasing slab dip which is corroborated by the k 2 o data calibrated against the presently active arc taupo volcanic zone with the exception of northern coromandel peninsula the andesitic magmas were generated at 85 100 km depth the interpretation of the dates adopted here indicates that the subducted slab originated at the ne sw trending kermadec hikurangi trench and implies a different and much simpler evolution of the australia pacific plate boundary in the vicinity of north island than other recent models subduction geometry has been found elsewhere to be a principal influence upon the state of stress and deformational style in an over riding plate the possibility is explored that the timing nature and pattern of the neogene to quaternary kaikoura orogeny in north island is due to this influence apart from the effect of oblique subduction in eastern north island there is an accord between the onset of deformation and the emplacement sequence of the shallow slab beneath north island and between the change in subduction geometry and a progressive north to south change in northern north island from compression to extension

the ordovician sierras pampeanas puna basin connection basement thinning and basin formation in the proto andean back arc article i nfo the ordovician sierras pampeanas located in a continental back arc position at the proto andean margin of southwest gondwana experienced substantial mantle heat transfer during the ordovician famatina orogeny converting neoproterozoic and early cambrian metasediments to migmatites and granites the high grade metamorphic basement underwent intense extensional shearing during the early and middle ordovician contemporaneously up to 7000 m marine sediments were deposited in extensional back arc basins covering the pre ordovician basement extensional ordovician tectonics were more effective in mid and lower crustal migmatites than in higher levels of the crust at a depth of about 13 km the separating boundary between low strain solid upper and high strain lower migmatitic crust evolved to an intra crustal detachment the detachment zone varies in thickness but does not exceed about 500 m the formation of anatectic melt at the metamorphic peak and the resulting drop in shear strength initiated extensional tectonics which continued along localized ductile shear zones until the migmatitic crust cooled to amphibolite facies p t conditions p t d t data in combination with field evidence suggest significant ca 52 crustal thinning below the detachment corresponding to a thinning factor of 2 1 ductile thinning of the upper crust is estimated to be less than that of the lower crust and might range between 25 and 44 constituting total crustal thinning factors of 1 7 2 0 while the migmatites experienced retrograde decompression during the ordovician rocks along and above the detachment show isobaric cooling this suggests that the magnitude of upper crustal extension controls the amount of space created for sediments deposited at the surface upper crustal extension and thinning is compensated by newly deposited sediments maintaining constant pressure at detachment level thinning of the migmatitic lower crust is compensated by elevation of the crust mantle boundary the degree of mechanical coupling between migmatitic lower and solid upper crust across the detachment zone is the main factor controlling upper crustal extension basin formation and sediment thickness in the back arc basin the initiation of crustal extension in the back arc however crucially depends on the presence of anatectic melt in the middle and lower crust consumption of melt and cooling of the lower crust correlate with decreasing deposition rates in the sedimentary basins and decreasing rates of crustal extension

tuffs tectonism and glacially related sea level changes carboniferous permian southern namibia abstract the carboniferous permian glacigene dwyka group in southern africa has been widely regarded as amagmatic only in the higher parts of the permian and subsequent younger strata has evidence of magmatism been recognised particularly in the form of pyroclastic fallout ash beds the source of these tuffs has been suggested to be patagonia and west antarctica forming part of an extensive volcanic arc which lay some 1500 to 2000 km to the south and west in their pre break up gondwana positions rhyolitic dacitic and basaltic andesitic tuff beds within the glaciomarine sediments of the dwyka group in namibia reveal new evidence for an early onset of proximal bimodal volcanic activity in southern africa contemporaneous tectonism is recorded by type 1 unconformities and systematic thickness changes across nw se trending extensional normal faults we suggest that this tectonomagmatic period marks initial extensional events in southern namibia and the carboniferous permian volcanic trend appears to coincide with the position of the eventual namibian continental margin significantly the stratigraphic positions of the tuff beds show a distinct association with enclosing marine transgressive depositional sequences most of them coincide with sharply defined flooding surfaces within relatively thick shaly successions reflecting their multiple transgressive nature widening the study context to the latest early permian in southern africa the main karoo basin of southern africa exhibits a similar coincidence of marine transgressions and extrusion of volcanics the relationship between tuff beds and transgressive depositional sequences is not only the result of enhanced preservational potential during transgression but of a genetic coupling between magmatism extensional tectonics and basin subsidence the latter accelerating a rise in relative sea level this interaction is significant not only for understanding potential controls on relative sea level change but also for understanding the early geodynamic evolution of the southern south atlantic continental margin

detailed three dimensional shear wave velocity structure of the northwestern united states from rayleigh wave tomography abstract since the mid miocene the northwestern united states has experienced extensive flood basalt volcanism followed by the formation of two time progressive tracks of silicic volcanism the yellowstone snake river plains ysrp and the high lava plains hlp the ysrp track progresses towards the northeast parallel to north american plate motion and has therefore often been attributed to a deep mantle plume source however the hlp track progresses to the northwest over the same time frame in a direction not consistent with any regional plate motion the causes of the mid miocene flood basalts and the tracks of the ysrp and hlp are a matter of ongoing debate we present results of rayleigh wave phase velocity inversions and inversions for 3 d shear wave velocity structure of the northwestern united states using data collected from the high lava plains seismic experiment and the earthscope usarray transportable array ta the large number of stations used in these inversions allows us to show an unprecedented level of detail in the seismic velocity structures of this tectonically complex area our velocity images indicate that low s wave velocities in the uppermost mantle do not well match the track of hlp volcanism while at the surface the newberry caldera appears to anchor the nw end of the hlp hotspot track the seismic results show that it lies in a separate north south trending low velocity band just east of the cascades that is distinct from the main hlp trace the ultra low s wave velocities beneath the ysrp track extend locally to at least 175 km depth and are by far the most prominent seismic anomalies in the region along axis the ysrp hotspot track is characterized by a discrete low velocity channel in the upper mantle that shallows narrows and intensifies to the northeast but then deepens rapidly to the north beneath yellowstone the shallowing of the low velocity anomaly to the northeast is consistent with a stationary heat source beneath a moving plate coupled with basin and range extension and heating

intraoceanic subduction of heterogeneous oceanic lithosphere in narrow basins 2d numerical modeling abstract we provide 2d thermomechanical numerical models to constrain the subduction process of narrow 600 km wide oceanic basins resulting in formation and exhumation of serpentinite bearing high pressure hp and ultrahigh pressure uhp complexes presently exposed as dismembered massifs in orogenic belts we simulated subduction of heterogeneous i e non layered oceanic lithosphere testing different serpentinite rheologies and varying other major parameters initial slab dip convergence rate age and structure of oceanic lithosphere distance subduction zone continental margin weak zone geometry we show that serpentinite rheology is the most important parameter affecting both the geometry and behavior of serpentinite channels that form in the mantle wedge as the result of upward migration of slab fluids a strong serpentinite power law rheology is responsible for the continuous accretion and underthrusting of slab slices at the base of mantle wedge to form planar or wedge shape serpentinitic channels the great part of sediments is buried only during continental crust subduction but is never exhumed moreover slab slices do not mix with mantle wedge rocks exhumation of mafic ultramafic rocks occurs from pressure not exceeding 15 kbar on the other hand a weak newtonian behavior of serpentine best describes subduction and exhumation dynamics this rheology produces channels that usually widen to depth slab and sediment slices are easily incorporated in mantle wedge serpentinite during ongoing convergence in this case slab and mantle wedge serpentinite are mixed together and can be closely associated in the channel the models running a newtonian rheology highlight that slices of the oceanic overriding plate can be eroded and dragged in the channel till great depth a weak power law serpentinite rheology produces a transitional behavior between that of strong power law and newtonian rheology the serpentinitic melange produced using a newtonian rheology is finally exhumed from pressure 20 kbar during the final stages of continental crust subduction and collision the rise to shallow depth of the hp melange is driven by concurrent serpentinite buoyancy slab roll back and asthenospheric mantle upwelling moreover we propose that the subduction exhumation dynamics observed in the model acted during closure of the narrow alpine tethys ocean in particular we suggest that some slices presently outcropping in hp ophiolitic massifs of the western alps voltri and monviso massifs were tectonically coupled at slab mantle interface and exhumed during collision in a soft low viscosity serpentinite channel we finally investigated intensity and dynamics of magmatic processes during subduction our numerical experiments suggest that even subduction of narrow oceanic basins may indeed allow the development of a magmatic arc in such relatively small short living arcs steeper faster and younger slabs produce earlier onset of the volcanism after the beginning of subduction

relationships between plate motions and late cretaceous to paleogene magmatism in southwestern alaska four major late cretaceous to paleogene magmatic belts have been recognized in southwestern alaska and tentative boundaries have been drawn for each the belts indicate minimum ages by which the tectonostratigraphic terranes they overprint must have been linked together though there has been significant postmagmatic tectonism the kuskokwim mountains belt 73 60 m y b p the alaska range belt 74 55 m y b p and the gulf of alaska belt 62 57 m y b p all may have been related to a period of rapid north northeastward motion of the kuia plate from about 74 to 56 m y b p the angle of convergence was low with respect to the present orientation of the belts suggesting the possibility of postmagmatic rotation of the belts a higher angle of convergence may have existed with respect to the alaska range belt a probable arc if postmagmatic rotation resulted from 300 400 km of right lateral displacement on the denali fault the gulf of alaska belt is in an anomalous forearc position with respect to the alaska range belt but probably formed by a convergence related mechanism the kuskokwim mountains belt is in an apparent back arc position but the mechanism of its origin is uncertain it may have been part of an unusually broad alaska range arc or less likely it may have been a separate arc magmatic quiescence occurred throughout south western alaska during a period of plate reorganization from about 56 to 43 m y b p a period which began with extremely rapid northward motion of the kuia plate magmatism resumed when the present period of rather slow northwestward pacific plate motion began at about 43 m y b p the aleutian arc probably formed at this time and has been episodically active in about the same location ever since

the role of mantle delamination in widespread late cretaceous extension and magmatism in the cordilleran orogen western united states extension during plate convergence and mountain building is widely recognized yet the causes of synconvergent extension remain controversial here we propose that delamination of lithospheric mantle aided by decoupling of the crust from the mantle via a reduction in the viscosity of the lower crust through heating incursion of fl uids and partial melting explains many enigmatic yet prevalent aspects of the metamorphic magmatic and kinematic history of the sevier laramide orogen of the western united states during the late cretaceous extension heating anatexis magmatism and perhaps rock uplift were widespread during a restricted time interval in the late cretaceous 75 67 ma along the axis of maximum crustal thickening within the mojave sector of the sevier orogen and to a lesser extent within the interior of the idaho utah wyoming sector to the north similar processes may have been active in the peninsular range sierran western mojave and salinian segments of the mesozoic cordilleran arc these processes are viewed as predictable consequences of the thermal rheological and dynamic state of the overlying crust following delamination of mantle lithosphere beneath isostatically compensated mountain belts the proposed delamination would have occurred immediately prior to eastward propagation of lowangle subduction of the farallon plate during the inception of the laramide orogeny following delamination extension and anatexis of the north american crust were aided locally by egress of slab derived fl uids from the low angle farallon slab we suggest that lithospheric delamination may have aided in the shallowing of the slab to achieve lowangle subduction geometry delamination has been proposed to be common in areas of thickened continental lithosphere in the terminal phase or late in orogenesis the late cretaceous delamination event proposed here for the sevier laramide orogen occurred during protracted plate convergence and was synchronous with and followed by continued shortening in the external part of the orogen

effects of synchronous uplift and intrusion during magmatic arc construction abstract individual areas in the southern patagonian batholith spb contain similar sequences of old mafic to young felsic plutons that crystallized during overlapping time intervals a positive correlation between age and initial 87 sr 86 sr sri suggests magmatism was not diachronous field relations and geobarometry show that mafic plutons crystallized at greater depths than felsic plutons an arc construction model involving continual uplift driven by continual magmatic inflation can explain the observed relationships in the model old mafic plutons are uplifted to shallower depths where younger more felsic plutons intrude adjacent to them model age lithology histograms show that exposed plutonic suites inaccurately record the timing and nature of arc magmatism the observed mafic to felsic sequence can result from uplift and erosion rather than from protracted changes in magma composition differences in plutonic ages crystallization depths and length of magmatic records can be used to interpret local and regional uplift histories during magmatic inflation a nodal plane must exist between the upper crust that experiences shallowing due to erosion and the lower crust that experiences deepening due to crustal thickening the depth of the nodal plane is located were the mass of magmatic additions above is balanced by the mass removed by erosion and is controlled by the relative rates of erosion and magmatic inflation and the vertical distribution of magmatic additions this depth may be buffered because increased magmatism shallows the nodal plane and accelerates erosion which lowers the nodal plane under steady state conditions rocks below the nodal plane will not be exposed by uplift and erosion

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

assessing the role of slab rheology in coupled plate mantle convection models reconstructing the 3d structure of the earth s mantle has been a challenge for geodynamicists for about 40 yr although numerical models and computational capabilities have substantially progressed parameterizations used for modeling convection forced by plate motions are far from being earth like among the set of parameters rheology is fundamental because it defines in a non linear way the dynamics of slabs and plumes and the organization of lithosphere deformation in this study we evaluate the role of the temperature dependence of viscosity variations up to 6 orders of magnitude and the importance of pseudo plasticity on reconstructing slab evolution in 3d spherical models of convection driven by plate history models pseudo plasticity which produces plate like behavior in convection models allows a consistent coupling between imposed plate motions and global convection which is not possible with temperature dependent viscosity alone using test case models we show that increasing temperature dependence of viscosity enhances vertical and lateral coherence of slabs but leads to unrealistic slab morphologies for large viscosity contrasts introducing pseudo plasticity partially solves this issue producing thin laterally and vertically more continuous slabs and flat subduction where trench retreat is fast we evaluate the differences between convection reconstructions employing different viscosity laws to be very large and similar to the differences between two models with the same rheology but using two different plate histories or initial conditions

a review of the geodynamic evolution of flat slab subduction in mexico peru and chile abstract subducting plates around the globe display a large variability in terms of slab geometry including regions where smooth and little variation in subduction parameters is observed while the vast majority of subduction slabs plunge into the mantle at different but positive dip angles the end member case of flat slab subduction seems to strongly defy this rule and move horizontally several hundreds of kilometers before diving into the surrounding hotter mantle by employing a comparative assessment for the mexican peruvian and chilean flat slab subduction zones we find a series of parameters that apparently facilitate slab flattening among them trench roll back as well as strong variations and discontinuities in the structure of oceanic and overriding plates seem to be the most important however we were not able to find the necessary and sufficient conditions that provide an explanation for the formation of flat slabs in all three subduction zones in order to unravel the origin of flat slab subduction it is probably necessary a numerical approach that considers also the influence of surrounding plates and their corresponding geometries on 3d subduction dynamics

impact of aseismic ridges on subduction systems insights from analog modeling the influence of aseismic ridges on subduction kinematics and dynamics and on the deformation of the overriding plate remains a topic of debate this study presents laboratory based geodynamic models that simulate the process of aseismic ridge subduction below an overriding plate the analog experiments show that depending on its size subduction of an aseismic ridge can impact on the kinematics of subduction overriding plate deformation and to a lower extent on the geometry of the slab and that of the trench specifically it can cause a reduction of the subducting plate and trench retreat velocities decreases overriding plate extension and even drives local forearc shortening for the widest ridge increases the slab dip angle next to the aseismic ridge but decreases the dip angle at the aseismic ridge in case the ridge is thick and wide and produces a local perturbation indentation of the trench geometry the magnitude of the described modifications relies to a comparable extent on aseismic ridge width and thickness the effects of aseismic ridge subduction observed in our models compare to those of the louisville seamount chain the cocos ridge and the d entrecasteaux ridge system in nature and provide an explanation for the local deformation observed at these subduction zones

plio quaternary geodynamic evolution of a segment of the peruvian andean cordillera located above the change in the subduction geometry the cuzco region abstract the cuzco region which is located above a change in subduction geometry appears to be characterized by a variable plio quaternary tectono sedimentary evolution essentially located along the major fault system that separates the high plateaux from the eastern cordillera after the higher surface formation of the high plateaux a set of neogene basins were filled by miocene fluvio torrential series and by plio pleistocene fluvio lacustrine deposits the neogene series have been affected by compressional tectonic forces attributed to the late miocene this compression is followed by roughly e w trending syn sedimentary extensional tectonics attributed to the pliocene it is related to reactivation of the pre existing major faults basin evolution and volcanic activity concentrated along the faults in the early pleistocene fluvio lacustrine deposits are affected by syn and post sedimentary compressional tectonism it is characterized by shortening that trends both n s and e w and produces folding and faulting of the sedimentary cover extensional tectonism trending roughly n s has been taking place from the middle pleistocene to the present it is coeval with shoshonitic volcanic activity and with sedimentation of fluvio lacustrine terraces torrential fans and moraines quaternary and active normal faults due to this tectonism are located in a narrow zone more than 100 km long between the high plateaux and the eastern cordillera and two 15 km long fault sectors in the eastern cordillera characteristic pleistocene scarps 400 m or more high are due to the cumulative normal offset and there are also little scarps with heights ranging between 2 and 20 m which are related to holocene fault reactivations recent fault reactivation on the cuzco fault system during the april 5 1986 earthquake m b 5 3 is due to the n s trending extension this state of stress located at a mean elevation of roughly 3730 m is generally homogeneous to different scales the active cuzco normal faults may be a consequence of adjustment between the compensated western cordillera and the undercompensated eastern cordillera this latter being uplifted higher than its isostatic equilibrium due to compression acting on its eastern edge the variation of the state of stress during the plio quaternary is in agreement with the variations of the compressional boundary forces it may be explained by variation of the convergence rate or by the variation of pull slab forces

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

accretion of japanese island arcs and implications for the origin of archean greenstone belts the present day region of arc arc collision between the izu bonin arc and mainland japan honshu arc in the western pacific called the izu collision zone or icz provides a useful kinematic model for the development of many archean greenstone belts the compositional assemblage thermal structure and structural style of the crust in the icz are very similar to those inferred for many archean greenstone belts in both cases thick young and flaky arc crust with a high heat flux is in collision with another crustal block the icz has resulted in the imbrication of very thick slices of crust and associated sediments juxtaposed against older accretionary prisms many archean greenstone belts show similar geological features and we suggest that the incremental imbrication of the crust was a common tectonic style in some archean convergent margins

cocos nazca slab window beneath central america abstract integration of petrologic and tectonic data favours a model of slab window formation beneath central america in the pliocene pleistocene central america has been the site of voluminous cenozoic arc volcanism the cocos and nazca plates which are subducting beneath central america are diverging along the east trending cocos nazca spreading ridge since 25 ma the americas have advanced about 1800 km west over the ridge transform system since at least 8 ma plate integrity and the ridge transform configuration have been preserved during convergence resulting in subduction of the spreading ridge and development of a slab window the panama fracture zone an active transform fault is the part of the ridge transform system currently being subducted the ridge transform system formerly adjoining the northern end of the panama fracture zone is likely to have been left stepping we use present day plate motions to design a slab window to fit known variations in igneous composition hypocentre distribution and mantle anisotropy the modeling demonstrates that subduction of ridge segments and resultant slab window development began between 6 and 10 ma cessation of ridge subduction occurred between 1 and 3 ma when subduction of the panama fracture zone is considered to have begun the slab window is continuing to expand and migrate northeastward below the central american volcanic arc the absence of a wadati benioff zone from southeastern costa rica through panama corresponds to the position of the slab window within this region dacitic and rhyolitic volcanic rocks have adakitic compositions and are thought to result from anatexis of the young buoyant crust which forms the trailing edges of the slabs bounding the window basalts in this area were derived from an enriched ocean island type mantle source whereas basalts from the rest of the arc in nicaragua el salvador and guatemala are mainly derived from slab modified depleted mantle characteristic of volcanic arcs the presence of ocean island type mantle beneath southern costa rica and panama is explained by eastward flow of enriched asthenosphere from the galapagos plume head through the slab window and into the volcano source region eastward transfer of asthenosphere is consistent with global plate motion studies and seismic anisotropy in the asthenosphere beneath the nazca and caribbean plates the flow of peridotite is a consequence of progressive shrinkage of the pacific mantle reservoir and concurrent expansion of the atlantic mantle reservoir

horizontal subduction zones convergence velocity and the building of the andes we discuss the relationships between andean shortening plate velocities at the trench and slab geometry beneath south america although some correlation exists between the convergence velocity and the westward motion of south america on the one hand and the shortening of the continental plate on the other hand plate kinematics neither gives a satisfactory explanation to the andean segmentation in general nor explains the development of the bolivian orocline in paleogene times we discuss the cenozoic history of horizontal slab segments below south america arguing that they result from the subduction of oceanic plateaus whose effect is to switch the buoyancy of the young subducting plate to positive we argue that the existence of horizontal slab segments below the central andes during eocene oligocene times and below peru and north central chile since pliocene resulted 1 in the shortening of the continental plate interiors at a large distance from the trench 2 in stronger interplate coupling and ultimately 3 in a decrease of the trenchward velocity of the oceanic plate present day horizontal slab segments may thus explain the diminution of the convergence velocity between the nazca and south american plates since late miocene

the geometry of the wadati benioff zone under southern central america and its tectonic significance results from a high resolution local seismographic network abstract we present here a detailed geometry of the wadati benioff zone under costa rica obtained from seismicity recorded by a dense local seismographic network jointly operated by the costa rica volcanological and seismological observatory national university and the charles f richter seismological laboratory university of california santa cruz underneath the nicaragua costa rica border the wadati benioff zone smoothly contorts from steep to shallow dip angles nw to se but does not show evidence of a brittle tear as postulated by others however further to the se ne of puerto caldera costa rica the wadati benioff zone does show a segmentation the quesada sharp contortion at intermediate depths h 70km nw of this sharp contortion the deepest portion of the seismically active slab dips at about 80 and reaches maximum depths ranging from 200 km near the nicaragua costa rica border to 135 km under ciudad quesada to the se the deeper portion of the wadati benioff zone dips at about 60 and the seismicity does not extend below depths ranging from 125 km behind the volcanic arc to 50 km east of quepos in southern costa rica east of 83 55 w we find no evidence of the wadati benioff zone deeper than 50 km the obtained geometry and other known tectonic features related to the subduction of the cocos plate under the caribbean plate along the southern terminus of the middle america trench nicaragua and costa rica correlate well with along trench variations in age of the subducted cocos plate some of these tectonic features are 1 the shallowing of middle america trench bathymetry from nw to se 2 variations in the energy release within the subducted slab 3 differences in coupling between cocos and caribbean plates 4 the termination of the central america volcanic chain in central costa rica 5 distinct stress field variations on the overriding caribbean plate the subduction of the cocos ridge under southern costa rica is partially responsible for some of these features shallowing of the middle america trench in southern costa rica and for the high uplift rates of the outer arc however as the presence of the panama fracture zone limits the subducted extension of the cocos ridge to less than 100 km from the trench we propose that the overall geometry of the wadati benioff zone is controlled by abrupt along trench changes in the age of the subducted cocos plate

late paleogene extensional collapse of the cordilleran foreland fold and thrust belt the cordilleran foreland fold and thrust belt collapsed and spread to the west during a middle eocene to early miocene ca 49 20 ma episode of crustal extension the sedimentary and structural record of this event is preserved in a network of half grabens that extends from southern canada to central utah extensional structures superposed on this allochthonous terrain are rooted to the physical stratigraphy structural relief and sole faults of preexisting thrust fold structures the sole faults dip 3 6 west above an undeformed precambrian crystalline basement and accommodated tectonic transport of a thick up to 20 km eastward tapering hanging wall during regimes of crustal shortening and extension the chronology of tectonism for the foreland fold and thrust belt is established here by dating latest thrusting and initial normal faulting and is best defined where thrusts and normal faults are linked by common detachment surfaces dated movement on two extensionally reactivated thrusts the lewis thrust of northwest montana and southeast british columbia and the medicine butte thrust of southwest wyoming and northeast utah suggests that the hiatus between the end of crustal shortening in the early or early middle eocene and the start of extension in the early middle eocene was brief lateral spreading and extensional basin formation in the cordilleran foreland fold and thrust belt were partly concurrent with formation of metamorphic core complexes and regional magmatism conceptually linking extensional processes that were simultaneously deforming both the hinterland and foreland of the late paleogene cordilleran orogenic wedge is accomplished by applying the extensional wedge coulomb critical taper model the rapid drop in north america pacific plate convergence rate and or steepening of the subducted oceanic slab at ca 50 ma resulted in a large reduction in east west horizontal compressive stress in the cordillera as a result the cordilleran orogenic wedge was left unsupported and it gravitationally collapsed and horizontally spread west until a new equilibrium was established at ca 20 ma subsequently crustal extension and magmatism during the basin and range event ca 17 0 ma overprinted much of this earlier phase of extension

incipient india eurasia collision and plutonism the lower cenozoic batura granites hunza karakorum north pakistan plutonism of the karakorum north of the india eurasia suture zone comprises two major intrusive stages the oldest mid cretaceous predates the collision the youngest miocene is post suturing in between lower cenozoic plutonism roughly syn suturing and represented by the batura complex is made up of a few small and usually undeformed plutons the complex has a subalkaline i e intermediate between alkaline and calc alkaline and ferriferous association dominated by light coloured metaluminous and slightly peraluminous granites and adamellites with biotite amphibole titanite whole rock rb sr isochron ages of 63 ma this study and 43 ma debon et al 1987b have been obtained for two plutons with rather low initial 87 sr 86 sr ratios 0 7050 0 7056 a comparison with granitoids of similar age from kohistan suggests that the complex relates to the subduction of the tethys ocean and originates from a mantle source with but small crustal contribution through a two stage melting process the apparently limited extent of lower cenozoic plutonism in the karakorum contrasts with its wider development in the transhimalaya and kohistan this contrast could be partly accounted for by changes in tectonic and magmatic suduction related processes after the initial contact between india and eurasia was established in the nw himalaya

shear wave velocities in the pampean flat slab region from rayleigh wave tomography implications for slab and upper mantle hydration 1 the pampean flat slab region located in central argentina and chile between 29 and 34 s is considered a modern analog for laramide flat slab subduction within western north america regionally flat slab subduction is characterized by the nazca slab descending to 100 km depth flattening out for 300 km laterally before resuming a more normal angle of subduction flat slab subduction correlates spatially with the track of the juan fernandez ridge and is associated with the inboard migration of deformation and the cessation of volcanism within the region to better understand flat slab subduction we combine ambient noise tomography and earthquake generated surface wave measurements to calculate a regional 3d shear velocity model for the region shear wave velocity variations largely relate to changes in lithology within the crust with basins and bedrock exposures clearly defined as low and high velocity regions respectively we argue that subduction related hydration plays a significant role in controlling shear wave velocities within the upper mantle in the southern part of the study area where normal angle subduction is occurring the slab is visible as a high velocity body with a low velocity mantle wedge above it extending eastward from the active arc where flat slab subduction is occurring slab velocities increase to the east while velocities in the overlying lithosphere decrease consistent with the slab dewatering and gradually hydrating the overlying mantle the hydration of the slab may be contributing to the excess buoyancy of the subducting oceanic lithosphere helping to drive flat slab subduction

a modelling study of vertical surface displacements at convergent plate margins summary during the evolution of a subduction zone system variations are likely to occur in for example surface plate velocities and buoyancy of the subducting lithosphere we quantify vertical surface displacements at convergent plate margins resulting from such imposed variations for this purpose we use a 2 d numerical model in which the lithospheric plates have an effective elastic thickness we first define a model in which the subducting plate is driven by its negative buoyancy and a velocity at its surface side boundary its equilibrium topography after around 2 myr is the reference level for examination of surface displacements resulting from variations in buoyancy velocity of the surface plates friction along the interplate contact and subduction zone roll back we find that a decrease increase in buoyancy of the subducting material leads to a deepening uplift of the plate margins an increase in friction along the subduction fault deepens the overriding plate margin subduction zone roll back due to sinking of the negatively buoyant subducting plate induces subsidence of the overriding plate margin this subsidence is reduced when roll back takes place in a land locked basin setting trench retreat forced by the motion of the overriding plate is characterized by higher topography of the overriding plate margin than the case of retreat due to the sinking of the negatively buoyant slab in the first case in plane stress in the back arc region is compressive while it is tensional for roll back due to the sinking of the slab we conclude that vertical surface displacements during ongoing subduction may reach a magnitude of a few kilometres on the overriding and subducting plate margins