space geodetic observations of nazca south america convergence across the central andes space geodetic data recorded rates and directions of motion across the convergent boundary zone between the oceanic nazca and continental south american plates in peru and bolivia roughly half of the overall convergence about 30 to 40 millimeters per year accumulated on the locked plate interface and can be released in future earthquakes about 10 to 15 millimeters per year of crustal shortening occurred inland at the sub andean foreland fold and thrust belt indicating that the andes are continuing to build little 5 to 10 millimeters per year along trench motion of coastal forearc slivers was observed despite the oblique convergence

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

molecular and stable carbon isotopic compositions of saturated fatty acids within one sedimentary profile in the shenhu northern south china sea source implications abstract this study examined the distributions and stable carbon isotopic compositions of saturated fatty acids safas in one 300 cm long sedimentary profile which was named as site4b in shenhu northern south china sea the concentrations of total safas in sediments ranged from 1 80 to 10 16 g g g fa g dry sediment and showed an even over odd predominance in the carbon chain of c12 to c32 mostly with n c16 and n c18 being the two major components the short chain fatty acids scfas n c12 to n c18 mainly from marine microorganisms had average 13c values of 26 7 to 28 2 whereas some terrigenous sourced long chain fatty acids lcfas n c21 to n c32 had average 13c values of 29 6 to 34 1 the other lcfas n c24 n c26 n c28 average 13c values are 26 1 to 28 0 as well as n c19 and n c20 safas average 13c values are 29 1 and 29 3 respectively showed a mixed signal of carbon isotope compositions the relative bioproductivity calculation marine vs terrigenous demonstrated that most of organic carbon accumulation throughout the sedimentary profile was contributed by marine organism the high marine productivity in shenhu south china sea may be related to the hydrocarbon seepage which evidenced by diapiric structures interestingly there is a sever fluctuation of terrigenous inputs around the depth of 97 cm below the seafloor bsf probably resulting from the influence of the dansgaard oeschger events and the younger dryas event as revealed by 14c age measurements

neotectonic deformation of the central chile margin deepwater forearc basin formation in response to hot spot ridge and seamount subduction 1 underlying the midslope valparaiso terrace the valparaiso basin constitutes a unique and prominent deepwater forearc basin on the central chile margin the basin contains a 3 3 5 km thick sediment fill of late cenozoic age the seismic stratigraphic architecture of the valparaiso basin fill shows that margin wide subsidence above the underthrusting juan fernandez ridge an aseismic volcanic swell crested with seamounts on the nazca plate characterized the basin evolution contemporaneously with margin subsidence and seaward rotation of the upper slope and shelf the outer part of the forearc was deformed by thrust faulting above the leading flank of underthrusting seamounts that rotated the middle slope to a more horizontal attitude and formed seaward vergent outer forearc highs subsidence of the margin above the underthrusting juan fernandez ridge is ascribed to processes of subduction erosion enhanced by factors related to subduction of the ridge that caused dismemberment of upper plate material along the shallow aseismic part of the plate interface and transfer of material to the downgoing nazca plate accelerated midslope subsidence 55 75 km landward of the trench coincides with the seaward limit of interplate seismicity suggesting that enhanced basal erosion is related to the transition from aseismic to stick slip subduction behavior since subduction of the juan fernandez ridge began near 33 s 10 ma ago the chile trench north of 33 s is inferred to have retreated 3 km myr and the midslope valparaiso basin subsided 0 3 0 5 km myr approximating analogous neogene rates for the peruvian margin

bathymetry of mariana trench arc system and formation of the challenger deep as a consequence of weak plate coupling 1 the challenger deep in the southernmost mariana trench western pacific ocean is the deepest point on the earth s surface 10 920 m below sea level its location within a subduction trench where one plate bends and descends below another is not surprising however why is it located in the southernmost mariana trench and not at its central part where the rate of subduction is higher where the lithosphere is the oldest and densest on the earth and where the subducted lithosphere pulling down is the longest in the earth 1000 km or more according to seismic tomography we suggest that although subduction rate and slab age generally control trench depth the width of the plate coupling zone is more important beneath the central marianas the subducted slab is attached to the upper plate along a 150 km wide surface that holds the shallow portion of the subducted plate nearly horizontal in spite of its great load and thus counters trench deepening in contrast along the south mariana trench the subducted length of the lithosphere is much shorter but its attachment to the upper plate is only along a relatively narrow 50 km wide surface in addition a tear in the slab beneath this region helps it to sink rapidly through the mantle and this combination of circumstances allows the slab to steepen and form the deepest trench on the earth in a wider perspective the interrelations shown here between trench deepening ridge shallowing slab steepening and forearc narrowing may shed light on other subduction zones located near edges of rapidly steepening slabs

