joint inversions of seismic and geodynamic data for models of three dimensional mantle heterogeneity the seismic models of three dimensional 3 d mantle heterogeneity may be interpreted in terms of the density perturbations which drive mantle flow and thus provide important predictions of flow related observables such as the nonhydrostatic geoid the current models of global scale shear velocity heterogeneity provide reasonably good fits to the long wavelength nonhydrostatic geoid data 60 70 variance reductions but rather poor fits to the corresponding free air gravity anomalies 30 40 variance reductions this major difference is due to two factors 1 the very different amplitude spectrum of the gravity anomalies which is nearly flat for degrees l 3 the largest mismatch between the pattern of predicted and observed gravity anomalies is in the southern hemisphere this observation suggests that one reason for the poor overall match between the two fields may be that the global seismic data are not resolving sufficiently well the heterogeneity in the southern hemisphere portion of the deep mantle in contrast the gravity data provide accurate and geographically uniform constraints on the vertically integrated heterogeneity in the mantle we therefore perform a series of experiments in which we simultaneously invert a large set of seismic data which includes long period waveforms ss s and scs s differential travel times and normal mode structure coefficients and the long wavelength gravity anomaly data we thus determine whether it is possible to derive new 3 d heterogeneity models which satisfy both data sets the gravity anomaly data are interpreted in the context of spherically symmetric viscous flow models of the mantle in these inversion experiments we test several radial viscosity and ln ln profiles and thereby assess their plausibility the joint seismic geodynamic inversions reveal that it is indeed possible to greatly improve the fit to the free air gravity anomalies with variance reductions of 80 90 readily accessible while preserving the fit to the seismic data this improvement is achieved with some adjustment to the heterogeneity in the depth range 1500 2500 km where the seismic data constraints appear to be weakest the joint inversions also reveal new structures in the southern hemisphere portion of the lower mantle which apparently are not resolved by the seismic data alone

mesozoic plate motion history below the northeast pacific ocean from seismic images of the subducted farallon slab the high resolution seismic imaging of subducted oceanic slabs1 2 has become a powerful tool for reconstructing palaeogeography3 the images can now be interpreted quantitatively by comparison with models of the general circulation of the earth s mantle4 here we use a three dimensional spherical computer model of mantle convection5 6 to show that seismic images of the subducted farallon plate provide strong evidence for a mesozoic period of low angle subduction under north america such a period of low angle subduction has been invoked independently to explain rocky mountain uplift far inland from the plate boundary during the laramide orogeny7 the computer simulations also allow us to locate the largely unknown kula farallon spreading plate boundary the location of which is important for inferring the trajectories of suspect terrain across the pacific basin8

segmented african lithosphere beneath the anatolian region inferred from teleseismic p wave tomography summary lithospheric deformation throughout anatolia a part of the alpine himalayan orogenic belt is controlled mainly by collision related tectonic escape of the anatolian plate and subduction roll back along the aegean subduction zone we study the deeper lithosphere and mantle structure of anatolia using teleseismic finite frequency p wave traveltime tomography we use data from several temporary and permanent seismic networks deployed in the region approximately 34 000 p wave relative traveltime residuals measured in multiple frequency bands are inverted using approximate finite frequency sensitivity kernels our tomograms reveal segmented fast seismic anomalies beneath anatolia that corresponds to the subducted portion of the african lithosphere along the cyprean and the aegean trenches we identify these anomalies as the subducted aegean and the cyprus slabs that are separated from each other by a gap as wide as 300 km beneath western anatolia this gap is occupied by slow velocity perturbations that we interpret as hot upwelling asthenosphere the eastern termination of the subducting african lithosphere is located near the transition from central anatolia to the eastern anatolian plateau or arabian eurasian collision front that is underlain by large volumes of hot underplating asthenosphere marked by slow velocity perturbations our tomograms also show fast velocity perturbations at shallow depths beneath northwestern anatolia that sharply terminates towards the south at the north anatolian fault zone nafz the associated velocity contrast across the nafz persists down to a depth of 100 150 km hence our study is the first to investigate and interpret the vertical extent of deformation along this nascent transform plate boundary overall the resolved upper mantle structure of anatolia is directly related with the geology and tectonic features observed at the surface of the anatolian plate and suggest that the segmented nature of the subducted african lithosphere plays an important role in the evolution of anatolia and distribution of its tectonic provinces

trench migration and subduction zone geometry seismic tomography results show a large variety in upper mantle structure along convergent plate boundaries we numerically investigate the effect of trench migration on the evolution of a slab with temperature dependent viscosity encountering a viscosity interface we find that subduction zone geometry is sensitive to even small rates 1 cm yr of retrograde motion increase in trench migration rate decreases the dip angle of the slab and its ability to penetrate the lower mantle upon including a background mantle flow it turns out that trench migration relative to the upper mantle flow is more decisive than the absolute plate velocities subduction zone geometry appears to be strongly time dependent we conclude that individual tectonic setting and time dependent slab behaviour can account for many different types of observable subduction zone anomalies our model temperature can adequately account for magnitudes and patterns of seismic anomalies as obtained from seismic tomography

hotspots mantle plumes flood basalts and true polar wander persistent long lived stationary sites of excessive mantle melting are called hotspots hotspots leave volcanic trails on lithospheric plates passing across them the global constellation of fixed hotspots thus forms a convenient frame of reference for plate motions through the orientations and age distributions of volcanic trails left by these melting anomalies hotspots appear to be maintained by whole mantle convection in the form of upward flow through narrow plumes evidence suggests that plumes are deflected little by horizontal flow of the upper mantle mantle plumes are largely thermal features and arise from a thermal boundary layer most likely the mantle layer just above the core mantle boundary experiments and theory show that gravitational instability drives flow beginning with the formation of diapirs such a diapir will grow as it rises fed by flow through the trailing conduit and entrainment of surrounding mantle the structure thus develops a large spherical plume head and a long narrow tail on arrival at the base of the lithosphere the plume head flattens and melts by decompression producing enormous quantities of magma which erupt in a short period these are flood basalt events that have occurred on continents and in ocean basins and that signal the beginning of major hotspot tracks the plume supported hotspot reference frame is fixed in the steady state convective flow of the mantle and is independent of the core generated axial dipole paleomagnetic reference frame comparison of plate motions measured in the two frames reveals small but systematic differences that indicate whole mantle motion relative to the earth s spin axis this is termed true polar wander and has amounted to some 12 since early tertiary time the direction and magnitude of true polar wander have varied sporadically through the mesozoic probably in response to major changes in plate motions particularly subduction zone location that change the planet s moments of inertia

signature of slab fragmentation beneath anatolia from full waveform tomography when oceanic basins close after a long period of convergence and subduction continental collision and mountain building is a common consequence slab segmentation is expected to have been relatively common just prior to closure of other oceans in the geological past and may explain some of the complexity that geologists have documented in the tibetan plateau also we focus on the eastern mediterranean basin which is the last remainder of a once hemispherical neo tethys ocean that has nearly disappeared due to convergence of the india and africa arabia plates with the eurasia plate we present new results of full waveform tomography that allow us to image both the crust and upper mantle in great detail we show that a major discontinuity exists between western anatolia lithosphere and the region to the east of it also the correlation of geological features and the crustal velocities is substantially stronger in the west than in the east we interpret these observations as the imprint in the overriding plate of fragmentation of the neo tethys slab below it this north dipping slab may have fragmented following the eocene about 35 million years ago arrival of a continental promontory central anatolian core complex at the subduction contact from the eocene through the miocene slab roll back ensued in the aegean and west anatolia while the cyprus bitlis slab subducted horizontally beneath central and east anatolia following collision of arabia about 16 million years ago the cyprus bitlis slab steepened exposing the crust of central and east anatolia to high temperature and resulting in the velocity structure that we image today slab fragmentation thus was a major driver of the evolution of the overriding plate as collision unfolded

