lithospheric structure and along strike segmentation of the central andean plateau seismic q magmatism flexure topography and tectonics the central andean plateau is a 300 km wide nearly 4 km high plateau which is situated above a 30 e dipping segment of the subducted nazca plate systematic changes in the topography the upper mantle seismic attenuation q structure the lithospheric flexural rigidity and the distribution and chemistry of back arc lavas indicate that the upper mantle structure changes along strike beneath the plateau south of 23 s the upper mantle beneath the plateau becomes hotter and the lithosphere becomes thinner and weaker this lateral variation in upper mantle structure coincides with variations in the tectonic style and timing of deformation in two physiographically distinct segments of the plateau and its adjacent foreland thrust belt to the east the bolivian altiplano and subandean ranges in the north and the argentine puna and santa barbara system in the south the segments different mantle lid thicknesses result in an altiplano supported mainly by crustal thickening and a puna supported by a combination of a crustal root and a thermal mantle root a thicker stronger foreland lithosphere in the north favors thin skinned underthrusting of the foreland beneath the plateau while to the south the thinner weaker lithosphere allows more basement involved shortening within the foreland and less underthrusting beneath the plateau a weaker lithosphere results in a longer period of compressional deformation and greater local and structural relief in the puna than in the relatively stable altiplano to the north ultimately this observed n s directed segmentation on the plateau and in the foreland may reflect the different responses of strong and weak continental lithosphere to an andean type subduction margin

global mantle heterogeneity and its influence on teleseismic regional tomography abstract a new global p wave tomography model is determined using a flexible grid parameterization this new model better reveals the mantle structure under the polar regions than the previous tomographic models the subducting slabs are generally imaged clearly as high velocity high v zones the young slabs are still subducting in the upper mantle and the mantle transition zone mtz whereas the old and ancient slabs are either stagnant in the mtz or have subducted down to the lower mantle even reaching the core mantle boundary low velocity low v anomalies are generally revealed in the mantle under the hotspot regions it seems that a variety of mantle upwelling plumes exist some strong plumes are visible in the whole mantle under the long living hotspots such as those in south central pacific africa hawaii and iceland whereas weak plumes are visible in only some depth range under the minor hotspots under the intraplate volcanoes in east eurasia bering sea and west alaska significant low v anomalies are revealed in the upper mantle which may reflect hot and wet upwelling associated with corner flows in the big mantle wedge bmw above the stagnant pacific slab in the mtz and perhaps deep slab dehydration as well the subduction triggered magmatism in the bmw may be a new class of mantle plumes we also used the new global model to investigate the influence of whole mantle heterogeneity on the determination of upper mantle tomography under japan with a teleseismic inversion method the results show that the mantle heterogeneities outside the target volume of regional tomography can cause significant changes 0 2 0 4 s to the observed relative travel time residuals of a distant earthquake the pattern of regional tomography remains the same even after correcting for the whole mantle heterogeneity but there are some changes in the amplitude of velocity anomalies in the regional tomography hence it is necessary to correct for the mantle heterogeneities outside the target volume so as to obtain a better regional tomography

evidence for lithospheric detachment in the central andes from local earthquake tomography abstract data from three temporary seismic networks were merged for tomographic inversion although the deployments did not coincide in time spatial overlap was achieved by re occupying existing sites travel times and t operators of about 1600 earthquakes were inverted for 3d models of p p s and p wave attenuation q p 1 all three attributes provide a consistent image of the entire subduction zone on a lithospheric scale the tomographic images reveal low velocities and high attenuation in the crust and mantle underlying the western cordillera and most of the puna plateau indicative of weak rheology and mostly asthenospheric mantle in contrast forearc and eastern foreland are characterized by high q p values corresponding to cold temperatures in accordance with thermal models in the backarc between 23 s and 24 s a high velocity high q p structure beneath the eastern cordillera and eastern puna is interpreted as detaching continental lithosphere that has been thickened in the orogenic process south of this structure the mantle is characterized by low velocities high p s ratios and low q p values here it is believed that lithosphere originally underlying andean crust has already been removed this is supported by new estimates of crustal thickness and volcanic activity

the lost inca plateau cause of flat subduction beneath peru abstract since flat subduction of the nazca plate beneath peru was first recognized in the 1970s and 1980s a satisfactory explanation has eluded researchers we present evidence that a lost oceanic plateau inca plateau has subducted beneath northern peru and propose that the combined buoyancy of inca plateau and nazca ridge in southern peru supports a 1500 km long segment of the downgoing slab and shuts off arc volcanism this conclusion is based on an analysis of the seismicity of the subducting nazca plate the structure and geochemistry of the marquesas plateau as well as tectonic reconstructions of the pacific farallon spreading center 34 to 43 ma these restore three sub parallel pacific oceanic plateaus the austral tuamotu and marquesas to two farallon plate counterparts the iquique and nazca ridges inca plateau is apparently the sixth and missing piece in an ensemble of v shaped hotspot tracks formed at on axis positions we argue the mirror image of the inca plateau the marquesas plateau is an ancient edifice overprinted by recent volcanism in disagreement with the widely accepted young

receiver function images from the moho and the slab beneath the altiplano and puna plateaus in the central andes summary teleseismic data recorded during one and a half years are investigated with the receiver function technique to determine the crustal and upper mantle structures underneath the highly elevated altiplano and puna plateaus in the central andes a series of converting interfaces are determined along two profiles at 21 s and 25 5 s respectively with a station spacing of approximately 10 km the data provide the highest resolution gained from a passive project in this area so far the oceanic nazca plate is detected down to 120 km beneath the altiplano whereas beneath the puna the slab can unexpectedly be traced down to 200 km depth at longer periods a shallow crustal low velocity zone is determined beneath both plateaus exhibiting segmentation in the case of the altiplano the segments present vertical offsets and are separated by inclined interfaces which coincide with major fault systems at the surface an average depth to moho of about 70 km is determined for the altiplano plateau a strong negative velocity anomaly located directly below the moho along with local crustal thinning is interpreted beneath the volcanic arc of the altiplano plateau between 67 w and 68 5 w a deep section of the puna profile reveals thinning of the mantle transition zone although poorly resolved the detected anomaly may suggest the presence of a mantle plume which may constitute the origin of the anomalous temperatures at the depth of the upper mantle discontinuities

development of the accretionary prism along peru and material flux after subduction of nazca ridge results from ocean drilling program leg 112 indicated tectonic erosion of the peruvian convergent margin during subduction of the nazca ridge followed by renewed accretion against the erosional scar seismic images and swath mapping morphology show a regional geology shaped as subduction of nazca ridge migrated 800 km along the peru trench the accretionary prism is imaged at two stages of development since upper miocene time with depth migrated seismic records in which the back stop is well defined from these data material flux can be quantified after the ridge crest subducted and accretion dominated over erosion the prism grew rapidly to 10 to 15 km width during rapid growth the margin taper was large and about 60 percent of the sediment supply accreted as the prism growth slowed only 30 percent accreted despite an increased trench sediment supply since the convergence rate changed little the inverse relation between sediment supply and prism growth rate suggests that other processes allow most of the sediment input to be subducted the decreased prism growth may be partially controlled by varying the structure of the back stop the force required to activate faults cutting across the back stop is probably greater than the force allowed by basal friction across the plate boundary thus the shear stress transmitted across the plate boundary is insufficient to drive out of sequence thrust faults that would thicken the back stop or raise it beyond a limited height a constant back stop height limits accretion and the subduction window appears to open and accommodate increased trench sediment beneath the back stop rather than inducing a transport of sediment over the top of the back stop nazca ridge subduction affected the convergent margin mass flux for about 8 ma mass balancing indicates an order of magnitude increase in the current sediment volume subducted in the north compared to that near the ridge crest when growth of the accretionary prism slowed a greater terrigenous input through the trench was accommodated by sediment subduction the resulting range of material flux is estimated to increase the sediment interlayer between the plates about 500 m which may affect coupling

buckling instabilities of subducted lithosphere beneath the transition zone a sheet of viscous fluid poured onto a surface buckles periodically to generate a pile of regular folds recent tomographic images beneath subduction zones together with quantitative fluid mechanical scaling laws suggest that a similar instability can occur when slabs of subducted oceanic lithosphere encounter enhanced resistance to penetration at 660 km depth beneath subduction zones such as central america and java the fast anomalies in elastic wavespeed are too broad 400 km to be explained by simple compressional thickening of subducted lithosphere that is initially 50 100 km thick by contrast a periodic buckling mechanism can account for both the characteristic wedge shapes of the anomalies and for their widths which agree closely with those predicted by the scaling laws in addition to providing a plausible physical interpretation of tomographic observations beneath some subduction zones buckling instabilities have important consequences for other aspects of global geodynamics such as tectonic plate reconstructions interpretation of geoid anomalies and modeling of secular changes in earth s moment of inertia

moho topography in the central andes and its geodynamic implications p to s converted waves at the continental moho together with waves multiply reflected between the earth s surface and the moho have been used to estimate the moho depth and average crustal vp vs variations in the central andes our analysis confirms and significantly complements the moho depth estimates previously obtained from wide angle seismic studies and receiver functions the resulting crustal thickness varies from about 35 km in the forearc region to more than 70 km beneath the plateau and thins 30 km further to the east in the chaco plains beneath the andean plateau the moho is deeper in the north altiplano and shallower in the south puna where the plateau attains its maximum elevation a non linear relation exists between crustal thickness and elevation and bouguer gravity suggesting that the crust shallower than 50 55 km is predominately felsic in contrast to a predominately mafic crust below such a relation also implies a 100 km thick thermal lithosphere beneath the altiplano and with a lithospheric thinning of a few tens of kilometers beneath the puna absence of expected increase in lithospheric thickness in regions of almost doubled crust strongly suggests partial removal of the mantle lithosphere beneath the entire plateau in the subandean ranges at 19 20 s the relation between altitude and crustal thickness indicates a thick lithosphere up to 130 150 km and lithospheric flexure beneath a relative topographic low at the salar de atacama a thick crust 67 km suggests that the lithosphere in this region is abnormally cold and dynamically subsided possibly due to coupling with the subducting plate this may be related to the strongest ms 8 0 known intra slab earthquake in the central andes that happened very close to this region in 1950 the average crustal vp vs ratio is about 1 77 for the altiplano puna and it reaches the highest values 1 80 1 85 beneath the volcanic arc indicating high ambient crustal temperatures and wide spread intra crustal melting 2002 elsevier science b v all rights reserved