tectonic erosion of the peruvian forearc lima basin by subduction and nazca ridge collision 1 subsidence of lima basin part of the peruvian forearc is controlled by tectonic erosion by the subducting nazca plate multichannel seismic reflection data coupled with age and paleowater depth constraints derived from ocean drilling program odp coring now allow the rates of erosion to be reconstructed through time in trenchward locations the forearc has experienced limited recent relative uplift 700 850 m likely due to preferential basal erosion under the center of lima basin long term subsidence driven by basal tectonic erosion dominates and is fastest closest to the trench since 47 ma eocene up to 148 km of the plate margin have been lost at an average rate of up to 3 1 km myr 1 appoximately 110 km of that total appears to be lost since 11 ma implying much faster average rates of trench retreat 10 km myr 1 since collision of the nazca ridge with the lima basin at 11 ma although there is no clear subsidence event at odp site 679 during the time at which nazca ridge was subducting beneath this part of the forearc 4 11 ma the more trenchward odp sites 682 and 688 show significant deepening after 11 ma indicating that subduction of the ridge accelerates tectonic erosion long term rates of crustal erosion in the region of lima basin are greater than estimates of regional arc magmatic productivity implying that such margins are net sinks of continental crust

tectonic erosion and consequent collapse of the pacific margin of costa rica combined implications from odp leg 170 seismic offshore data and regional geology of the nicoya peninsula the convergent margin off the pacific coast of the nicoya peninsula of costa rica exhibits evidence for subduction erosion caused by the underthrusting cocos plate critical evidence for efficacy of this process was recovered at the ocean drilling program odp drilling site 1042 leg 170 positioned 7 km landward of the middle america trench axis off the nicoya peninsula the primary drilling objective at this site was to identify the age and origin of a regionally extensive and prominent seismic discontinuity the so called base of slope sediment boss horizon or surface the boss horizon which can be traced landward from near the trench to the nicoya coastal area and parallel to it for hundreds of kilometers separates a low velocity 2 0 2 5 km s 1 sequence of slope sediment from an underlying sequence of much higher velocity 4 4 5 km s 1 rock site 1042 reached the acoustically defined boss horizon at a below sea level depth of 3900 m and yielded a carbonate cemented calcarenitic breccia of early middle miocene age sedimentological geochemical paleontological and cement paragenesis data document that the breccia accumulated in a shallow water depositional environment on the basis of coastal exposures the boss horizon as a margin wide geologic interface can be temporally and lithostratigraphically correlated to a regional angular unconformity this unconformity known as the mal pais unconformity separates neogene and younger shelf to littoral beds from the underlying mafic units of the mesozoic nicoya complex and cretaceous and early tertiary sedimentary sequences at site 1042 it is inferred that tectonism caused the vertical subsidence of the early neogene breccia from a shallow to a deep water setting the mal pais unconformity of the boss horizon thus connects the rock fabric of the outermost part of margin to that of coastal nicoya and implies that in the early neogene the nicoya shelf extended seaward to near the present trench axis this circumstance requires that the early neogene trench axis was at least 50 km seaward of where it is now located the long term effects of subduction erosion similar to those described for the scientifically drilled japan tonga and peru margins best account for offshore and onshore evidence for a post paleogene history of crustal thinning and landward trench migration of costa rica s pacific margin during the past 16 17 myr the calculated mass removal and landward migration rates are 34 36 km myr 1 km 1 of margin and 3 km myr 1 respectively these values are similar to those found for other pacific margins dominated by nonaccretionary subduction zone processes