slab temperature and thickness from seismic tomography 1 method and application to tonga delay times of compressional body wave phases from both teleseismic and local events are used to invert for a high resolution p wave velocity model in the tonga subduction zone the images obtained show a high velocity subducting slab with velocity deviations of the order of 3 4 assuming to first order that the positive velocity anomalies within the slab are caused by a temperature effect a theoretical slab temperature model based on the diffusion equation is used to explain velocity anomalies within the tomographic slab temperature differences between the interior of the slab and the ambient mantle are converted to velocity perturbations using the scaling parameter dvp dt 4 8 10 4 km s 1 c 1 for lithosphere material the optimal values for the parameters in the temperature model are found using a nonlinear optimization that compares the integrated velocity anomalies in the tomographic slab region to integral of high velocities in a synthetic slab derived from a temperature model the parameters for slab thickness and mantle potential temperature are not uniquely determined therefore a fixed value for the mantle potential temperature based on laboratory values for the temperature of the spinel to perovskite transition at 660 km is used using 1180 c as the potential temperature the theoretical temperature model gives an optimal slab thickness of 82 km for a region near 29 s in tonga the uncertainty in the thickness is dominated by the uncertainty in the mantle temperature and would be 8 km for an uncertainty of 100 in mantle temperature but nonsystematic errors are less in order to enhance the tomographic result the velocity model is biased towards the theoretical slab model however a posteriori changes made to the tomogram will most likely violate the fit to the delay time data to prevent this the difference between the tomogram and the predicted slab model is projected onto the null space of the inversion to remove components which do not satisfy the seismic data using only null space components to modify the minimum norm solution an enhanced model is obtained which has been biased toward the theoretical solution but has the same data misfit as the minimum norm solution the final image shows a very narrow and continuous slab with maximum velocity anomalies of the order of 6 7 many of the gaps within the slab as well as artifacts around the slab which were present in the minimum norm solution are absent in the biased image

sprite initiation altitude measured by triangulation 1 high time resolution 10 000 frames per second images of sprites combined with multistation concurrent video recordings have provided data for triangulation of the altitude of the initial sprite onset the high speed images were obtained from the langmuir laboratory new mexico during summer campaigns in 2007 and 2008 with video observations from sites at portales new mexico and las vegas new mexico sprites start with one or more downward propagating streamer heads the triangulated onset altitudes of this initial downward streamer vary between 66 and 89 km in some sprites the downward streamers are followed a little later by upward propagating streamers the upward streamers start from a lower altitude and existing luminous sprite structures and their triangulated altitudes vary from 64 to 78 km the downward streamers create c sprite characteristics while the upward streamers form the broad diffuse tops of carrot sprites in the sprites analyzed the higher onset altitudes for the downward propagating initial streamers were associated with c sprites and the lower with carrot sprites but our larger data set indicates that this is not generally the case it appears that the dominant sprite types vary from year to year indicating that some longer lasting environmental parameter such as mesospheric conductivity and composition or thunderstorm cloud dynamics may play an important role in determining the types of sprites observed

relative hotspot motions versus true polar wander abstract the fixity of hotspots and mantle plume locations has long been axiomatic if the assumption of fixed hotspots is granted absolute plate motions and movements of the spin axis with respect to the hotspot framework defined by some as true polar wander tpw can be determined however this assumption can be tested by paleomagnetic data and such tests are gradually raising some doubts about the fixity of hotspots the result is that discrepancies between cretaceous and tertiary hotspot and paleomagnetic reference frames are now beginning to be interpreted as the result of plume drift within a convective mantle in the indo atlantic hotspots have remained relatively stationary with respect to the spin axis for the last 95 million yr however the pacific hotspots notably hawaii appear to have undergone large scale southward drift with respect to the spin axis during the early tertiary global paleomagnetic data do not indicate that any tpw occurred during the late cretaceous or tertiary although the early cretaceous paleomagnetic and hotspot frames for the indo atlantic realm can be interpreted as slow tpw direct estimates of paleolatitude and hotspot motion in particular the kerguelen hotspot challenge tpw as a global phenomenon at present we consider that the large early cretaceous discrepancy between hotspot and paleomagnetic data is best explained by southward drift of the atlantic hotspots prior to 95 ma

sensitivities of seismic velocities to temperature pressure and composition in the lower mantle we calculated temperature pressure and compositional sensitivities of seismic velocities in the lower mantle using latest mineral physics data the compositional variable refers to the volume proportion of perovskite in a simplified perovskite magnesiowustite mantle assemblage the novelty of our approach is the exploration of a reasonable range of input parameters which enter the lower mantle extrapolations this leads to realistic error bars on the sensitivities temperature variations can be inferred throughout the lower mantle within a good degree of precision contrary to the uppermost mantle modest compositional changes in the lower mantle can be detected by seismic tomography with a larger uncertainty though a likely trade off between temperature and composition will be largely determined by uncertainties in tomography itself given current sources of uncertainties on recent data anelastic contributions to the temperature sensitivities calculated using karato s approach appear less significant than previously thought recent seismological determinations of the ratio of relative s to p velocity heterogeneity can be entirely explain by thermal effects although isolated spots beneath africa and the central pacific in the lowermost mantle may ask for a compositional origin 2001 elsevier science b v all rights reserved

magma generation by slab steepening and breakoff beneath a subduction accretion complex an alternative model for collision related volcanism in eastern anatolia turkey 1 eastern anatolia is regarded as one of the best examples of continental collision represented by a 2 km high plateau it displays shallow and diffuse seismicity indicating that the crust is still being actively deformed perhaps the most striking aspect of eastern anatolia is the volume and compositional variability of collision related volcanic products erupted during the neogene and quaternary time a number of models have been proposed for the genesis of this volcanism in the region however recent geophysical studies have revealed that a mantle lithosphere is almost completely absent beneath a greater portion of the region which makes us question the validity of the previous models in this paper we propose a new model for magma genesis slab steepening and breakoff beneath a large subduction accretion complex this model holds that a northward subducting oceanic lithosphere beneath the eastern anatolia accretionary prism gets steepened and eventually detached from the continental lithosphere of the bitlis poturge massif following the continent accretionary complex collision this brings the asthenospheric mantle in contact with the accretionary complex at shallow depths e g around 50 km resulting in rapid regional uplift and extensive melting which leads to collisional volcanism

large scale mantle convection and the history of subduction a large scale pattern of mantle convention is evident in the shape of the earth s gravity field1 lateral variations in seismic velocity2 3 the distribution of hotspots4 and long wavelength bathv metric anomalies5 paradoxically these signals of mantle convection bear little resemblance to present day plate motions1 4 although they seem to be related to past subduction zone configurations6 here we quantify this relationship by correlating a time integrated history of subduction calculated in the hotspot reference frame for the past 180 myr with regions of the deep mantle inferred from other geophysical observations to be colder than average our results support suggestions6 9 that a direct relationship exists between large scale thermal heterogeneity in the lower mantle and subduction during the cenozoic and mesozoic