structure of the subduction system in southern peru from seismic array data the subduction zone in southern peru is imaged using converted phases from teleseismic p pp and pkp waves and pwave tomography using local and teleseismic events with a linear array of 50 broadband seismic stations spanning 300 km from the coast to near lake titicaca the slab dips at 30 and can be observed to a depth of over 200 km the moho is seen as a continuous interface along the profile and the crustal thickness in the back arc region the altiplano is 75 km thick which is sufficient to isostatically support the andes as evidenced by the gravity the shallow crust has zones of negative impedance at a depth of 20 km which is likely the result of volcanism at the midcrustal level of 40 km there is a continuous structure with a positive impedance contrast which we interpret as the western extent of the brazilian craton as it underthrusts to the west v p v s ratios estimated from receiver function stacks show average values for this region with a few areas of elevated v p v s near the volcanic arc and at a few points in the altiplano the results support a model of crustal thickening in which the margin crust is underthrust by the brazilian shield

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

the role of ridges in the formation and longevity of flat slabs flat slab subduction occurs when the descending plate becomes horizontal at some depth before resuming its descent into the mantle it is often proposed as a mechanism for the uplifting of deep crustal rocks thick skinned deformation far from plate boundaries and for causing unusual patterns of volcanism as far back as the proterozoic eon for example the formation of the expansive rocky mountains and the subsequent voluminous volcanism across much of the western usa has been attributed to a broad region of flat slab subduction beneath north america that occurred during the laramide orogeny 80 55 million years ago here we study the largest modern flat slab located in peru to better understand the processes controlling the formation and extent of flat slabs we present new data that indicate that the subducting nazca ridge is necessary for the development and continued support of the horizontal plate at a depth of about 90 kilometres by combining constraints from rayleigh wave phase velocities with improved earthquake locations we find that the flat slab is shallowest along the ridge while to the northwest of the ridge the slab is sagging tearing and re initiating normal subduction on the basis of our observations we propose a conceptual model for the temporal evolution of the peruvian flat slab in which the flat slab forms because of the combined effects of trench retreat along the peruvian plate boundary suction and ridge subduction we find that while the ridge is necessary but not sufficient for the formation of the flat slab its removal is sufficient for the flat slab to fail this provides new constraints on our understanding of the processes controlling the beginning and end of the laramide orogeny and other putative episodes of flat slab subduction

anisotropy and mantle flow in the chile argentina subduction zone from shear wave splitting analysis 1 we examine shear wave splitting in teleseismic phases to observe seismic anisotropy in the south american subduction zone data is from the charge network which traversed chile and western argentina across two transects between 30 s and 36 s beneath the southern and northwestern parts of the network fast polarization direction j is consistently trench parallel while in the northeast j is trench normal the transition between these two zones is gradual we infer that anisotropy sampled by teleseismic phases is localized within or below the subducting slab we explain our observations with a model in which eastward nazca entrained asthenospheric flow is deflected by retrograde motion of the subducting nazca plate resulting southward flow through this area produces n s j observed in the south and northwest e w j result from interaction of this flow with the local slab geometry producing eastward mantle flow under the actively flattening part of the slab index terms 7203 seismology body wave propagation 7218 seismology lithosphere and upper mantle 8123 tectonophysics dynamics seismotectonics 8150 tectonophysics plate boundary general 3040 9360 information related to geographic region south america citation anderson m l g zandt e triep m fouch and s beck 2004 anisotropy and mantle flow in the chile argentina subduction zone from shear wave splitting analysis geophys res lett 31 l23608 doi 10 1029 2004gl020906

p and s velocity images of the lithosphere asthenosphere system in the central andes from local source tomographic inversion summary about 50 000 p and s arrival times and 25 000 values of t recorded at seismic arrays operated in the central andes between 20 s and 25 s in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3 d p s velocity and qp models the study volume in the reference model is subdivided into three domains slab continental crust and mantle wedge a starting velocity distribution in each domain is set from a priori information in the crust it is based on the controlled sources seismic studies in slab and mantle wedge it is defined using relations between p and s velocities temperature and composition given by mineral physics each iteration of tomographic inversion consists of the following steps 1 absolute location of sources in 3 d velocity model using p and s arrival times 2 double difference relocation of the sources and 3 simultaneous determination of p and s velocity anomalies p and s station corrections and source parameters by inverting one matrix velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double sided nodes at the domain boundaries the next iteration is repeated with the updated velocity model and source parameters obtained at the previous step different tests aimed at checking the reliability of the obtained velocity models are presented in addition we present the results of inversion for vp and vp vs parameters which appear to be practically equivalent to vp and vs inversion a separate inversion for qp has been performed using the ray paths and source locations in the final velocity model the resulting vp vs and qp distributions show complicated essentially 3 d structure in the lithosphere and asthenosphere p and s velocities appear to be well correlated suggesting the important role of variations of composition temperature water content and degree of partial melting

lithospheric scale structure across the bolivian andes from tomographic images of velocity and attenuation for p and s waves we have developed a three dimensional lithospheric scale model across the bolivian andes at 20 s based on tomographic images of velocity and attenuation for both p and s waves observations of travel time and attenuation for this study are from regional mantle earthquakes in the subducted nazca plate recorded on a portable broadband seismic array broadband andean joint experiment and seismic exploration of the deep andes in bolivia and chile the shallow mantle under the altiplano from 18 s to 21 s is high velocity and moderately high q vp 8 3 vs 4 7 qp 500 and qs 200 suggesting lithospheric mantle high velocity material in the altiplano extends to a depth of 125 150 km the shallow mantle of the western cordillera is characterized by high vp vs 1 83 suggesting a correlation between vp vs and arc volcanism seismic velocity in the western cordillera mantle is on average only slightly reduced from global averages however velocity and attenuation anomalies are locally strong vp 7 8 vs 4 3 qp 200 and qs 100 consistent with partial melt conditions under the los frailes volcanic field in the eastern cordillera shallow mantle velocity and q decrease drastically from the neighboring altiplano vp 7 8 vs 4 3 qp 300 qs 100 however high vp vs is not as pervasive as it is in the western cordillera we believe that slab derived water and perhaps other volatiles strongly influence the western cordillera while the eastern cordillera low velocity region is more affected by partial melt and or compositional changes average velocity and q in the shallow mantle across the bolivian andes where the tomographic images are best resolved are significantly higher than in most mantle wedge environments where corresponding images are available this is likely the result of a compressional back arc setting in the andes this implies that lithospheric shortening and thickening associated with the formation of the andes has profoundly influenced the shallow mantle structure across the range shallow mantle structure is locally influenced by the subduction processes particularly under the western cordillera however the differing volcanism and seismic character under the two cordilleras suggest that the volcanic process in the eastern cordillera may be distinct from arc volcanism tertiary volcanism in the eastern cordillera is located in the region where mantle shortening is suspected to be greatest both the timing and location of volcanism are consistent with upward migration of mantle wedge asthenosphere following the removal of over thickened lithosphere

imaging of subducted lithosphere beneath south america tomographic images are produced for the deep structure of the andean subduction zone beneath western south america the data used in the imaging are the delay times of p pp and pwp phases from relocated teleseismic earthquakes in the region regionally structural features larger than about 150 km are resolved by the data presentations of layer anomaly maps and cross sections reveal 1 the nazca slab is probably continuous laterally and at depth over most regions 2 the offset between the north and south deep earthquake zones containing the 1994 deep bolivia main shock and its aftershocks can be modelled by a northwest striking and steeply northeast dipping slab structure and 3 the nazca slab clearly penetrates the lower mantle beneath central south america but is partly deflected in the southern deep zone

elastic thickness structure of south america estimated using wavelets and satellite derived gravity data abstract we used a wavelet formulation of the classical spectral isostatic analysis to invert satellite derived gravity and topography bathymetry for elastic thickness te over south america and its surrounding plates to provide a homogeneous representation of the gravity field for this vast region we corrected free air anomalies derived from a combination of terrestrial marine gravity data with data from the grace and champ satellite missions model eigen cg03c by a simple bouguer slab using a smoothed representation of surface relief wavelengths 125 km the resulting bouguer anomaly compares well with terrestrial data acquired in the central andes and allows te to be confidently estimated for values greater than 10 km the te map resolves regional scale features that are well correlated with known surface structures and shows maximum values of 100 15 km over the archean neoproterozoic core of the continent decreasing to less than 30 km around continental margins several regions of the oceanic plates and continental margins have an elastic thickness less than 10 km we performed a quantitative analysis by comparing the elastic thickness with the thermal structure predicted from the age of oceanic crust and igneous metamorphic rocks this demonstrates that oceanic plates have been weakened by thermal interaction with hotspots and locally by fracturing and hydration near the trench we observe that only the nucleus of the continent has resisted the thermomechanical weakening induced by the rifting of africa and south america along the passive margin and the andean orogeny along the active margin this latter region shows along strike variations in te that correlate with the geotectonic segmentation of the margin and with the pattern of crustal seismicity our results reveal that the rigidity structure follows the segmentation of the seismogenic zone along the subduction fault suggesting a causal relationship that should be investigated in order to improve the understanding and predictability of great earthquakes and tsunamis