global tectonic significance of the solomon islands and ontong java plateau convergent zone abstract oceanic plateaus areas of anomalously thick oceanic crust cover about 3 of the earth s seafloor and are thought to mark the surface location of mantle plume heads hotspot tracks represent continuing magmatism associated with the remaining plume conduit or tail it is presently controversial whether voluminous and mafic oceanic plateau lithosphere is eventually accreted at subduction zones and therefore 1 influences the eventual composition of continental crust and 2 is responsible for significantly higher rates of continental growth than growth only by accretion of island arcs the ontong java plateau ojp of the southwestern pacific ocean is the largest and thickest oceanic plateau on earth and the largest plateau currently converging on an island arc solomon islands for this reason this convergent zone is a key area for understanding the fate of large and thick plateaus on reaching subduction zones this volume consists of a series of four papers that summarize the results of joint us japan marine geophysical studies in 1995 and 1998 of the solomon islands ontong java plateau convergent zone marine geophysical data include single and multi channel seismic reflection ocean bottom seismometer obs refraction gravity magnetic sidescan sonar and earthquake studies objectives of this introductory paper include 1 review of the significance of oceanic plateaus as potential contributors to continental crust 2 review of the current theories on the fate of oceanic plateaus at subduction zones 3 establish the present day and neogene tectonic setting of the solomon islands ontong java plateau convergent zone 4 discuss the controversial sequence and timing of tectonic events surrounding ontong java plateau solomon arc convergence 5 present a series of tectonic reconstructions for the period 20 ma early miocene to the present day in support of our proposed timing of major tectonic events affecting the ontong java plateau solomon islands convergent zone and 6 compare the structural and deformational pattern observed in the solomon islands to ancient oceanic plateaus preserved in precambrian and phanerozoic orogenic belts our main conclusion of this study is that 80 of the crustal thickness of the ontong java plateau is subducted beneath the solomon island arc only the uppermost basaltic and sedimentary part of the crust 7 km is preserved on the overriding plate by subduction accretion processes this observation is consistent with the observed imbricate structural style of plateaus and seamount chains preserved in both precambrian and phanerozoic orogenic belts

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

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

influence of the subducting plate velocity on the geometry of the slab and migration of the subduction hinge abstract geological observations indicate that along two active continental margins east asia and mediterranean major phases of overriding plate extension resulting from subduction hinge retreat occurred synchronously with a reduction in subducting plate velocity in this paper results of fluid dynamical experiments are presented to test the influence of the velocity of the subducting plate on the hinge migration velocity and on the geometry of the slab results show that hinge retreat decreases with increasing subducting plate velocity in addition phases of hinge retreat alternate with phases of hinge advance for relatively high subducting plate velocities due to interaction of the slab with the bottom of the box simulating the upper lower mantle discontinuity such slab kinematics could explain the episodic behaviour of back arc opening observed in convergent settings the geometry of the slab and the kinematics of subduction are significantly affected by the velocity of the subducting plate three subduction modes with accompanying slab geometry can be recognized a relatively low subducting plate velocity is accompanied by relatively fast hinge retreat with backward sinking of the slab and a backward draping slab geometry with increasing subducting plate velocity hinge migration is relatively small resulting in subvertical sinking of the slab and a folded slab piling geometry for a very high velocity the hinge migrates forward resulting in forward oriented subduction vectors and a forward draping slab geometry

terrestrial and gps measurements of deformation across the taupo back arc and hikurangi forearc regions in new zealand geodetic measurements of deformation across the taupo back arc and hikurangi forearc regions of new zealand are derived from global positioning system gps measurements made in 1990 and 1991 and triangulation observations made in the 1920s 1950s and 1970s the gps horizontal coordinate differences have precisions of 5 6 mm and the triangulation observations have precisions of 0 7 1 2 arc sec these different kinds of observation allow simultaneous estimation of strain parameters and single epoch coordinates for a total of 184 stations at 68 of which gps observations were made under the assumption that no shift of the scale or orientation of the gps reference frame is embodied in the data from 13 gps stations within the taupo volcanic zone tvz common to the 1990 and 1991 surveys the rotational dilatational and shear components of the deformation rate tensor within the tvz can be estimated for that 1 year interval the principal extension rate from 1990 to 1991 is 0 21 0 09 10 6 yr 68 confidence at an azimuth of 124 13 corresponding to 8 4 mm yr extension over 40 km neither the orthogonal principal extension rate 0 02 0 03 10 6 yr nor the rotation rate 0 06 0 04 rad yr is significant the dilatation rate of 0 23 0 09 10 6 yr is therefore produced by uniaxial extension the data from 75 stations distributed from behind the back arc region and well into the forearc region which were surveyed in at least two of the triangulation or gps epochs between the 1920s and 1991 allow the spatial variation of the shear components of the deformation rate tensor to be estimated the maximum engineering shear rates lie in the range 0 1 0 2 10 6 yr within the tvz and maintain these values with similar orientations in the forearc region to the east the western and southern margins of active extension are reasonably well determined by changes in orientation and magnitude of the shear component these results provide confirmation of the previously less well determined deformation field and in comparison with subduction models indicate that the surface deformation is reflecting variation from north to south of coupling of the subduction plate interface and that both trench suction and gravitational collapse probably contribute to the extension in the forearc region

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