the role of slab detachment processes in the opening of the western central mediterranean basins some geological and geophysical evidence abstract a review of the geological and geophysical data from the central and western mediterranean region and the present day upper mantle structure derived from tomographic studies are utilized in order to define the oligocene recent geodynamic evolution for the area in line with previous work we suggest that the miocene quaternary opening of the western and central mediterranean basins is the result of back arc extension due to the roll back toward the southeast of a northwestward subducting african slab in a geodynamic setting pinned between the alpine and betic collisional zones we find however that this general pattern is complicated by four different detachment events which occurred beneath the alps early oligocene the betic chain aquitanian northern africa langhian and the apennines late miocene pliocene we show that each of these events determines a major tectonic reorganization within the european plate

slabs in the lower mantle results of dynamic modelling compared with tomographic images and the geoid abstract the recent increase in resolution of tomographic images of the earth s interior has enabled us to see slabs in the lower mantle on the other hand the distribution of slabs can be inferred from subduction history using a dynamical model of mantle flow comparison of tomographic images with model results can help to distinguish between alternative models of subduction history tomographic models and mantle flow models and thus improve our understanding of the earth here a simple dynamic model of subduction is used and some results are presented amounts and locations of subduction are inferred from published models of plate motion and plate boundaries for the past 120 ma the latter have been interpolated on 2 ma time intervals mantle flow driven by the density anomalies corresponding to subducted slabs is calculated with the method of hager and o connell for a viscosity structure that features an increase in viscosity from 4 10 20 pa below the lithosphere to 4 10 22 pa in the lower part of the lower mantle slabs are advected in the flow field in the model slabs sink on average about 1700 km through the mantle during 100 ma however the calculated average horizontal motion of slabs throughout most of the mantle is only about 600 700 km for slabs that have sunk to the lowermost mantle it is about 1000 km this is in accordance with the previously observed good correlation between subduction locations and lowermost mantle heterogeneities slabs sink fastest in regions with abundant subduction and little trench migration e g below east asia in the lowermost mantle predicted slab locations tend to lie in areas with high seismic velocities however there are large areas of fast seismic velocities where no slabs are predicted and which therefore presumably correspond to earlier subduction in some regions e g below japan our model approximately reproduces the observed shape of the slab in other regions where the agreement is poor the results can help to further constrain models of subduction history and associated mantle flow additionally the predicted geoid is computed and compared with the actual geoid predicted and observed geoid have similar magnitude and a number of features in common however an optimization of the fit is not attempted here since in order to do this also the hot upwellings from the lower boundary layer presumably in the form of plumes would have to be included

4 d evolution of se asia s mantle from geological reconstructions and seismic tomography abstract how the collision between india and asia is related to processes deeper in the mantle is unclear here we compare geological reconstructions of block motions within asia since 50 ma with the tomographically imaged three dimensional 3 d morphology of subducted lithosphere to obtain insight into the spatiotemporal evolution of mantle structure past positions of the convergent margin show remarkable similarities with slab geometry at specific depths the striking change in slab geometry from a linear structure beneath 1100 km to an increasingly distorted shape at depths of less than 700 km results from collision the slab contours match the progressive deformation of asia s margin including india s indentation and sundaland s extrusion ever since the onset of collision the indian plate appears to have overridden its own sinking mantle and it does not seem at present to underthrust tibet significantly north of the zangbo suture if correct this observation would provide further evidence against models of plateau build up involving indian lithosphere the tomographic images beneath india confirm that asian deformation has absorbed at least 1500 km of convergence since collision began from the match between the southeastward motion of sundaland between 40 and 20 ma and the principal change in slab structure between 700 and 1100 km depths we infer that lateral advection in the mantle is small and that the sinking rate beneath sunda was 2 cm yr in the lower mantle and 5 cm yr above the transition zone

large scale motion between pacific and atlantic hotspots studies of true polar wander tpw the rotation of the solid earth with respect to the spin axis1 have suggested that there has been 10 15 of relative motion over the past 130 myr refs 2 4 in such studies the orientation of the spin axis is recovered from continental palaeomagnetic poles corrected for relative plate motions and compared with a deep mantle reference frame defined by hotspot locations but deducing relative plate motions becomes increasingly difficult for older mesozoic time periods hindering tests of tpw on timescales comparable to those of large scale mantle convection moreover the assumption of hotspot fixity is controversial5 6 we examine here a more direct approach7 8 using palaeolatitudes derived from pacific guyots contrary to predictions from tpw models these data suggest only minor latitudinal shifts of pacific hotspots during the cretaceous period instead of tpw relative motion between the atlantic and pacific hotspot groups9 is required at a velocity of approximately 30 mm yr 1 more than 50 larger than previously proposed5

geodynamic implications of moving indian ocean hotspots abstract the reconstructed positions of many of the indian ocean intraplate volcanic features are inconsistent with the present positions of the hotspots to which they have been ascribed the discrepancy can be explained by the past motion of the hotspots themselves and in the errors introduced in an absolute plate model based on fixed hotspots previous models for the motion of hotspots in the indian ocean e g steinberger and o connell in the history and dynamics of global plate motions 2000 pp 377 398 and geophys j int 132 1998 412 434 antretter et al earth planet sci lett 203 2002 635 650 have had success in matching paleolatitudes of the kerguelen hotspot but not the marion or reunion hotspots we calculate the motion of the indian ocean hotspots using a method described in steinberger and o connell op cit employing a plate motion model revised using identification of mesozoic magnetic anomalies in the enderby basin gaina et al in the evolution and dynamics of the australian plate in press we find a motion of 7 south for the kerguelen plume similar to previous models antretter et al op cit this motion alleviates the discrepancies in the position of the kerguelen plume and its purported involvement in the formation of the ninetyeast ridge rajmahal traps bunbury basalts and assorted volcanics of the western australian margin we find a motion of 5 north and 7 south for the reunion and marion plumes respectively the motion of the reunion hotspot is consistent with paleolatitude estimates the apparent fixity of the marion hotspot torsvik et al earth planet sci lett 164 1998 221 232 is not well constrained but may be due to the combination of hotspot motion and true polar wander acting in opposing directions we present revised absolute finite rotations for the indian plate for the last 65 ma based on the motion of these hotspots which place india farther north in the past than for fixed hotspot models

a slab breakoff model for the neogene thermal evolution of south karakorum and south tibet on the south karakorum margin neogene high temperature medium pressure ht mp gneisses define an east west trending thermal anomaly these rocks have been heated from 600 to 750 c during a slight pressure drop from 0 7 to 0 5 gpa their retrogressive path cross cuts the relaxed geotherm of tectonically thickened crust such a p t evolution occurs only if an advective source of heat is involved involvement of an advective heat source is also implied by the occurrence of neogene granitoids and lamprophyres within the ht mp gneiss area these rocks are strongly enriched in large ion lithophile elements relative to primitive mantle and show negative high field strength element anomalies we interpret these geochemical characteristics to be the result of melting of metasomatized asian lithospheric mantle the nd and sr isotopic compositions of the south karakorum neogene magmatic rocks nd 12 to 7 and 87sr 86sr 0 705 0 725 further suggest they could have originated from mixing between asian variously metasomatized mantle and precambrian crust by contrast the origin of the youngest magmatic rocks 10 myr here exemplified by the hemasil syenite and associated lamprophyres requires involvement of a depleted mantle the combined hf nd signature of these rocks hf 10 4 11 5 and nd 3 4 4 3 suggests that the source of the hemasil syenite could have been depleted mantle contaminated by oceanic sediments likely during the earlier subduction of the tethyan ocean neogene magmatic rocks with the same geochemical characteristics and evolution as those of south karakorum have previously been described in south tibet based on their location and the geochemical evolution of their source region we here propose that the neogene magmatic and metamorphic evolution of the south asian margin was controlled by slab breakoff of the subducting indian continental margin starting at about 25 ma this model is supported by available geophysical data from south karakorum and south tibet