the migration history of the nazca ridge along the peruvian active margin a re evaluation abstract the collision zone of the 200 km wide and 1 5 km high nazca ridge and the peruvian segment of the convergent south american margin between 14 s and 17 s is characterized by deformation of the upper plate and several hundred meters of uplift of the forearc this is evident by a narrowing of the shelf a westward shift of the coastline and the presence of marine terraces as the nazca ridge is oblique with respect to both trench and convergence direction of the nazca plate it migrates southward along the active plate boundary for reconstructing the migration history of the nazca ridge this study uses updated plate motion data resulting from a revision of the geomagnetic time scale the new model suggests that the ridge crest moved laterally parallel to the margin at a decreasing velocity of 75 mm a before 10 8 ma 61 mm a 10 8 4 9 ma and 43 mm a 4 9 ma to present intra plate deformation associated with mountain building in the peruvian andes since the miocene reduces the relative convergence rate between nazca plate and peruvian forearc taking an intra plate deformation at a rate of 10 mm a estimated from space geodetic and geological data into account does not significantly reduce these lateral migration velocities constraining the length of the original nazca ridge by its conjugate feature on the pacific plate yields a length of 900 km for the subducted portion of the ridge using this constraint ridge subduction began 11 2 ma ago at 11 s therefore the nazca ridge did not affect the northern sites of ocean drilling program odp leg 112 located at 9 s this is supported by benthic foraminiferal assemblages in odp leg 112 cores indicating more than 1000 m of subsidence since at least middle miocene time and by continuous shale deposition on the shelf from 18 to 7 ma recorded in the ballena industrial well at 11 5 s the model predicts the passage of the ridge crest 9 5 ma ago this agrees with the sedimentary facies and benthic foraminiferal stratigraphy of odp leg 112 cores which argue for deposition on the shelf in the middle and late miocene with subsequent subsidence of a minimum of several hundred meters onshore at 12 s the sedimentary record shows at least 500 m uplift prior to the end of the miocene also in agreement with the model

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

response of the mantle to flat slab evolution insights from local s splitting beneath peru the dynamics of flat subduction particularly the interaction between a flat slab and the overriding plate are poorly understood here we study the seismically anisotropic properties and deformational regime of the mantle directly above the peruvian flat slab we analyze shear wave splitting from 370 local s events at 49 stations across southern peru we find that the mantle above the flat slab appears to be anisotropic with modest average delay times 0 28 s that are consistent with 4 anisotropy in a 30 km thick mantle layer the most likely mechanism is the lattice preferred orientation of olivine which suggests that the observed splitting pattern preserves information about the mantle deformation we observe a pronounced change in anisotropy along strike with predominately trench parallel fast directions in the north and more variable orientations in the south which we attribute to the ongoing migration of the nazca ridge through the flat slab system

searching for travertines calcretes and speleothems in deep time processes appearances predictions and the impact of plants common models for modern calcite precipitation in and around caves soils springs and streams involve co 2 supplied by thick high p co 2 biogenic soils which were probably thin or non existent before vascular plants indeed plant influenced chemical weathering might have caused accelerated terrestrial carbonate production from the devonian onwards however terrestrial carbonates have also been documented from the archaean proterozoic cambrian ordovician and silurian mechanisms which could have caused non marine carbonates to precipitate without organic rich soils are described and some geological events likely to have influenced non marine carbonate precipitation up to the origin of vascular plants are highlighted as organisms have evolved so have the petrographic characteristics of non marine carbonates some examples of this are also given here

the lack of correlation between flat slabs and bathymetric impactors in south america flat slab subduction has been attributed to various causes including mantle wedge dynamics overriding by the upper plate age of the subducting plate and subduction of anomalously thick oceanic crust one often favored explanation for flat slabs is the subduction of buoyant features on the oceanic plate in the form of an aseismic ridge or oceanic plateau we show through plate tectonic reconstructions of the marquesas tuamotu and austral plateau assuming that features on the conjugate plate can be used as proxies for subducted bathymetric anomalies that there is very little correlation between the subduction of such anomalies and historic zones of flat subduction in south america it is apparent that subduction of a bathymetric anomaly need not lead to a flat slab and not all flat slabs are associated with the subduction of a bathymetric anomaly

composition of the crust and upper mantle in the central andes 19 30 s inferred from p wave velocity and poisson s ratio the present day crustal structure of the andean chain results from deformation processes and magmatic activity since the onset of subduction one line of investigation to resolve questions about the mechanism of crustal thickening is to get information about the composition of the crust and upper mantle by imaging the structures in velocity and poisson s ratio or vp vs a high density seismic network was operated during six months in 1994 across the whole andean chain the p and s travel times of the best located local earthquakes mainly situated in the subduction zone have been inverted for vp and vp vs using thurber s 3d iterative simultaneous inversion method the resulting models giving short scale information about the crustal and upper mantle composition provide new limits on the respective involvement of crustal shortening and magmatism in the thickening of the andean crust beneath the western cordillera the seismic properties are well explained by typical active volcanic arc processes and the thickening of the crust is explained by magmatic addition conversely beneath the central altiplano basin and the eastern cordillera our results preclude significant volume of magmatic addition from the mantle and reinforce crustal shortening as the main thickening process moreover cratonic lower crust and subcontinental lithospheric upper mantle are observed under the eastern cordillera up to 67 w evidence of the underthrusting of the brazilian craton which drives the crustal shortening beneath the central altiplano basin a thin lithospheric mantle is observed above a well marked asthenospheric wedge which extends in depth below the western part of the eastern cordillera thus partial melt is observed in the shallow mantle west of the western cordillera up to the brazilian craton

three dimensional structure of vp vs and vp vs beneath northeastern japan implications for arc magmatism and fluids we estimated both p v p and s wave velocity v s structures beneath northeastern japan by applying a tomographic method to 169 712 p and 103 993 s wave arrival time data from 4338 local events the average value of v p v s ratio is 1 69 in the upper crust 1 75 in the lower crust and 1 77 in the uppermost mantle these differences in v p v s ratio may mainly reflect the lithological variations with depth low v p low v s and high v p v s zones are extensively distributed along the volcanic front in the uppermost mantle and are downward to the back arc side in the mantle wedge here the low v p low v s and high v p v s values are interpreted because of partial melting materials which imply the presence of a vast amount of melt in the uppermost mantle in the lower crust low v p low v s and high v p v s zones are not continuously visible along the volcanic front but are confined to individual volcanic areas which suggests that melting of lower crustal materials occurs just beneath each active volcano in contrast the upper crust beneath active volcanoes exhibits low v p low v s and low v p v s these results suggest that the partial melting zones under northeastern japan spread out from the uppermost mantle along the volcanic front up to the midcrust right beneath active volcanoes the low v p low v s and low v p v s features in the upper crust suggest the existence of h 2 o rather than melt right beneath active volcanoes

stress in the contorted nazca plate beneath southern peru from local earthquakes we study earthquake focal mechanisms in a region of highly contorted subducting lithosphere to identify dominant sources of stress in the subduction process we observe a stress pattern in the contorted nazca plate beneath southern peru from an analysis of hypocentral trend and focal mechanisms of intermediate depth earthquakes expanding on previous studies we examine the hypocentral trend using 1673 of 2178 well located local events from the nine station arequipa network the dip of the plate beneath southern peru averages 25 30 from 25 to 100 km depth below this depth there is an 80 to 100 km wide contortion between a zone of increasing dip convex to the southeast and a flat lying concave zone to the northwest using more than 6000 p wave first motions of events deeper than 50 km we derive stress orientations from a moving average of composite focal mechanisms across a 200 by 350 km region including the contortion the in plate distribution of tension t and compression p axes reveals a coherent stress pattern the trend is most clear beneath south central peru nw section and below 100 km depth in southernmost peru se section both t and p axes tend to be dominantly in plate especially below 100 km depth t axes orient toward the contortion in a fan shaped trend which suggests that the deepest part of the seismic zone within the convex se section is sinking and pulling the more buoyant nw section we conclude that from 50 to 200 km depth slab pull forces are dominant in the observed stress our results suggest that a significant amount of plate extension occurs in this region of intermediate depth subduction

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

imaging the transition from flat to normal subduction variations in the structure of the nazca slab and upper mantle under southern peru and northwestern bolivia

late cretaceous to recent plate motions in western south america revisited abstract the andean cordillera has evolved since the late cretaceous in the context of subduction of oceanic lithosphere beneath continental lithosphere making the kinematics between south america and its adjacent oceanic plates in the pacific basin valuable to analyze the development of the andean orogen the latest cretaceous cenozoic convergence history in western south america may be divided into three stages the youngest stage 1 25 0 ma is characterized by ene directed convergence of the nazca plate toward most of south america and by e w subduction of the antarctic plate beneath southern patagonia the nazca south america convergence rate in stage 1 shows a continuous decrease from the highest values in the cenozoic 15 cm yr to the present day values from gps measurements 7 cm yr stage 2 47 28 ma is characterized by ne directed subduction of farallon with the convergence rate remaining almost constant during the entire interval in those times obliquity was dextral in chile sinistral in southern peru while almost head on convergence occurred in central and northern peru during latest cretaceous to early eocene times stage 3 the farallon plate was subducted beneath peru and the phoenix plate was subducted farther south where a triple junction migrated southward along the chilean margin the subduction of the farallon plate was rather slow with variable direction imposed by the position of the triple junction whereas subduction of the phoenix plate was rapid 10 cm yr and ese directed we present a working hypothesis suggesting no major changes in the age of subducted lithosphere in the chile trench from middle eocene to late oligocene followed by subduction of progressively older oceanic lithosphere in the early neogene and progressively younger lithosphere during the late neogene and the quaternary in addition it is shown that south american motion as predicted by available hotspot models has insufficient resolution to be applied to the analysis of cenozoic andean deformation

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

nonuniform surface uplift of the andean plateau revealed by deuterium isotopes in miocene volcanic glass from southern peru abstract proposals for rapid late miocene surface uplift driven by large scale lithospheric removal beneath the central andean plateau have been based largely on temperature sensitive paleoaltimeters both the magnitude and mechanism of this proposed pulse of uplift have been challenged first climatic general circulation models support protracted uplift with predicted temperature and isotopic shifts enhanced by attainment of threshold elevations second tectonic models in which surface elevations are compensated by regional contraction and crustal thickening question the need for lithospheric removal and predict broadly coeval uplift of the entire plateau we present hydrogen isotope data using a novel temperature insensitive volcanic glass proxy from continuous well dated lower to middle miocene basin fill in the western cordillera of the northern plateau that show a rapid decrease in d values 62 8 at 19 16 ma with extremely negative values continuing into the pliocene we propose that the basin reached its current elevation by 16 ma 6 myr earlier than proposed for the central plateau the rapid decrease in d values is consistent with punctuated surface uplift of 2 2 3 7 km between 19 and 16 ma whereas the 3 7 km upper estimate assumes a static climate similar to modern the 2 2 km lower estimate conservatively incorporates modeled changes in the isotopic composition of precipitation associated with elevation change comparison of these results to existing paleoelevation estimates from the andean plateau facilitates a tentative reconstruction of earliest middle miocene paleotopography showing a central depression with flanking hinterland and thrust belt highlands this apparent pattern of nonuniform plateau uplift contradicts shortening proportional topographic growth over tens of myr we propose that temporally and spatially irregular surface uplift may be linked to pronounced local variability in crustal shortening and or piecemeal removal of dense mantle lithosphere insofar as these irregularities reflect the scale of heterogeneity for surface uplift processes available paleoelevation estimates point to a principal uplift mechanism that varies over relatively short