hotspots in the southern oceans an absolute frame of reference for motion of the gondwana continents abstract the geometry and geochronology of aseismic ridges and oceanic islands in the southern oceans provide a good test of the proposition that hotspots remain fixed over long periods of time that is motion of an order of magnitude less than the relative motion between plate pairs in most cases it is concluded that inter hotspot movement cannot be discerned for the period 100 m y to present and that widely distributed hotspots in the atlantic and indian oceans provide a frame of reference for plate motions following the disintegration of gondwanaland which is independent of paleomagnetism this frame of reference is absolute in that it gives the motion of the lithosphere with respect to the mantle hotspots the absolute motion model indicates that africa and antarctica are now moving only very slowly that there has been significant relative movement between east and west antarctica since the cretaceous and prescribes the relative motion between the somali and african plates

a three plate kinematic model for lau basin opening 1 we present revised compilations of bathymetry magnetization acoustic imagery and seismicity data for the lau basin and surrounds we interpret these data to more precisely locate the plate boundaries in the region and to derive a three plate kinematic model for the opening of the lau basin during the brunhes chron our tectonic model includes the niuafo ou microplate which is separated from the australian plate along the peggy ridge lau extensional transform zone central lau spreading center and from the tongan plate along the newly discovered fonualei rift and spreading center our model shows that australia tonga plate motion between 15 5 s and 19 s is partitioned across the microplate and it resolves the former apparent conflict between geodetic versus seafloor spreading velocities regionally this implies that the angular rate of opening of the lau basin has been faster than that of the havre trough to the south and that the rapid current rates of pacific subduction at the tonga trench 240 mm yr at 16 s have been sustained for at least 0 78 m y

interseismic and coseismic motions in gps series related to the ms 7 3 july 13 1994 malekula earthquake central new hebrides subduction zone on july 13 1994 an earthquake ms 7 3 occurred at malekula in the new hebrides archipelago the gps data collected across the new hebrides trench between 1990 and 1996 were processed in order to separate the interseismic and coseismic motions from the drifts related to the convergence with the australian plate the gps derived coseismic displacements at the gps site in malekula are 49 15 mm southward 230 30 mm westward and 170 37 mm downward when the cmt derived displacements are 50 mm southward 210 mm westward and 150 mm downward taking into account the interseismic strain accumulation 25 mm yr at the source established from historical seismicity 7 5 mm yr at the gps site the strain free convergence rate at malekula is 49 3 mm yr other gps derived convergence rates are 95 1 mm yr at efate and 37 2 mm yr at santo these rates imply a regional right lateral motion between the efate and the santo malekula segments in contrast the focal mechanism of the earthquake mostly indicates a left lateral motion therefore we hypothesize that the earthquake is related to variations in the interplate coupling along the converging boundary of the santo malekula segment

transient and permanent deformation of central japan estimated by gps 2 strain partitioning and arc arc collision abstract deformation of central japan can be divided into two main components elastic deformation due to the interseismic locking of subduction thrusts and deformation associated with intra arc tectonics in a companion paper henry et al earth planet sci lett 2000 this issue we estimated an elastic deformation field based on the modeling of global positioning system gps strain rates in this study we show that the corrected field gps minus predicted elastic deformation is consistent with geological and seismological observations and should represent long term plastic deformation we use this corrected gps deformation field to describe the kinematics of partitioning of the philippine sea plate amurian plate amu convergence in relation with the izu bonin arc collision about one third of the relative philippine sea amur motion is transferred to the central japan crust north of the collision front this transfer induces high strain rates in southern kanto and in the northern japanese alps the latter region corresponds to the eastern boundary of the amu we propose that the independent motion of central japan results from the delamination and the underthrusting of the izu bonin arc lower crust and the upward bent oceanic slab

underthrusting at the hjort trench australian pacific plate boundary incipient subduction 1 the southernmost portion of the macquarie ridge complex mrc comprises the hjort ridge trench and plateau the topographic expression of the australian pacific plate boundary south of new zealand on the basis of marine geophysical swath bathymetry reflectivity seismic reflection gravity magnetic data teleseismic data and gravity modeling we argue that the australian plate is actively underthrusting the pacific plate along the hjort trench but self sustaining subduction does not yet appear to have commenced we interpret a crustal discontinuity in the trench as a shallowly dipping thrust splay off a sub vertical transform fault in the crest of the hjort ridge hjort fault the hjort fault separates oceanic crust generated at the southeast indian ridge seir crust located to the south from oceanic crust generated at the extinct australian pacific spreading center located to the north mrc crust oblique convergence on the transform has initiated a significant plate boundary decollement in the trench but our estimate of 50 km of underthrusting does not support the existence of an eastwardly dipping australian slab below 20 km depth for the length of the hjort ridge and trench system oblique convergence is partitioned between the decollement in the trench and the hjort fault south of 57 5 s the trench decollement accommodates thrusting north of 57 5 s the boundary parallel component of convergence becomes dominant the trench gradually shallows and the trench decollement evolves into a low angle oblique slip fault the trench decollement and strike slip system in the hjort ridge are structural boundaries that contain tectonic slivers of inferred seir crust of the western flank of the hjort ridge that currently belong to neither the australian nor the pacific plates