the post eocene evolution of the doruneh fault region central iran the intraplate response to the reorganization of the arabia eurasia collision zone

mantle conveyor beneath the tethyan collisional belt abstract collisional belts are generated by the arrival of continental lithosphere into a subduction zone the tethyan suture from the bitlis to the himalayas is a prime example where the arabian and indian plates collided with eurasia during the cenozoic while the kinematics of this process are well established its dynamics are more uncertain india and arabia intriguingly keep advancing in spite of large collisional resisting forces and in the absence of a substantial upper mantle slab driving force at present day we perform global mantle circulation computations to test the role of deep mantle flow as a driving force for the kinematics of the tethyan collisional belt evaluating different boundary conditions and mantle density distributions as inferred from seismic tomography or slab models our results show that mantle drag exerted on the base of the lithosphere by a large scale convective conveyor belt with an active upwelling component is likely the main cause for the ongoing indentation of the indian and arabian plates into eurasia

paleomagnetic results along the bitlis zagros suture zone in se anatolia turkey implications for the activation of the dead sea fault zone abstract the tectonic evolution of se anatolia is associated with the collision between the taurides and the arabian platform in the early cenozoic after the final closure of the southern branch of the neotethys ocean during late cretaceous to early eocene times the ongoing deformation was characterized by northward movement of arabia followed by westward extrusion of anatolia having been displaced along the northern and eastern anatolian transform faults the neotectonic deformation history between the anatolian plate and the arabian platform since the closure of the southern neotethyan oceanic basin from late cretaceous to present has not been studied in detail paleomagnetically therefore we carried out a paleomagnetic study of mesozoic and cenozoic rocks at 34 localities along the bitlis zagros suture zone in se anatolia located in the boundary region between the taurides and the arabian platform and delimited by the dead sea fault paleomagnetic results from the area show clockwise rotations of 35 6 14 1 and 26 7 8 0 with respect to the african plate from upper jurassic lower cretaceous and upper cretaceous limestones respectively the paleomagnetic results indicate also an earlier deformation opening due to the dead sea fault zone dsfz farther north on the tauride block counterclockwise rotation of 34 5 12 5 is recorded by upper cretaceous limestones in the se of the study area on the northern margin of arabia however smaller counterclockwise rotations of 12 3 8 3 and 2 5 10 0 in late cretaceous and 12 7 13 6 in middle eocene rocks show relative tectonic stability middle miocene rocks close to the suture zone show counterclockwise rotation by about 39 4 10 9 while farther away from the suture zone counterclockwise rotations of 21 9 9 9 and 7 5 5 8 are observed paleolatitudinal motion from south to north from the late jurassic to the present is in concordance with the apparent polar wander path of eurasia and africa

a laboratory investigation of effects of trench migration on the descent of subducted slabs a laboratory investigation of viscous slabs subducted from a migrating trench reveals a range of possible behaviours the slab dip in a uniform mantle is found to be steady or oscillatory depending on the rates of descent and trench migration in addition density and viscosity interfaces are used to model the increased resistance to sinking of stiff slabs through the seismic discontinuity at a depth of about 660 km in earth s mantle if the slab is denser than the lower layer and its dip in the upper layer is steady it can continue to descend through the lower layer in a tabular form and in either a steady or oscillatory manner or it can be laid horizontally on the interface for a distance before descending in a chain of diapirs resulting from gravitational instability at the interface if the slab dip in the upper layer is unstable the slab sinks into the lower layer in a chain of large diapirs at a spacing determined by the frequency of oscillations set by instability of the slab within the upper layer the style of penetration depends on the trench migration speed and the ratio of sinking velocities in the two layers estimates for mantle slabs indicate that they may range across the major regime transitions the experimental system provides a gross simplification of mantle conditions but the results are a testimony to the possibility of a range of complex behaviour of subducted lithosphere they also indicate that the relationship between the structural features in the lower mantle revealed by seismic imaging and present day tectonic processes at the surface may not be obvious tomographic images of western pacific subduction zones and data for the migration rates of associated trenches suggest a dependence of slab behaviour on migration rate similar to that seen in the experiments

evolution of the downgoing lithosphere and the mechanisms of deep focus earthquakes summary the thermal evolution of the lithospheric slab at subduction zones and its geophysical effects are numerically calculated an alternating direction implicit finite difference scheme is used to compute the thermal models taking into account all heating sources and phase boundaries these models with the appropriate spreading rates and dip angles are compared with different island arc systems temperatures inside the slab are strongly controlled by the conductivity and by the time elapsed since the initiation of descent the depth to which temperature anomalies persist is generally about 700 km or less the thermal results are used to construct the seismic velocities and ray paths density anomalies and the resulting stress distribution comparing the theoretical stress distribution and the focal mechanism studies of intermediate and deep focus earthquakes indicates the importance of both the mantle s rheology and the temperature dependence of the slab s elastic properties the intermediate and deep focus earthquakes are located along the coolest region of the slab the theoretical results explain the source mechanisms and the orientation of principal stresses under major island arcs

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

tearing of stagnant slab subducted slabs of oceanic lithosphere below the western pacific tend to be stagnant in the transition zone with poorly known mechanical properties typical examples are the izu bonin and japan slabs that meet each other to form a cusplike junction beneath southwest japan here we show that these two slabs are torn apart at their junction when they bend to flatten over the 660 kilometer discontinuity as is expected from a simple geometric argument we present three lines of evidence for this ongoing slab tear

reconciling the geological history of western turkey with plate circuits and mantle tomography abstract we place the geological history since cretaceous times in western turkey in a context of convergence subduction collision and slab break off to this end we compare the west anatolian geological history with amounts of africa europe convergence calculated from the atlantic plate circuit and the seismic tomography images of the west anatolian mantle structure western turkish geology reflects the convergence between the sakarya continent here treated as eurasia in the north and africa in the south with the anatolide tauride block atb between two strands of the neotethyan ocean convergence between the sakarya and the atb started at least 95 90 myr ago marked by ages of metamorphic soles of ophiolites that form the highest structural unit below sakarya these are underlain by high pressure low temperature metamorphic rocks of the tavsanli and afyon zones and the oren unit which in turn are underlain by the menderes massif derived from the atb underthrusting of the atb below sakarya was since 50 ma associated with high temperature metamorphism and widespread granitic magmatism thrusting in the menderes massif continued until 35 ma after which there is no record of accretion in western turkey plate circuits show that since 90 ma 1400 km of africa europe convergence occurred of which 700 km since 50 ma and 450 km since 35 ma seismic tomography shows that the african slab under western turkey is decoupled from the african plate this detached slab is a single coherent body representing the lithosphere consumed since 90 ma there was no subduction re initiation after slab break off atb collision with europe therefore did not immediately lead to slab break off but instead to delamination of subducting lithospheric mantle from accreting atb crust while staying attached to the african plate this led to asthenospheric inflow below the atb crust high temperature metamorphism and felsic magmatism slab break off in western turkey probably occurred 15 myr ago after which overriding plate compression and rotation accommodated ongoing africa europe convergence slab break off was accommodated along a vertical ne trending subduction transform edge propagator step fault zone accelerating southwestward slab retreat of the aegean slab the se aegean slab edge may have existed already since early miocene times or before but started to rapidly roll back along the southeastern aegean step in middle miocene times penetrating the aegean region in the pliocene