high resolution image of the geometry and thickness of the subducting nazca lithosphere beneath northern chile 1 results obtained from s and p receiver functions produced a clear image of the top and bottom of the subducting nazca lithosphere beneath northern chile using data from the teleseismic events recorded at 15 permanent integrated plate boundary observatory chile ipoc stations we obtained new constraints on the geometry and thickness of the descending nazca lithosphere we observed the subducted crust of the nazca plate at depths ranging from 50 km beneath the coastal cordillera down to 110 km beneath the western cordillera we found significant along strike variations in the geometry of the nazca plate beneath northern chile on closer inspection it appears that the oceanic nazca plate is divided into two distinct segments as it descends beneath the continental south american plate the transition from the relatively steeper 23 and deeper slab to the north of 21 s to the flatter southern segment 19 is shown reasonably clearly by our data this feature could well be associated with variations in the curvature of the plate margin and the geometry of the chile trench which is mainly curved to the north of 21 s we have also mapped the continental moho of the south american plate at depths ranging between 60 and 70 km to the east of the longitudinal valley beneath the coastal cordillera this boundary becomes invisible probably due to the serpentinization of the forearc mantle wedge that reduces the velocity in the uppermost mantle the base of the subducted nazca plate was clearly identified as a sharp boundary in the results obtained from the p and s receiver functions the thickness of the subducted oceanic nazca plate which has an age of 50 my is estimated to be 50 km although this thickness is consistent with that predicted by thermal gradients the explanation of the sharpness of the lithosphere asthenosphere boundary may require another mechanism such as hydration or melting

subduction of the nazca ridge and the inca plateau insights into the formation of ore deposits in peru a large number of ore deposits that formed in the peruvian andes during the miocene 15 5 ma are related to the subduction of the nazea plate beneath the south american plate here we show that the spatial and temporal distribution of these deposits correspond with the arrival of relatively buoyant topographic anomalies namely the nazca ridge in central peru and the now consumed inca plateau in northern peru at the subduction zone plate reconstruction shows a rapid metallogenic response to the arrival of the topographic anomalies at the subduction trench this is indicated by clusters of ore deposits situated within the proximity of the laterally migrating zones of ridge subduction it is accordingly suggested that tectonic changes associated with impingement of the aseismic ridge into the subduction zone may trigger the formation of ore deposits in metallogenically fertile suprasubduction environments c 2005 elsevier b v all rights reserved

teleseismic tomography of the southern puna plateau in argentina and adjacent regions abstract we performed a teleseismic p wave tomography study using seismic events at both teleseismic and regional distances recorded by a temporary seismic array in the argentine puna plateau and adjacent regions the tomographic images show the presence of a number of positive and negative anomalies in a depth range of 20 300 km beneath the array the most prominent of these anomalies corresponds to a low velocity body located in the crust most clearly seen in the center of the array 27 s 67 w between the cerro peinado volcano the cerro blanco caldera and the farallon negro in the east this anomaly southern puna magmatic body extends from the northern most part of the array and follows the line with the highest density of stations towards the south where it becomes smaller it is flanked by high velocities on the west and the east respectively on the west the high velocities might be related to the subducted nazca plate on the northeast the high velocity block coincides with the position of the hombre muerto basin in the crust and could be indicating an area of lithospheric delamination where we detected a high velocity block at 100 km depth on the eastern border of the puna plateau north of galan this block might be related to a delamination event in an area with a thick crust of paleozoic metamorphic rocks at the border between puna and eastern cordillera in the center of the array the southern puna magmatic body is also flanked by high velocities but the most prominent region is located on the east and is interpreted as part of the sierras pampeanas lithosphere with high velocities the position of the sierras pampeanas geological province is key in this area as it appears to limit the extension of the plateau towards the south

uplift of the central andean plateau and bending of the bolivian orocline the topography of the central andes can be considered the primary tectonic signal of late cenozoic mountain building in an arid region where the effects of uplift and magmatism are little obscured by denudation the spatial coverage of the topographic signal is more complete than that for sparsely sampled geological and geophysical data a color coded image of digitized topography between 12 s and 37 s highlights the altiplano puna one of the world s most remarkable plateaus and reveals important physiographic clues about the formation of that major feature the topographic data combined with information on structure magmatism seismicity and paleomagnetism support a simple kinematical model for the late cenozoic evolution of the central andes the model does not require collisional effects or enormous volumes of intrusive additions to the crust but instead calls upon plausible amounts of crustal shortening and lithospheric thinning the model interrelates andean uplift a changing geometry of the subducted nazca plate and a changing outline in map view of the leading edge of the south american plate crustal shortening has accommodated convergence between the chilean peruvian forearc and the south american foreland the altiplano puna plateau can be constructed by a combination of crustal shortening and thickening and lithospheric thinning above a shallow dipping 20 30 subducted plate the seawardly concave bend of western south america the bolivian orocline was enhanced but not completely produced by an along strike variation in the amount of late cenozoic shortening maximum shortening in bolivia both produced the widest part of the plateau and increased the seaward concavity of the bolivian orocline the along strike variations of shortening are hypothesized to result from corresponding along strike variations in the width of a weakened zone in the overriding plate weakening occurs above the wedge of asthenosphere located between the subducted and overriding plates hence the width of the zone of weakening depends upon the dip of the subducted plate two types of shortening are recognized 1 a widespread basin and range laramide like shortening that characterizes modem activity in the sierras pampeanas and late miocene deformation of the altiplano puna and 2 on the eastern side of the cordilleras and plateau an east verging foreland fold thrust belt in which the underthrust foreland compresses and thickens the ductile lower crust and produces a plateau uplift of the upper crust the second type of shortening can be applied to plio quaternary deformations throughout the central andes but with a substantial narrowing of the region of plateau uplift in peru and south of 28 s a proposed monoclinal flexure of the upper crust on the western side of the plateau uplift explains the remarkably simple and regular morphology of the main western slope of the central andes the monocline is located above the tip of the asthenospheric wedge between the converging plates it is postulated to occur above the western limit of lower crustal thickening in the regions of horizontal subduction the monocline can be associated with a late miocene asthenospheric wedge tip

shallowing and steepening subduction zones continental lithospheric loss magmatism and crustal flow under the central andean altiplano puna plateau

flat subduction dynamics and deformation of the south american plate insights from analog modeling received 14 june 2007 revised 13 january 2008 accepted 12 march 2008 published 21 june 2008 1 we present lithospheric scale analog models investigating how the absolute plates motion and subduction of buoyant oceanic plateaus can affect both the kinematics and the geometry of subduction possibly resulting in the appearance of flat slab segments and how it changes the overriding plate tectonic regime experiments suggest that flat subductions only occur if a large amount of a buoyant slab segment is forced into subduction by kinematic boundary conditions part of the buoyant plateau being incorporated in the steep part of the slab to balance the negative buoyancy of the dense oceanic slab slab flattening is a long term process 10 ma which requires the subduction of hundreds of kilometers of buoyant plateau the overriding plate shortening rate increases if the oceanic plateau is large enough to decrease the slab pull effect slab flattening increases the interplate friction force and results in migration of the shortening zone within the interior of the overriding plate the increase of the overriding plate topography close to the trench results from 1 the buoyancy of the plate subducting at trench and 2 the overriding plate shortening experiments are compared to the south american active margin where two major horizontal slab segments had formed since the pliocene along the south american subduction zone flat slab segments below peru and central chile nw argentina appeared at 7 ma following the beginning of buoyant slab segments subduction in northern ecuador and northern chile the process of slab flattening resulting from the carnegie and iquique ridges subductions respectively seems to be active but not completed the formation of flat slab segments below south america from the pliocene may explain the deceleration of the nazca plate trenchward velocity citation espurt n f funiciello j martinod b guillaume v regard c faccenna and s brusset 2008 flat subduction dynamics and deformation of the south american plate insights

comment and reply on weathering before the advent of land plants evidence from detrital k feldspars in cambrian ordovician arenites reply

magnetostratigraphy of the kavur tepe section southwestern turkey a magnetic polarity time scale for the norian we present the magnetostratigraphy of the late triassic pelagic kavur tepe section from the taurides southwestern turkey the section investigated is about 30 m thick and consists of fine grained whitish limestones the biostratigraphy based on the conodont zonation gives an age from the lacian 2 zone lower norian to the sevatian 2 zone upper norian the paleomagnetic analyses indicate very favourable behaviour with two magnetic components which can be isolated both by thermal and af demagnetization the magnetization is carried by magnetite and goethite which are clearly evident in irm experiments the first component of normal magnetic polarity has roughly the direction of the present day field the high unblocking temperature component shows directions of both polarities in total fifteen magnetic polarity intervals are observed within the kavur tepe section the characteristic calculated fisherian mean direction is particularly tightly grouped d 220 2 and i 32 1 with k 40 9 95 1 7 and nsamp 179 after bedding correction polarity correlations with two other sections in turkey and in austria northern calcareous alps show that the kavur tepe sediments were deposited in the southern hemisphere by combining the magnetostratigraphic results from the bolucektasi tepe 1 and the kavur tepe sections an almost complete norian magnetic polarity sequence is established showing 30 magnetic intervals that are well calibrated with the conodont and ammonoid zonations a magnetic reversal frequency close to about 2 7 reversals myr is derived for the norian this medium frequency is in agreement with the hypothesis of a long time constant 150 200 myr in the changes of the magnetic reversal frequency since the kiaman reversed superchron