evolving force balance during incipient subduction nearly half of all active subduction zones initiated during the cenozoic all subduction zones associated with active back arc extension have initiated since the eocene hinting that back arc extension may be intimately associated with an interval several tens of myr following subduction initiation that such a large proportion of subduction zones are young indicates that subduction initiation is a continuous process in which the net resisting force associated with forming a new subduction zone can be overcome during the normal evolution of plates subduction initiation is known to have occurred in a variety of tectonic settings old fracture zones transform faults and extinct spreading centers and through polarity reversal behind active subduction zones although occurring within different tectonic settings four known subduction initiation events izu bonin mariana ibm along a fracture zone tonga kermadec along an extinct subduction boundary new hebrides within a back arc and puysegur fiordland along a spreading center were typified by rapid uplift within the forearc followed by sudden subsidence other constraints corroborate the compressive nature of ibm and tonga kermadec during initiation using an explicit finite element method within a two dimensional domain we explore the evolving force balance during initiation in which elastic flexure viscous flow plastic failure and heat transport are all considered in order to tie theory with observation known tectonic settings of subduction initiation are used as initial and boundary conditions we systematically explore incipient compression of a homogeneous plate a former spreading center and a fracture zone the force balance is typified by a rapid growth in resisting force as the plate begins bending reaching a maximum value dependent on plate thickness but typically ranging from 2 to 3 1012 n m for cases that become self sustaining this is followed by a drop in stress once a shear zone extends through the plate the formation of a throughgoing fault is associated with rapid uplift on the hanging wall and subsidence on the footwall cumulative convergence not the rate of convergence is the dominant control on the force balance viscous tractions influence the force balance only if the viscosity of the asthenosphere is 1020 pa s and then only after plate failure following plate failure buoyancy of the oceanic crust leads to a linear increase with crustal thickness in the work required to initiate subduction the total work done is also influenced by the rate of lithospheric failure a self sustaining subduction zone does not form from a homogeneous plate a ridge placed under compression localizes subduction initiation but the resisting ridge push force is not nearly as large as the force required to bend the subducting plate the large initial bending resistance can be entirely eliminated in ridge models explaining the propensity for new subduction zones to form through polarity reversals a fracture zone fz placed in compression leads to subduction initiation with rapid extension of the overriding plate a fz must be underthrust by the older plate for 100 150 km before a transition from forced to self sustaining states is reached in fz models the change in force during transition is reflected by a shift from forearc uplift to subsidence subduction initiation is followed by trench retreat and back arc extension moderate resisting forces associated with modeled subduction initiation are consistent with the observed youth of pacific subduction zones the models provide an explanation for the compressive state of western pacific margins before and during subduction initiation including ibm and tonga kermadec in the eocene and the association of active back arcs with young subduction zones on the basis of our dynamic models and the relative poles of rotation between pacific and australia during the eocene we predict that the northern segment of the tonga kermadec convergent margin would have initiated earlier with a progressive southern migration of the transition between forced and self sustaining states

