earthquakes related to active salt doming near kulyab tadjikistan ussr in the region near kulyab tadjikistan hundreds of shallow earthquakes with magnitudes greater than 2 have been reported in soviet yearly catalogs since 1964 this area appears as a well defined cluster of activity distant from the line of epicenters that defines the gissar kokshal seismic zone to the north of the pamir ranges the geology of this region is dominated by the presence of numerous salt domes surrounded by neogene and quaternary continental deposits the spatial relationship between these earthquakes and the salt domes suggests that the two phenomena may be related moderate earthquakes m 5 occurred in 1972 and 1973 and intensities of surface shaking greater than mm 6 were reported from earthquakes in 1937 1952 1969 1972 1973 and 1978 the earthquake on 2 april 1973 and its aftershocks were located in a region where no salt domes have been mapped at the surface however a buried salt diapir has been mapped at depth by geophysical means these earthquakes may result from active salt diapirism at depth the mechanism for producing this seismicity could be either the active fracturing of the cap rock by the rising diapir or the concentration of tectonic stresses in the thinned section above and adjacent to the diapir the salt related earthquakes may produce lower frequency radiation than other events of the same size

active thrusting and folding along the northern tien shan and late cenozoic rotation of the tarim relative to dzungaria and kazakhstan we have studied geometries and rates of late cenozoic thrust faulting and folding along the northern piedmont of the tien shan mountain belt west of urumqi where the m 8 3 manas earthquake occurred on december 23 1906 the northern range of the tien shan rising above 5000 m overthrusts a flexural foredeep filled with up to 11 000 m of sediment of the dzungarian basement our fieldwork reveals that the active thrust reaches the surface 30 km north of the range front within a 200 km long zone of neogene quaternary anticlines fault scarps are clearest across inset terraces within narrow valleys incised through the anticlines by large rivers flowing down from the range in all the valleys the scarps offset vertically the highest terrace surface by the same amount 10 2 0 7 m inferring an early holocene age 10 2 kyr for this terrace which is continuous with the largest recent fans of the piedmont yields a rate of vertical throw of 1 0 0 3mm yr on the main active thrust at the surface a quantitative morphological analysis of the degradation of terrace edges that are offset by the thrust corroborates such a rate and yields a mass diffusivity of 5 5 2 5 m 2 kyr a rather fresh surface scarp 0 8 0 15 m high that is unlikely to result from shallow earthquakes with 6 m 7 in the last 230 years is visible at the extremities of the main fold zone we associate this scarp with the 1906 manas earthquake and infer that a structure comprising a deep basement ramp under the range gently dipping flats in the foreland and shallow ramps responsible for the formation of the active fault propagation anticlines could have been activated by that earthquake if so the return period of a 1906 type event would be 850 380 years the small size of the scarp for an earthquake of this magnitude suggests that a large fraction of the slip at depth 2 3 is taken up by incremental folding near the surface comparable earthquakes might activate flat detachments and ramp anticlines at a distance from the front of other rising quaternary ranges such as the san gabriel mountains in california or the mont blanc aar massifs in the alps we estimate the finite cenozoic shortening of the folded dzungarian sediments to be of the order of 30 km and the cenozoic shortening rate to have been 3 1 5 mm yr assuming comparable shortening along the tarim piedmont and minor additional active thrusting within the mountain belt we infer the rate of shortening across the tien shan to be at least 6 3 mm yr at the longitude of manas 85 5 e a total shortening of 125 30 km is estimated from crustal thickening assuming local airy isostatic equilibrium under the same assumption serial n s sections imply that cenozoic shortening across the belt increases westwards to 203 50 km at the longitude of kashgar 76 e as reflected by the westward increase of the width of the belt this strain gradient implies a clockwise rotation of tarim relative to dzungaria and kazakhstan of 7 2 5 around a pole located near the eastern extremity of the tien shan west of hami 96 e 43 5 n comparable to that revealed by paleomagnetism between tarim and dzungaria 8 6 8 7 a 6 mm yr rate of shortening at the longitude of manas would imply a rate of rotation of 0 45 m y and would be consistent with a shortening rate of 12 mm yr north of kashgar taking such values to be representative of late cenozoic rates would place the onset of reactivation of the tien shan by the india asia collision in the early to middle miocene 16 22 9 m y in accord with the existence of particularly thick late neogene and quaternary deposits such reactivation would thus have started much later than the collision roughly at the time of the great mid miocene changes in tectonic regimes denudation and sedimentation rates observed in southeast asia the himalayas and the bay of bengal and of the correlative rapid change in seawater sr isotopic ratio 20 to 15 ma like these other changes the rise of the tien shan might be a distant consequence of the end of indochina s escape

exhumational history of the north central pamir 1 the pamir plateau forms a prominent tectonic salient along the western end of the tibet tarim margin despite its tectonic significance relatively little is known about the timing of major cenozoic tectonic events in the pamir here we present new apatite and zircon u th he ages bulk rock geochemistry and al in hornblende barometry results from the karakul graben a prominent north south oriented rift basin located 50 km south of the main pamir thrust although cooling ages do not record the onset of extension graben bounding normal faults provide exposures of otherwise slowly eroding rocks which record two cenozoic thermal events existing geochronology and new results suggest that granitic rocks in the karakul region were shallowly emplaced cooled very quickly through 300 c and have experienced less than 10 km of exhumation since the late triassic a long period of relatively slow exhumation throughout much of the late mesozoic and cenozoic was punctuated by two periods of accelerated exhumation during the middle eocene 50 40 ma and early miocene 25 16 ma we interpret the first period of accelerated exhumation as a result of tectonic uplift and subsequent erosion due to the northward propagation of the india asia collision we attribute the second period of rapid exhumation to a renewed phase of tectonism and plateau uplift in the pamir perhaps related to a break off event along the down going indian plate at 25 ma or to the onset of slip along the nascent karakoram fault

tsr origin of sulfur in permian and triassic reservoir bitumen east sichuan basin china solid bitumen h2s pyrite and anhydrite samples from the lower triassic feixianguan formation and upper permian changxing formation and source rock kerogen samples from the cambrian ordovician silurian and permian in the east sichuan basin were analyzed for their sulfur content and d 34 s values to determine if the measurements were consistent with a thermochemical sulfate reduction tsr origin of the produced h2s and solid bitumen the results show that the solid bitumen samples have sulfur contents from 2 6 to 12 5 wt and d 34 s values from 13 4 to 22 1 the d 34 s values are significantly heavier than those of potential upper permian source rock kerogens that average 1 3 n 3 suggesting that most of the sulfur in the solid bitumen is derived from isotopically heavy h2s of tsr origin from 12 0 to 13 7 interestingly two gases produced from gas to water transition zones during drill stem testing dst show significantly lighter d

convergence rate across the nepal himalaya and interseismic coupling on the main himalayan thrust implications for seismic hazard we document geodetic strain across the nepal himalaya using gps times series from 30 stations in nepal and southern tibet in addition to previously published campaign gps points and leveling data and determine the pattern of interseismic coupling on the main himalayan thrust fault mht the noise on the daily gps positions is modeled as a combination of white and colored noise in order to infer secular velocities at the stations with consistent uncertainties we then locate the pole of rotation of the indian plate in the itrf 2005 reference frame at longitude 1 34 3 31 latitude 51 4 0 3 with an angular velocity of 0 5029 0 0072 myr the pattern of coupling on the mht is computed on a fault dipping 10 to the north and whose strike roughly follows the arcuate shape of the himalaya the model indicates that the mht is locked from the surface to a distance of approximately 100 km down dip corresponding to a depth of 15 to 20 km in map view the transition zone between the locked portion of the mht and the portion which is creeping at the long term slip rate seems to be at the most a few tens of kilometers wide and coincides with the belt of midcrustal microseismicity underneath the himalaya according to a previous study based on thermokinematic modeling of thermochronological and thermobarometric data this transition seems to happen in a zone where the temperature reaches 350 c the convergence between india and south tibet proceeds at a rate of 17 8 0 5 mm yr in central and eastern nepal and 20 5 1 mm yr in western nepal the moment deficit due to locking of the mht in the interseismic period accrues at a rate of 6 6 0 4 10 19 nm yr on the mht underneath nepal for comparison the moment released by the seismicity over the past 500 years including 14 m w 7 earthquakes with moment magnitudes up to 8 5 amounts to only 0 9 10 19 nm yr indicating a large deficit of seismic slip over that period or very infrequent large slow slip events no large slow slip event has been observed however over the 20 years covered by geodetic measurements in the nepal himalaya we discuss the magnitude and return period of m 8 earthquakes required to balance the long term slip budget on the mht

near surface diagnostics of dripping or delaminating lithosphere 1 in various geological regions it has been postulated that the mantle lithosphere has been thinned or completely removed two of the primary removal mechanisms that have been put forward include 1 delamination a wholesale peeling away of a coherent block of the mantle lithosphere and 2 lithospheric dripping viscous rayleigh taylor instability of the mantle lithosphere using computational models we investigate several near surface observables to determine if these may be diagnostic of either often ambiguous removal mechanism surface topography associated with delamination has a broad region of uplift above the lithospheric gap and a localized and mobile zone of subsidence at the delaminating hinge with dripping lithosphere the topographic expression is symmetric and fixed above the downwelling delamination of mantle lithosphere is more efficient than dripping for thermal heating of the crust the onset is more rapid and the elevated temperatures persist longer the resultant crustal p t t paths show modest pressure variations and high temperature increases with large scale delamination or dripping delamination also causes contraction directly above the migrating hinge and distal extension dripping lithosphere induces superimposed contraction and extension above and symmetric about the viscous instability for all the observables if only a portion of the mantle lithosphere is removed by viscous instability delamination inherently removes all of the mantle lithosphere the differences between the two removal mechanisms are even more pronounced with only partial removal of the mantle lithosphere uppermost mantle lithosphere remains well coupled to the crust leading to lower surface temperature variations and broad zones of crustal deformation thickening

the case for dynamic subsidence of the u s east coast since the eocene the dynamic subsidence of the united states east coast is addressed using the discrepancy between regional and global estimates of sea level elevation of paleoshorelines and adjoint models of mantle convection that assimilate plate motions and seismic tomography the positions of eocene and miocene paleoshorelines are lower than predicted by global sea levels suggesting at least 50 m and possibly as much as 200 m of subsidence since the end of the eocene dynamic models predict subsidence of the east coast since the end of eocene although the exact magnitude is uncertain this subsidence has been occurring during an overall global sea level fall with the eustatic change being larger than the dynamic subsidence this results in a regional sea level fall in the absence of land subsidence dynamic subsidence is consistent with the difference between eustasy and regional sea level at the new jersey coastal plain