subduction history of the tethyan region derived from seismic tomography and tectonic reconstructions 1 in the mantle underneath the tethyan suture zone large volumes of positive velocity anomalies have been imaged by seismic tomography and interpreted as the present day signature of subducted tethyan lithosphere we investigate the mesozoic cenozoic subduction history of the region by integrating independent information from mantle tomography and tectonic reconstructions three different subduction scenarios for the tethyan oceanic lithosphere representative for the available tectonic reconstructions are used to predict the present thermally anomalous volumes associated with the lithospheric surface subducted since the late mesozoic next these predicted thermal volumes and their expected positions are compared to the relevant anomalous volumes derived from seismic tomographic images in this analysis we include among others the possible effects of ridge subduction and slab detachment after the cenozoic continental collisions absolute plate motion and slab thickening in the mantle our preferred subduction model comprises the opening of large back arc oceanic basins within the eurasian margin the model points to slab thickening by a factor of 3 in the mantle in which case the estimated volumes allow for active oceanic spreading 1 2 5 cm yr in the tethyan lithosphere during convergence our results further indicate the occurrence of early oligocene slab detachment underneath the northern zagros suture zone followed by both westward and eastward propagation of the slab tear and diachronous eocene to miocene slab detachment below the eastern to western himalayas free sinking rates of the detached material of 2 cm yr in the lower mantle provide the best fit between the tomographic mantle structure and our tethyan subduction model

slab breakoff a model of lithosphere detachment and its test in the magmatism and deformation of collisional orogens we present a model proposing that oceanic lithosphere detaches from continental lithosphere during continental collision slab breakoff allowing an explanation of syn to post collisional magmatism and metamorphism continental collisions are preceded by subduction of dense oceanic lithosphere and followed by attempted subduction of buoyant continental lithosphere this situation of opposing buoyancy forces leads to extensional deformation in the subducted slab a narrow rifting mode of deformation will result if there is strain localization slab breakoff results we have assessed the plausiblity of this process by quantitatively evaluating an upper bound for the strength of the lithosphere and have compared it with the change in buoyancy force during continental sbduction whether breakoff will occur and the depth at which it will occur is a strong function of temperature and hence the subduction velocity for a subduction velocity of 1 cm yr breakoff could occur at depths of between 50 and 120 km while at higher velocities it is still likely to occur but at deeper depths as a result of the rifting during breakoff the asthenosphere upwells into the narrow rift and following breakoff it impinges on the mechanical lithosphere of the overriding plate the resulting conducted thermal perturbation leads to melting of the metasomatised overriding mantle lithosphere producing basaltic magmatism that leads to granitic magmatism in the crust dry asthenospheric mantle will melt only if breakoff occurs at a depth shallower than ca 50 km breakoff removes the force at the downdip side of the continental crust while the enhanced heating leads to a reduction of the strength of the underlying crust both effects facilitate the freeing of buoyant crustal sheets which can then rise towards the surface leading to the rapid exhumation of eclogite facies continental crust the cessation of subduction and replacement of the cold oceanic lithosphere by asthenosphere leads to rapid uplift of the orogen we have tested the variety of predicted expressions of slab breakoff in the geological record of the alps the aegean islands and the dabie shan a comparison of the various belts highlights 1 that the magmatism and metamorphism are found near the suture and the centre of the orogen demonstrating the general steepening of suture during collision 2 that a vague inverse correlation exists between maximum depth of metamorphism and volume of syn orogenic magmatism and 3 that the melts can be emplaced in both compressive and extensional environments we suggest that slab breakoff is an important step in the evolution of many orogenic belts and it allows an explanation for the combined presence in collisional orogens of magmatism with a mantle parentage and ultra high pressure metamorphics

miocene unroofing of the canyon range during extension along the sevier desert detachment west central utah this is the published version copyright 2001 american geophysical union all rights reserved

geomorphological analysis of the drainage system on the growing makran accretionary wedge abstract the morphology of six adjacent major catchments draining the onshore makran accretionary wedge in southeast iran and southwest pakistan was studied to examine how the channel pattern and the length profiles may reflect the recent and active growth of the wedge qualitative field surveys were combined with the quantitative analysis of channel steepness and concavity measured from digital elevation models these profiles were compared with modelled profiles using a stream power approach assuming homogeneously uplifting uniform rock substratum results show a distinct difference between the studied western and eastern catchments the three western rivers are in morphological equilibrium whereas the two eastern rivers exhibit profiles with prominent convexities and knickpoints thus notably diverging from equilibrium concave up shapes all the studied catchments share the same base level flow on similar lithologies and developed under uniform climate conditions therefore we interpret the morphometric differences in terms of differential rock uplift rate as a response to local tectonic activity this interpretation is consistent with both uplift rates of marine terraces along the coast of makran and the recorded seismicity the geomorphological analysis presented here extends coastal information to wide inland areas and documents longer term tectonic behaviour than the seismotectonic record hence the steeper surface slope and faster surface uplift rates recorded by the eastern catchments compared to equilibrium of the western catchments are regional long term signals relevant for a wide part of the makran accretionary wedge we attribute the regional geomorphic difference to quaternary variations in tectonic regimes that forced differential uplift rates of the wedge surface we hypothesize that the different tectonic regimes are related to different subduction rates

laboratory experiments on the structure of subducted lithosphere abstract mantle seismic tomography images indicate a variety of shapes for slabs of oceanic lithosphere subducted into the transition zone in the vicinity of the 660 km discontinuity both horizontally deflected slabs and nearly vertical slabs have been inferred at subduction zones on the pacific and indian plates here we determine experimentally the behavior of a cold negatively buoyant and highly viscous slab incident on a fluid interface with a viscosity and a density increase the experiments are scaled to mantle conditions with realistic subduction peclet number ratio of trench migration to subduction velocities and physical properties of the interface a wide range of slab deformation styles is observed when slight changes in the parameters are imposed sinking slab stagnant slab spreading slab sinking pile and stagnant pile represent the five different modes identified we show that these deformation modes are governed by two velocity ratios which characterize the horizontal and the vertical components of the slab velocity near the interface slab penetration into the lower layer takes the form of folded slab piles for large values of the velocity ratios and stagnant behavior is promoted by rapid retrograde trench motion the sensitivity of slab behavior to transition zone density and viscosity changes and to trench migration can account for the coexistence of different modes of deformation in subducted lithosphere

tectonic triggers for postsubduction magmatic hydrothermal gold metallogeny in the late cenozoic anatolian metallogenic trend turkey