global trench migration velocities and slab migration induced upper mantle volume fluxes constraints to find an earth reference frame based on minimizing viscous dissipation sincetheadventofplatetectonicsdifferentglobalreferenceframeshavebeenusedtodescribethemotionofplatesandtrenches thedifferencein plate motion and trench migration between different reference frames can be substantial up to 4 cm yr this study presents an overview of trench migration velocities for all the mature and incipient subduction zones on earth as calculated in eight different global reference frames calculations show that irrespective of the reference frame 1 trench retreat always dominates over trench advance with 62 78 of the 244 trench segments retreating 2 the mean and median trench velocity are always positive retreating and within the range 1 3 1 5 cm yr and 0 9 1 3 cm yr respectively 3 rapid trench retreat is only observed close to lateral slab edges b1500 km and 4 trench retreat is always slow far from slab edges n2000km thesecalculationsarepredictedbygeodynamicmodelswithavaryingslabwidth inwhichplatemotion trenchmotionandmantleflow result from subduction ofdense slabs suggesting that trenchmotion is indeedprimarily drivenby slab buoyancyforces andthat proximity to alateral slab edge exerts a dominant control on the trench migration velocity despite these four general conclusions significant differences in velocities between such reference frames remain it is therefore important to determine which reference frame most likely describes the true absolute velocities togetanunderstandingoftheforcesdrivingplatetectonicsandmantleconvection itishereproposedthat basedonfluiddynamicconsiderationsand predictions from geodynamic modelling the best candidate is the one which optimises the number of trench segments that retreat minimizes the trench perpendicular trench migration velocity vt in the centre of wide n4000 km subduction zones maximizes the number of retreating trench segments located within 2000 km of the closest lateral slab edge minimizes the average of the absolute of the trench perpendicular trench migration velocity vt forallsubductionzonesonearth andminimizestheglobaluppermantletoroidalvolumeflux to thatresultsfromtrenchmigration and associated lateral slab migration i e slab rollback or slab roll forward calculations show that these conditions are best met in one particular indo atlantic hotspot reference frame where 75 of the subduction zones retreat vt in the centre of wide subduction zones ranges between 3 5 and 1 8 cm yr 83 of the trench segments located within 2000 km of the closest lateral slab edge retreat the average of vt is 2 1 cm yr and to 456 km 3 yr lower limit and 539 km 3 yr upper limit inclusion of all the incipient subduction zones on earth results in slightly greater fluxes of 465 km 3 yr lower limit and 569 km 3 yr upper limit it is also found that this reference frame is close to minimizing the total sub lithospheric uppermantlevolumeflux k associatedwithmotionofcontinentalkeelslocatedbelowthemajorcratons itisstressed however that kisanorder of magnitude smaller than to and thus of subordinate importance in conclusion the indo atlantic hotspot reference frame appears preferable for calculating plate velocities and plate boundary velocities

influence of trench width on subduction hinge retreat rates in 3 d models of slab rollback subduction of tectonic plates limited in lateral extent and with a free trailing tail i e free subduction is modeled in a three dimensional 3 d geometry the models use a nonlinear viscoplastic rheology for the subducting plate and exhibit a wide range of behaviors depending on such plate characteristics as strength width and thickness we investigate the time evolution of this progressive rollback subduction measure the accompanying return flow in the upper mantle and quantify the plate kinematics due to the 3 d geometry flow is allowed to accompany slab rollback around the lateral edges of the slab the toroidal component as opposed to 2 d geometry where material is forced to flow underneath the slab tip the poloidal component a simple force balance is provided which relates the speed of backward trench migration to the resistive forces of generating flow and weakening the plate our results indicate most of the gravitational energy of the system i e the negative buoyancy of the subducting slab is converted into a toroidal flow 69 a much smaller amount goes into weakening the plate 18 and the remaining amount goes into driving flow parallel to displacement of the slab 13 for the trench widths w we investigate 1500 km a maximum trench retreat rate occurs for trenches 600 km wide which is attributed to the interaction between a plate of finite width and the induced flow which has a lengthscale in the horizontal direction these numerical results quantitatively agree with comparable 3 d laboratory experiments using analogue models with a purely viscous plate material schellart 2004a 2004b including correlations between increasing retreat rate with increasing plate thickness trench width for maximum retreat rate 500 km and estimated amount of slab buoyancy used to drive rollback induced flow 70 several implications for plate tectonics on earth result from these models such as rollback subduction providing a physical mechanism for ephemeral slab graveyards situated above the more viscous lower mantle and endothermic phase transition prior to a flushing event into the lower mantle mantle avalanche

genetic relations between the central and southern philippine trench and the sangihe trench we surveyed the junction between the central and southern philippine trench and the sangihe trench near 6on using swath bathymetry gravity and magnetics these data along with seismicity allow us to discuss the genetic relations between these trenches and the forces acting on converging plates our final model favors the northern extension of the halmahera arc up to 8on with three segments offset left laterally along nw se transform faults accretion of the northern segment to mindanao island 4 to 5 m y ago resulted in the failure within the philippine sea plate east of the arc initiation of the philippine trench between 7on and 10on agrees with the maximum recorded depth of the philippine trench floor 10 000 m below sea level and philippine sea slab 200 km south of 6on trench junction another segment of the arc is being subducted beneath the sangihe margin while south of 3on the southern segment of the halmahera arc is still active the rapid southward shallowing of the trench floor along the southern philippine trench the type of faulting affecting both sides of the trench the lack of significant interplate seismicity and the concentration of the seismicity beneath the miangas talaud ridge are interpreted as a slowing down of the subduction along this branch of the philippine trench compared with the rest of the subduction zone the sangihe deformation front has been recognized up to 7on but seems active only south of 6on