lithologic and climatologic controls of river chemistry abstract the chemistries of rivers draining a variety of lithologic and climatic regions have been surveyed for the purpose of quantifying the fluxes of bicarbonate and silica from rivers with respect to bedrock lithology and runoff in all 101 different rivers each draining a primary lithology were examined across the united states puerto rico and iceland to minimize seasonal effects only rivers with at least two years of data were used basaltic catchments were examined in the greatest detail in a survey of hawaiian island watersheds the average river chemistries could be related to the distribution of soil associations within each catchment an analysis of cation activity relationships among rivers draining basaltic catchments shows that the river compositions define slopes which are consistent with an equilibrium ion exchange control on cation ratios among different lithologies unique weathering rate relationships were developed with yields at typical present day runoff rates 1 100 cm y increasing in the order sandstones granites basalts shales and carbonates the bicarbonate and silica fluxes for each of these lithologies have been quantified for use in global studies of chemical denudation our study confirms that the dissolved yield of a given drainage basin is determined by a balance between physical and chemical weathering thus a warm wet climate or the presence of abundant vegetation cannot guarantee high rates of chemical denudation unless accompanied by high rates of physical removal

interaction between the nazca and south american plates and formation of the altiplano puna plateau review of a flexural analysis along the andean margin 15 34 s abstract the results of a two dimensional flexural analysis applied to the andean margin which is based on the correlation between topography and bouguer anomaly are here reviewed in order to characterize rigidity variations across and along the forearc arc transition of the central andes and to understand the role of the forearc in the formation of the altiplano plateau the forearc has maximum rigidities between 15 and 23 s forearc rigidity decreases gradually southward and sharply toward the plateau the main orogen elevations higher than 3000 m is very weak along the entire central andes a semi quantitative interpretation of these trends based on the relationship between flexural rigidity and the thermo mechanically and compositionally controlled strength of the lithosphere allows the following conclusions to be made 1 across strike rigidity variations are dominated by the thermal structure derived from the subduction process 2 the forearc constitutes a strong cold and rigid geotectonic element 3 southward weakening of the forearc is directly related to the decreasing thermal age of the subducted slab 4 very low rigidities along the main orogen are caused by the existence of a thick quartz rich crust with a low strain rate to heat flow ratio 5 the strength of the plateau lithosphere is localized in an upper crustal layer whose base at 15 km could be correlated with a p to s seismic wave converter trac1 of yuan et al 2000 yuan x sobolev s kind r oncken o et al 2000 subduction and collision processes in the central andes constrained by converted seismic phases nature v 408 21 28 diciembre p 958 961 6 the forearc plateau rigidity boundary corresponds to a zone of changing thermal conditions eastward increasing crustal thickness and felsic component in the crust and low strain rate deformation which correlates with a west verging structural system at the surface these conclusions suggest that the rigid forearc acts as a pseudo indenter against the weak plateau and allows the accumulation of ductile crustal material that moves westward from the eastern foreland this pseudo indenter is geometrically represented by a crustal scale triangular zone rooted at trac1 this model allows the integration of existing contradictory ideas on the dynamics of forearc plateau interaction that are related to the relative importance of upper crustal compressive structures and lower crustal accumulation below the forearc

the impact of the south american plate motion and the nazca ridge subduction on the flat subduction below south peru 1 flat subduction near peru occurs only where the thickened crust of the nazca ridge subducts furthermore the south america continent shows a westward absolute plate motion both the overriding motion of south america and the subduction of the nazca ridge have been proposed to explain the flat slab segment below south peru we have conducted a series of numerical model experiments to investigate the relative importance of both mechanisms results suggest that the average upper mantle viscosity should be about 3 5 1020 pa s or less and basaltic crust should be able to survive 600 to 800 c ambient temperature before transforming into eclogite to explain the slab geometry below peru the effect of the overriding plate is estimated to be as large or twice as large as that of the plateau subduction

seismological constraints on the mechanism of deep earthquakes temperature dependence of deep earthquake source properties abstract seismological observations are of vital importance for understanding the mechanism of deep earthquakes most of the seismological observables for earthquakes deeper than 300 km are similar to shallow earthquakes deep earthquakes clearly represent shear failure on a planar surface as shown by their double couple mechanisms deep earthquake aftershock sequences show temporal decay rates and magnitude frequency relations b values that are similar to shallow earthquakes with the aftershocks occurring preferentially along the mainshock fault planes in addition deep earthquake rupture velocities are similar to those of shallow earthquakes however there are also some observations that are clearly distinctive relative to shallow earthquakes stress drops of deep earthquakes show a large variation but are larger on average than shallow earthquakes different deep seismic zones show very different b values in contrast to shallow earthquakes which show similar b values worldwide in addition deep earthquakes show fewer aftershocks than shallow earthquakes a variety of observations suggest that deep earthquakes are highly sensitive to the temperature of the slab both deep earthquake b values and the rate of deep earthquake aftershock occurrence are inversely correlated with the temperature of the deep slabs suggesting that these factors are temperature controlled to an extent much greater than with shallow earthquakes large deep earthquakes in warm slabs show slower rupture velocities larger stress drops and lower seismic efficiencies than similar earthquakes in cold slabs the width of deep seismic zones is also probably temperature controlled but deep earthquake rupture can propagate outside the normal limits of benioff zone seismicity simple thermal models for the tonga slab near the 1994 deep earthquake suggest that the temperature at the rupture termination point was at least 200 c warmer than the temperature that limits smaller earthquakes in the slab these observations can be used to evaluate physical models for deep earthquakes including brittle slip along fluid weakened faults transformational faulting and thermal and perhaps melt lubricated shear instabilities it is difficult to explain the large fault widths of some deep earthquakes using a fluid weakening model since hydrated materials are expected in only a narrow depth zone at the top of the slab and the fault planes do not have the expected orientations for reactivated faults the lateral extent of the largest deep earthquakes cast doubt on the transformational faulting model in which events should be confined within a narrow metastable wedge seismological studies have also failed to find evidence for the existence of metastable olivine in slabs the temperature dependence of deep earthquakes argues in favor of a temperature activated phenomenon such as thermal shear instabilities and perhaps fault zone melting

plateau style accumulation of deformation southern altiplano 1 employing surface mapping of syntectonic sediments interpretation of industry reflection seismic profiles gravity data and isotopic age dating we reconstruct the tectonic evolution of the southern altiplano 20 22 s between the cordilleras defining its margins the southern altiplano crust was deformed between the late oligocene and the late miocene with two main shortening stages in the oligocene 33 27 ma and middle late miocene 19 8 ma that succeeded eocene onset of shortening at the protoplateau margins shortening rates in the southern altiplano ranged between 1 and 4 7 mm yr with maximum rates in the late miocene summing rates for the southern altiplano and the eastern cordillera we observe an increase from eocene times to the late oligocene to some 8 mm yr followed by fluctuation around this value during the miocene prior to shutoff of deformation at 7 8 ma and transfer of active shortening to the sub andean fold and thrust belt shortening inverted early tertiary graben and half graben systems and was partitioned in three fault systems in the western central and eastern altiplano respectively the east vergent fault systems of the western and central altiplano were synchronously active with the west vergent altiplano west flank fault system from these data and from section balancing we infer a kinematically linked western altiplano thrust belt that accumulated a minimum of 65 km shortening evolution of this belt contrasts with the eastern cordillera which reached peak shortening rates 8 mm yr in between the above two stages hence local shortening rates fluctuated across the plateau superimposed on a general trend of increasing bulk rate with no trend of lateral propagation this observation is repeated at the shorter length and time scales of individual growth structures that show evidence for periods of enhanced local rates at a timescale of 1 3 myr we interpret this irregular pattern of deformation to reflect a plateau style of shortening related to a self organized state of a weak crust in the central south american back arc with a fault network that fluctuated around the critical state of mechanical failure tuning of this state may have occurred by changes in plate kinematics during the paleogene initially reactivating crustal weak zones and by thermal weakening of the crust with active magmatism mainly in the neogene stage

tectonic segmentation of the north andean margin impact of the carnegie ridge collision the north andean convergent margin is a region of intense crustal deformation with six great subduction earthquakes mw 7 8 this century the regional pattern of seismicity and volcanism shows a high degree of segmentation along strike of the andes segments of steep slab subduction alternate with aseismic regions and segments of flat slab subduction this segmentation is related to heterogeneity on the subducting nazca plate in particular the influence of the carnegie ridge collision is investigated four distinct seismotectonic regions can be distinguished region 1 from 6on to 2 5on with steep ese dipping subduction and a narrow volcanic arc region 2 from 2 5on to 1os showing an intermediate depth seismic gap and a broad volcanic arc region 3 from 1os to 2os with steep ne dipping subduction and a narrow volcanic arc region 4 south of 2os with flat subduction and no modern volcanic arc the carnegie ridge has been colliding with the margin since at least 2 ma based on examination of the basement uplift signal along trench parallel transects the subducted prolongation of carnegie ridge may extend up to 500 km from the trench as suggested by the seismic gap and the perturbed broad volcanic arc these findings conflict with previous tectonic models suggesting that the carnegie ridge entered the trench at 1 ma furthermore the anomalous geochemical adakitic signature of the volcanoes in the broad ecuador volcanic arc and the seismicity pattern are proposed to be caused by lithospheric tears separating the buoyant shallowly subducting carnegie ridge from segments of steep subduction in regions 1 and 3 it is further suggested that carnegie ridge supports a local flat slab segment similar to that observed in peru the impact of the carnegie ridge collision on the upper plate causes transpressional deformation extending inboard to beyond the volcanic arc with a modern level of seismicity comparable to the san andreas fault system the pattern of instrumental and historical seismicity indicates 1 great earthquakes on the northern and southern flanks of the colliding ridge 2 a slight reduction in observed seismicity at the trench ridge intersection 3 increased stress far into the continent and 4 a nne displacement of the n andes block to be further effects of the collision 1999 elsevier science b v all rights reserved