tethyan collision forces and the stress field of the eurasian plate resistive forces along convergent plate boundaries have a major impact on surface deformation most visibly at collisional plate boundaries although quantification of these forces is key to understanding the evolution and present state of mountain belts they remain highly uncertain due to the complexity of plate boundary structures and rheologies in previous studies of the eurasian plate we have analysed the balance of plate boundary forces tractions resulting from lithosphere mantle coupling and intraplate variations in topography and density structure this yielded a range of acceptable force distributions in this study we investigate to which extent the observed present day stress field provides further constraints on the distribution of forces we address the dynamics of the eurasian plate as a whole this enables us to base our analysis on mechanical equilibrium of a tectonic plate and to evaluate all forces as part of an internally consistent set of forces driving and deforming eurasia we incorporate tractions from convective mantle flow modelling in a lithospheric model in which edge and lithospheric body forces are modelled explicitly and compute resulting stresses in a homogeneous elastic thin shell intraplate stress observations used are from the world stress map project eurasia s stress field turns out to be particularly sensitive to the distribution of collision forces on the plate s southern margin and to a much lesser extent to lithospheric density structure and tractions from mantle flow stress observations require collision forces on the india eurasia boundary of 7 0 10 5 tn m 1 and on the arabia eurasia boundary of 1 3 2 7 tn m 1 implication of mechanical equilibrium of the plate is that forces on the contacts with the african and australian plates amount to 1 0 2 5 and 0 1 3 tn m 1 respectively we use our results to assess the validity of the classical view that the mean elevation of an orogenic plateau can be taken as a measure of the magnitude of the compressive in this case collision related forces involved for both the tibetan and the iranian plateaus two plateaus with significantly different average elevations we find that the horizontal force derived from the excess gravitational potential energy collapse force is in balance with the collision force

regional scale tectonic evolution and the seismic velocity structure of the lithosphere and upper mantle the mediterranean region we present a numerical modeling approach designed for the verification of tectonic reconstructions and for a quantitative analysis of the seismic velocity structure of the lower lithosphere and upper mantle in this forward modeling procedure we use published reconstructions of the tectonic evolution of a region to predict the thermal structure from which we calculate physical quantities like p wave velocity basin subsidence and basement heat flow we show an application of this approach to two reconstructions of the late mesozoic and cenozoic evolution of the mediterranean region the resulting synthetic velocity models are tested against seismological models of the mantle structure obtained independently by delay time tomography subsidence and heat flow predictions are compared with actual observations to further constrain the geological assumptions underlying the forward models the mantle models predicted from reconstructions of the evolution of the mediterranean region correlate with major features in the tomographic images from our model results we conclude that the present distribution of seismic velocities in the mediterranean region is largely determined by the tectonic evolution during at least the last 40 m y this means we can verify processes represented in different tectonic scenarios with the seismologically determined velocity structure seen from a different perspective the model results contribute to understanding the causes of the velocity structure imaged by delay time tomography by showing how the tectonic evolution results in this structure

the fate of slabs inferred from seismic tomography and 130 million years of subduction the volume and location of subducted plate since 130 ma are reconstructed from magnetic anomalies and updated finite rotation parameters describing relative motions in the ocean basins and between plates and hotspots the area of subducted plate is calculated from the relative motions between overriding and overridden plates along the length of convergence in the hotspot reference frame in 5 ma intervals ridge locations are reconstructed by rotating magnetic anomalies from their present to their former positions the distances between trenches and ridges at a certain point of time in the spreading direction are thus estimated and are used to determine the ages of trenches according to the spreading history recorded in present day oceans the thickness of the oceanic plate is obtained from the age based on the half space cooling model about 3 45 10 10 km 3 of oceanic plate has been subducted during the past 130 ma with maximum accumulation beneath southeast eurasia at spherical harmonic degree l 2 where subduction flux peaks excellent correlation is found between four seismic tomographic models at the top of the lower mantle about 800 1100 km deep and predicted slab locations for some seismic models excellent correlation is also found at other degrees in that depth range some models also have good correlations with predicted slab locations in the deep mantle at l 2 direct comparison between subduction and seismic tomographic patterns at 800 1100 km shows that subduction history correlates with tomography very well in terms of location and amplitude of anomalies cold downwellings which may be related to slabs appear to be trapped in the mesosphere or middle mantle correlations between seismic tomography and subduction over different time periods support this conclusion there are also some correlations in the upper mantle and in the deep lower mantle although they are generally not as significant as those in the 800 1100 km depth range there may be a significant boundary in the mantle near the 800 1100 km depth from a geodynamic and chemical point of view the lower mantle may start at a depth closer to 1000 km than to 670 km

the fate of descending slabs oceanic lithosphere the stiff plate formed by the cooling of chemically diffcrentiated material that rises and mclts under mid ocean ridges and spreads laterally in the process called sea floor spreading has an ephemeral residence at the surface of the earth typically within 100 200 ma densi fication of the old oceanic lithosphere leads to a dynamic instability causing it to sink into the interior descending as a layered slab with highly anom alous thermal and chemical characteristics relative to the surrounding ambient mantle the slab heats as it sinks but the warming is sluggish slow to overcome the effect of more than 100 ma of cooling at the surface and the negative buoyancy of the slab continues to drive its descent increasing pressure and temperature with depth drive various chemical reactions in the slab including expulsion of water and other volatiles accumulated at the surface and destabilization of associated hydrous phases suites of mineralogical phase transformations occur within the slab s crustal and depleted mantle components because the slab has been processed through the melting zone at the ridge partially hydrated and cooled significantly the phase equilibria within the slab differ from those of the surrounding mantle near 400 km depth an exothermic phase transformation occurs with mg fe 2si04 olivine in the slab converting to the high pressure modified j spinel structure the spinel transition helps to drive the slab onward although a central core of metastable olivine may persist to greater depths with eventual transformational faulting

petrology and geochemistry of post collisional middle eocene volcanic units in north central turkey evidence for magma generation by slab breakoff following the closure of the northern neotethys ocean the eocene volcano sedimentary units of northern anatolia are confined into a narrow zone trending parallel to the intra pontide and izmir ankara erzincan sutures along which the northern branch of the neotethys ocean was closed during a period between late maastrichtian and paleocene the middle eocene formations overlie both the imbricated and highly deformed units of the suture zone which are paleocene or older in age as well as the formations of adjacent continental blocks with a regional disconformity therefore they can be regarded to be post collisional these units are composed of subaerial to shallow marine sedimentary beds i e the orencik formation at the base and a subaerial volcanic unit i e the hamamozu formation in the middle and at the top this sudden facies change from marine to subaerial environment in the middle eocene is a common phenomenon across northern turkey implying that a regional uplift event occurred possibly across the suture zone before the initiation of the volcanism during lutetian the middle eocene lavas span the whole compositional range from basalts to rhyolites and display a calc alkaline character except for alkaline to mildly alkaline lavas from the top of the sequence all lavas display a distinct subduction signature our geochemical data indicate that calc alkaline lavas were derived from a subduction modified source whereas alkaline to mildly alkaline lavas of the late stage were possibly sourced by an enriched mantle domain magmas evolved in magma chambers emplaced possibly at two different crustal levels magmas in deeper 13 km and possibly larger chambers fractionated hydrous mafic minerals e g amphibole and biotite two pyroxenes and plagioclase and assimilated a significant amount of crustal material intermediate to acid calc alkaline lavas and pyroclastics were derived from these chambers magmas in the shallower chambers on the other hand 12 km crystallized anhydrous mineral assemblages assimilated little or no crustal material and fed basic to intermediate lavas in the region both deep and shallow chambers were periodically replenished by mafic magmas we argue that a slab breakoff model explains better than any alternative model i why the volcanism during the middle eocene was confined into a rather narrow belt along the suture zone ii why it initiated almost contemporaneous with a regional uplift after the continental collision event iii why it postdated arc volcanism along the pontides in the north by 15 20 my iv why it assimilated significant amount of crustal material and v why alkalinity of lavas increased in time

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