lateral extrusion in the eastern alps part 1 boundary conditions and experiments scaled for gravity lateral extrusion encompasses extensional collapse gravitational spreading away from a topographic high in an orogenic belt and tectonic escape plane strain horizontal motion of wedges driven by forces applied to their boundaries in the eastern alps it resulted from 1 an overall northerly compression apulia against eurasia 2 a strong foreland bohemian massif 3 lack of constraint along a lateral boundary carpathian region and 4 a previously thickened gravitationally unstable thermally weakened crust eastern alpine orogenic belt six indentation experiments reproduce lateral extrusion at lithospheric scale the models have two to four lithospheric layers with a mohr coulomb rheology for the upper and a viscous rheology for the lower crust the lithosphere rests upon a low viscosity asthenosphere a broad indenter a narrow deformable area and a weakly constrained eastern margin fullfill as closely as possible conditions in the eastern alps indentation produces both thickening in front of the indenter and escape of triangular wedges lateral variations in crustal thickness become attenuated by gravitational spreading the overall fault pattern includes domains of reverse strike slip oblique normal and pure normal faults strike slip faults in conjugate sets develop serially the narrow width of the deformable area and the strength of the foreland determine the angles between the sets gravitational spreading produces a rhombohedral pattern of oblique and pure normal faults along the unconstrained margin opposite the unconstrained margin the indenter front shows thrusts and folds intersecting with the conjugate strike slip sets a triangular indenter favors spreading high velocity of indentation favors escape high confinement limits lateral motion inhibits spreading and favors thickening lateral extrusion in the eastern alps is best modeled by 1 a weak lateral confinement 2 a broad and straight indenter 3 a narrow width of the deformable area and 4 a rigid foreland crustal thickening lateral escape and gravitational spreading all contribute to the overall deformation

deep burial of asian continental crust beneath the pamir imaged with local earthquake tomography an inclined zone of intermediate depth seismicity beneath the pamir orogen in central asia has been interpreted as southward subduction of a slab of asian lithosphere however it is not known whether asian lithosphere subducts intact or only partially we used arrival times of shallow and intermediate depth earthquakes recorded with a temporary 2008 2010 seismic network in this region to invert for 3d models of seismic velocities in an attempt to answer this question with local seismicity reaching depths of up to 240 km the deep structure of the pamir could be illuminated with high resolution the resulting velocity models show a north south contrast in crustal seismic velocities in the pamir with very low p velocities 5 7 5 9 km s at 15 30 km depth coupled with relatively low vp vsvp vs 1 70 at mid crustal levels in the southern part of the orogen at sub moho depths we image an arcuate high velocity 8 2 8 6 km s slab dipping south in the eastern pamir and east in the pamir s southwest underlying the intermediate depth earthquakes on top of the high velocity slab and just above the onset of deep seismicity between a depth of 60 to 100 km very low compressional wavespeeds around 7 1 km s and high vp vsvp vs ratios 1 80 attest to subducted crustal rocks additionally we carried out 2d numerical thermomechanical modeling of the continental collision in the pamir focusing on the fate of the crust and mantle lithosphere of the asian and indian plates seismic velocities were computed from the modeling results and the resulting images were compared with the velocity distributions obtained from seismic traveltimes combining tomography and modeling results we infer that a substantial amount of crustal material is pulled down beneath the pamir by cold mantle lithosphere to depths of at least 80 100 km from there on only lower crust and mantle lithosphere continue their subduction and earthquakes occur inside the lower crustal layer probably due to metamorphic reactions the buoyant asian upper and middle crust does not penetrate deeper into the mantle but pools at this depth level from where it might eventually exhume or relaminate

on the temporal evolution of long wavelength mantle structure of the earth since the early paleozoic the seismic structure of the earth s lower mantle is characterized by a dominantly degree 2 pattern with the african and pacific large low shear velocity provinces i e llsvp that are separated by circum pacific seismically fast anomalies it is important to understand the origin of such a degree 2 mantle structure and its temporal evolution in this study we investigated the effects of plate motion history and mantle viscosity on the temporal evolution of the lower mantle structure since the early paleozoic by formulating 3 d spherical shell models of thermochemical convection for convection models with realistic mantle viscosity and no initial structure it takes about 50 myr to develop dominantly degree 2 lower mantle structure using the published plate motion models for the last either 120 ma or 250 ma however it takes longer time to develop the mantle structure for more viscous mantle while the circum pangea subduction in plate motion history models promotes the formation of degree 2 mantle structure the published pre pangea plate motions before 330 ma produce relatively cold lower mantle in the african hemisphere and significant degree 1 structure in the early pangea 300 ma or later times even if the lower mantle has an initially degree 2 structure and a viscosity as high as 1023 pas this suggests that the african llsvp may not be stationary since the early paleozoic with the published plate motion models and lower mantle viscosity of 1022 pas our mantle convection models suggest that the present day degree 2 mantle structure may have largely been formed by 200 ma

thermoelasticity of a mineral composite and a reconsideration of lower mantle properties a re examination of thermal and elastic properties of a mineral composite emphasising derivative properties especially pressure and temperature dependences of incompressibility resolves residual doubts about the application to the lower mantle of thermodynamic estimates of these properties the only counter intuitive conclusion is that ks p s 4 21 at p 0 which is significantly higher than the values for component minerals this leads to a revised lower mantle equation of state and also to a new one which requires as input only k0 4 21 and k 1 425 k k0 k k0 k p k these equations and three other finite strain theories fit the lower mantle prem tabulation very well and agree on the extrapolation to p 0 k0 203 1 1 6 gpa 0 128 2 1 2 gpa 0 3972 5 4 0 kg m 3 a further six finite strain equations are rejected as incompatible with homogeneity of the lower mantle application of the theory to thermal properties gives s 1 ks ks t p decreasing with depth from 2 7 to 1 6 over the lower mantle depth range conversely e 1 t p increases with depth from 7 1 to 8 7 corresponding values of lnvs lnvp p vary from 1 7 to 2 3 this conclusion is not consistent with a simple thermal interpretation of tomographic observations of lower mantle velocity anomalies compositional variations must be invoked

physical and chemical signatures of a developing anticyclonic eddy in the leeuwin current eastern indian ocean a multidisciplinary cruise aboard the r v southern surveyor was conducted in may 2006 to sample a developing anticyclonic eddy of the leeuwin current off western australia the eddy formed from a meander of the leeuwin current in mid april 2006 and remained attached to the current until mid august in this study a combination of satellite data altimeter sea surface temperature and chlorophyll a and shipboard measurements acoustic doppler current profiler and conductivity temperature depth were used to characterize the physical and chemical signatures of the eddy the temperature salinity properties of the mixed layer waters within the anticyclonic eddy and on the shelf were both connected to that of the leeuwin current indicating the water mass in the eddy is mainly derived from the leeuwin current and the modified leeuwin current water on the shelf above the salinity maximum near the eddy center there was a regionally significant concentration of nitrate 0 9 mol l 1 and the maximum 2 mol l 1 was at 150 in depth below the photic zone nitrification within the eddy and or local upwelling due to the forming eddy could be responsible for this high concentration of nitrate near the eddy center which potentially makes the eddy a relatively productive feature in the leeuwin current

late cenozoic slip on the talas ferghana fault the tien shan central asia although cenozoic crustal shortening and thickening by thrust faulting have built the present tien shan active right lateral shear on the northwest trending talas ferghana fault appears to be the most rapid localized deformation in the belt ephemeral stream valleys have been offset right laterally tens of metres new and published radiocarbon dates of organic material deposited in depressions blocked by offset ridges place upper bounds on the average holocene slip rate at 18 localities uncertainties allow 14 upper bounds to overlap the range of 8 16 mm yr and 95 confidence limits on such bounds at 11 sites are entirely within this range we infer that the rate of 10 mm yr is not simply an upper bound but applies to the late holocene epoch although the bounds on rates permit more rapid slip in the northwest than the southeast they do not place a useful constraint on variations in slip rate along the fault offsets of paleozoic facies boundaries imply a total right lateral shear of 180 250 km but early cretaceous sedimentary rock appears to have been offset only 60 10 km published paleomagnetic declinations of cretaceous miocene rock demonstrate 20 30 of counterclockwise rotation of the ferghana valley which lies just west of the talas ferghana fault with respect to stable parts of eurasia and 20 11 with respect to the central tien shan east of the fault these declinations are consistent with a maximum northwestward translation of 70 210 km of the ferghana valley at the talas ferghana fault and therefore with a similar maximum horizontal shortening across the chatkal ranges which lie between the ferghana valley and the kazakh platform estimates of crustal thickness beneath the chatkal ranges however permit only 60 100 km of cenozoic shortening if 100 km of slip on the talas ferghana fault accumulated at a constant rate of 10 mm yr it would imply an initiation of slip more recently than ca 10 ma long after india collided with the rest of eurasia

on the deep mantle origin of the deccan traps the deccan traps in west central india constitute one of earth s largest continental flood basalt provinces whose eruption played a role in the cretaceous paleogene extinction event the unknown mantle structure under the indian ocean at the start of the cenozoic presents a challenge for connecting the event to a deep mantle origin we used a back and forth iterative method for time reversed convection modeling which incorporates tomography based present day mantle heterogeneity to reconstruct mantle structure at the start of the cenozoic we show a very low density deep seated upwelling that ascends beneath the reunion hot spot at the time of the deccan eruptions we found a second active upwelling below the comores hot spot that likely contributed to the region of partial melt feeding the massive eruption

joint seismic geodynamic mineral physical modelling of african geodynamics a reconciliation of deep mantle convection with surface geophysical constraints recent progress in seismic tomography provides the first complete 3 d images of the combined thermal and chemical anomalies that characterise the unique deep mantle structure below the african continent we present a tomography based model of mantle convection that provides an excellent match to fundamental surface geodynamic constraints on internal density heterogeneity that includes both compositional and thermal contributions where the latter are constrained by mineral physics the application of this thermochemical convection model to the problem of african mantle dynamics yields a reconciliation of both surface gravity and topography anomalies to deep seated mantle flow under the african plate over a wider range of wavelengths than has been possible to date on the basis of these results we predict flow in the african asthenosphere characterised by a clear pattern of focussed upwellings below the major centres of late cenozoic volcanism including the kenya domes hoggar massif cameroon volcanicline cape verde and canary islands the flow predictions also reveal a deep seated large scale active hot upwelling below the western margin of africa under the cape verde islands that extends down to the core mantle boundary the scale and dynamical intensity of this west african superplume is comparable to the south african superplume that has long been assumed to dominate the large scale flow dynamics in the deep mantle under africa we evaluate the plausibility of the predicted asthenospheric flow patterns through a comparison with seismic azimuthal anisotropy derived from independent analyses of african shear wave splitting data