three dimensional models of p wave velocity and p to s velocity ratio in the southern central andes by simultaneous inversion of local earthquake data the pisco 94 proyecto de investigation sismologica de la cordillera occidental 1994 seismological network of 31 digital broad band and short period three component seismometers was deployed in northern chile between the coastal cordillera and the western cordillera more than 5300 local seismic events were observed in a 100 day period a subset of high quality p and s arrival time data was used to invert simultaneously for hypocenters and velocity structure additional data from two other networks in the region could be included the velocity models show a number of prominent anomalies outlining an extremely thickened crust about 70 km beneath the forearc region an anomalous crustal structure beneath the recent magmatic arc western cordillera characterized by very low velocities and a high velocity slab a region of an increased vp vs ratio has been found directly above the wadati benioff zone which might be caused by hydration processes a zone of lower than average velocities and a high vp vs ratio might correspond to the asthenospheric wedge the upper edge of the wadati benioff zone is sharply defined by intermediate depth hypocenters while evidence for a double seismic zone can hardly be seen crustal events between the precordillera and the western cordillera have been observed for the first time and are mainly located in the vicinity of the salar de atacama down to depths of about 40 km

seismic structure in southern peru evidence for a smooth contortion between flat and normal subduction of the nazca plate rapid changes in slab geometry are typically associated with fragmentation of the subducted plate however continuous curvature of the slab is also possible the transition from flat to normal subduction in southern peru is one such geometrical change the morphology of the subducted nazca plate along this transition is explored using intraslab earthquakes recorded by temporary regional seismic arrays observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest warping of the slab concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry the presence of a thin ultra slow velocity layer usl atop the horizontal nazca slab is identified and located the lateral extent of this usl is coincident with the margin of the projected linear continuation of the subducting nazca ridge implying a causal relationship wherein increased hydration of the ridge results in the formation of the usl downdip waveform modelling of the 2 d structure in southern peru using a finite difference algorithm provides constraints on the velocity and geometry of the slab s seismic structure and confirms the absence of any tears in the slab the seismicity and structural evidence suggests smooth contortion of the nazca plate along the transition from flat to normal subduction the slab is estimated to have experienced 10 per cent strain in the along strike direction across this transition

subduction of the nazca plate beneath peru as determined from seismic observations a detailed reinvestigation of seismicity in the southern half of peru using data from a local seismic network and a focal mechanism study using p wave first motions reveal that for the first 100 km of descent the slab as defined by this seismicity enters at a normal dip angle near 30 and that below this depth it is bent to a nearly horizontal angle this horizontal slab extends eastward for about 300 km and then dips steeply below the thick continental lithosphere scsp observations are found to be consistent with the profile of the subducting plate as determined from the local seismicity beneath the region from southern peru to northern chile the nazca plate descends with an almost constant dip angle 30 down to at least 300 km depth between the 30 dipping plate in southern peru and the horizontal profile beneath central peru there is a contortion over a 100 km wide continuous lateral section of the subducting plate we infer that at least in the upper 150 km depth range the descending nazca plate is contorted rather than torn

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

local earthquake tomography of the altiplano and the eastern cordillera of northern bolivia the velocity structure in the crest beneath the central andes of northern bolivia is obtained by the simultaneous inversion of local earthquake arrival time data for velocity and hypocentral parameters the data were recorded between october 1990 and february 1991 by the french lithoscope network the 34 vertical short period seismic stations were installed along a 320 km long profile crossing two main morphological units of the andean chain the altiplano and the eastern cordillera the arrival times of seismic waves generated by earthquakes located in the subduction zone of the nazca plate underneath the south american plate were inverted using a modified version of thurber s 3 d iterative simultaneous inversion code the results are consistent with the inversion of teleseismic travel time residuals from the same network and closely correlate with the mapped surface geological features the altiplano basin is characterized by a low velocity zone down to 20 km enclosed by the san andres and laurani faults at the surface the cordillera real fault zone crfz bounds to the west the high p wave velocity anomaly lying under the axial zone of the eastern cordillera the velocity discontinuity dips steeply to the southwest to the bottom of the crust from the velocity model the depth of the moho is estimated to be 66 km under the altiplano moreover the different behaviour of p and s waves found below the western cordillera and above the slab is consistent with the presence of partially melted asthenosphere that may be responsible for the magmatic activity finally the seismicity defines between 90 and 225 kin a part of the slab dipping to the north northeast with a low angle of about 30

crustal structure of the altiplano from broadband regional waveform modeling implications for the composition of thick continental crust we have modeled the full waveforms from six intermediate depth and two shallow earthquakes recorded at regional distances by the banjo broadband andean joint experiment banjo and seismic exploration of the deep altiplano seda portable seismic networks in the central andes in this study we utilize data from those banjo and seda stations located within the altiplano and eastern cordillera we used reflectivity synthetic seismograms and a grid search to constrain four parameters of the altiplano eastern cordillera lithosphere crustal thickness average crustal velocity vp and crustal and upper mantle poisson s ratios crust and mantle using our grid search we investigated the crustal and upper mantle structure along 36 individual event station paths and applied forward modeling to 56 event station paths robust models for the altiplano that provide the best overall fit between the data and synthetic seismograms are characterized by an average vp of 5 75 6 25 km s crustal thicknesses of 60 65 km crust 0 25 and mantle 0 27 0 29 we find a north south variation in the structure of the altiplano with the crust south of the banjo transect characterized by either lower than average crustal p wave velocities or a slightly higher crust relative to crust north of the banjo transect these results are consistent with a model of crustal thickening caused predominantly by tectonic shortening of felsic crust rather than by underplating or magmatic intrusion from the mantle

spatial gaps in arc volcanism the effect of collision or subduction of oceanic plateaus abstract one of the major processes in the formation and deformation of continental lithosphere is the process of arc volcanism the plate tectonic theory predicts that a continuous chain of arc volcanoes lies parallel to any continuous subduction zone however the map pattern of active volcanoes shows at least 24 areas where there are major spatial gaps in the volcanic chains 200 km a significant proportion 30 of oceanic crust is subducted at these gaps all but three of these gaps coincide with the collision or subduction of a large aseismic plateau or ridge the idea that the collision of such features may have a major tectonic impact on the arc lithosphere including cessation of volcanism is not new however it is not clear how the collision or subduction of an oceanic plateau perturbs the system to the extent of inhibiting arc volcanism three main factors necessary for arc volcanism are 1 source materials for the volcanics either volatiles or melt from the subducting slab and or melt from the overlying asthenospheric wedge 2 a heat source either for the dehydration or the melting of the slab or the melting within the asthenosphere and 3 a favorable state of stress in the overlying lithosphere the absence of any one of these features may cause a volcanic gap to form there are several ways in which the collision or subduction of an oceanic plateau may affect arc volcanism the clearest and most common cases considered are those where the feature completely resists subduction causing local plate boundaries to reorganize this includes the formation of new plate bounding transform faults or a flip in subduction polarity in these cases subduction has slowed down or stopped and the lack of source material has created a volcanic gap there are a few cases most notably in peru chile and the nankai trough where the dip of subduction is so shallow that effectively no asthenospheric wedge exists to produce source material for volcanism the shallow dip of the slab may be a buoyant effect of the plateau imbedded in the oceanic lithosphere the cases which are the most enigmatic are those where subduction is continuous the oceanic plateau is subducted along with the slab and the dip of the slab is clearly steep enough to allow arc volcanism yet a volcanic gap exists in these areas the subducted plateau may have a fundamental effect on the physical process of arc volcanism itself the presence of a large topographic feature on the subducting plate may affect the stress state in the are by increasing the amount of decoupling between the two plates alternatively the subduction of the plateau may change the chemical processes at depth if either the water rich top of the plateau with accompanying sediments are scraped off during subduction or if the ridge is compositionally different

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

rise of the andes the surface uplift of mountain belts is generally assumed to reflect progressive shortening and crustal thickening leading to their gradual rise recent studies of the andes indicate that their elevation remained relatively stable for long periods tens of millions of years separated by rapid 1 to 4 million years changes of 1 5 kilometers or more periodic punctuated surface uplift of mountain belts probably reflects the rapid removal of unstable dense lower lithosphere after long term thickening of the crust and lithospheric mantle

slab1 0 a three dimensional model of global subduction zone geometries 1 we describe and present a new model of global subduction zone geometries called slab1 0 an extension of previous efforts to constrain the two dimensional non planar geometry of subduction zones around the focus of large earthquakes slab1 0 describes the detailed non planar three dimensional geometry of approximately 85 of subduction zones worldwide while the model focuses on the detailed form of each slab from their trenches through the seismogenic zone where it combines data sets from active source and passive seismology it also continues to the limits of their seismic extent in the upper mid mantle providing a uniform approach to the definition of the entire seismically active slab geometry examples are shown for two well constrained global locations models for many other regions are available and can be freely downloaded in several formats from our new slab1 0 website http on doi gov d9arbs we describe improvements in our two dimensional geometry constraint inversion including the use of average active source seismic data profiles in the shallow trench regions where data are otherwise lacking derived from the interpolation between other active source seismic data along strike in the same subduction zone we include several analyses of the uncertainty and robustness of our three dimensional interpolation methods in addition we use the filtered subduction related earthquake data sets compiled to build slab1 0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far concluding that the width of these seismogenic zones is on average 30 larger than previous studies have suggested

ridge subduction at an erosive margin the collision zone of the nazca ridge in southern peru the 1 5 km high obliquely subducting nazca ridge and its collision zone with the peruvian margin have been imaged by wide angle and reflection seismic profiles swath bathymetry and gravity surveying these data reveal that the crust of the ridge at its northeastern tip is 17 km thick and exhibits seismic velocities and densities similar to layers 2 and 3 of typical oceanic crust the lowermost layer contributes 10 12 km to the total crustal thickness of the ridge the sedimentary cover is 300 400 m thick on most parts of the ridge but less than 100 m thick on seamounts and small volcanic ridges at the collision zone of ridge and margin the following observations indicate intense tectonic erosion related to the passage of the ridge the thin sediment layer on the ridge is completely subducted the lower continental slope is steep dipping at 9 and the continental wedge has a high taper of 18 tentative correlation of model layers with stratigraphy derived from ocean drilling program leg 112 cores suggests the presence of eocene shelf deposits near the trench continental basement is located 15 km landward of the trench normal faults on the upper slope and shelf indicate extension a comparison with the peruvian and northern chilean forearc systems currently not affected by ridge subduction suggests that the passage of the nazca ridge along the continental margin induces a temporarily limited phase of enhanced tectonic erosion superposed on a long term erosive regime