tethyan subducted slabs under india abstract tomographic imaging of the mantle under tibet india and the adjacent indian ocean reveals several zones of relatively high p wave velocities at various depths under the hindu kush region in northeastern afghanistan and southern tajikistan a regional northward dipping slab is seen in the entire upper 600 km of the mantle and is apparently still attached to the lithosphere of the indian plate under northern pakistan this same slab shows a roll over structure with the deeper portion overturned and dipping southward as can also be seen in the distribution of earthquake hypocenters farther east southeast e g in the vicinity of nepal a well resolved anomaly below 450 km depth is connected to the slab under the hindu kush but seems to be separated from the lithosphere above 350 km these upper mantle anomalies are interpreted as the remnants of delaminated sub continental lithosphere that went down when greater india continued to converge northward with asia after 45 ma the deeper high velocity anomalies under the indian sub continent appear clearly separated from the shallower ones as well as from each other and are inferred to represent remnants of oceanic lithospheric slabs that have sunk into the lower mantle and were subsequently overridden by the indian plate they occur at depths between 1000 and 2300 km and occasionally descend down to the core mantle boundary the anomalies form three parallel wnw ese striking zones we interpret the two southern zones as remnants of oceanic lithosphere that was subducted when the neo tethys ocean closed between india and tibet in the cretaceous and earliest tertiary the northern deep mantle zone under northern afghanistan the himalayas and the lhasa block in southern tibet may represent the last subducted remnant of the paleo tethys ocean which is thought to have closed before the hauterivian stage of the early cretaceous the middle zone continues southeastward as a rather straight high velocity zone towards sumatra where it becomes convex southward and parallel to the subduction zone under the sunda arc comparison of this straight middle zone near india with the shallower upper 600 1000 km northern zone which displays a cusp like shape near the yunnan sw china syntaxis of the eastern himalayas supports the notion that the shallow northern zone represents later subduction than the deeper middle zone the suggestion of a counterclockwise rotation 20 of the indian plate during tertiary indentation of asia is supported by these features the present day latitudes of 5 35 n of the deep slabs under india and adjacent areas correspond to the approximate paleolatitudes of the cretaceous subduction zones the slab remnants in the middle mantle occur therefore near the ancient locations where they started their downward journey which implies that lateral movements in the deeper mantle were not large

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

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

neogene and quaternary rollback evolution of the tyrrhenian sea the apennines and the sicilian maghrebides reconstruction of the evolution of the tyrrhenian sea shows that the major stage of rifting associated with the opening of this basin began at similar to10 ma it involved two episodes of back arc extension which were induced by the rollback of a west dipping subducting slab the first period of extension 10 6 ma was prominent in the northern tyrrhenian sea and in the western part of the southern tyrrhenian sea the second period of extension mainly affected the southern tyrrhenian sea began in the latest messinian 6 5 ma and has been accompanied by subduction rollback at rates of 60 100 km myr 1 slab reconstruction combined with paleomagnetic and paleogeographic constraints indicates that in the central apennines the latest messinian 6 5 ma arrival of a carbonate platform at the subduction zone impeded subduction and initiated a slab tear and major strike slip faults these processes resulted in the formation of a narrow subducting slab beneath the ionian sea that has undergone faster subduction rollback and induced extreme rates of back arc extension

scotia arc kinematics from gps geodesy 1 gps crustal velocity data from the scotia and south sandwich plates transform azimuths spreading data and an updated earthquake slip vector catalog provide the first scotia and south sandwich plate euler vector estimates not dependent on closure as the gps data tie them to the global plate circuit neither the gps data which sample limited portions of the plates nor the geologic data which are not tied to the global spreading circuit are sufficient individually to define the euler vectors as scotia plate gps measurements do not sample the stable plate interior plate boundary deformation field modeling is necessary for euler vector estimation our south america antarctic and scotia south sandwich euler pole estimates agree with previous estimates from either gps or geologic data our south america scotia euler vector however is significantly different and near the south america antarctic euler vector producing an approximately coaxial motion of scotia between south america and antarctica