seismic domains within the gissar kokshal seismic zone soviet central asia by evaluating seismological data from the belt of shallow intraplate seismicity known as the gissar kokshal seismic zone in soviet central asia in conjunction with recent data on the shallow crustal structure and longer term geologic deformation we show that certain size classes of earthquakes are characteristic of adjacent and in some places overlapping geologic tectonic settings we call these settings seismic domains they are both geologically and seismologically distinct i e they contain different rocks and structures and have different frequency magnitude relationships for example in terms of number of events gissar kokshal seismicity is dominated by the microearthquakes that have occurred within the sedimentary rocks of the vakhsh fold thrust belt these are all shallow and have occurred mostly within the frontal portions of the thrust sheets more rarely moderate size m 5 5 earthquakes have ruptured the thrust ramps on the other hand the major strain release in this region has accompanied the large 5 5 m 8 1 shallow earthquakes on basement faults in the southern gissar range adjacent to the fold thrust belt within these ranges there is relatively little microearthquake activity that dichotomy indicates that a significantly different seismic hazard is associated with each of these areas but the seismological distinction is obscured along the front of the fold thrust belt where the domain of small shocks in the sedimentary rocks overlaps the domain of large earthquakes on basement faults to explain the localization of the large basement events along the southern flank of the fold thrust belt adjacent to the fold thrust belt we hypothesize a mechanism of thrust sheet loading thrust loads are greatest in this region since thrust sheets have overridden large basement block uplifts beneath the illiac and peter the first range faults

role of dynamic topography in sustaining the nile river over 30 million years the nile is the longest river on earth and has persisted for millions of years it has been suggested that the nile in its present path is 6 million years old whereas others argue that it may have formed much earlier in geological history here we present geological evidence and geodynamic model results that suggest that the nile drainage has been stable for 30 million years we suggest that the nile s longevity in essentially the same path is sustained by the persistence of a stable topographic gradient which in turn is controlled by deeper mantle processes we propose that a large mantle convection cell beneath the nile region has controlled topography over the last 30 million years inducing uplift in the ethiopian yemen dome and subsidence in the levant sea and northern egypt we conclude that the drainage system of large rivers and their evolution over time can be sustained by a dynamic topographic gradient the path of the river nile has been stable for as long as 30 million years sustained by mantle convection according to geophysical and geological evidence and geodynamic model simulations

new space geodetic constraints on the distribution of deformation in central asia abstract results from repeated measurements across a 90 sites global positioning system network in central asia provide direct evidence of current high rates of tectonic deformation far north of the india eurasia suture zone the data help to quantify both the partitioning of deformation within the seismically active tien shan and northern pamirs as well as the ongoing rotation of virtually undeformed blocks such as the tarim with respect to stable eurasia the nnw to sse shortening rates derived from the computations reach values of 23 3 mm yr 1 over the ranges studied suggesting that regional deformation rates have increased considerably within the study area since the onset of cenozoic shortening 20 25 myr ago

evidence for mechanical coupling and strong indian lower crust beneath southern tibet how surface deformation within mountain ranges relates to tectonic processes at depth is not well understood the upper crust of the tibetan plateau is generally thought to be poorly coupled to the underthrusting indian crust because of an intervening low viscosity channel here however we show that the contrast in tectonic regime between primarily strike slip faulting in northern tibet and dominantly normal faulting in southern tibet requires mechanical coupling between the upper crust of southern tibet and the underthrusting indian crust such coupling is inconsistent with the presence of active channel flow beneath southern tibet and suggests that the indian crust retains its strength as it underthrusts the plateau these results shed new light on the debates regarding the mechanical properties of the continental lithosphere and the deformation of tibet

time dependent convection models of mantle thermal structure constrained by seismic tomography and geodynamics implications for mantle plume dynamics and cmb heat flux summary one of the outstanding problems in modern geodynamics is the development of thermal convection models that are consistent with the present day flow dynamics in the earth s mantle in accord with seismic tomographic images of 3 d earth structure and that are also capable of providing a time dependent evolution of the mantle thermal structure that is as realistic earth like as possible a successful realization of this objective would provide a realistic model of 3 d mantle convection that has optimal consistency with a wide suite of seismic geodynamic and mineral physical constraints on mantle structure and thermodynamic properties to address this challenge we have constructed a time dependent compressible convection model in 3 d spherical geometry that is consistent with tomography based instantaneous flow dynamics using an updated and revised pseudo spectral numerical method the novel feature of our numerical solutions is that the equations of conservation of mass and momentum are solved only once in terms of spectral green s functions we initially focus on the theory and numerical methods employed to solve the equation of thermal energy conservation using the green s function solutions for the equation of motion with special attention placed on the numerical accuracy and stability of the convection solutions a particular concern is the verification of the global energy balance in the dissipative compressible mantle formulation we adopt such validation is essential because we then present geodynamically constrained convection solutions over billion year timescales starting from present day seismically constrained thermal images of the mantle the use of geodynamically constrained spectral green s functions facilitates the modelling of the dynamic impact on the mantle evolution of 1 depth dependent thermal conductivity profiles 2 extreme variations of viscosity over depth and 3 different surface boundary conditions in this case mobile surface plates and a rigid surface the thermal interpretation of seismic tomography models does not provide a radial profile of the horizontally averaged temperature i e the geotherm in the mantle one important goal of this study is to obtain a steady state geotherm with boundary layers which satisfies energy balance of the system and provides the starting point for more realistic numerical simulations of the earth s evolution we obtain surface heat flux in the range of earth like values 37 tw for a rigid surface and 44 tw for a surface with tectonic plates coupled to the mantle flow also our convection simulations deliver cmb heat flux that is on the high end of previously estimated values namely 13 tw and 20 tw for rigid and plate like surface boundary conditions respectively we finally employ these two end member surface boundary conditions to explore the very long time scale evolution of convection over billion year time windows these billion year scale simulations will allow us to determine the extent to which a memory of the starting tomography based thermal structure is preserved and hence to explore the longevity of the structures in the present day mantle the two surface boundary conditions along with the geodynamically inferred radial viscosity profiles yield steady state convective flows that are dominated by long wavelengths throughout the lower mantle the rigid surface condition yields a spectrum of mantle heterogeneity dominated by spherical harmonic degree 3 and 4 and the plate like surface condition yields a pattern dominated by degree 1 our exploration of the time dependence of the spatial heterogeneity shows that for both types of surface boundary condition deep mantle hot upwellings resolved in the present day tomography model are durable and stable features these deeply rooted mantle plumes show remarkable longevity over very long geological time spans mainly owing to the geodynamically inferred high viscosity in the lower mantle

late cenozoic extension and crustal doming in the india eurasia collision zone new thermochronologic constraints from the ne chinese pamir 1 the northward motion of the pamir indenter with respect to eurasia has resulted in coeval thrusting strike slip faulting and normal faulting the eastern pamir is currently deformed by east west oriented extension accompanied by uplift and exhumation of the kongur shan 7719 m and muztagh ata 7546 m gneiss domes both domes are an integral part of the footwall of the kongur shan extensional fault system kes a 250 km long north south oriented graben why active normal faulting within the pamir is primarily localized along the kes and not distributed more widely throughout the orogen has remained unclear in addition relatively little is known about how deformation has evolved throughout the cenozoic despite refined estimates on present day crustal deformation rates and microseismicity which indicate where crustal deformation is presently being accommodated to better constrain the spatiotemporal evolution of faulting along the kes we present 39 new apatite fission track zircon u th sm he and 40ar 39ar cooling ages from a series of footwall transects along the kes graben shoulder combining these data with present day topographic relief 1 d thermokinematic and exhumational modeling documents successive stages rather than synchronous deformation and gneiss dome exhumation while the exhumation of the kongur shan commenced during the late miocene extensional processes in the muztagh ata massif began earlier and have slowed down since the late miocene we present a new model of synorogenic extension suggesting that thermal and density effects associated with a lithospheric tear fault along the eastern margin of the subducting alai slab localize extensional upper plate deformation along the kes and decouple crustal motion between the central western pamir and eastern pamir tarim basin

an enhanced image of the pamir hindu kush seismic zone from relocated earthquake hypocentres summary we determine the shape of the seismic zone in the pamir hindu kush region de ned by 30 42 0 n 68 78 0 e by obtaining improved hypocentral locations with 90 per cent con dence limits of less than 30 km the depth error bar for most of the earthquakes of about 6000 shallow and intermediate depth earthquakes available s and depth phase arrival times are also used together with the p wave arrival times in the joint hypocentre determination technique to obtain the best possible hypocentral locations the study region is divided into three depth ranges 0 60 60 160 and 160 km the 0 60 km depth zone is then subdivided laterally into 19 blocks with the deeper regions divided into two blocks each the improved delineation of the seismic zone obtained by using the relocated hypocentres implies that the intermediate depth seismicity in the pamir hindu kush region is most simply explained by a single s shaped seismic zone 700 km long and no more than 30 km wide and with most activity concentrated at 100 300 km depth the main features observed are 1 the eastward steepening of the north dipping hindu kush seismic zone through to its overturning at its eastern end beneath the pamirs where it dips to the southeast 2 the curvature and forking of the subducting slab at depths greater than 200 km within the eastern part of the hindu kush seismic zone 3 the very abrupt cut oi in intermediate depth seismicity at 90 110 km depth with no extension to shallower depths under the pamirs and with a persistent gap between the intermediate and shallow seismicity in the northern pamirs and 4 the unusual horizontal t axes for intermediate depth earthquakes of the pamir seismic zone which align with its curvature this study shows that the seismic zone under the hindu kush has stress axes which follow the classical pattern for subducting slabs controlled by gravity whereas the pamir region has horizontal t axes that follow the trend of the contorted seismic zone this suggests that the pamir seismic zone is a slab deformed due to ow in the upper mantle

gypsum a joint tomographic model of mantle density and seismic wave speeds 1 gypsum is a 3 d model of mantle shear wave s speeds compressional wave p speeds and density the model is developed through simultaneous inversion of seismic body wave travel times p and s and geodynamic observations while using realistic mineral physics parameters linking wave speeds and density geodynamic observations include the global free air gravity field divergence of the tectonic plates dynamic topography of the free surface and the flow induced excess ellipticity of the core mantle boundary gypsum is built with the philosophy that heterogeneity that most closely resembles thermal variations is the simplest possible solution models of the density field from earth s free oscillations have provided great insight into the density configuration of the mantle but are limited to very long wavelength solutions alternatively scaling higher resolution seismic images to obtain density anomalies generates density fields that do not satisfy geodynamic observations the current study provides a 3 d density model for the mantle that directly satisfies geodynamic and seismic observations through a joint seismic geodynamic inversion process notable density field observations include high density piles at the base of superplume structures supporting the general results of past normal mode studies however we find that these features are more localized and have lower amplitude than past studies would suggest when we consider both fast and slow seismic anomalies in gypsum we find that p and s wave speeds are strongly correlated throughout the mantle however we find a low correlation of fast s wave zones in the deep mantle 1500 km depth with the corresponding p wave anomalies suggesting a systematic divergence from simplified thermal effects in ancient subducted slab anomalies the cratonic lithosphere and d regions are shown to have strong compositional signatures however we argue that temperature variations are the primary cause of p wave speed s wave speed and density anomalies throughout most of the mantle