building the zagros collisional orogen timing strain distribution and the dynamics of arabia eurasia plate convergence he zagros mountains are the result of the arabia eurasia collision initiated at 35 ma as the rifted arabian lithosphere was underthrusted beneath the iranian plate due to its negative buoyancy the onset of crustal thickening started at 25 ma as recorded by the hinterland exhumation and foreland clastic deposition deformation throughout the arabia eurasia collision zone and the uplift of the iranian plateau occurred after 15 12 ma as a result of shortening thickening of the thin iranian crust we emphasize that only 42 of the post 35 ma convergence is partitioned by shortening within central iran tomographic constraints show ongoing slab steepening or breakoff in the nw zagros whereas underthrusting of the arabian plate is observed beneath central zagros the current subduction dynamics can be explained by the original lateral difference in the buoyancy of the distal margin that promoted slab sinking in nw zagros and underthrusting in central zagros critical wedge approach applied to the zagros favors the hypothesis of strong brittle crust detached above a viscous lower crust in contrast the weak sedimentary cover deforms by buckling of a thick multilayered cover thrust faulting associated with folding occurs in the competent layers and is responsible for most of the earthquakes there is evidence that the role of the slab pull force in driving the arabian plate motion was reduced after 12 ma large scale mantle flow induced by mantle upwelling at the afar plume appears to be the main driver of the arabia plate motion we stress that the main kinematic change in the zagros region occurred at 15 12 ma as the zagros uplifted before the arabian slab detached the zagros appears key to investigate coupling between continental rheology plate driving forces and mountain building in which the role of rift inheritance appears to be central

large scale subduction of continental crust implied by india asia mass balance calculation continental crust is buoyant compared with its oceanic counterpart and resists subduction into the mantle when two continents collide the mass balance for the continental crust is therefore assumed to be maintained here we use estimates of pre collisional crustal thickness and convergence history derived from plate kinematic models to calculate the crustal mass balance in the india asia collisional system using the current best estimates for the timing of the diachronous onset of collision between india and eurasia we find that about 50 of the pre collisional continental crustal mass cannot be accounted for in the crustal reservoir preserved at earth s surface today represented by the mass preserved in the thickened crust that makes up the himalaya tibet and much of adjacent asia as well as southeast asian tectonic escape and exported eroded sediments this implies large scale subduction of continental crust during the collision with a mass equivalent to about 15 of the total oceanic crustal subduction flux since 56 million years ago we suggest that similar contamination of the mantle by direct input of radiogenic continental crustal materials during past continent continent collisions is reflected in some ocean crust and ocean island basalt geochemistry the subduction of continental crust may therefore contribute significantly to the evolution of mantle geochemistry buoyant continental crust is thought to resist subduction calculation of the crustal mass balance during the collision between india and eurasia indicates that about 50 of pre collisional continental crust has been subducted into the mantle

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

imaging composite dike propagation etna 2002 case 1 late on the night of 26 october 2002 a dike intrusion started suddenly at mount etna producing intense explosive activity and lava effusion on the southern flank five to six hours afterward a long field of eruptive fractures propagated radially along the northeastern flank of the volcano producing marked variations at the permanent tilt network the dike propagation velocity was inferred by the associated seismicity we modeled the temporal evolution of the continuously recorded tilt data both during the vertical dike propagation on the high south flank on 26 october and during the radial propagation along the northeast flank between 27 and 28 october the reproduction of the recorded tilt signal allowed us to describe the geometry and characteristics of the two dikes in greater detail than the previous static inversion we deduced that the eruption was characterized by an unusual composite mechanism clearly showing a transition from a nearly pure opening mode displacement to a mechanism characterized by an equally strong normal dip slip component and a smaller left lateral strike slip component in this study we demonstrate the interaction between the final segment of the dike and a preexisting structure that was reactivated in response to the intrusion we show that tilt and its modeling represent a powerful tool to verify and constrain dike intrusions in detail

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

restoration of cenozoic deformation in asia and the size of greater india a long standing problem in the geological evolution of the india asia collision zone is how and where convergence between india and asia was accommodated since collision proposed collision ages vary from 65 to 35 ma although most data sets are consistent with collision being underway by 50 ma plate reconstructions show that since 50 ma 2400 3200 km west to east of india asia convergence occurred much more than the 450 900 km of documented himalayan shortening current models therefore suggest that most post 50 ma convergence was accommodated north of the indus yarlung suture zone we review kinematic data and construct an updated restoration of cenozoic asian deformation to test this assumption we show that geologic studies have documented 600 750 km of n s cenozoic shortening across and north of the tibetan plateau the pamir hindu kush region accommodated 1050 km of n s convergence geological evidence from tibet is inconsistent with models that propose 750 1250 km of eastward extrusion of indochina approximately 250 km of indochinese extrusion from 30 to 20 ma of indochina suggested by se asia reconstructions can be reconciled by dextral transpression in eastern tibet we use our reconstruction to calculate the required size of greater india as a function of collision age even with a 35 ma collision age the size of greater india is 2 3 times larger than himalayan shortening for a 50 ma collision the size of greater india from west to east is 1350 2600 km consistent with robust paleomagnetic data from upper cretaceous paleocene tethyan himalayan strata these estimates for the size of greater india far exceed documented shortening in the himalaya we conclude that most of greater india was consumed by subduction or underthrusting without leaving a geological record that has been recognized at the surface

petrologic case for eocene slab breakoff during the indo asian collision the greater himalayan sequence is the metamorphic core of the himalaya and has been a focus of considerable study yet its petrologic evolution remains controversial pre oligocene metamorphism was nearly obliterated by miocene metamorphism and melting which many workers ascribe to shear heating unusually high concentrations of radioactive elements accompanying burial and or several kilometers of miocene exhumation due to extensional faulting sparse oligocene to eocene ages are often assigned to a tectonically unspecified event an alternative slab breakoff model subduction of greater himalayan rocks to 100 km depth in the eocene followed by buoyant extrusion due to decoupling of the oceanic lithosphere has also been proposed based on rare eocene eclogites but without explanation of miocene melting or metamorphic petrogenesis we argue that slab breakoff readily explains the eocene eclogites miocene partial melts and late eocene k rich magmas in southeastern tibet and that metamorphic and plutonic ages help define the timing and rates of breakoff and extrusion this model implies that 1 much of the greater himalayan sequence was subducted to depths greater than commonly considered 2 fluid absent decompression melting at 30 35 km depth was the consequence of as much as 100 km of extrusion rather than radioactive heating or a smaller crustal level extensional or erosional event and 3 eruption of eocene k rich tibetan plateau lavas has no implication for topography of the tibetan plateau