a teleseismic study of the altiplano and the eastern cordillera in northern bolivia new constraints on a lithospheric model the altiplano one of the world s great plateaus is the most significant geomorphological unit of the central andes knowledge of its lithospheric structure is essential to the understanding of the mountain building in this region in order to study this we performed a teleseismic field experiment in northern bolivia during a 4 month period 34 short period seismic stations were installed along a 320 km long profile crossing the altiplano and the eastern cordillera from the volcanic arc zone to the sub andean zone perpendicular to the main structural direction a careful study of relative travel time residuals led us to propose a qualitative model related to the structural patterns the maximum amplitude of these residuals along the profile reaches 3 s for p phases and have a strong azimuthal dependence the residuals show a noticeable sudden increase clearly associated with the fault system bordering the eastern cordillera to the west the cordillera real fault system the high quality of the data allows us to perform a velocity inversion and to calculate at what depth lateral variations occur the origin of the perturbations in the upper crust is reasonably accounted for by the depth variations of the sedimentary fill with the maximum thickness found under the altiplano basin the velocity perturbations in the lower crust are interpreted as variations of the moho depth which decreases from about 60 km below the altiplano to 50 km below the eastern cordillera in the upper mantle a high velocity zone is observed below the eastern cordillera the major feature of this model is a subvertical boundary dipping slightly to the southwest which separates two strongly contrasting velocity units throughout the model from the surface down to 140 km this boundary coincides with the cordillera real fault system and is interpreted as an old suture as the altiplano is characterized by low velocities in the crust the high velocities beneath the eastern cordillera down to 120 km may correspond to the brazilian craton in the region studied the western limit of the underthrusting of the craton corresponds to the western border of the eastern cordillera

seismic evidence for continuity of the deep slab beneath central and eastern peru seismograms of peru brazil deep focus earthquakes recorded on the nearby carnegie broadband station at cuzco peru cus are characterized by an anomalous large amplitude p wave arrival that follows about 1 5 s after direct p the anomalous arrival is interpreted to be an underside wide angle reflection from the upper surface of the descending nazca plate somewhere in the depth range 150 400 km and in a region where the slab is wholly aseismic particle motion analysis of the cus seismograms shows that the direct p wave has arrival angles as predicted from simple ray tracing the large impulsive anomalous phase termed pr also arrives as a p wave but is of opposite polarity and about 1 5 times the amplitude of the direct p arrival the back azimuth of pr is about 20 clockwise from the path of direct p and its angle of emergence is significantly larger forward modeling calculations indicate that the seismic observations are consistent with a model in which pr originates as a p wave that travels up inside the slab from the source and is then reflected at wide angle from the steeply dipping upper plate boundary results from modeling lead to the following conclusions 1 the descending nazca plate exists in and is probably continuous through the aseismic region between 150 and 500 km depth 2 the plate is very steeply dipping 70 below depths of about 150 or 200 km and 3 the velocity discontinuity at the reflection point must be relatively sharp probably less than about 10 km in thickness

seismic heterogeneity and small scale convection in the southern california upper mantle we invert traveltime residuals of teleseismic p and s phases for 3 d perturbations in vp vs and vp vs structures in the southern california upper mantle the tomographic inversion uses frequency dependent 3 d sensitivity kernels to interpret traveltime residuals measured in multiple frequency bands and recent advances in regional crustal thickness and velocity models to better isolate the mantle component of traveltime residuals in addition to separate vp and vs tomography we jointly invert the p and s data sets for vp vs perturbations by imposing a smoothness constraint on the lnvs lnvp field the regional upper mantle is very heterogeneous with the greatest amplitude of velocity and vp vs perturbations in the upper 200 km and mantle structures that correlate spatially with the major physiographic provinces of southern california our imaging improves the resolution of the geometry and amplitude of these features in addition to the major structures imaged and discussed in previous papers we find a large high velocity anomaly at depths between 340 and 500 km beneath the borderland transverse ranges and northern peninsular ranges this anomaly is separate from the well known uppermost mantle transverse ranges anomaly which extends no deeper than 175 200 km the strongest low velocity high vp vs anomaly is found beneath and west of the salton trough and is attributed to relatively high partial melt in the asthenosphere as a result of lithospheric thinning the magnitudes of vp vs and vp vs perturbations and knowledge of regional 1 d velocity and attenuation lead us to conclude that the asthenosphere contains up to 1 partial melt extending to depths as great as 150 200 km

the structure of the nazca ridge and sala y gomez seamount chain from the dispersion of rayleigh waves summary we report single station measurements of group velocity for fundamental mode rayleigh waves that propagated along two aseismic ridges in the eastern pacific ocean basin the nazca ridge and the sala y gomez chain we find that the nazca ridge dispersion curves when regionalized to isolate the ridge contribution indicate group velocities of about 3 5 km s for short periods these values are slower than predicted for the geometry of the propagation paths interpreting this velocity heterogeneity is complicated by the recognition that the nazca ridge is about 1000 km longer than previously thought but models of vyv for the ridge derived by inverting the dispersion curves reveal abnormally thick 18 f 3 km oceanic crust along its longitudinal axis this value compares favourably with previously published results over the tuamotu plateau thus strengthening the hypothesis that the nazca ridge and tuamotu plateau are mirror images formed by a hotspot coincident with the farallon pacific spreading centre before 20 ma in contrast the measurements along the sala y gomez chain are consistent with typical oceanic crust and mantle when these values are inverted they reveal no thick crust 6 km along the sala y gomez chain we conclude from this result that the chain was likely formed as a leaky fracture zone not as a simple hotspot track

anomalous crust of the bolivian altiplano central andes constraints from broadband regional seismic waveforms a one year deployment of broadband seismographs in the bolivian altiplano recorded numerous intermediate depth earthquakes at near regional distances we modeled the associated broadband waveforms of two earthquakes to estimate an average crustal structure for the altiplano the resulting model is characterized by an anomalously low mean p velocity of 6 0 km s a low poisson s ratio of 0 25 and a crustal thickness of 65 km the combination of the low mean velocity and low poisson s ratio can be explained only by a predominantly quartz rich felsic bulk composition this constraint precludes significant volumes of magmatic addition from the mantle contributing to the great thickness of the altiplano crust but is consistent with thickening by compressive shortening concentrated in a weak felsic layer

the nature of orogenic crust in the central andes 1 the central andes 16 22 s are part of an active continental margin mountain belt and the result of shortening of the weak western edge of south america between the strong lithospheres of the subducting nazca plate and the underthrusting brazilian shield we have combined receiver function and surface wave dispersion results from the banjo seda project with other geophysical studies to characterize the nature of the continental crust and mantle lithospheric structure the major results are as follows 1 the crust supporting the high elevations is thick and has a felsic to intermediate bulk composition 2 the relatively strong brazilian lithosphere is underthrusting as far west 65 5 w as the high elevations of the western part of the eastern cordillera ec but does not underthrust the entire altiplano 3 the subcrustal lithosphere is delaminating piecemeal under the altiplano ec boundary but is not completely removed beneath the central altiplano the altiplano crust is characterized by a brittle upper crust decoupled from a very weak lower crust that is dominated by ductile deformation leading to lower crustal flow and flat topography in contrast in the high relief inland sloping regions of the ec and sub andean zone the upper crust is still strongly coupled across the basal thrust of the fold thrust belt to the underthrusting brazilian shield lithosphere subcrustal shortening between the altiplano and brazilian lithosphere appears to be accommodated by delamination near the altiplano ec boundary our study suggests that orogenic reworking may be an important part of the felsification of continental crust

lithospheric evolution of the andean fold thrust belt bolivia and the origin of the central andean plateau we combine geological and geophysical data to develop a generalized model for the lithospheric evolution of the central andean plateau between 188 and 208 s from late cretaceous to present by integrating geophysical results of upper mantle structure crustal thickness and composition with recently published structural stratigraphic and thermochronologic data we emphasize the importance of both the crust and upper mantle in the evolution of the central andean plateau four key steps in the evolution of the andean plateau are as follows 1 initiation of mountain building by 70 ma suggested by the associated foreland basin depositional history 2 eastward jump of a narrow early fold thrust belt at 40 ma through the eastward propagation of a 200 400 km long basement thrust sheet 3 continued shortening within the eastern cordillera from 40 to 15 ma which thickened the crust and mantle and established the eastern boundary of the modern central andean plateau removal of excess mantle through lithospheric delamination at the eastern cordillera altiplano boundary during the early miocene appears necessary to accommodate underthrusting of the brazilian shield replacement of mantle lithosphere by hot asthenosphere may have provided the heat source for a pulse of mafic volcanism in the eastern cordillera and altiplano at 24 23 ma and further volcanism recorded by 12 7 ma crustal ignimbrites 4 after 20 ma deformation waned in the eastern cordillera and interandean zone and began to be transferred into the subandean zone long term rates of shortening in the fold thrust belt indicate that the average shortening rate has remained fairly constant 8 10 mm year through time with possible slowing 5 7 mm year in the last 15 20 myr we suggest that cenozoic deformation within the mantle lithosphere has been focused at the eastern cordillera altiplano boundary where the mantle most likely continues to be removed through piecemeal delamination d 2005 elsevier b v all rights reserved

electrical conductivity of the pampean shallow subduction region of argentina near 33 s evidence for a slab window 1 we present a three dimensional 3 d interpretation of 117 long period 20 4096 s magnetotelluric mt sites between 31 s and 35 s in western argentina they cover the most horizontal part of the pampean shallow angle subduction of the nazca plate and extend south into the more steeply dipping region sixty two 3 d inversions using various smoothing parameters and data misfit goals were done with a nonlinear conjugate gradient nlcg algorithm a dominant feature of the mantle structure east of the horizontal slab is a conductive plume rising from near the top of the mantle transition zone at 410 km to the probable base of the lithosphere at 100 km depth the subducted slab is known to descend to 190 km just west of the plume but the wadati benioff zone cannot be traced deeper if the slab is extrapolated downdip it slices through the plume at 250 km depth removal of portions of the plume or blocking vertical current flow at 250 km depth significantly changes the predicted responses this argues that the plume is not an artifact and that it is continuous the simplest explanation is that there is a wedge shaped slab window that has torn laterally and opens down to the east with its apex at the plume location stress within the slab and seismic tomography support this shape its northern edge likely explains why there is no deep seismicity south of 29 s