scaling relations of earthquake source parameter estimates with special focus on subduction environment abstract earthquake rupture length and width estimates are in demand in many seismological applications earthquake magnitude estimates are often available whereas the geometrical extensions of the rupture fault mostly are lacking therefore scaling relations are needed to derive length and width from magnitude most frequently used are the relationships of wells and coppersmith 1994 derived on the basis of a large dataset including all slip types with the exception of thrust faulting events in subduction environments however there are many applications dealing with earthquakes in subduction zones because of their high seismic and tsunamigenic potential there are no well established scaling relations for moment magnitude and length width for subduction events within this study we compiled a large database of source parameter estimates of 283 earthquakes all focal mechanisms are represented but special focus is set on large subduction zone events in particular scaling relations were fitted with linear least square as well as orthogonal regression and analyzed regarding the difference between continental and subduction zone oceanic relationships additionally the effect of technical progress in earthquake parameter estimation on scaling relations was tested as well as the influence of different fault mechanisms for a given moment magnitude we found shorter but wider rupture areas of thrust events compared to wells and coppersmith 1994 the thrust event relationships for pure continental and pure subduction zone rupture areas were found to be almost identical the scaling relations differ significantly for slip types the exclusion of events prior to 1964 when the worldwide standard seismic network was established resulted in a remarkable effect on strike slip scaling relations the data do not show any saturation of rupture width of strike slip earthquakes generally rupture area seems to scale with mean slip independent of magnitude the aspect ratio l w however depends on moment and differs for each slip type

a new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data we present new profiles of mantle viscosity derived on the basis of non linear occam style joint inversions of an extensive set of data associated with mantle convection and glacial isostatic adjustment gia the convection related observables include satellite derived free air gravity harmonics the geodetically inferred excess ellipticity of the cmb and tectonic plate motions the gia constraints involve two classes of observables previously shown to be relatively insensitive to errors in the late pleistocene ice history the so called fennoscandian relaxation spectrum frs and a set of site specific decay times determined from the postglacial sea level history in hudson bay and sweden the inverted viscosity profiles show a significant three orders of magnitude increase from the upper mantle mean value of f4 10 20 pa s to a high viscosity 10 23 pa s peak at

tectonic and climatic forcing on the panj river system during the quaternary surface processes involve complex feedback effects between tectonic and climatic influences in the high mountains of pamir the ongoing india asia collision provokes the development of east west trending mountain ranges that impose structural control on flow directions of the pamir rivers the evolving relief is further controlled by strong moisture gradients the decreasing precipitations from the southern and western margins of the pamir plateau to its center in their turn control the emplacement of glaciers chronologies of glacial records from the pamir plateau attest for strong climatic variability during the quaternary corresponding remnants of glacial advances suggest glacial morphodynamic restricted to 4 000 m a s l since marine isotope stage 4 the panj the trunk river of pamir deflects from the predominant westward drainage connecting its main tributaries at the western margin of the drainage basin the geometry of the river network and the pattern of incision characterize the panj as a composite river river reaches of indicated low incision coincide with west trending valleys parallel to domes and their bounding faults valley shape ratios reflect increased incision in north trending sections but do not match with changes in the catchment geometry or erodibility of rock types modelled riverbed profiles distinguish three panj reaches the upstream increase in convexity suggests successive river captures in response to local base level changes the northward deflected river reaches link the local base levels which coincide with the southern boundaries of the shakhdara and yazgulom dome and darvaz range we argue that tectonics plays a large role controlling the drainage system of the panj and hence surface processes in the pamir mountains

erosional processes as a control on the structural evolution of an actively deforming fold and thrust belt an example from the pamir tien shan region central asia the foreland of the pamir mountains displays an unusual orogenic topography in which an active north vergent thin skinned fold and thrust belt shows no topographic foredeep in advance of the thrust sheets instead the thrusts are impinging on a topographic barrier formed by the southern edge of the tien shan the deformation front of this thrust belt is coincident with a deeply incised river valley occupied by a large river the surkhob we suggest that this river played a key role in the evolution of this orogenic system by acting as a lateral conveyor belt for the transport of mass delivered to the deformation front by thrusting on the basis of reconnaissance geologic mapping it is clear that the present system developed in the late miocene or pliocene when the tien shan were uplifted within an existing foreland basin as the tien shan highland developed a longitudinal river within the foreland basin was trapped between the north directed thrusts and the tien shan highland with continued deformation the river was squeezed into a narrow canyon and since that time the system has been in an unusual equilibrium where thrusting acts to push the deformation front northward toward the tien shan but fluvial erosion and transport act to keep it in a constant position geologic relationships indicate that the vakhsh thrust the frontal thrust of the north directed system is now emerging on an erosion surface beneath a thin quaternary cover throughout the canyon suggesting that the river has been unable to keep pace with the tectonic mass flux and the thrust front as well as the river have been pushed northward across the beveled edge of what was originally the southern edge of the tien shan we suggest that the original deformation front of the tien shan now lies 10 km south of the present day topographic front of the tien shan beneath the north vergent vakhsh thrust as indicated by 1 a near vertical zone of intense basement seismicity the petrovsky seismic zone 2 an en echelon array of backthrusts along the south flank of the peter the first range and 3 footwall cutoff relationships of the vakhsh thrust system order of magnitude mass balance considerations based on sediment budgets support this conclusion but suggest that the mass imbalance may be small thus large volumes of material may have passed through the thrust fron and been carried out of the system by erosion a relationship that may account for the narrow 30 50 km wide band of deformed sedimentary rock between the pamir salient and the tien shan the present near equilibrium is maintained by long term tectonic transport rates and fluvial transport rates the latter closely linked to climate particularly precipitation thus although there is apparently a mass pile up within this system other systems of this type might show grossly different structural and stratigraphic relationships depending on the amount of rainfall and corresponding variations in rates of fluvial sediment transport similar settings can be recognized in the active thick skinned thrust systems of the tien shan and ancient systems like the laramide belt of western north america

two stage subduction history under north america inferred from multiple frequency tomography eastward subduction of oceanic tectonic plates has shaped the geologic history of western north america over the past 150 million years the mountain building and volcanism that brought forth the spectacular landscapes of the west are credited to the vast ancient farallon plate which interacted mechanically and chemically with the overlying continent as it plunged back into the mantle here we use finite frequency travel time and amplitude measurements of teleseismic p waves in seven frequency bands to obtain a high resolution tomographic image to approx1 800 km depth we discover several large previously unknown pieces of the plate which show that two distinct stages of whole mantle subduction are present under north america the currently active one descends from the pacific northwest coast to 1 500 km depth beneath the great plains whereas its stalled predecessor occupies the transition zone and lower mantle beneath the eastern half of the continent we argue that the separation between them is linked to the laramide era 70 50 myr ago a time of unusual volcanism and mountain building far inland generally explained by an episode of extremely flat subduction

kongur shan normal fault type example of mountain building assisted by extension karakoram fault eastern pamir the northernmost segment of the karakoram fault is an active normal dextral wrench fault that bounds the muji tashgurgan pliocene quaternary basin eastern pamir the kongurshan 7719 m and mustaghata 7546 m mountains are interpreted as twin crustal ramp anticlines formed en echelon on the eastern side of the karakoram fault by subduction of the tarim lithosphere under the pamir west dipping s c mylonitic gneisses that have downdip stretching lineations form the western flank of kongur mountain kinematic criteria consistently indicate a top to the west sense of shear the mylonitic foliation is cut by a steeper active normal fault which is marked by prominent triangular facets on west facing spurs the ductile normal shearing deformation occurred during the interval 5 1 ma therefore like the himalayas the pamir mountain range shows the coeval development of frontal thrusting and rear normal faulting the kongur fault thus appears to be an excellent example of synorogenic extension where normal faulting contributes to relief

oceanic style subduction controls late cenozoic deformation of the northern pamir orogen abstract the northern part of the pamir orogen is the preeminent example of an active intracontinental subduction zone in the early stages of continent continent collision such zones are the least understood type of plate boundaries because modern examples are few and of limited access and ancient analogs have been extensively overprinted by subsequent tectonic and erosion processes in the pamir it has been assumed that most of the plate convergence was accommodated by overthrusting along the plate bounding main pamir thrust mpt which forms the principal northern mountain and deformation front of the pamir however the synopsis of our new and previously published thermochronologic data from this region shows that the hanging wall of the mpt experienced relatively minor amounts of late cenozoic exhumation the pamir orogen as a whole is an integral part of the overriding plate in a subduction system while the remnant basin to the north constitutes the downgoing plate with the bulk of the convergence accommodated by underthrusting herein we demonstrate that the observed deformation of the upper and lower plates within the pamir alai convergence zone resembles highly arcuate oceanic subduction systems characterized by slab rollback subduction erosion subduction accretion and marginal slab tear faults we suggest that the curvature of the north pamir is genetically linked to the short width and rollback of the south dipping alai slab northward motion indentation of the pamir is accommodated by crustal processes related to this rollback the onset of south dipping subduction is tentatively linked to intense pamir contraction following break off of the north dipping indian slab beneath the karakoram

active tectonics of the pamirs and karakorum the source parameters of moderate to large sized earthquakes and the spatial distribution of earthquake hypocenters are used to investigate the active tectonics and three dimensional configuration of the subducted lithosphere in the pamir karakorum intracontinental convergence zone the continental collision has resulted in intense deformation in the upper crust the deformation is distributed over a broad area and is absorbed by thrust and strike slip faulting along the northern margin of the pamirs where the deformation due to continental collision is most intense the northern pamir thrust and several major fault zones mark the present boundary most of the earthquakes that occurred in these fault zones are shallow crustal events with focal mechanism solutions correlated with the tectonic features of the region along the pamir front thrust faulting dominates while on the western and eastern edges of the pamirs the style of deformation is characterized by oblique thrusting with a component of strike slip motion thrust type events beneath the tadjik depression indicate that both the sedimentary rocks and the basement are involved in shortening right lateral strike slip motion is observed on the eastern edge of the pamirs where they border the tarim basin and in the talas fergana kun lun fault zone a southward dipping seismic zone beneath the pamirs and karakorum indicates that the asian lithosphere has been subducted along the pamir front to a depth of at least 200 km this interpretation is consistent with published geological and geophysical data focal mechanism solutions of some intermediate depth events beneath the pamirs show strike slip faulting with approximately n s horizontally oriented p axes indicating that nearly horizontal compression at the intermediate depth is the predominant mode of deformation a 90 km deep event beneath the karakorum is interpreted as occurring in the subducted asian lithosphere the fact that the p axis of this event is oriented parallel to the descending plate suggests downdip compression within the asian lithosphere the overall tectonics of the region is interpreted as a consequence of the underthrusting asian lithosphere being impinged upon by the shallow northward underthrusting indian lithosphere