the relation between grain size and 40ar 39ar date for alpine white mica from the siviez mischabel nappe switzerland abstract variations in alpine white mica 40 ar 39 ar dates from the cover units of the siviez mischabel nappe relate to regional variations in the thermal history of the nappe we focus on three regions within the nappe the central siviez mischabel csm the southern siviez mischabel ssm and the eastern siviez mischabel esm our approach weaves together observations of quartz and mica textures in thin section the variation of 40 ar 39 ar date with grain size considerations of the effective diffusion dimension edd of argon in white mica and a comparison of dates with diffusion model results in the csm pressure solution of quartz and dislocation glide in mica accommodated alpine deformation dates record mica growth during nappe emplacement from 40 to 36 ma and do not vary with grain size in the ssm and esm both mica and quartz show textures associated with dynamic recrystallization and dates decrease with grain size in the ssm dates also agree with the timing of nappe emplacement but in the esm dates significantly post date the timing of nappe emplacement a comparison of dates with diffusion model results supports inferences from rock fabrics that the ssm experienced higher peak temperatures than the csm even though dates from both units approximate the timing of mica growth dates obtained from the esm however do not compare well with simple models and the thermal evolution of this region of the nappe in the neighborhood of the simplon fault zone is not well understood

kinematics and flow patterns in deep mantle and upper mantle subduction models influence of the mantle depth and slab to mantle viscosity ratio 1 three dimensional fluid dynamic laboratory simulations are presented that investigate the subduction process in two mantle models an upper mantle model and a deep mantle model and for various subducting plate mantle viscosity ratios sp m 59 1375 the models investigate the mantle flow field geometrical evolution of the slab sinking kinematics and relative contributions of subducting plate motion and trench migration to the total rate of subduction all models show that the subducting plate is always moving trenchward resulting from slab pull furthermore all deep mantle models show trench retreat as do upper mantle models in the initial stage of subduction before slab tip transition zone interaction upper mantle models with a low sp m 66 217 continue to show trench retreat after interaction upper mantle models with a high sp m 378 709 show a period of trench advance after interaction followed by trench retreat upper mantle models with a very high sp m 1375 show continued trench advance after interaction the difference in trench migration behavior and associated slab geometries is attributed to both sp m and the mantle depth to plate thickness ratio tm tsp which both affect the slab bending radius to mantle thickness ratio rb tm four subduction regimes can be defined regime i with rb tm 0 3 trench retreat slab draping and a concave trench regime ii with 0 3 rb tm 0 5 episodic trench migration slab folding and a concave trench regime iii with rb tm 0 5 trench advance slab rollover geometries and minor trench curvature and regime iv with rb tm 0 8 trench retreat slab draping and a rectilinear trench in all models slab parallel downdip motion induces poloidal mantle flow structures in addition trench retreat and rollback motion of the slab induce quasi toroidal return flow around the lateral slab edges toward the mantle wedge rollback induced poloidal flow around the slab tip is not observed in any of the experiments finally comparison between the slab geometries observed in the upper mantle models and slab geometries observed in nature imply that the effective viscosity ratio between slab and ambient upper mantle in nature is less than 103 and of the order 1 7 102 with a best estimate of 1 3 102

arc parallel extension and localization of volcanic complexes in guadeloupe lesser antilles 1 subduction of atlantic seafloor under the caribbean plate causes shallow earthquakes within the lesser antilles volcanic arc such earthquakes above the subduction interface show strike slip or normal fault plane solutions the latter with e w striking nodal planes to better assess seismic hazard and the coupling between volcanism and tectonics we investigated faulting related to overriding plate deformation in the guadeloupe archipelago using aerial photographs satellite spot images and topographic maps 1 25000 scale we mapped active and middle to late pleistocene fissures and normal fault systems that cut the uplifted coral platforms grande terre and marie galante and the volcanic rocks of basse terre the available marine geophysical data show that the faults extend offshore to bound submarine rifts the e w striking 1500 m deep v shaped marie galante rift separates the two islands of marie galante and grande terre normal faults in the north of grande terre appear to mark the similarly v shaped western termination of the 5000 m deep n 50e to n130 e striking desirade graben three shallow m 5 5 earthquakes 6 may 1851 29 april 1897 3 august 1992 appear to have ruptured segments of the marie galante rift boundary faults the young la grande decouverte volcanic complex of basse terre including the 1440 a d soufriere dome lies within the western termination of the marie galante rift the ancient volcanic shoulders of the rift buttress the active dome to the north and south which may explain why major prehistoric sector collapses and pyroclastic avalanches have been directed southwestward into the caribbean sea or southeastward into the atlantic ocean the marie galante rift is typical of other troughs transverse to the northeastern edge of the caribbean plate we interpret such troughs which are roughly orthogonal to the arc to result from slip partitioning and extension perpendicular to plate convergence that they disappear southward implies that they result from interaction between the caribbean and north american plates