complex anisotropy beneath the peruvian flat slab from frequency dependent multiple phase shear wave splitting analysis flat or shallow subduction is a relatively widespread global occurrence but the dynamics remain poorly understood in particular the interaction between flat slabs and the surrounding mantle flow has yet to be studied in detail here we present measurements of seismic anisotropy to investigate mantle flow beneath the peruvian flat slab segment the largest present day region of flat subduction we conduct a detailed shear wave splitting analysis at a long running seismic station nna located near lima peru we present measurements of apparent splitting parameters fast direction and delay time t for sks scs and local s phases from 80 events we observe well defined frequency dependence and backazimuthal variability indicating the likely presence of complex anisotropy forward modeling the observations with two or three layers of anisotropy reveals a likely layer with a trench normal fast direction underlying a layer with a more trench oblique to trench subparallel fast direction in order to further constrain the anisotropic geometry we analyzed the source side splitting from events originating within the slab measured at distant stations beneath the flat slab segment we found trench normal fast splitting directions in the subslab mantle while within the dipping portion of the slab further to the east likely trench subparallel anisotropy within the slab itself this subslab pattern contradicts observations from elsewhere in south america for normal i e more steeply dipping slab conditions it is similar however to inferences from other shallowly dipping subduction zones around the world while there is an apparent link between slab dip and the surrounding mantle flow at least beneath peru the precise nature of the relationship remains to be clarified

continental and oceanic crustal structure of the pampean flat slab region western argentina using receiver function analysis new high resolution results summary the pampean flat slab of central chile and argentina 30 32 s has strongly influenced cenozoic tectonics in western argentina which contains both the thick skinned basement cored uplifts of the sierras pampeanas and the thin skinned andean precordillera fold and thrust belt in this region of south america the nazca plate is subducting nearly horizontally beneath the south american plate at 100 km depth to gain a better understanding of the deeper structure of this region including the transition from flat to normal subduction to the south three iris passcal arrays of broad band seismic stations have been deployed in central argentina using the dense siembra array combined with the broader charge and esp arrays the flat slab is imaged for the first time in 3 d detail using receiver function rf analysis a distinct pair of rf arrivals consisting of a negative pulse that marks the top of the oceanic crust followed by a positive pulse which indicates the base of the oceanic crust can be used to map the slab s structure depths to moho and oceanic crustal thicknesses estimated from rf results provide new more detailed regional maps an improved depth to continental moho map shows depths of more than 70 km in the main cordillera and 50 km in the western sierras pampeanas that shallow to 35 km in the eastern sierras pampeanas depth to moho contours roughly follow terrane boundaries offshore the hotspot seamount chain of the juan fernandez ridge jfr is thought to create overthickened oceanic crust providing a mechanism for flat slab subduction by comparing synthetic rfs based on various structures to the observed rf signal we determine that the thickness of the oceanic crust at the top of the slab averages at least 13 19 km supporting the idea of a moderately overthickened crust to provide the additional buoyancy for the slab to remain flat the overthickened region is broader than the area directly aligned with the path of the jfr however and indicates along with the slab earthquake locations that the flat slab area is wider than the jfr volcanic chain observed in the offshore bathymetry further rfs indicate that the subducted oceanic crust in the region directly along the path of the subducted ridge is broken by trench parallel faults one explanation for these faults is that they are older structures within the oceanic crust that were created when the slab subducted alternatively it is possible that faults formed recently from tectonic underplating caused by increased interplate coupling in the flat slab region

rapid uplift of the altiplano revealed through 13c 18o bonds in paleosol carbonates the elevation of earth s surface is among the most difficult environmental variables to reconstruct from the geological record here we describe an approach to paleoaltimetry based on independent and simultaneous determinations of soil temperatures and the oxygen isotope compositions of soil waters constrained by measurements of abundances of 13 c 18 o bonds in soil carbonates we use this approach to show that the altiplano plateau in the bolivian andes rose at an average rate of 1 03 0 12 millimeters per year between 10 3 and 6 7 million years ago this rate is consistent with the removal of dense lower crust and or lithospheric mantle as the cause of elevation gain

metastable mantle phase transformations and deep earthquakes in subducting oceanic lithosphere earth s deepest earthquakes occur as a population in subducting or previously subducted lithosphere at depths ranging from about 325 to 690 km this depth interval closely brackets the mantle transition zone characterized by rapid seismic velocity increases resulting from the transformation of upper mantle minerals to higher pressure phases deep earthquakes thus provide the primary direct evidence for subduction of the lithosphere to these depths and allow us to investigate the deep thermal thermodynamic and mechanical ferment inside slabs numerical simulations of reaction rates show that the olivine spinel transformation should be kinetically hindered in old cold slabs descending into the transition zone thus wedge shaped zones of metastable peridotite probably persist to depths of more than 600 km laboratory deformation experiments on some metastable minerals display a shear instability called transformational faulting this instability involves sudden failure by localized superplasticity in thin shear zones where the metastable host mineral transforms to a denser finer grained phase hence in cold slabs such faulting is expected for the polymorphic reactions in which olivine transforms to the spinel structure and clinoenstatite transforms to ilmenite it is thus natural to hypothesize that deep earthquakes result from transformational faulting in metastable peridotite wedges within cold slabs this consideration of the mineralogical states of slabs augments the traditional largely thermal view of slab processes and explains some previously enigmatic slab features it explains why deep seismicity occurs only in the approximate depth range of the mantle transition zone where minerals in downgoing slabs should transform to spinel and ilmenite structures the onset of deep shocks at about 325 km is consistent with the onset of metastability near the equilibrium phase boundary in the slab even if a slab penetrates into the lower mantle earthquakes should cease at depths near 700 km because the seismogenic phase transformations in the slab are completed or can no longer occur substantial metastability is expected only in old cold slabs consistent with the observed restriction of deep earthquakes to those settings earthquakes should be restricted to the cold cores of slabs as in any model in which the seismicity is temperature controlled via the distribution of metastability however the geometries of recent large deep earthquakes pose a challenge for any such models transformational faulting may give insight into why deep shocks lack appreciable aftershocks and why their source characteristics including focal mechanisms indicating localized shear failure rather than implosive deformation are so similar to those of shallow earthquakes finally metastable phase changes in slabs would produce an internal source of stress in addition to those due to the weight of the sinking slab such internal stresses may explain the occurrence of earthquakes in portions of lithosphere which have foundered to the bottom of the transition zone and or are detached from subducting slabs metastability in downgoing slabs could have considerable geodynamic significance metastable wedges would reduce the negative buoyancy of slabs decrease the driving force for subduction and influence the state of stress in slabs heat released by metastable phase changes would raise temperatures within slabs and facilitate the transformation of spinel to the lower mantle mineral assemblage causing slabs to equilibrate more rapidly with the ambient mantle and thus contribute to the cessation of deep seismicity because wedge formation should occur only for fast subducting slabs it may act as a parachute and contribute to regulating plate speeds wedge formation would also have consequences for mantle evolution because the density of a slab stagnated near the bottom of the transition zone would increase as it heats up and the wedge transforms to denser spinel favoring the subsequent sinking of the slab into the lower mantle

did shortening in thick crust cause rapid late cenozoic uplift in the northern bolivian andes abstract this study explores the role of crustal shortening as a mechanism for late cenozoic surface uplift in the northern bolivian andes between c 16 s and c 19 s based on new geological cross sections across the corque corocoro basin in the altiplano c 60 km thick crust today and published radiometric ages of volcanic rocks in the context of a simple 2d crustal balancing model deformation since c 10 ma mainly between 10 and 6 ma involved 28 7 km of shortening across the corque corocoro basin with a wider zone of lower crustal strain and a total of c 70 km of underthrusting in the subandes the uplift model is supported by dynamical analyses and is identical within 2 error to published palaeo elevation proxies with 2 5 0 5 km of surface uplift in the corque corocoro basin since c 10 ma 1 8 0 5 km between 10 and 6 ma and 1 6 0 5 km outside the basin in this way the late cenozoic palaeo elevation record in the northern bolivian andes can be explained within error entirely in terms of the late cenozoic crustal thickening expected for the observed history of crustal shortening

structure of the subduction transition region from seismic array data in southern peru data from three seismic arrays installed in southern peru were analysed using receiver functions from p pp and pkp wave phases in order to image the subducted nazca slab the arrays cover the transition region from flat slab subduction in central peru to normal subduction with an angle of about 30 further south a previous study used data from the first array in the normal subduction region to image the moho depth and slab and showed the existence of a mid crustal structure at 40 km depth that is suggested to be a possible underthrusting of the brazilian shield here we discuss new observations from the other two arrays that span the transition between the two subduction regimes and the flat subduction region the results provide an image of the flattened slab from the coast to approximately 300 km inland and also across the transition region from flat to 30 subduction which appears to be a bend rather than a tear in the slab in the flat slab region the slab is well defined near the coast and flattens out at 100 km depth beneath the altiplano the slab appears to start flattening some 400 km in advance of the subduction of the nazca ridge and the flattening is maintained for 1300 km after its passage the moho begins at a depth of around 30 km near the coast and has a maximum depth of 75 km beneath the altiplano consistent with the results of the other arrays both arrays also show a positive impedance mid crustal structure at 40 km depth which if explainable by underthrusting of the brazilian shield would add further support to the observations from the normal subduction region and show the northward and westward extent of the signal the underthrusting hypothesis would explain the missing crust from the shortening budget needed to support the altiplano the vp vs ratios for both arrays exhibit average values between 1 73 and 1 75 indicating a lesser likelihood that there is a high degree of partial melting or magma bodies at depth in this region the receiver function results provide new imaging of the flat slab and transition from normal to flat slab subduction which allows for comparison of different subduction regimes