earthquake focal mechanisms deformation state and seismotectonics of the pamir tien shan region central asia this paper presents a review of the seismotectonics of the pamir tien shan collision zone in the garm region tajikistan based on geological structure seismicity and focal mechanism solutions the region is dominated by horizontal compression manifested by imbricate low angle thrust faults that separate the upper crust into a series of tectonic sheets these thrust systems verge northward from the tajik depression toward the southern tien shan and southward toward the northern pamir the pattern of seismicity across the region suggests that similar low angle thrust faults exist within the crystalline basement as well in order to reconstruct the present day stress strain state of the region we used data from over 15 000 earthquake focal mechanism solutions for small earthquakes m 1 0 gathered over the 27 year period 1963 1989 the method of reconstruction of the seismotectonic deformation std field involves dividing the data set into small spatial windows and summing the individual focal mechanism solutions to form an average earthquake mechanism tensor for each cell the std state of that cell can then be represented by the orientation of the principal axes of compression and tension the intensity or relative uniformity of the average mechanism and the lode nadai coefficient which defines the relation between the magnitudes of the principal strain components we observe throughout the study area a general predominance of subhorizontal compression manifested in a mixture of thrust and strike slip deformation normal fault deformation is observed in only two small portions of the study area the orientation of principal compression varies from nearly n s in the northern pamir and the southern tien shan to nw se in the peter the first range in general the std field has a clearly developed mosaic structure defined by spatial groupings of nearly uniform orientations of compression and tension axes the boundaries of these groups cannot in most cases be directly correlated with known geological contacts the std structure also shows some variation with depth with more or less uniform std orientation in the upper crust 0 14 km depth changing significantly at greater depths examination of various subsets of the focal mechanism catalog shows the spatial structure of the std field to be largely stable with respect to both time and magnitude the reliability of the reconstructed std field was verified using three approaches 1 comparison of individual focal mechanism determinations for a large subset of the data using both soviet and u s algorithms 2 comparison of std reconstruction results using the two independent focal mechanism catalogs and 3 examination of the effect of spatial sampling on the results the dominance of generally north oriented subhorizontal compression is interpreted to be primarily the result of convergence between the pamir and tien shan ranges in turn caused by the ongoing collision of india and eurasia

quasi reversibility method for data assimilation in models of mantle dynamics summary rapid progress in imaging deep earth structures using seismic tomography and in studies of physical and chemical properties of mantle rocks facilitates research in assimilation of data related to mantle dynamics in this paper we present a new numerical approach for data assimilation which allows for incorporating observations at present and unknown initial conditions in the past for mantle temperature and flow into a 3 d dynamic model in order to determine the initial conditions the dynamic model is described by the backward heat motion and continuity equations the use of the quasi reversibility qrv method implies the introduction into the backward heat equation of the additional term involving the product of a small regularization parameter and a higher order temperature derivative the data assimilation in this case is based on a search of the best fit between the forecast model state and the observations by minimizing the regularization parameter we apply the qrv data assimilation method to restore the evolution of i mantle plumes a synthetic case study and ii the lithospheric slab imaged by teleseismic body wave tomography in the southeastern carpathians for both models the present temperature and mantle flow are assimilated to the geological past and the prominent features of mantle structures are recovered we then model the evolution of the mantle structures forward in time starting from the restored state to the present state and estimate the accuracy of the model predictions the results of the qrv data assimilation are compared to that obtained by the variational var and backward advection data assimilation although the accuracy of the qrv data assimilation is lower than that of the var data assimilation the qrv method does not require any additional smoothing of the input data or filtering of temperature noise as the var method does based on the results and the comparison of the methods we consider the qrv method to be a highly promising approach to assimilation of data related to mantle dynamics

thermochemical structures beneath africa and the pacific ocean large low velocity seismic anomalies have been detected in the earth s lower mantle beneath africa and the pacific ocean that are not easily explained by temperature variations alone the african anomaly has been interpreted to be a northwest southeast trending structure with a sharp edged linear ridge like morphology the pacific anomaly on the other hand appears to be more rounded in shape mantle models with heterogeneous composition have related these structures to dense thermochemical piles or superplumes it has not been shown however that such models can lead to thermochemical structures that satisfy the geometrical constraints as inferred from seismological observations here we present numerical models of thermochemical convection in a three dimensional spherical geometry using plate velocities inferred for the past 119 million years we show that earth s subduction history can lead to thermochemical structures similar in shape to the observed large lower mantle velocity anomalies we find that subduction history tends to focus dense material into a ridge like pile beneath africa and a relatively more rounded pile under the pacific ocean consistent with seismic observations

slip partitioning in the northeast pamir tian shan convergence zone abstract based on a detailed analysis of satellite imagery combined with field geologic and geomorphic observations we have mapped late cenozoic folds and faults in the northeastern pamir tian shan convergence zone it is a unique example to understand intracontinental ongoing mountain building within india eurasia collision system in the front of northeastern pamir our investigations reveal that the nw wnw trending folds display a right stepping en echelon pattern and nw wnw striking faults are mainly characterized by south dipping thrusts with an extensive dextral strike slip component drainage systems across the active faults show a systematic right lateral offset in contrast structural style of the ene trending fold and thrust belts are predominated by south north directed shortening southwest of the tian shan our results also infer that oblique thrusting accommodates as long term dextral slip rate of ca 4 0 mm yr during the late cenozoic time north of the pamir topographic front tectono stratigraphic evidence suggests that the tectonic deformation was initiated at ca 3 5 ma in the study area we suggest that intracontinental mountain building in the pamir tian shan convergence zone should be attributed to the crustal shortening caused by folding and thrusting as well as block rotation related to strike slip faulting within the india eurasia collision system

variations in grain size and viscosity based on vacancy diffusion in minerals seismic tomography and geodynamically inferred mantle rheology the effect of grain size on mantle viscosity is comparable to that of temperature and pressure the current 3 d distribution of grain size in the mantle is however unknown to explore the possible variability of grain size we use the following geodynamic inferences of effective viscosity vacancy diffusion rates in upper mantle minerals and perovskite in the lower mantle lateral variations in temperature derived from seismic images and different geotherms an important outcome of this modeling is a new mapping of lateral viscosity variations throughout the mantle the corresponding 3 d variations in grain size are characterized by 2 orders of magnitude changes we find a correlation between grain size variability in the mantle and absolute viscosity changes with depth our findings suggest that the traditional assumption of arrhenius temperature dependence for vacancy diffusion in the lower mantle is not sufficient to constrain the deformation mechanisms that determine its effective bulk viscosity

dynamic topography and long term sea level variations there is no such thing as a stable continental platform abstract globally uniform i e eustatic sea level trends with time scales approaching 100 myr have been inferred from both seismic and backstripping stratigraphic analysis at a small set of geographic sites that are presumed to lie on stable continental platforms and passive continental margins characterized by simple thermal subsidence histories e g the new jersey margin the western african margin we demonstrate using mantle flow simulations based on high resolution seismic tomography that both the new jersey margin and the conjugate western african margin have been subject to orders of 100 m of dynamic i e flow induced topography change over the last 30 myr we also show that the changing pattern of downwelling mantle flow associated with plate subduction is a significant contributor to the background eustatic sea level trend which is also of order 100 m during the 30 myr time window therefore late cenozoic variations of dynamic topography on these passive margin sites are comparable to the eustatic sea level changes and can partially mask the latter furthermore even if the observed trend could be accurately corrected for local dynamic topography variation the residual eustatic signal does not merely reflect changes in mean spreading rates at mid ocean ridges we conclude that the observed long term sea level variations at so called stable sites cannot be interpreted as eustatic moreover previous analyses that have used long term sea level trends as a proxy for spreading rates and geochemical fluxes must be revisited

plate motions slab dynamics and back arc deformation back arc extension or compression is often related to a particular behavior of the slab dip change and or forward backward migration with respect to the upper plate slabs act either passively when anchored in the deep mantle whereas back arc deformation accommodates the backward or forward motion of the upper plate with respect to the trench or actively when slab pull increasing with slab age forces the trench to rollback in addition to these two mechanisms i e slab anchor and slab rollback numerous observations support the existence of dynamic mantle flow that can exert an overpressure on one side of the slab causing its forward or rearward migration with respect to the arc based on a compilation of upper plate absolute motion trench absolute motion back arc deformation rate upper plate strain regime and slab age for all oceanic subduction zones excluding any kind of collision with continents arcs or plateaus we have examined how the combined effects of these parameters can account for the observed back arc deformations our main results are 1 a global correlation exists between upper plate absolute motion and back arc deformation i e back arc extension when upper plate retreats and vice versa 2 there are as many advancing trenches as retreating ones with trench motion globally limited to 50 mm y 1 furthermore there is no positive correlation between trench retreat and slab age and 3 upper plate absolute motion often fails to explain the back arc deformation rates and the trench motions observed at several subduction zones as tonga new hebrides sandwich or ryukyu we propose that the trench migration of these subduction zones are on the influence of mantle flows we conclude that spontaneous trench rollback related to slab pull is negligible with respect to upper plate motion back arc deformation regime is mostly controlled by the upper plate absolute motion relative to a partly anchored slab trench motion is limited to 50 mm y 1 even if locally mantle flows force slab to move

stick slip advance of the kohat plateau in pakistan in the westernmost himalaya the indian plate is thought to slip beneath the potwar and kohat plateaux on a layer of viscous material in an entirely aseismic manner analysis of insar data from 1992 shows that slip occurred during a rare mw 6 0 earthquake implying that the kohat plateau is locally grounded

least squares restoration of tertiary thrust sheets in map view tajik depression central asia the tajik depression located west of the pamirs and south of the tien shan is a compressional intermontane basin bounded by basement overthrusts and filled with mesozoic and cenozoic sediments the internal structure is typical of a thin skinned tbld and thrust belt kinematic data available in the literature suggest that indentation of the pamirs into asia durino the cenozoic collision of india and asia has been accommodated in various ways within the depression including westward extrusion thickening combined with wrenching along n s tblds and thrusts and counterclockwise rotations these various detbrmation processes can be analyzed and quantified by reconstruction of the predeformed state of the depression however the combina tion of thrusting wrenching and block rotations implies a nonplane deformation which cannot be restored properly using balanced cross sections alone we have therefore developed a numerical method for restoration of stratigraphic surfaces designed for regions of nonplane compressional tectonics the deformed region is represented in map view as a mosaic of fault bounded blocks overlapping each other along the faults blocks are separately unfolded and then numerically packed together by least squares minimization of overlaps yielding fields of finite horizontal translations and rotations about vertical axes to analyze the deformation postdating the collision of india and asia we have restored a stratigraphic surface at the base of the cenozoic first in order to test the numerical method we restored a map that had previously been restored by a purely manual method restoration of a second map drawn from newly available subsurface data leads to geometrical inconsistencies overlaps and gaps which cannot be reduced remain in the restored state after correction of these inconsistencies restoration yields a complex mode of deformation for the depression individual thrust slices have undergone counterclockwise rotations about verti cal axes the magnitudes increasing from west to east up to a maximum of 40 o near the pamirs horizontal shortening is 150 km 35 in the center of the depression and 240 km 85 in the northeastern part between the pamirs and the tien shan strike slip motions commonly not revealed by balancing cross sections are associated with thrusting on faults striking n s inde pendent paleomagnetic data and slip directions on small scale faults provide positive checks on the results of our restoration

reconstruction of the main parameters of old large earthquakes in soviet central asia using the paleoseismogeological method abstract the paleoseismogeological method has been improved on the basis of the frequent occurrence of recent seismogenic dislocations in soviet central asia the possibilities are studied and methods for reconstructing epicentral zones and determining the magnitude and age of large earthquakes in the past are proposed worldwide and regional scales of magnitude determination by dislocations are compared on the basis of concrete data new methods of age determination are cited examples are given in which parameters of previously unknown large paleoearthquakes at the pamirs tien shan boundary zone are studied the first for soviet central asia experience in map making for paleoseismodislocations and paleoearthquakes in one of the seismically active and dangerous areas in soviet central asia is described the existing limitations of the method are analysed