how valid are dynamic models of subduction and convection when plate motions are prescribed a detailed comparison between fully dynamic and kinematic plate formulations has been made in models of mantle convection plate velocity is computed self consistently from fully dynamic plate models with temperature and stress dependent viscosity and preexisting mobile faults in fully dynamic models the flow is driven solely by internal buoyancy while in kinematic models the flow is driven by a combination of the prescribed surface velocity and internal buoyancy only a temperature dependent viscosity close to the effective viscosity determined from the fully dynamic models is used in the kinematic models the two types of models give very similar temperature structures and slab evolutionary histories when the effective viscosity and surface velocity are nearly identical in kinematic plate models the additional work introduced by the prescribed velocity boundary condition is apparently dissipated within the lithosphere and has little influence on the convection under the lithosphere in models with periodic lateral boundary conditions slabs sink into the lower mantle at an oblique angle and this contrasts with the vertical sinking which occurs with reflecting boundary conditions models show that we can simulate fully dynamic models with kinematic models under either periodic boundary conditions or reflecting boundary conditions

constraining central neo tethys ocean reconstructions with mantle convection models a striking feature of the indian ocean is a distinct geoid low south of india pointing to a regionally anomalous mantle density structure equally prominent are rapid plate convergence rate variations between india and se asia particularly in late cretaceous paleocene times both observations are linked to the central neo tethys ocean subduction history for which competing scenarios have been proposed here we evaluate three alternative reconstructions by assimilating their associated time dependent velocity fields in global high resolution geodynamic earth models allowing us to predict the resulting seismic mantle heterogeneity and geoid signal our analysis reveals that a geoid low similar to the one observed develops naturally when a long lived back arc basin south of eurasia s paleo margin is assumed a quantitative comparison to seismic tomography further supports this model in contrast reconstructions assuming a single northward dipping subduction zone along eurasia s margin or models incorporating a temporary southward dipping intra oceanic subduction zone cannot sufficiently reproduce geoid and seismic observations

subducted slabs beneath the eastern indonesia tonga region insights from tomography tomographic images of mantle structure beneath the region north and northeast of australia show a number of anomalously fast regions these are interpreted using a recent plate tectonic reconstruction in terms of current and former subduction systems several strong anomalies are related to current subduction the inferred slab lengths and positions are consistent with neogene subduction beneath the new britain and halmahera arcs and at the tonga and the new hebrides trenches where there has been rapid rollback of subduction hinges since about 10 ma there are several deeper flat lying anomalies which are not related to present subduction and we interpret them as former subduction zones overridden by australia since 25ma beneath the bird s head and arafura sea is an anomaly interpreted to be due to north dipping subduction beneath the philippines halmahera arc between 45and 25ma a very large anomaly extending from the papuan peninsula to the new hebrides and from the solomon islands to the east australian margin is interpreted to be the remnant of south dipping subduction beneath the melanesian arc between 45and 25ma this interpretation implies that a flat lying slab can survive for many tens of millions of years at the bottom of the upper mantle in the lower mantle there is a huge anomaly beneath the gulf of carpentaria and east papua new guinea this is located above the position where the tectonic model interprets a change in polarity of subduction from north dipping to south dipping between 45and 25ma we suggest this deep anomaly may be a slab subducted beneath eastern australian during the cretaceous or subducted north of australia during the cenozoic before 45ma the tomography also supports the tectonic interpretation which suggests little neogene subduction beneath western new guinea since no slab is imaged south of the new guinea trench however one subduction zone in the tectonic model and many others that associated with the trobriand trough east of papua new guinea and the miocene maramuni arc is not seen in the tomographic images and may require reconsideration of currently accepted tectonic interpretations 4 2002 elsevier science b v all rights reserved

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

global positioning system constraints on plate kinematics and dynamics in the eastern mediterranean and caucasus we present and interpret global positioning system gps measurements of crustal motions for the period 1988 1997 at 189 sites extending east west from the caucasus mountains to the adriatic sea and north south from the southern edge of the eurasian plate to the northern edge of the african plate sites on the northern arabian platform move 18 2 mm yr at n25 5 w relative to eurasia less than the nuvel 1a circuit closure rate 25 1 mm yr at n21 7 w preliminary motion estimates 1994 1997 for stations located in egypt on the northeastern part of africa show northward motion at 5 6 2 mm yr also slower than nuvel ia estimates 10 1 mm yr at n2 4 e eastern turkey is characterized by distributed deformation while central turkey is characterized by coherent plate motion internal deformation of 2 mm yr involving westward displacement and counterclockwise rotation of the anatolian plate the anatolian plate is de coupled from eurasia along the right lateral strike slip north anatolian fault naf we derive a best fitting euler vector for anatolia eurasia motion of 30 7 0 8 n 32 6 0 4 e 1 2 0 1 myr the euler vector gives an upper bound for naf slip rate of 24 1 mm yr we determine a preliminary gps arabia anatolia euler vector of 32 9 1 2 n 40 3 1 1 e 0 8 0 2 myr and an upper bound on left lateral slip on the east anatolian fault eaf of 9 1 mm yr the central and southern aegean is characterized by coherent motion internal deformation of 2 mm yr toward the sw at 30 1 mm yr relative to eurasia stations in the se aegean deviate significantly from the overall motion of the southern aegean showing increasing velocities toward the trench and reaching 10 1 mm yr relative to the southern aegean as a whole

geology of the onshore makran accretionary wedge synthesis and tectonic interpretation abstract this contribution is a synthesis of studies carried out since 2004 in the makran accretionary wedge exposed in se iran these studies were intended to improve our understanding of the behaviour of accretionary wedges and of the significance of the makran wedge in the frame of tethys tectonics reassessing the stratigraphy and structure of the iranian makran has allowed distinguishing four main litho tectonic units the lithological content of these four units is briefly described to discuss the sedimentary environments which record the general filling of an oceanic basin with clastic sediments derived from the north the petrology of ophiolites exposed in the backstop region to the north of the accretionary wedge itself demonstrates the existence of an oceanic basin formed in the jurassic late cretaceous to the north of continental slivers which had separated from central iran there is petrological evidence for a subduction related magmatic arc during the late cretaceous which sourced much of the clastic material in the cretaceous lower miocene turbidites present throughout the wedge the literature survey shows that the late cretaceous to present evolution of the pakistan part of the makran fold and trust belt is in its major attributes similar to that of the iranian makran these observations and results are discussed in a tectonic model that integrates knowledge on the oman mountains on the other side of the current oman subduction

thermal and dynamic modeling of deep subduction of a spreading center implications for the fate of the subducted chile rise southern chile our objective is to determine whether the chile rise spreading center currently subducting south of the chile triple junction in southern chile continues or ceases to separate after subduction a two dimensional kinematic thermal model is used to predict the temperature history of the subduction zone during spreading center subduction density anomalies calculated from the temperature fields of the thermal model are then used to define body forces that drive the motion of a dynamic viscous fluid model of spreading center subduction viscous models that include only local in plane driving forces predict that the chile rise spreading center ceases to separate after ridge subduction if global in plane plate forces are also included in the viscous model continuing separation of the subducted chile rise is predicted but at a greatly diminished rate fully continuing separation of the subducted chile rise only occurs if the subducted nazca and antarctic plates remain rigid allowing the lateral transmission of global plate forces from along strike however this situation is unlikely as young subducted plates are believed to fragment continuing separation but at a greatly diminished rate is believed to be the most plausible hypothesis for spreading center subduction in southern chile dynamic viscous modeling of subduction in the central andes suggests a viscosity of 2 5 x 10 20 pa s for upper mantle above the 400 km phase transition