reconstructing the cenozoic evolution of the mantle implications for mantle plume dynamics under the pacific and indian plates article history the lack of knowledge of the initial thermal state of the mantle in the geological past is an outstanding problem in mantle convection the resolution of this problem also requires the modelling of 3 d mantle evolution that yields maximum consistency with a wide suite of geophysical constraints quantifying the robustness of the reconstructed thermal evolution is another major concern to solve and estimate the robustness of the time reversed inverse problem of mantle convection we analyse two different numerical techniques the quasi reversible qrv and the backward advection bad methods our investigation extends over the 65 myr interval encompassing the cenozoic era using a pseudo spectral solution for compressible flow thermal convection in 3 d spherical geometry we find that the two dominant issues for solving the inverse problem of mantle convection are the choice of horizontally averaged temperature i e geotherm and mechanical surface boundary conditions we find in particular that the inclusion of thermal boundary layers that yield earth like heat flux at the top and bottom of the mantle has a critical impact on the reconstruction of mantle evolution we have developed a new regularisation scheme for the qrv method using a time dependent regularisation function this revised implementation of the qrv method delivers time dependent reconstructions of mantle heterogeneity that reveal 1 the stability of pacific and african large low shear velocity provinces llsvp over the last 65 myr 2 strong upward deflections of the cmb topography at 65 ma beneath the north atlantic the south central pacific the east pacific rise epr and the eastern antarctica 3 an anchored deep mantle plume ascending directly under the epr easter and pitcairn hotspots throughout the cenozoic era and 4 the appearance of the transient reunion plume head beneath the western edge of the deccan plateau at 65 ma our reconstructions of cenozoic mantle evolution thus suggest that mantle plumes play a key role in driving surface tectonic processes and large scale volcanism 2014 elsevier b v all rights reserved

gps velocity field for the tien shan and surrounding regions 1 measurements at 400 campaign style gps points and another 14 continuously recording stations in central asia define variations in their velocities both along and across the kyrgyz and neighboring parts of tien shan they show that at the longitude of kyrgyzstan the tarim basin converges with eurasia at 20 2 mm yr nearly two thirds of the total convergence rate between india and eurasia at this longitude this high rate suggests that the tien shan has grown into a major mountain range only late in the evolution of the india eurasia collision most of the convergence between tarim and eurasia within the upper crust of the tien shan presumably occurs by slip on faults on the edges of and within the belt but 1 3 mm yr of convergence is absorbed farther north at the dzungarian alatau and at a lower rate with the kazakh platform to the west the tarim basin is thrust beneath the tien shan at 4 7 mm yr with respect to eurasia the ferghana valley rotates counterclockwise at 0 7 myr 1 about an axis at the southwest end of the valley thus gps data place a bound of 4 mm yr on the rate of crustal shortening across the chatkal and neighboring ranges on the northwest margin of the ferghana valley and they limit the present day slip rate on the right lateral talas ferghana fault to less than 2 mm yr gps measurements corroborate geologic evidence indicating that the northern margin of the pamir overthrusts the alay valley and require a rate of at least 10 and possibly 15 mm yr

extrapolating oceanic age distributions lessons from the pacific region abstractextrapolation of the age distribution of oceanic lithosphere has played a significant role in assessments of variations in global mean spreading rate global mean ocean basin depth and implications for global mean sea level subduction has already removed 50 of oceanic lithosphere younger than 55 7 ma making some level of extrapolation a necessary part of global plate reconstructions an area equal in size to the pacific basin oceanic lithosphere must be extrapolated for ages older than 29 1 ma three modes of extrapolation are identified mode 1 extrapolation uses the preserved history as recorded on one plate to infer the history of the previously adjacent plate this mode of extrapolation is exemplified by the inferred history of the farallon vancouver nazca and cocos plates relative to the pacific plate on which this record is preserved mode 2 involves extrapolation beyond the preserved age extent of a given ridge system no observable data exist that directly constrain the motions bey

landslides triggered by the 1949 khait earthquake tajikistan and associated loss of life abstract earthquake triggered landslides are a major geological hazard in central asia in july 1949 the m7 4 khait earthquake triggered many hundreds of landslides in a mountainous region near the southern limit of the tien shan mountains central tajikistan these landslides involved widespread rock slope failure as well as large numbers of flowslides in loess that mantles the steep slopes of the region in the yasman valley hundreds of loess landslides coalesced to form a massive loess flow est vol 245 mm 3 that travelled up to 20 km on a slope of only 2 in an adjacent valley the khait landslide involved transformation of an earthquake triggered rockslide into a very rapid flow by the entrainment of saturated loess into its movement it travelled 7 41 km over a vertical distance of 1421 m with an estimated average velocity of 30 m s we estimate its volume as 75 mm 3 an order of magnitude less that previously published estimates the khait landslide was simulated using dan the number of casualties due to earthquake triggered landslides in the epicentral region was considerable approximately 4000 people were killed in the yasman valley loess flow as 20 villages kishlaks were overwhelmed in the khait landslide alone we estimate ca 800 people lost their lives when the villages of khait and khisorak were overrun by rapidly moving debris our data indicates that a total of approximately 7200 people were killed by earthquake triggered landslides in the epicentral region of the khait earthquake and that in terms of loss of life the 1949 yasman valley loess flow was one of the most destructive landslides in recent history

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

the giant shakhdara migmatitic gneiss dome pamir india asia collision zone 1 geometry and kinematics 1 cenozoic gneiss domes comprise one third of the surface exposure of the pamir and provide a window into the deep crustal processes of the india asia collision the largest of these are the doubly vergent composite shakhdara alichur domes of the southwestern pamir tajikistan and afghanistan they are separated by a low strain horst top to sse noncoaxial pervasive flow over the up to 4 km thick south pamir shear zone exhumed crust from 30 40 km depth in the 250 80 km shakhdara dome the top to nne alichur shear zone exposed upper crustal rocks in the 125 25 km alichur dome the gunt shear zone bounds the shakhdara dome in the north and records alternations of normal shear and dextral transpression it contributed little to bulk exhumation footwall exhumation along two low angle normal sense detachments resulted in up to 90 km syn orogenic n s extension extension in the southwestern pamir opposes shortening in a fold thrust belt north of the domes and in particular in the tajik depression where an evaporitic decollement facilitated upper crustal shortening gravitational collapse of the pamir plateau margin drove core complex formation in the southwestern pamir and shortening of the weak foreland adjacent to the plateau overall this geometry defines a vertical extrusion scenario comprising frontal and basal underthrusting and thickening and hanging gravitationally driven normal shear in contrast to the himalayan vertical extrusion scenario erosion in the pamir was minor preserving most of the extruded deep crust including the top of the south pamir shear zone at peak elevations throughout the dome

crustal and uppermost mantle velocity structure along a profile across the pamir and southern tien shan as derived from project tipage wide angle seismic data summary utilizing seismic refraction wide angle reflection data from 11 approximately in line earthquakes 2 d p and s velocity models and a poisson s ratio model of the crust and uppermost mantle beneath the southern tien shan and the pamir have been derived along the 400 km long main profile of the tipage tien shan pamir geodynamic program project these models show that the crustal thickness varies from about 65 5 km close to the southern end of the profile beneath the south pamir through about 73 6 km under lake karakul in the north pamir to about 57 7 km 50 km south of the northern end of the profile in the southern tien shan average crustal p velocities are low with respect to the global average varying from 6 26 to 6 30 km s 1 the average crustal s velocity varies from 3 54 to 3 70 km s 1 along the profile and thus average crustal poisson s ratio varies from 0 23 beneath the central pamir in the south central part of the profile to 0 265 towards the northern end of the profile beneath the southern tien shan the main layer of the upper crust extending from about 2 km below the earth s surface to 27 km depth below sea level b s l has average p velocities of about 6 05 6 1 km s 1 except beneath the south central part of the profile where they decrease to around 5 95 km s 1 this is in contrast to the s velocities which range from 3 4 to 3 6 km s 1 and exhibit the highest values of 3 55 3 6 km s 1 where the p velocity is lowest thus for the main layer of the upper crust is 0 26 beneath the profile except beneath the south central part of the profile where it decreases to 0 22 the low value of 0 22 for under the central pamir the along strike equivalent of the qiangtang terrane in tibet is similar to that within the corresponding layer beneath the northern lhasa and southern qiangtang terranes in central tibet and is indicative of felsic rocks rich in quartz in the state the lower crust below 27 km b s l has p velocities ranging from 6 1 km s 1 at the top to 7 1 km s 1 at the base further for this layer is 0 27 0 28 towards the northern end of the profile but is low at about 0 24 beneath the central and southern parts of the profile which is similar to the situation found in the northeast tibetan plateau the low values can be explained by felsic schists and gneisses in the upper part of the lower crust transitioning to granulite facies and possibly also eclogite facies metapelites in the lower part within the uppermost mantle the average p velocity is about 8 10 8 15 km s 1 and is about 0 26 assuming an isotropic situation then a relatively cool 700 800 c uppermost mantle beneath the profile is indicated this would in turn indicate an intact mantle lid beneath the profile an upper mantle reflector dipping from 104 km b s l 120 km from the southern end of the profile to 86 km b s l 155 km from the northern end of the profile has also been identified the preferred model presented here for the crustal and lithospheric mantle structure beneath the pamir calls for nearly horizontal underthrusting of relatively cool indian mantle lithosphere the leading edge of which is outlined by the pamir seismic zone this cool indian mantle lithosphere is overlain by significantly shortening warm asian crust the moho trough that is a feature seen beneath some other orogenic belts for example the alps and the urals beneath the northern pamir may mark the southern tip of the actively underthrusting tien shan crust along the main pamir thrust

late miocene pliocene deceleration of dextral slip between pamir and tarim implications for pamir orogenesis abstract the timing of the late cenozoic collision between the pamir salient and the tien shan as well as changes in the relative motion between the pamir and tarim are poorly constrained the northern margin of the pamir salient indented northward by 300 km during the late cenozoic accommodated by south dipping intracontinental subduction along the main pamir thrust mpt coupled to strike slip faults on the eastern flank of the orogen and both strike slip and thrust faults on the western margin the kashgar yecheng transfer system kyts is the main dextral slip shear zone separating tarim from the eastern pamir with an estimated cumulative offset of 280 km at an average late cenozoic dextral slip rate of 11 15 mm a cowgill 2010 in order to better constrain the slip history of the kyts we collected thermochronologic samples along the eastward flowing deeply incised antecedent tashkorgan yarkand river which crosses the fault system on the eastern flank of the orogen we present 29 new biotite 40 ar 39 ar ages apatite and zircon u th sm he ages and apatite fission track aft analysis combined with published muscovite and biotite 40 ar 39 ar and aft data to create a unique thermochronologic dataset in this poorly studied and remote region we constrain the timing of four major n trending faults the latter three are strands of the kyts the westernmost the kuke fault experienced significant dip slip west side up displacement between 12 and 6 ma to the east within the kyts our new thermochronologic data and geomorphic observations suggest that the kumtag and kusilaf dextral slip faults have been inactive since at least 3 5 ma long term incision rates across the aertashi dextral slip fault the easternmost strand of the kyts are compatible with slow horizontal slip rates of 1 7 5 3 mm a over the past 3 to 5 ma in summary these data show that the slip rate of the kyts decreased substantially during the late miocene or pliocene furthermore miocene present regional kinematic reconstructions suggest that this deceleration reflects the substantial increase of northward motion of tarim rather than a significant decrease of the northward velocity of the pamir

plate motions and stresses from global dynamic models delineating the driving forces behind plate motions is important for understanding the processes that have shaped earth throughout its history however the accurate prediction of plate motions boundary zone deformation rigidity and stresses remains a difficult frontier in numerical modeling we present a global dynamic model that produces a good fit to such parameters by accounting for lateral viscosity variations in the top 200 kilometers of earth together with forces associated with topography and lithosphere structure as well as coupling with mantle flow the relative importance of shallow structure versus deeper mantle flow varies over earth s surface our model reveals where mantle flow contributes toward driving or resisting plate motions furthermore subducted slabs need not act as strong stress guides to satisfy global observations of plate motions and stress

the angular velocities of the plates and the velocity of earth s centre from space geodesy summary using space geodetic observations from four techniques gps vlbi slr and doris we simultaneously estimate the angular velocities of 11 major plates and the velocity of earth s centre we call this set of relative plate angular velocities geodvel for geodesy velocity plate angular velocities depend on the estimate of the velocity of earth s centre and on the assignment of sites to plates most geodetic estimates of the angular velocities of the plates are determined assuming that earth s centre is fixed in an international terrestrial reference frame itrf and are therefore subject to errors in the estimate of the velocity of earth s centre in itrf2005 and itrf2000 earth s centre is the centre of mass of earth oceans and atmosphere cm the velocity of cm is estimated by slr observation of lageos s orbit herein we define earth s centre to be the centre of mass of solid earth ce we determine the velocity of ce by assuming that the portions of plate interiors not near the late pleistocene ice sheets move laterally as if they were part of a rigid spherical cap the geodvel estimate of the velocity of ce is likely nearer the true velocity of cm than are the itrf2005 and itrf2000 estimates because 1 no phenomena can sustain a significant velocity between cm and ce 2 the plates are indeed nearly rigid aside from vertical motion and 3 the velocity of cm differs between itrf2005 and itrf2000 by an unacceptably large speed of 1 8 mm yr 1 the velocity of earth s centre in geodvel lies between that of itrf2000 and that of itrf2005 with the distance from itrf2005 being about twice that from itrf2000 because the geodvel estimates of uncertainties in plate angular velocities account for uncertainty in the velocity of earth s centre they are more realistic than prior estimates of uncertainties geodvel differs significantly from all prior global sets of relative plate angular velocities determined from space geodesy for example the 95 per cent confidence limits for the angular velocities of geodvel exclude those of revel sella et al for 34 of the 36 plate pairs that can be formed between any two of the nine plates with the best constrained motion the median angular velocity vector difference between geodvel and revel is 0 028 myr 1 which is up to 3 1 mm yr 1 on earth s surface geodvel differs the least from the geodetic angular velocities that altamimi et al determine from itrf2005 geodvel s 95 per cent confidence limits exclude 11 of 36 angular velocities of altamimi et al and the median difference is 0 015 myr 1 geodvel differs significantly from nearly all relative plate angular velocities averaged over the past few million years including those of nuvel 1a the difference of geodvel from updated 3 2 myr angular velocities is statistically significant for all but two of 36 angular velocities with a median difference of 0 063 myr 1 across spreading centres eight have slowed down while only two have sped up we conclude that plate angular velocities over the past few decades differ significantly from the corresponding angular velocity averaged over the past 3 2 myr

seismotectonic range front segmentation and mountain belt growth in the pamir alai region kyrgyzstan india eurasia collision zone the trans alai range along the northern perimeter of the pamir region of kyrgyzstan is at the northern edge of the india eurasia collision zone the range defines a section of main pamir thrust which is divided into eastern central and western segments that record differential absorption of plate convergence the 50 km long east west striking central fault segment of the main pamir thrust has dip slip thrust fault offsets with a minimum holocene dip slip rate of 6 mm yr and is linked to the other segments via northwest striking dextral transfer the western transfer faults are well characterized geomorphically and the westernmost records the predominant transition to north vergent thrusting along the western segment via north vergent low angle thrust faults rooted in a steeply south dipping dextral shear zone in contrast the eastern segment deformation is widely distributed and geomorphically less evident but the transfer also takes place in a structurally complex zone asymmetric offset of a regionally correlated terrace geomorphic zonation and the distribution of historic earthquakes suggest that the central segment is mechanically linked with the western segment but not with the eastern segment the progressive closure of the alai valley by the northward advance of thrusting exemplifies the annihilation of an intramontane basin if the high surface uplift rates implied by these geologic observations have been maintained for the last few million years they attest to the high level of tectonic activity and the recency of construction of the trans alai range

seismic imaging of subducting continental lower crust beneath the pamir abstract exhumation of ultra high pressure metamorphic rocks testifies that the continental crust can subduct to significant depth into the mantle despite its buoyancy however direct observation of ongoing subduction of continental crust is rare the pamir is regarded as a possible place of active continental subduction because of the intermediate depth seismicity crustal xenoliths and estimates of crustal shortening versus convergence rates here we present for the first time receiver function images from a passive source seismic array traversing the tien shan and the pamir plateau showing southward subduction of eurasian continental crust in the eastern pamir we observe a southerly dipping 10 15 km thick low velocity zone lvz that extends from 50 km depth near the base of the crust to more than 150 km depth with a dip angle increasing to subvertical while the upper and mid crustal material seems to be shortened and incorporated into the pamir the lower eurasian crust detaches and subducts in its deeper part 80 km the lvz envelopes the intermediate depth earthquakes our observations imply that the complete arcuate intermediate depth seismic zone beneath the pamir traces a slab of subducting eurasian continental lower crust

the relationship between surface kinematics and deformation of the whole lithosphere the variation of mechanical properties with depth in the lithosphere determines the relationship between surface deformation and whole lithosphere deformation hence between surface deformation and whole lithosphere dynamics where viscosity or elastic strength is a continuous function with depth surface deformation can be used to constrain both force balance and rheological parameters where viscosity is discontinuous but the upper crust and mantle lithosphere have comparable maximum values surface deformation can be used to approximate force balance and rheological parameters but tradeoffs mean that estimates of stress and viscosity are effective equivalent values rather than actual values where viscosity is both discontinuous and differs by much more than an order of magnitude between the upper crust and mantle lithosphere information about both force balance and rheology are absent from the surface deformation so surface observations alone are insufficient to estimate either the dynamic or mechanical state of the lithosphere

himalayan tectonics explained by extrusion of a low viscosity crustal channel coupled to focused surface denudation recent interpretations of himalayan tibetan tectonics have proposed that channel flow in the middle to lower crust can explain outward growth of the tibetan plateau1 2 3 and that ductile extrusion of high grade metamorphic rocks between coeval normal and thrust sense shear zones can explain exhumation of the greater himalayan sequence4 5 6 7 here we use coupled thermal mechanical numerical models to show that these two processes channel flow and ductile extrusion may be dynamically linked through the effects of surface denudation focused at the edge of a plateau that is underlain by low viscosity material our models provide an internally self consistent explanation for many observed features of the himalayan tibetan system8 9 10

grain size variations in the earth s mantle and the evolution of primordial chemical heterogeneities 1 chemical heterogeneities with relatively high viscosity are known to be poorly mixed by mantle convection one of the factors which can cause viscosity variations is the grain size analysis of various physical processes in the earth s mantle suggests that variation in the grain size can be extremely large even outside of subducting slabs the largest grains are formed in the early molten earth where ostwald ripening is controlled by liquid state diffusion the smallest grains are formed due to recrystallization by the spinel perovskite phase transformation in downwelling flow and subsequent ostwald ripening the viscosity contrast between primordial large grain size material and recrystallized material is likely to be many orders of magnitude preliminary numerical simulations show that strongly time dependent mantle convection can process primordial large grain size heterogeneities through the spinel perovskite phase boundary only after a few overturns this reduces the viscosity contrasts and leads to rapid mixing by mantle convection a more steady flow in the lower mantle may lead to the formation of stagnant poorly mixed regions which may never have had a chance to cross the 660 km discontinuity and be recrystallized by the spinel perovskite phase transformation heterogeneities can also avoid grain size reduction and mixing if their density is slightly larger than the surrounding mantle in this case they accumulate at the core mantle boundary

descent of the ancient farallon slab drives localized mantle flow below the new madrid seismic zone 1 the great earthquake sequence that occurred in the central mississippi river valley during the winter of 1811 1812 is unprecedented in the historical record of seismicity within stable continental plate interiors we show using viscous flow models based on high resolution seismic tomography that the descent of the ancient farallon slab into the deep mantle beneath central north america induces a highly localized flow directly below the new madrid seismic zone nmsz this localization arises because of structural variability in the farallon slab and the low viscosity of the sub lithospheric upper mantle and it represents a heretofore unrecognized and possibly significant driving mechanism for the enigmatic intraplate seismicity

a new back and forth iterative method for time reversed convection modeling implications for the cenozoic evolution of 3 d structure and dynamics of the mantle

quaternary deformation in the eastern pamirs tadzhikistan and kyrgyzstan active deformation in the eastern pamir of central asia is concentrated on the margins of the orogen with minor deformation within the high terrain along the trans alai mountain front at the northern perimeter of the orogen quaternary thrusting is documented by uplifted pediments now at 500 m above the piedmont holocene fault scarps and large earthquakes with n to nw oriented p axes seismicity in the interior of the orogen outlines a n s belt that includes normal faulting events with e w oriented t axes n s striking active normal faults in the interior lake karakul region are compatible with these earthquakes they define an asymmetric graben with a master fault at the western basin margin in the southern pamirs dextral strike slip faults root in the dextral karakorum fault which bounds the pamirs to the east a mixture of dextral and reverse offsets totalling 135 m in pleistocene terraces and 8 m in late pleistocene holocene deposits demonstrates contemporary transpression indicating average displacement rates of 1 mm yr the concentration of young thrusts along the trans alai the northward migration of thrusting and the scarcity of other large scale shortening features within the eastern pamirs suggest that this part of the orogen moves northward en bloc and causes the progressive annihilation of the intermontane alai valley widespread dextral shear in the eastern pamirs both to the south and north of the extensional karakul depression and combined dextral strike slip and normal faulting in the muji tashgorgan graben of the chinese pamirs are interpreted as localized space accommodation phenomena formed during progressive transfer of compressional deformation along a dextral strike slip deformation zone with extensional stepovers

late cenozoic tectonic development of the intramontane alai valley pamir tien shan region central asia an example of intracontinental deformation due to the indo eurasia collision 1 the pamir indentor of the northwestern himalayan syntaxis is a first order feature demonstrating partly the northward extent of deformation due to the cenozoic indo eurasia collision the alai valley of kyrgyzstan at the northern end of the indentor is a strategically positioned e w trending intramontane basin that constrains the timing and extent of crustal deformation in this area of the collision zone to quantify the convergence accommodated across the alai valley during the late cenozoic we collected structural and stratigraphic field data reviewed existing soviet literature and analyzed migrated seismic reflection profiles to construct and restore two regional sections crosscutting the basin our study suggests that the development progressive closure and conversion of this formerly marine basin into a terrestrial intramontane basin result from two main deformation events 1 distributed north south contraction took place during the late oligocene early miocene accommodated one third to half of the total shortening and was followed by the formation of a regional erosion surface and 2 n s shortening resumed in the mid miocene and continues today during this second episode the thrust front migrated southward localized along the trans alai ranges and failed to reactivate earlier neogene structures horizontal shortening of about 35 across the alai valley implies relatively low strain rates and displacement rates of about 4 18 4 69 10 16 s 1 and 0 66 0 78 mm yr 1 respectively for the last 25 myr our study confirms other regional observations indicating that contractional deformation occurred far in the interior of the asian continent as early as the late oligocene

plate tectonics and aspherical earth structure the importance of poloidal toroidal coupling we derive poloidal flow internal loading green functions for incompressible fluid models of the mantle consisting either of a single constant viscosity spherical shell or of two adjacent spherical shells having different viscosities from the green function we obtain kernels connecting the surface divergence the geoid and the surface topography fields to the lateral density heterogeneity inferred on the basis of the application of seismic tomographic imaging techniques and with these kernels we argue that both the surface divergence and geoid fields consisting of harmonic degrees 2 5 may be reasonably fit with only a factor of 8 viscosity increase at a depth of 1200 km we point out however that the coupling from poloidal to toroidal flow which is required to understand surface velocity spectra may allow these observations to be understood in terms of a viscosity increase at depth which is smaller than required by nonhydrostatic geoid data in the context of pure poloidal models using the observed kinetic energy in the surface plate motions as a constraint allows us to infer a value for the steady state upper mantle viscosity of 2 0 0 5 1021 pa s we comment on the implication of the difference between this number and the somewhat lower value 1 1021 pa s which has been derived on the basis of analyses of signatures of the glacial isostatic adjustment process

geometry of the pamir hindu kush intermediate depth earthquake zone from local seismic data 1 we present new seismicity images based on a two year seismic deployment in the pamir and sw tien shan a total of 9532 earthquakes were detected located and rigorously assessed in a multistage automatic procedure utilizing state of the art picking algorithms waveform cross correlation and multi event relocation the obtained catalog provides new information on crustal seismicity and reveals the geometry and internal structure of the pamir hindu kush intermediate depth seismic zone with improved detail and resolution the relocated seismicity clearly defines at least two distinct planes one beneath the pamir and the other beneath the hindu kush separated by a gap across which strike and dip directions change abruptly the pamir seismic zone forms a thin approximately 10 km width curviplanar arc that strikes east west and dips south at its eastern end and then progressively turns by 90 to reach a north south strike and a due eastward dip at its southwestern termination pamir deep seismicity outlines several streaks at depths between 70 and 240 km with the deepest events occurring at its southwestern end intermediate depth earthquakes are clearly separated from shallow crustal seismicity which is confined to the uppermost 20 25 km the hindu kush seismic zone extends from 40 to 250 km depth and generally strikes east west yet bends northeast toward the pamir at its eastern end it may be divided vertically into upper and lower parts separated by a gap at approximately 150 km depth in the upper part events form a plane that is 15 25 km thick in cross section and dips sub vertically north to northwest seismic activity is more virile in the lower part where several distinct clusters form a complex pattern of sub parallel planes the observed geometry could be reconciled either with a model of two sided subduction of eurasian and previously underthrusted indian continental lithosphere or by a purely eurasian origin of both pamir and hindu kush seismic zones which necessitates a contortion and oversteepening of the latter

assembly of the pamirs age and origin of magmatic belts from the southern tien shan to the southern pamirs and their relation to tibet 1 magmatic rocks and depositional setting of associated volcaniclastic strata along a north south traverse spanning the southern tien shan and eastern pamirs of kyrgyzstan and tajikistan constrain the tectonics of the pamirs and tibet the northern pamirs and northwestern tibet contain the north facing kunlun suture the south facing jinsha suture and the intervening carboniferous to triassic karakul mazar subduction accretion system the latter is correlated with the songpan garze hoh xi system of tibet the kunlun arc is a composite early paleozoic to late paleozoic triassic arc arc formation in the pamirs is characterized by 370 320 ma volcanism that probably continued until the triassic the cryptic tanymas suture of the southern northern pamirs is part of the jinsha suture a massive 227 ma batholith stitches the karakul mazar complex in the pamirs there are striking similarities between the qiangtang block in the pamirs and tibet like tibet the regional structure of the pamirs is an anticlinorium that includes the muskol and sares domes like tibet the metamorphic rocks in these domes are equivalents to the karakul mazar songpan garze system granitoids intruding the qiangtang block yield 200 230 ma ages in the pamirs and in central tibet the stratigraphy of the eastern pshart area in the pamirs is similar to the bangong nujiang suture zone in the amdo region of eastern central tibet but a triassic ocean basin sequence is preserved in the pamirs arc type granitoids that intruded into the eastern pshart oceanic basin arc sequence 190 160 ma and granitoids that cut the southern qiangtang block 170 160 ma constitute the rushan pshart arc cretaceous plutons that intruded the central and southern pamirs record a long lasting magmatic history their zircons and those from late miocene xenoliths show that the most distinct magmatic events were cambro ordovician 410 575 ma triassic 210 250 ma likely due to subduction along the jinsha suture middle jurassic 147 195 ma subduction along rushan pshart suture and mainly cretaceous middle and late cretaceous magmatism may reflect arc activity in asia prior to the accretion of the karakoram block and flat slab subduction along the shyok suture north of the kohistan ladakh arc respectively before india and asia collided the pamir region from the indus yarlung to the jinsha suture was an andean style plate margin our analysis suggests a relatively simple crustal structure for the pamirs and tibet from the kunlun arc in the north to the southern qiangtang block in the south the pamirs and tibet likely have a dominantly sedimentary crust characterized by karakul mazar songpan garze accretionary wedge rocks the crust south of the southern qiangtang block is likely of granodioritic composition reflecting long lived subduction arc formation and cretaceous cenozoic underthrusting

constraining the extent of crust mantle coupling in central asia using gps geologic and shear wave splitting data abstract we have obtained constraints on mechanical crust mantle coupling for tibet and yunnan indo china by comparing the observed surface deformation field inferred from gps and quaternary fault slip rate data with the mantle deformation field inferred from several sks shear wave splitting data sets we first determined whether the anisotropy is dominantly asthenospheric or lithospheric by testing simple models of both types against the observed values of the fast polarization direction which is assumed to be parallel to the horizontal projection of the direction of maximum shear for asthenospheric flow we solved for a best fitting uniform sub asthenospheric velocity model for eurasia the fit however was not satisfactory rms misfit 20 solving for separate uniform flow fields in each region improves the fit 15 for tibet 11 for yunnan although the resulting flow fields are inconsistent with several geophysical and geological constraints and thus considered unlikely we then considered lithospheric models for tibet vertically coherent deformation i e maximum shear direction from surface deformation is parallel to yields an improved match 11 for left lateral shear both the goodness of fit and the dominance of left lateral surface faulting in tibet argue for a lithospheric source of anisotropy the misfit for yunnan is large for either right 53 or left lateral 49 shear which argues for complete crust mantle decoupling in yunnan we show that the fast polarization directions throughout both the tibet and yunnan region can be fit by a single lithospheric dynamic model in which there is strong coupling between crust and mantle beneath tibet but a complete decoupling between crust and mantle beneath yunnan crust this dynamic model predicts left lateral maximum shear directions within the mantle that align with fast polarization directions in both regions 9 these maximum shear directions within the mantle align with the left lateral maximum shear directions in the crustal deformation field in tibet but are not in yunnan our results have the following implications first the coherence between crust and mantle deformation in tibet implies strong crust mantle mechanical coupling since this is the only way that crustal buoyancy forces required to account for the surface deformation field can be transmitted into the mantle this behavior is consistent with a uniform strength lithosphere or strong crust but not with a substantially weaker crust for example it is inconsistent with the popular jelly sandwich rheology and thus precludes behavior such as large scale lower crustal flow in tibet the crust mantle coherence is also incompatible with mantle delamination second crust mantle decoupling within the yunnan lithosphere argues that mantle deformation there is controlled only by boundary conditions crustal buoyancy forces are not transmitted into the mantle beneath yunnan moreover the dynamic model for the mantle shows that the yunnan crust is moving south with respect to the mantle at rates as high as 30 mm yr third there is a fundamental rheological lithospheric transition between tibet and yunnan that may provide a key to understanding this significant orogen

two deep mantle sources for paleocene doming and volcanism in the north atlantic the north atlantic igneous province naip erupted in two major pulses that coincide with the continental breakup and the opening of the north atlantic ocean over a period from 62 to 54 ma the unknown mantle structure under the north atlantic during the paleocene represents a major missing link in deciphering the geodynamic causes of this event to address this outstanding challenge we use a back and forth iterative method for time reversed global convection modeling over the cenozoic era which incorporates models of present day tomography based mantle heterogeneity we find that the paleocene mantle under the north atlantic is characterized by two major low density plumes in the lower mantle one beneath greenland and another beneath the azores these strong lower mantle upwellings generate small scale hot upwellings and cold downwellings in the upper mantle the upwellings are dispersed sources of magmatism and topographic uplift that were active on the rifted margins of the north atlantic during the formation of the naip while most studies of the paleocene evolution of the north atlantic have focused on the proto icelandic plume our cenozoic reconstructions reveal the equally important dynamics of a hot buoyant mantle wide upwelling below the azores

equivalency of geologic and geodetic rates in contractional orogens new insights from the pamir frontal thrust across contractional orogens the equivalency between decadal convergence rates from geodetic gps data and geologic shortening rates at time scales of thousands or millions of years has rarely been documented here we present an example from the northern margin of chinese pamir where the main pamir thrust is tectonically quiescent and recent deformation is concentrated on the pamir frontal thrust pft based on dated and faulted fluvial terraces magnetostratigraphy and mapping the horizontal shortening rate of the pft is 6 7 mm a at time scales of both 18 4 ka and 0 35 ma comparable to the geodetic rate of 6 9 mm a across the same zone implying that modern geodetic rates are a reasonable proxy for geologic rates since 0 35 ma comparing this example with studies in other contractional orogens we conjecture that a match or mismatch of geologic geodetic rates typically depends on the time scale of observation fault geometry and fault mechanics 2012 by the american geophysical union