a bound on the rheology of continental lithosphere using very long baseline interferometry the velocity of south china with respect to eurasia very long baseline interferometric vlbi measurements using the antenna near shanghai china between 1987 and 1994 show that this site moves at 8 0 5 mm yr 1 in the direction n116 5 e 4 1 with respect to vlbi sites in europe on the eurasia plate independent measurements of the site at shanghai with respect to sites on nearly all major plates eurasia north america pacific australia and africa determine this velocity south china translates east southeast relative to eurasia more slowly than 20 of india s roughly 50 mm yr convergence rate with eurasia and hence accommodates only a small fraction of india s penetration into eurasia apparently crustal thickening absorbs most of india s convergence the relatively low speed of south china supports treatments of continental lithosphere as a thin viscous sheet whose strength is governed primarily by power law creep in the upper mantle and not by friction in the crust

lateral extrusion of lower crust from under high topography in the isostatic limit where there is isostasy the rocks of the lower continental crust are subject to an effective lateral pressure gradient equal to the gradient of topographic load whether compensation is in crustal roots or in the mantle the result is a poiseuille flow planar channel flow in the weak lower crust which removes crust from under mountains and smooths and levels the topography assuming cubic power law creep the flux of crust is proportional to the third power of the topographic gradient to the usual arrhenius term to the tenth power of the absolute temperature of the lower crust and to the negative fifth power of the geothermal gradient the result is that any initial condition with an isolated high tends eventually toward a state where the topography is pancake shaped with a flat central plateau and steeper flanks spreading outward when adjacent pancakes merge most of the relief is eliminated flow may either roughen or smooth the moho depending on whether there are lateral density contrasts in the mantle lithosphere or not although it is difficult to find analytic solutions for this evolution it is easily simulated numerically using various published flow laws for plausible lower crustal rocks lateral extrusion is shown to be insignificant under the oceans marginally significant under shields and platforms important under elevated plains and dominant beneath high plateaus and or hot delaminated regions in particular the basin and range province of the western united states cannot maintain short wavelength 100 km moho relief of more than 1 km for times greater than 10 20 m y at most if the tibetan plateau of china has been delaminated as recently suggested then it may flatten even faster reducing 100 km wavelength topographic features to 200 m relief in no more than 0 03 0 13 m y and leaving only features with wavelengths over 400 km at the present therefore the lower crust of tibet may resemble a hydraulic reservoir as previously suggested because the flatness of the moho is self maintained in these regions information on past tectonics is continually lost and the present moho shape cannot be used to balance cross sections for times in the past the flatness of the moho today is not a constraint on the mechanisms of extension or shortening in these regions nor is it evidence for underplating by intrusions it also follows that very intense and inhomogeneous dilational strains may have occurred locally as in the metamorphic core complexes of the basin and range province

constraints from gps measurements on the dynamics of the zone of convergence between arabia and eurasia we investigate the dynamics of deformation in the zone of convergence between the arabian and eurasian plates using a physical model that treats the lithosphere as a thin fluid sheet deforming in response to lateral variations in gravitational potential energy gpe this model has two free material parameters the power law exponent n of the vertically averaged rheology of the lithosphere and the argand number ar which expresses the relative importance of gpe and stresses required to deform the lithosphere with boundary conditions described by three free parameters the model fits the observed deformation as measured by 367 gps velocities with a root mean square residual of 2 4 mm yr we find negligible improvement when variations in material properties are introduced that represent increases or decreases in the lithospheric strength of the central iranian and turkish iranian plateaux and the zagros mountains effective viscosity of the lithosphere ranges from 5 1022 pa s at 10 nanostrain yr to 1022 pa s at 100 nanostrain yr as well as matching the decadal timescale predominantly interseismic strain rates determined from gps the computed distribution of strain rates is also consistent with the distribution and types of earthquake focal mechanism significant historical earthquake activity is seen in regions with strain rates lower than 20 nanostrain yr implying that it is prudent to base assessments of seismic hazard on regional strain rates

calculating velocity fields across plate boundaries from observed shear rates summary a general horizontal velocity field can be expressed uniquely in terms of the velocities at points along a line and an integral function of the 2 d shear rate consequently the relative motion across a plate boundary can be calculated from observed shear rates by taking the velocities along a line within one of the plates to be zero

vertical averages of rheology of the continental lithosphere relation to thin sheet parameters abstract we calculate vertically averaged rheologies for continental lithosphere subject to biaxial states of stress corresponding to extensional compressional and strike slip regimes we assume that brittle deformation occurs by byerlee s law and that ductile deformation occurs by steady state creep mechanisms under these conditions the vertically averaged rheology may be approximated by a single power law over the range of strain rates from 10 13 to 10 17 s 1 this is true not only when either brittle or ductile behavior dominates the strength of the lithospere but also when the contributions of brittle and ductile behavior are comparable the value of the power law exponent n and the strength of the lithosphere depend on the depth and stress difference at the brittle ductile transition the temperature at the moho and the geothermal gradient in the range of vertically integrated stress between 1012 and 1013 n m 1n reaches a minimum of between 5 and 10 under other circumstancesn increases rapidly as brittle behavior dominates over ductile behavior the strength of the lithosphere depends most strongly on the thermal state of the upper mantle but also is influenced by the stress supported in brittle portions of the lithosphere

contemporary crustal deformation in east asia constrained by global positioning system measurements global positioning system gps measurements collected since the early 90 s allow us to derive geodetic velocities at 16 permanent stations in east asia and 68 campaign mode sites in north china the resulting velocity field shows the following 1 contrary to the early inferences that the shanxi rift has accommodated significant right slip motion our results suggest that the rift system at least in its northern part in north china is under ese wnw extension at a rate of 4 2 mm yr the velocity field also suggests that an ese wnw trending left lateral shear zone deforming at a rate of 2 1 mm yr may exist along the north rim of north china at the latitude of 40 n separating actively extending north china in the south from relatively stable mongolia in the north 2 central and east china move at a rate of 8 11 mm yr east southeast with respect to siberia implying that the overall east southeastward motion is the dominant mode of deformation in east china 3 the india plate moves at a rate of 6 1 mm yr slower than the nuvel 1a model prediction relative to the eurasia plate 4 significant eastward motion 20 2 mm yr with respect to siberia occurs in southeastern tibet about half of this eastward motion 11 mm yr is absorbed by structures along the eastern boundary of the tibetan plateau

kinematic model of active deformation in central asia the velocity field of present day deformation in central asia is modelled using a set of four rotating blocks siberia tarim tibet india on a spherical earth a best fit is inverted on the basis of estimated shortening rates across the main thrust zones himalayas tien shan and measured slip rates along the principal strike slip faults altyn tagh and karakorum separating those blocks the fit to the data implies that nearly all the present convergence between india and asia can be accounted for by slip partitioning on these four zones with as much as 50 absorbed by northeastwards extrusion of tibet this suggests that localised deformation governs the present mechanical behaviour of the central asian lithosphere

surface deformation and lower crustal flow in eastern tibet field observations and satellite geodesy indicate that little crustal shortening has occurred along the central to southern margin of the eastern tibetan plateau since about 4 million years ago instead central eastern tibet has been nearly stationary relative to southeastern china southeastern tibet has rotated clockwise without major crustal shortening and the crust along portions of the eastern plateau margin has been extended modeling suggests that these phenomena are the result of continental convergence where the lower crust is so weak that upper crustal deformation is decoupled from the motion of the underlying mantle this model also predicts east west extension on the high plateau without convective removal of tibetan lithosphere and without eastward movement of the crust east of the plateau

block modeling with connected fault network geometries and a linear elastic coupling estimator in spherical coordinates geodetic observations of interseismic deformation provide constraints on the partitioning of fault slip across plate boundary zones the spatial distribution of both elastic and inelastic strain accumulation and the nature of the fault system evolution here we describe linear block theory which decomposes surface velocity fields into four components 1 plate rotations 2 elastic deformation from faults with kinematically consistent slip rates 3 elastic deformation from faults with spatially variable coupling and 4 homogeneous intrablock strain elastic deformation rates are computed for each fault segment in a homogeneous elastic half space using multiple optimal planar cartesian coordinate systems to minimize areal distortion and triangular dislocation elements to accurately represent complex fault system geometry block motions fault slip rates elastic coupling and internal block strain rates are determined simultaneously using a linear estimator with constraints from both geodetically determined velocity fields and geologic fault slip rate estimates we also introduce algorithms for efficiently implementing alternative fault network geometries to quantify parameter sensitivity to nonlinear perturbations in model geometry

dynamics of the india eurasia collision zone we present simple new dynamic calculations of a vertically averaged deviatoric stress field over a depth average of 100 km for asia from geodetic geologic topographic and seismic data a first estimate of the minimum absolute magnitudes and directions of vertically averaged deviatoric stress is obtained by solving force balance equations for deviatoric stresses associated with gravitational potential energy differences within the lithosphere plus a first order contribution of deviatoric stresses associated with stress boundary conditions this initial estimate of the vertically averaged deviatoric stress field is obtained independent of assumptions about the rheology of the lithosphere absolute magnitudes of vertically averaged deviatoric stresses vary between 5 and 40 mpa assuming bulk viscous behavior for the lithosphere the magnitudes of deviatoric stresses together with the magnitudes of strain rates inferred from quaternary fault slip rate and gps data yield vertically averaged effective viscosities for tibet of 0 5 5 1022 pa s compared with 1 2 5 1023 pa s in more rigid areas elsewhere in the region a forward modeling method that solves force balance equations using velocity boundary conditions allows us to refine our estimates of the vertically averaged effective viscosity distribution and deviatoric stress field the total vertically averaged deviatoric stress and effective viscosity field are consistent with a weak lower crust in tibet they are consistent with some eastward motion of tibet and south china lithosphere relative to eurasia and they confirm that gravitational potential energy differences have a profound effect on the spatially varying style and magnitude of strain rate around the tibetan plateau our results for the vertically averaged deviatoric stress argue for a large portion of the strength of the lithosphere to reside within the seismogenic upper crust to get deviatoric stress magnitudes there to be as high as 100 300 mpa in accord with laboratory and theoretical friction experiments indicating that stress drops in earthquakes are small fractions of the total deviatoric stress

evidence for mechanically coupled lithosphere in central asia and resulting implications the recent dramatic increase in seismic anisotropy and surface global positioning system gps data for central asia permits a comprehensive examination of the mantle s role in mountain building a joint analysis of 178 shear wave splitting and 2000 gps observations using a new technique reveals that the crust and lithospheric mantle deform coherently arguing for crust mantle mechanical coupling during deformation the observed spatial variations in anisotropy refl ect the large scale pattern of lithospheric deformation as well as a change in deformational style from simple shear on the tibetan plateau transitioning to pure shear in surrounding regions

tectonics and topography for a lithosphere containing density heterogeneities the purpose of this paper is to clarify the dynamic role of lithospheric density heterogeneities in particular with respect to mountain building and other processes of intraplate deformation density anomalies within or just beneath the lithosphere constitute major sources for tectonic stress fields the product of their magnitude by their depth is shown to characterize their ability to induce deformation this rule of the density moment directly yields the lithospheric thickening or thinning rate when applied to structures of large lateral extent for anomalies of lateral extent that is small in comparison with their depth the deformation is vertically inhomogeneous and has been computed with the help of simple physical models of a stratified viscous newtonian lithosphere the analytical treatment is based on fourier transform continent continent collision thickens not only the crust but the entire lithosphere the cold root underlying a mountain chain induces strong regional compressive stresses able to sustain the mountain bulding process without further help from forces transmitted from far away thus the continental lithosphere is in a somewhat metastable mechanical state adiabatic i e rapid thickening or thinning tends to grow further once initiated tectonic phases of strong compression correspond to the climax of such instabilities the response of models with cold lithospheric roots of various intensities has been computed both in two and three dimensions they yield velocity distributions and stress fields instructive comparisons are made with earthquake focal mechanisms and in situ stress measurements in the alpine and appalachian regions in the presence of lateral variations of the mechanical properties of the lithosphere the tectonic style is not only function of the local topography and of the nature of its compensation deformations in neighbouring provinces are coupled as shown by 3 dimensional models for example thickening sustained by a cold lithospheric root may generate extension in peripheral zones of weakness these last results illustrate the point that the computation of regional tectonic stresses requires the knowledge of the density anomalies within the lithosphere on the one hand and of geometrical constraints related to lateral mechanical heterogeneities on the other

late quaternary to decadal velocity fields in asia 1 fields of crustal velocity and strain rate within asia averaged over horizontal distances of 300 km are calculated from estimates of quaternary slip rates on faults and rates of change of lengths of baselines between gps sites the calculated fields fit strain rate components estimated from quaternary slip rates within 95 confidence limits and fit the gps data with an rms residual of 2 5 mm yr 1 the strain rate field constrained only by quaternary slip rates of faults differs from the combined solution at the level of 20 nstrain yr 1 and the differences are explicable by the limited accuracy of geological estimates of slip rate therefore the strain rate fields determined by geodetic measurements taken on the 10 year timescale are consistent with geological observations on the 104 year timescale the combined solution for velocities shows plate or block like motion within deforming asia only in tarim south china and the amurian region elsewhere the kinematics of rigid blocks does not provide a useful description of the motion in particular velocities relative to eurasia point at high angles to contours of their magnitudes not parallel to them as they would if large regions moved as rigid blocks the axes of principal horizontal compressional strain align with topographic gradients in agreement with previous suggestions that asia deforms as a continuum under the influence of gravity

source parameters for 11 earthquakes in the tien shan central asia determined by p and sh waveform inversion the tien shan mountain belt of central asia was formed in late paleozoic time and reactivated in tertiary time following the collision of india with eurasia more than 1500 km to the south to study the style and distribution of faulting occurring today in the tien shan we digitized long period world wide standard seismograph network p and sh waveforms of 11 of the largest tien shan earthquakes between 1965 and 1982 and then used a least squares inversion routine to constrain their fault plane solutions and depths four of the earthquakes occurred near the southern edge of the tien shan and two occurred in the intermontane fergana basin these earthquakes occurred at depths of 10 20 km and are associated with thrust faulting on east west to southwest northeast striking fault planes that dip moderately 35 55 the other five earthquakes in the northern tien shan were mostly deeper 15 44 km and the fault plane solutions of four of them are similar to the events farther south the exception an event on june 2 1973 probably occurred on a gently northward dipping east west striking fault plane all events occur within the basement at depths of 10 km or greater waveforms for earthquakes near the edges of the dzungarian fergana and kucha basins are fit best using velocity structures with thick sediment layers implying that sediments from these basins have been and are being underthrust beneath the neighboring mountains we examined only the relatively small mb 5 5 6 2 events after 1965 and did not study any of the largest events ms 8 that occurred in the tien shan earlier in this century still we can conclude that north south shortening is presently occurring in the tien shan with the formation of basement uplifts flanked by moderately dipping thrust faults the present day tectonics of the tien shan seem to be analogous to those of the rocky mountains in colorado wyoming and utah during the laramide orogeny in late cretaceous and early tertiary time

mantle dynamics uplift of the tibetan plateau and the indian monsoon convective removal of lower lithosphere beneath the tibetan plateau can account for a rapid increase in the mean elevation of the tibetan plateau of 1000 m or more in a few million years such uplift seems to be required by abrupt tectonic and environmental changes in asia and the indian ocean in late cenozoic time the composition of basaltic volcanism in northern tibet which apparently began at about 13 ma implies melting of lithosphere not asthenosphere the most plausible mechanism for rapid heat transfer to the midlithosphere is by convective removal of deeper lithosphere and its replacement by hotter asthenosphere the initiation of normal faulting in tibet at about 8 3 ma suggests that the plateau underwent an appreciable increase in elevation at that time an increase due solely to the isostatic response to crustal thickening caused by india s penetration into eurasia should have been slow and could not have triggered normal faulting another process such as removal of relatively cold dense lower lithosphere must have caused a supplemental uplift of the surface folding and faulting of the indo australian plate south of india the most prominent oceanic intraplate deformation on earth began between about 7 5 and 8 ma and indicates an increased north south compressional stress within the indo australian plate a tibetan uplift of only 1000 m if the result of removal of lower lithosphere should have increased the compressional stress that the plateau applies to india and that resists india s northward movement from an amount too small to fold oceanic lithosphere to one sufficient to do so the climate of the equatorial indian ocean and southern asia changed at about 6 9 ma monsoonal winds apparently strengthened northern pakistan became more arid but weathering of rock in the eastern himalaya apparently increased because of its high altitude and lateral extent the tibetan plateau provides a heat source at midlatitudes that should oppose classical symmetric hadley circulation between the equator and temperate latitudes and that should help to drive an essentially opposite circulation characteristic of summer monsoons for the simple case of axisymmetric heating no dependence on longitude of an atmosphere without dissipation theoretical analyses by hou lindzen and plumb show that an axisymmetric heat source displaced from the equator can drive a much stronger meridianal monsoonlike circulation than such a source centered on the equator but only if heating exceeds a threshold whose level increases with the latitude of the heat source because heating of the atmosphere over tibet should increase monotonically with elevation of the plateau a modest uplift 1000 2500 m of tibet already of substantial extent and height might have been sufficient to exceed a threshold necessary for a strong monsoon the virtual simultaneity of these phenomena suggests that uplift was rapid approximately 1000 m to 2500 m in a few million years moreover nearly simultaneously with the late miocene strengthening of the monsoon the calcite compensation depth in the oceans dropped plants using the relatively efficient c4 pathway for photosynthesis evolved rapidly and atmospheric co2 seems to have decreased suggesting causal relationships and positive feedbacks among these phenomena both a supplemental uplift of the himalaya the southern edge of tibet and a strengthened monsoon may have accelerated erosion and weathering of silicate rock in the himalaya that in turn enhanced extraction of co2 from the atmosphere thus these correlations offer some support for links between plateau uplift a downdrawing of co2 from the atmosphere and global climate change as proposed by raymo ruddiman and froehlich mantle dynamics beneath mountain belts not only may profoundly affect tectonic processes near and far from the belts but might also play an important role in altering regional and global climates

fault plane solutions of earthquakes and active tectonics of the tibetan plateau and its margins summary fault plane solutions of earthquakes within and on the margins of the tibetan plateau show diverse styles of faulting and deformation with thrust faulting and crustal shortening normal to the margins of the plateau and with normal and strike slip faulting resulting in roughly east west crustal extension within the plateau the direction of overthrusting of the himalaya onto the indian shield is radially outward varying from southwest in the western himalaya to south southeast in the east assuming that the indian shield behaves rigidly this requires a west northwest divergence of western tibet from southeastern tibet at a rate of 18 9 mm yr 1 comparable with the rate of convergence at the himalaya fault plane solutions of earthquakes in the southern portion of the tibet plateau consistently show large components of normal faulting on roughly north striking planes and corroborate such extension within the high plateau where elevations exceed 5000 m normal and strike slip faulting occur so that an overall east southeast west northwest extension of the region at about 10 mm yr 1 is partitioned into roughly equal parts of crustal thinning and north northeast south southwest crustal shortening about 5 mm yr 1 in general strike slip faulting characterizes solutions for earthquakes within eastern tibet where mean elevations drop below 4500 5000 m but the orientations of the strike slip faults vary across the region in central tibet left lateral slip occurs on planes trending roughly northeast but for earthquakes farther east the orientations of that plane become progressively east west and then southeast this variation in orientation implies a rotation of material along curved left lateral shear zones thus the eastward extrusion of tibet appears to be facilitated not only by rapid left lateral shear but also by large clockwise rotations of the material in eastern tibet the rate of eastward extrusion of material in eastern tibet relative to the tarim basin to its north is roughly 30 40 mm yr 1 fault plane solutions of earthquakes in the northern and eastern margins of tibet show large components of thrust faulting with the p axes oriented radially outward from the plateau and approximately perpendicular to the regional topographic contours of the plateau the orientation of this crustal shortening is northeast southwest on the northeastern margin east west on the eastern margin and northwest southeast in the longmenshan on the southeastern margin thus at least some of the extrusion of eastern tibet out of india s northward path into asia is absorbed by crustal shortening on the margins of the plateau the variation from normal faulting in the high tibetan plateau where elevations exceed 5000 m to dominantly strike slip faulting farther east where elevations are lower and then to thrust faulting on the margins of the plateau where elevations drop below 3000 m surely results at least in part from a decrease in the value of the vertical stress the magnitude of the east west compressive stress need not vary across the plateau

contemporary crustal deformation of northeast tibet from geodetic investigations and a comparison between the seismic and geodetic moment release rates abstract we processed the data recorded by 32 continuous gps stations from the crustal movement observation network of china and combined them with published gps observations in a common geodetic reference frame to better depict the ongoing crustal velocities and deformation field of northeast tibet the derived velocity field indicates that the deformation of northeast tibet is characterized mainly by crustal shortening coupled with the eastward extrusion of crustal material and clockwise block rotation we further used the gps velocity field to evaluate the strain rate tensors by the method of sandwell and wessel 2016 the strain rate field shows an obvious transition of the contraction rate from ne sw to nee sww and a decrease in the maximum shear strain rate from west to east was observed along the east kunlun fault and haiyuan fault indicating that the block rotation and lateral extrusion therein may be partitioned by ne sw oriented compression to some extent the inhomogeneity of the rotation strain rate implies that northeast tibet does not behave like a single rigid block but instead resembles a left lateral shear zone supporting the hypothesis that the northward transfer of motion along the east kunlun fault may play a significant role in the deformation style of this area the dilatation rate verified that most of northeast tibet experiences crustal thickening which coincides with the average uplift rate of 1 mm y revealed by gps velocities and channel steepness indices we also estimated the geodetic moment rates and compared them with the seismic moment rates estimated using the truncated gutenberg richter earthquake distribution the ratio of the seismic strain rate to the geodetic strain rate in the regions near the altyn tagh fault the qaidam basin adjoining the qilian mountains the haiyuan fault and the central part of the east kunlun fault are near unity indicating that the deformation in these areas is dominated by seismic release and may represent fully coupled seismogenic zones the occurrence of large earthquakes in the past estimated magnitudes m 6 coupled with a high geodetic deformation rate in these regions could indicate overdue m 6 earthquakes however other factors including the occurrence of creep and the long recurrence time of mw 8 0 earthquakes should be taken into consideration along the haiyuan fault by contrast in other regions including the southern and eastern parts of the haiyuan fault along the elashan fault and the northern qilian mountains the geodetic moment rate is higher than the seismic moment rate by a factor between 3 and 9 we propose that this phenomenon is related to the undersampling of long term earthquake rates within the seismic catalogue long term aseismic crustal deformation or a combination of these factors our study has implications for better comprehending the current seismotectonic pattern of northeast tibet

formation and evolution of strike slip faults rifts and basins during the india asia collision an experimental approach the processes which have governed the formation and evolution of large tertiary strike slip faults during the penetration of india into eastern asia are investigated by means of plane strain indentation experiments on layered plasticine models the steady state deformation of plasticine obeys a power creep flow law e c t n the stress exponent n is between 6 and 9 at 25 c uniaxial plane strain tests on cubic specimens show that the growth of faults in layered plasticine results from strain softening a process observed for strain rates ranging from 3 5 10 5 to 3 6 10 3 s 1 fault or shear zones form after only 7 10 bulk strain subsequent deformation is controlled by the geometry of the fault pattern rather than the physical properties of the plasticine a series of nine plane strain indentation experiments shows the influence of boundary conditions as well as that of the internal structure of the plasticine model on the faulting sequence the ubiquity of strain softening in experimental deformation of a variety of rocks as well as the widespread occurrence of shear zones in nature suggest that long term deformation of the continental lithosphere may also be primarily influenced by the geometry of large faults which rapidly develop with increasing strain the deformation and faulting sequence observed in the plasticine indentation experiments may thus be compared to collision induced strikeslip faulting in asia particularly to total offsets and rates of movements on the faults the experiments simulate the evolution of the western ends of the strike slip faults which have probably been analogous to trench fault fault triple junctions the experiments also illustrate mechanisms for the formation of extensional basins such as the south china sea north china basin and andaman sea near active continental margins the basins which appear to absorb terminal offsets along major strike slip faults near such margins may result from mismatch between the sharply angular shape of the deformed continental edge and the more regularly curved trench along which the smoothly flexed oceanic lithosphere subducts the existence of distinct phases of strike slip extrusion corroborates the idea that the discontinuities in time which typify intracontinental tectonics and orogenic cycles may often result from strain localization and the ensuing discontinuous non steady state deformation of the continental lithosphere

velocity field in asia inferred from quaternary fault slip rates and global positioning system observations we perform a joint inversion of quaternary strain rates and 238 global positioning system gps velocities in asia for a self consistent velocity field the reference frames for all geodetic velocity observations are determined in our inversion procedure india in moves relative to eurasia eu about a pole of rotation at 29 78 n 7 51 e 0 353 myr 1 which yields a velocity along the himalaya within india that is 73 76 of the magnitude of the in eu nuvel 1a velocity and a vector azimuth that is 8 10 clockwise of nuvel 1a in eu vector azimuth relative to eurasia south china moves at 9 11 mm yr in the direction 110 120 with a pole position 64 84 n 156 74 e 0 12 myr 1 amurian block motion has a pole position in a similar location but at a slower rate 64 61 n 158 23 e 0 077 myr 1 and most of the amurian eurasia motion is accommodated by extension across lake baikal tarim basin moves relative to eurasia about a pole of rotation at 39 24 n 98 2 e 0 539 myr 1 and 16 18 mm yr of shortening is accommodated across the west central tien shan there is distributed e w extension throughout both southern and north central tibet within southern tibet between the longitudes of 77 e to 92 e the deformation field accommodates 16 19 mm yr of e w extension we compare predicted seismic moment rates with those observed in this century in asia total observed seismic moment rates within the entire area of central and east asia 2 2 107 km2 in this century are 2 26 0 7 1020 n m yr 1 as compared with a predicted total rate of 2 03 0 066 1020 n m yr 1 comparisons between observed and predicted moment rates within 42 subregions reflect the generally unstable process of inferring long term seismic moment rates from a catalog of limited duration 94 years an observation period of 10 000 years would be required to reduce uncertainties in observed seismic moment rate to the same size as the uncertainties in model tectonic moment rates inferred from the joint inversion of gps and quaternary rates of strain we show that in general a better correlation with model tectonic moment rate is inferred from the seismicity catalog by considering the numbers of earthquakes above a cutoff magnitude mb 5 0 for the period january 1 1965 to january 1 1999

modelling waves and currents at the lstf and other laboratory facilities the main objective of the paper is to use the detailed measurements of wave and current motions in the large scale sediment transport facility lstf at the coastal and hydraulics laboratory chl erdc vicksburg to compare with model simulations and particularly to investigate the accuracy of two wave models used for driving nearshore circulation computations the major questions asked are how accurate are the models and how do the inevitable inaccuracies in the prediction of the wave quantities influence the current predictions the lstf is one of a few large scale facilities around the world designed to accurately reproduce the 2 d horizontal flow and sediment conditions in the nearshore region of a littoral coast the measurements include all relevant wave data such as wave heights setup variations 2 d current distributions across the basin and a detailed cross shore array of measurements of the 3 d current structures the two wave model formulations used for the comparisons are versions of the kinematic wave model one using sinusoidal the other using non sinusoidal wave shapes both before and after wave breaking as wave drivers in the quasi 3 d model shorecirc sc the comparisons with the measurements are aimed at providing an overall picture of both the flow conditions and the model performance the discussions provide insight into the mechanisms behind the complex nearshore dynamics and demonstrate the strengths and weaknesses of the two methods in a comprehensive comparison of available wave and current data the comparisons also reveal that an essential requirement for a wave driver is its capability to represent the phase motion of the wave and its capability to correctly predict the phase speed of waves

geodetic evidence for a low slip rate in the altyn tagh fault system the collision between india and asia has been simulated with a variety of computational models that describe or predict the motions of the main faults of east asia geological slip rate estimates of 20 30 mm yr 1 suggest that the largest of these faults the 2 000 km long altyn tagh fault system on the northern edge of the tibetan plateau absorbs as much of the indo asian convergence signal as do the himalayas partly by oblique slip and partly by contraction and mountain growth however the predictions of dynamic models for asian deformation and the lower bounds of some geological slip rates estimates 3 9 mm yr 1 refs 7 8 suggest that the altyn tagh system is less active here we report geodetic data from 89 91 degrees e that indicate left lateral shear of 9 5 mm yr 1 and contraction of 3 1 mm yr 1 across the altyn tagh system this result combined with our finding that at 90 degrees e tibet contracts north south at 9 1 mm yr 1 supports the predictions of dynamic models of asian deformation

how the continents deform the evidence from tectonic geodesy space geodesy now provides quantitative maps of the surface velocity field within tectonically active regions supplying constraints on the spatial distribution of deformation the forces that drive it and the brittle and ductile properties of continental lithosphere deformation is usefully described as relative motions among elastic blocks and is block like because major faults are weaker than adjacent intact crust despite similarities continental block kinematics differs from global plate tectonics blocks are much smaller typically 100 1000 km in size departures from block rigidity are sometimes measurable and blocks evolve over 1 10 ma timescales particularly near their often geometrically irregular boundaries quantitatively relating deformation to the forces that drive it requires simplifying assumptions about the strength distribution in the lithosphere if brittle elastic crust is strongest interactions among blocks control the deformation if ductile lithosphere is the stronger its fl

how flat is tibet high resolution digital topography three arc second grid for most of tibet provides new information to characterize the relief of the highest and largest plateau on earth the arid to semiarid central and northern part of the plateau interior has low relief average slopes of over 250 m windows and a mean elevation of 5023 m above sea level at moderate wavelengths of m relief is or less for most of tibet as opposed to the much higher relief of up to 6 km on the plateau edges where glacial and fluvial dissection is greater because of higher levels of precipitation the only faults manifesting significant topographic relief are the relatively small scale generally north trending graben systems primarily in southern tibet and several large scale fault systems near the edges of tibet the flatness of tibet implies that 1 there has been little deformation especially shortening of the uppermost crust north of the graben systems during the late cenozoic and 2 shallow crustal isostatic compensation has been acting to level the surface of the plateau

the amurian plate motion and current plate kinematics in eastern asia we use global positioning system gps velocity data to model eastern asian plate kinematics out of 15 stations in korea russia china and japan studied here three sites considered to be on the stable interior of the hypothetical amurian plate showed eastward velocities as fast as 9 10 mm yr with respect to the eurasian plate they were stationary relative to each other to within 1 mm yr and these velocity vectors together with those of a few additional sites were used to accurately determine the instantaneous angular velocity euler vector of the amurian plate the predicted movement between the amurian and the north american plates is consistent with slip vectors along the eastern margin of the japan sea and sakhalin which reduces the necessity to postulate the existence of the okhotsk plate the euler vector of the amurian plate predicts left lateral movement along its boundary with the south china block consistent with neotectonic estimates of the displacement at the qinling fault possibly the southern boundary of the amurian plate the amurian plate offers a platform for models of interseismic strain buildup in southwest japan by the philippine sea plate subduction at the nankai trough slip vectors along the baikal rift the boundary between the amurian and the eurasian plates are largely inconsistent with the gps based euler vector suggesting an intrinsic difficulty in using earthquake slip vectors in continental rift zones for such studies

present day crustal motion within the tibetan plateau inferred from gps measurements 1 using the measurements of 726 gps stations around the tibetan plateau we determine the rigid rotation of the entire plateau in a eurasia fixed reference frame which can be best described by an euler vector of 24 38 0 42 n 102 37 0 42 e 0 7096 0 0206 ma the rigid rotational component accommodates at least 50 of the northeastward thrust from india and dominates the eastward extrusion of the northern plateau after removing the rigid rotation to highlight the interior deformation within the plateau we find that the most remarkable interior deformation of the plateau is a glacier like flow zone which starts at somewhere between the middle and western plateau goes clockwise around the eastern himalayan syntaxis ehs and ends at the southeast corner of the plateau with a fan like front the deformation feature of the southern plateau especially the emergence of the flow zone could be attributed to an eastward escape of highly plastic upper crustal material driven by a lower crust viscous channel flow generated by lateral compression and gravitational buoyancy at the later developmental stage of the plateau the first order feature of crustal deformation of the northeastern plateau can be well explained by a three dimensional elastic half space dislocation model with rates of dislocation segments comparable to the ones from geological observations in the eastern plateau although gps data show no significant convergence between the eastern margin of the plateau and the sichuan basin a small but significant compressional strain rate component of 10 5 2 8 nstrain yr exists in a relatively narrow region around the eastern margin in addition a large part of the eastern plateau northeast of the ehs is not undergoing shortening along the northeastward convergence direction of the ehs but is stretching

tertiary strike slip faulting in southeastern mongolia and implications for asian tectonics abstract geologic maps have long portrayed the late cretaceous recent geologic history of southeastern mongolia as tectonically quiescent we present new data based on outcrop observations that indicate the northeast trending east gobi fault zone egfz was reactivated in the cenozoic as a sinistral strike slip fault system inversions of cenozoic fault slip data imply that faulting was associated with north northwest subhorizontal shortening and east northeast subhorizontal extension we propose that faulting is tertiary in age based on published interpretations of seismic reflection data which reveal that the mid cretaceous 100 95 ma unconformity is deformed by strike slip faults and based on field observation of strike slip faults and fracture sets that cut upper cretaceous and cenozoic strata but lack evidence for neotectonic activity published seismicity maps also appear to argue against significant quaternary faulting within the egfz these new data may lend credence to published models proposing a middle miocene or older kinematic linkage between the egfz and the altyn tagh fault in china the recognition that the egfz has a history of left lateral displacement in both the early mesozoic and cenozoic means that currently available estimates of offset based on displaced paleozoic rocks constrain total offset only this reactivation history supports the notion that inherited lithospheric structures are important in controlling the location and thus modes of intracontinental deformation in asia as a function of collisional far field effects and evolving boundary conditions of the pacific margin

active deformation of asia from kinematics to dynamics using estimated strain rates across the actively deforming region of asia to infer variations in viscous stress it is shown that gradients of stress point uphill toward the center of tibet where gravitational potential energy per unit surface area reaches a maximum thus the dynamics of the deformation seem to obey the equation of creeping flow which expresses a balance between gradients in stress and the gravitational body force this balance in the region of the tibetan plateau yields an estimate of 10 22 pascal second for the average viscosity of the tibetan lithosphere which is only about 10 to 100 times greater than the viscosity of the convecting upper mantle

relatively recent construction of the tien shan inferred from gps measurements of present day crustal deformation rates the tien shan a high seismically active intracontinental mountain belt 1 000 2 000 km north of the himalaya has grown as a result of india s collision with asia1 the crustal shortening 200 50 km refs 2 3 and thickening that gave rise to the tien shan accommodates only a small fraction of india s total penetration into asia 2 000 3 000km and the temporal relationship of deformation in this belt to the india asia collision remains unclear here we report geodetic measurements of the tien shan using the global positioning system gps that indicate that the current crustal shortening rate is nearly half of india s convergence rate with eurasia in this area4 we infer a total shortening rate for the tien shan of 20 mm yr 1 which is approximately twice that inferred previously from the extrapolation of slip rates in the holocene3 and earthquake induced displacements during this century5 suggesting that the rate of mountain building in this region has accelerated several fold since the onset of collision 50 55 myr ago6 7 if as we argue the current shortening rate can be extrapolated to geological timescales then our results suggest that most of the tien shan has been constructed during the past 10 myr perhaps in response to an increased horizontal force following an abrupt rise of the tibetan plateau8 9

implications of crustal property variations for models of tibetan plateau evolution shear coupled teleseismic p waves sampling the interior of the tibetan plateau provide evidence of systematic variations in crustal structure the crust thins by up to 20 km from south to north with a concomitant increase in poisson s ratio from normal values in the south to unusually high values in the north this suggests that the crust of the northern plateau is partially melted due to high temperatures these changes imply spatial and perhaps temporal variations in the way the elevation of the high plateau is created and maintained

quaternary extension in southern tibet field observations and tectonic implications we summarize evidence for quaternary and active faulting collected in the field during three sino french expeditions to southeastern tibet 1980 1982 detailed mapping of quaternary and active faults as well as microtectonic measurements indicate that normal faulting has been the dominant tectonic regime north of the himalayas in the last 2 0 5 m y the maximum horizontal principal stress in south tibet appears to be only the intermediate principal stress 2 1 being vertical south of the chord joining the eastern and western syntaxes of the himalayan arc extensional strains are principally localized within seven regularly spaced rift zones three of which have been studied in some detail the extension direction is determined to be n96 e 7 mainly from statistical averaging of strikes of newly formed normal faults throw rates on normal faults are evaluated for different time spans 2 0 5 m y 60 40 kyr and 10 2 kyr b p using structural and topographic reliefs as well as synglacial and postglacial vertical offsets the rate of quaternary extension is about 1 m y 1 along an 1100 km long ese traverse across south tibet this corresponds to a spreading rate of 1 0 6 cm yr this rate and the divergent horizontal projections of slip vectors of earthquakes along the himalayan front constrain the rate at which rigid india underthrusts southern tibet to be 2 1 cm yr although most of the normal faults appear to be independent of and nearly orthogonal to the e w mesozoic tertiary tectonic fabric the yadong gulu rift appears to be guided for over 130 km by the older oblique ne sw nyainqentanglha range and fault zones along it this reactivated zone is the most prominent left lateral strike slip fault system in se tibet excepting this zone and the vicinity of the se extremity of the karakorum fault quaternary strike slip faulting is rare in south tibet i e south of the chord between the syntaxes of the himalayan arc north of the chord the tectonic style is different there minor conjugate strike slip faulting is widespread and appears to control quaternary normal faulting which is more diffuse and subdued than in south tibet along the chord the presence of a major zone of active right lateral en echelon strike slip faults karakorum jiali fault zone probably reflects the greater facility of eastward extrusion in north central tibet in response to the northward push of india the eastern himalayan syntaxis may be an obstacle to such extrusion movements south of the chord

observing plate motions in s e asia geodetic results of the geodyssea project this paper presents the final geodetic results of the geodyssea project the gps data from a 42 station network observed during two field campaigns 1994 1996 were analyzed by four groups using different software packages and analysis strategies the precision of both campaign coordinate solutions was found to be 4 7 mm for the horizontal and 1 cm for the vertical component the campaign solutions were merged into one unique solution which was accurately mapped into the itrf 96 reference frame the global accuracy of this solution with respect to itrf 96 is 1 cm while the resolution of the relative horizontal velocities is estimated to be at the level of 2 3 mm yr this solution was used as the basis for all scientific interpretations which are published in separate papers the velocity estimates of a part of the network provided the first direct measurement of a relative motion of the sundaland block with respect to eurasian plate

global positioning system measurements from eastern tibet and their implications for india eurasia intercontinental deformation we present and interpret global positioning system gps measurements of crustal motions for the period 1991 1998 for a network encompassing the eastern part of the tibetan plateau and its foreland relative to a eurasian frame defined by minimizing the velocities of 16 gps stations in europe central asia and siberia stations within all parts of the plateau foreland in south china move 6 10 mm yr east southeast indicating that the eastward movement within the plateau is part of a broader eastward movement that involves the plateau and its eastern and northern foreland north of the plateau foreland stations move northeastward at 10 mm yr indicating that the northern boundary of the deformation zone lies north of the plateau with this realization of a eurasian frame the velocity of the gps station at bangalore in southern india implies that the northward motion of india is 5 12 mm yr slower than that predicted from the nuvel 1a plate reconstruction viewed relative to the south china block stations of the northeast plateau bounded on the north by the qilian shan and the altyn tagh fault move nne to ne with velocities ranging from 19 mm yr within the plateau to 5 11 mm yr in its foreland the altyn tagh fault shows left lateral slip of 10 mm yr at 95 e and shortening across the fault of 5 mm yr stations south and west of the xianshuihe xiaojiang fault system define a crustal fragment rotating clockwise at 10 mm yr relative to the south china block around the eastern himalayan syntaxis the gps measurements indicate no significant shortening 3 mm yr within the longmen shan of the central eastern plateau and its adjacent foreland although the longmen shan rise over 6 km in 100 km horizontal distance geological studies indicate that the deformational field was established diachronously in late miocene to pliocene time was characterized by no east west shortening of tibetan crust and has an inhomogeneous style of deformation resulting from a balance of different tectonic processes

numerical modeling of intraplate deformation simple mechanical models of continental collision we study a model of continental collision in which one of the continents acts as a rigid die indenting the other plate which flows as an incompressible viscoplastic medium we consider two extreme cases of plane deformation 1 plane strain which corresponds to an infinitely thick lithospheric plate and 2 plane stress corresponding to a very thin plate deformation of the lithosphere a thick plate should be intermediate between those extremes we found that the flow in the plane strain case is quite similar to that obtained by slip line theory the plane stress results are quite different since in this case most of the plate shortening is taken up by the thickening of the lithosphere we also explored the role of boundary conditions on the flow in particular the role of the side walls containing the flow of the lithosphere in the case of a free lateral boundary the main feature is a flow of matter toward this free wall and a s like pattern for the horizontal stress field for a rigid wall on the other hand the plane strain and the plane stress results are quite different in the first case there is a large return flow on the sides of the punch the material being extruded along the only free surface available in the plane stress case the return flow disappears and the material displaced by the penetration of the die tends to thicken the plate the role of a nonlinear constitutive relation is studied for power law creep as the power of the flow limit increases the flow retains its general features but the deformation localizes creating sharper contrasts between high and low strain rate areas in plane stress the effect of nonlinearity is to increase the contrasts in vertical motion available data for asia are discussed in the light of the new results

a thin viscous sheet model for continental deformation summary for the purposes of describing its large scale and long term deformation the continental lithosphere is regarded as a continuum obeying a newtonian or a power law rheology the flow of a thin sheet of power law material overlying an inviscid substrate is studied under the assumption that vertical gradients of the horizontal velocity are negligible a numerical model is used to investigate the deformation of such a sheet under conditions approximating those of continent continent collision the material flows in response to forces applied to its boundaries for example the indenting of one continent by another and to forces in its interior arising from gradients in crustal thickness the horizontal divergence of the flow produces changes in the crustal thickness and hence a time dependent form to the flow itself for a given set of boundary conditions the flow depends on the stress exponent in the power law rheology n and on the argand number which is a measure of the ratio between the stress arising from crustal thickness contrasts and the stress required to deform the material at the ambient strain rates when the effective viscosity of the medium is very high ar o crustal thickness variations do not influence the flow if the material is newtonian n l the deformation associated with an influx of material approximating an indenter is of much greater lateral dimension than the width of the indenter whereas when material has a power law rheology n 3 5 are used the deformation is confined to a region of lateral extent comparable to that of the indenter as ar increases the forces arising from crustal thickness contrasts exert more influence on the flow and the maximum crustal thickness that can be sustained by a given influx of material is related by a simple expression to the effective viscosity of the medium at the ambient strain rates in the limit of a very weak medium ar 10 the lithosphere is unable to sustain appreciable crustal elevation contrasts the results of these numerical experiments show that systems in which the effective viscosities are such that the maximum deviatoric stresses are between 1 kbar

tomographic evidence for localized lithospheric shear along the altyn tagh fault seismic tomography across the altyn tagh fault at the north edge of the tibetan plateau reveals a low p wave velocity anomaly below the fault down to 140 kilometers this anomaly probably reflects strike slip shear in the lithosphere slip partitioning may also induce a wedge of crust from the tarim basin to plunge into the mantle

parallel thrust and normal faulting in peru and constraints on the state of stress abstract we present evidence for recent normal faulting in the cordillera blanca of peru along planes trending approximately n40 w parallel to the chain this type of faulting suggests that there is an important component of extension perpendicular to the chain and contrasts markedly with the underthrusting of the nazca plate on the west side of peru and with the folding and thrusting observed in the sub andes on the eastern flank of the chain these contrasting styles of deformation and associated stress fields can be explained by buoyancy forces due to the high mountains and crustal roots in so far as this explanation is correct the difference in elevation between zones of normal and thrust faulting can be used to place an upper bound of 500 bars on the average stress difference in the crust

gps observation of crustal deformation in the taiwan luzon region fifteen gps stations in the batan luzon area were surveyed three times during 1996 to 1998 in order to study crustal deformation in the taiwan luzon region the data from these annual surveys are processed together with data from 20 permanent stations in the taiwan area and igs stations in the asia pacific region we found that the chinese continental margin is moving in directions of 102 112 with rates of 11 12 mm yr relative to stable eurasia there is a significant southward increase in crustal motion with respect to eurasia along the n s trend of the luzon arc in the coastal range of taiwan the velocity is 53 mm yr in the direction 313 further south at two volcanic islands lanhsu and batan the velocities increase to 68 mm yr in the direction 306 and 81 mm yr in the direction 299 respectively in northern luzon the velocity reaches 86 mm yr in the direction 301 remarkable fault slips of 17 31 mm yr across the philippine fault system are detected

finite strain calculations of continental deformation 2 comparison with the india asia collision zone numerical experiments on a thin viscous sheet model for deformation of continental lithosphere subjected to an indenting boundary condition yield distributions of crustal thickness of stress and strain rate and of latitudinal displacements that may be compared with observations in the india asia collision zone a simple indenting boundary condition applied to initially laterally homogeneous sheets obeying a power law rheology produces results that are in broad agreement with the observations provided that the power law exponent is three or greater and the sheet can support vertically integrated stress differences of 2 1013 5 1012 n m 1 in regions in front of the indenter under these conditions the calculated deformation shows accommodation of convergence primarily by crustal thickening to produce a plateau in front of the indenter palaeomagnetic data from india and tibet and the observed distribution of topography suggest that much of the post eocene convergence of india with asia has been taken up by deformation within asia that involved crustal thickening the principal difference between calculation and observation is the absence from the calculated strain rate fields of east west extension of the plateau in front of the indenting boundary the calculations show that once such a plateau is formed the buoyancy force associated with the crustal thickness contrast inhibits further thickening and the plateau strains at less than half the rate of its immediate surroundings seismically determined regional strain rates exhibit a similar distribution with the tibetan plateau straining at about one quarter the rate of the tien shan and ningxia gansu regions calculated principal compressive stress orientations and regional strain rates agree with the seismically determined quantities in the mongolia baikal tien shan tibet and ningxia gansu regions of asia to within the uncertainty of the latter the vertically integrated stresses that are calculated for the viscous sheet are comparable with those that can be supported by a theologically stratified continental lithosphere obeying laboratory determined flow laws we suggest that the thin viscous sheet model described in this paper and its companion gives a simple and physically plausible description of the observed deformation in central asia in this description the predominant mechanism of accommodation of continental convergence is diffuse crustal thickening with shear on vertical planes playing a subsidiary role once large crustal thickness contrasts have been established

earthquake strain rates and instantaneous relative motions within central and eastern asia summary we use the spatial distribution of moment tensors of earthquakes in this century to estimate the velocity field in asia within a eurasian reference frame in a least squares inversion strain rates on the surface of the earth are matched with continuous spline functions in order to recover the velocity gradient tensor associated with the seismic moment release in asia earthquakes account for 40 60 per cent of the expected motion of india relative to eurasia with the missing component of strain rate equivalent to about 20 mm yr 1 of n s shortening between siberia and india in this solution south china rotates counterclockwise and moves eastwards relative to siberia using rigid plate constraints we next investigate the characteristics of the complete horizontal strain field in asia that accommodates plate motions our strain rate solutions are analogous to the response of a newtonian thin viscous sheet in which the rate of work done by the straining medium in accommodating the velocity boundary conditions is a minimum in these solutions the euler pole for india relative to eurasia is constrained nuvel 1a demets et al 1994 but in the process of fitting the vlbi velocity at shanghai china ward 1994 heki et al 1995 the euler pole for south china is determined in the inversion a solution that both fits the velocity at shanghai china and yields a strain rate field consistent with the earthquake mechanisms is one where the south china block has a motion relative to siberia described by the pole at 51 n 131 e 0 3 deg myr 1 comparison of the complete strain field that accommodates plate motion with the seismic strains indicates that earthquake moment release rates in this century within mongolia are about a factor of 4 larger than the long term rate within gansu ningxia the earthquake moment rates have been about a factor of 2 higher than the long term rate the strike slip faulting within mongolia gansu ningxia western sichuan and yunnan is possibly a direct result of velocity boundary conditions imposed on the south china block by forces unrelated to continental collision of india and eurasia such as forces associated with subduction along the margins of south eastern asia verification of this requires a better understanding of the role of pre existing zones of weakness within the asian continental lithosphere

evidence from gravity and topography data for folding of tibet bouguer gravity and topography data from tibet suggest that the surface of the plateau and the subsurface density interfaces are warped into two series of ridges and troughs trending parallel to the collision zone with wavelengths of 150 and 500 km these folds are superimposed on an overall state of isostatic compensation owing to crustal thickening such folding is predicted from models of compressional shortening of a rheologically layered plate the results thus suggest the presence of a weak decoupling zone between the tibetan crust and upper mantle

dust fall and biological soil crust distribution as indicators of the aeolian environment in china s shapotou railway protective system abstract spatial and temporal distributions of dust fall along with the spatial variation of the biological crust were measured along a cross section in china s shapotou railway protective system the variation of the aeolian environment was evident from the spatial variation of near surface wind flow and both grain size and accumulation rate of dust fall these varied greatly along the selected cross section that extended from the mobile dune through the straw checkerboards and vegetation zone to the gravel platform of the railway on the temporal scale the dust fall concentrated from march through may the grain size composition changed monthly these temporal changes were related to the wind force variation and surface condition the surface conditions along the cross section were characterized by crust type crust thickness and ground 0 5 cm compaction crust could only be found in a relatively stable environment the distribution of dust fall and crust indicated the intensity of aeolian activity and the stability of the aeolian environment further improvement of the aeolian environment by natural processes or by human construction would promote additional dust fall and accelerate the development of the crust

a no lid zone in the central chang thang platform of tibet evidence from pure path phase velocity measurements of long period rayleigh waves pure path phase velocities of fundamental mode rayleigh waves in the period range 30 80s are obtained for different paths across the tibet plateau using a two event method with events located in tibet and aligned with world wide standard seismograph network stations at teleseismic distance evidence is found for significant lateral variations of structure of the crust and uppermost mantle in tibet in particular phase velocities obtained on paths contained in the chang thang platform are higher than average in tibet at periods shorter than 60s while remaining very low at longer periods when interpreting these phase velocities in terms of simple layered s velocity models it is found that contrary to other areas in tibet they are incompatible with the existence of a high velocity lid s velocities of 4 4 4 5 km s are found immediately beneath the moho and the crust is somewhat thinner 50 60 km than average in tibet this indicates a possible upwelling of hot asthenospheric material compatible with landsat evidence for subrecent basaltic volcanism and independent seismic evidence such as inefficient transmission of sn waves on the other hand a 70 km thick crust and very high upper mantle s velocities are found on the western border of tibet karakorum hindu kush compatible with previous finer scale results concerning the associated zone of intermediate depth seismicity

upper mantle velocity structure beneath the tibetan plateau from pn travel time tomography we inverted 1510 p arrival times from regional distances 333 1600 km in and around the tibetan plateau to map the lateral velocity variation within the uppermost mantle previous studies have placed first order constraints on upper mantle velocities but relied on data recorded almost exclusively at stations outside of the plateau we improve resolution by using 40 events recorded at stations within the tibetan plateau we combine these data with observations obtained from the international seismological centre isc to extend our coverage by including pn arrivals from 85 additional plateau events relocated in previous studies and recorded at stations in and around the tibetan plateau we use synthetic travel time data to evaluate the resolution of our data set the observations provide good resolution to about 1 over most of the plateau and surrounding regions our results show average pn velocities that are about 3 lower in the northern plateau relative to the southern plateau these variations correlate well with major tectonic features and previous geophysical observations in the qiangtang terrane of the northern plateau an area known to be inefficient for sn propagation pn is slow relative to both the plateau south of the banggong nujiang suture and the tectonically stable tarim basin north of the plateau this is strong evidence for the existence of partial melt within the uppermost mantle beneath the northern tibetan plateau however when laboratory estimates of relationships between temperature velocity and attenuation are applied a relatively small temperature variation 240 to 370 c is required to explain our pn velocity observations when combined with geochemical constraints from volcanics in the northern plateau our results strongly suggest that the mantle lid is intact beneath the northern plateau this result would preclude tectonic models involving wholesale delamination of the mantle lithosphere in the northern tibetan plateau

a stress interpretation scheme applied to lunar gravity and topography data 1 we present an approach of the inverse gravimetric problem that allows the gravity to be directly related to the deviatoric stresses without any rheological assumptions in this approach a new set of parameters is considered 1 the density variations over equipotential surfaces and the height of interfaces above the corresponding equipotential surfaces and 2 the stress difference the method is applied to lunar topographic and gravimetric data that are interpreted in term of transversally isotropic deviatoric stress within the moon it also provides inference on density and crustal thickness variations the estimated lateral variation in deviatoric stress is about 500 bars within the crust and upper mantle in the crust because of topography the strongest stress differences take place on the far side with large lateral compressions beneath the south pole aitken basin vertical compression under the mascons of the nearside is the main feature within the upper mantle

the active tectonics of the eastern himalayan syntaxis and surrounding regions source parameters of 53 moderate sized earthquakes obtained from the joint inversion of regional and teleseismic distance long period body waves provide the data set for an analysis of the style of deformation and kinematics in the region of the eastern himalayan syntaxis focal mechanisms of eastern himalayan events show oblique thrust consistent with the n ne directed movement of the indian plate as it underthrusts a boundary that strikes at an oblique angle to the direction of convergence earthquakes near the sagaing fault show strike slip mechanisms with right lateral slip earthquakes on its northern splays however indicate predominant thrusting evidence that the dextral motion on the sagaing fault which accommodates a portion of the lateral motion between india and southeast asia terminates in a zone of thrust faulting at the eastern himalayan syntaxis remaining motion between india and southeast asia is accommodated in a zone of distributed shear in east burma and yunnan manifested by strike slip and oblique normal faulting east west extension crustal thinning and clockwise rotation of crustal blocks we determined strain rates throughout the region with a moment tensor summation using 25 years modern and 85 years modern and historic of earthquake data we matched the observed strains with a fifth order polynomial function and from this we determined both the velocity field and rotations with respect to a specified region velocities calculated relative to south china stationary show that the entire area extending from 20 n 36 n within deforming asia yunnan western sichuan and east tibet constitutes a distributed dextral shear zone with clockwise rotations up to 1 7 m y maximum in the region of the eastern syntaxis proper integrated strains across this zone relative to south china stationary show 38 mm yr 12mm yr of north directed motion at the himalaya remaining plate motion relative to south china fixed must be taken up by the underthrusting of india beneath the lesser himalaya strike slip motion on the sagaing fault and intraplate ne directed shortening within ne india as well as ne directed shortening within the eastern syntaxis proper 10 mm yr 2 mm yr of relative right lateral motion between india and southeast asia is absorbed in the region between the sagaing and red river faults 94 e 100 e it is the clockwise vorticity relative to south china associated with the deformation in yunnan east tibet and western sichuan that provides the relative north directed motion of 38 12 mm yr at the himalaya not all of the deformation is accommodated in right lateral shear between india and south china and between east tibet and south china velocity gradients exist that are parallel to the trend of the shear zone relative to a point within western sichuan 32 n 100 e the velocity field shows that the yunnan crust is moving s se at rates of 8 10 mm yr relative to south china there is no eastward expulsion of crustal material beyond the eastern margin of the tibetan plateau

present day kinematics of asia derived from geologic fault rates active fault geometry available quaternary rates on major faults and the far field plate motion are used to determine crustal kinematics in the collision zone between india and asia using a finite element formalism to construct a spherical shell model of the asian continent with embedded faults we seek a velocity solution approaching the motion of rigid blocks by minimizing the elastic strain energy in the fault bounded blocks by doing so we test the assumption that long term deformation within continents is mostly localized into major faults in the solution fault motion accounts for more than 80 of the deformation allowing us to describe our velocity model in terms of quasi rigid block rotations on the sphere south china is rotating clockwise about a pole located southwest of borneo implying an e ese velocity vector of 11 mm yr for a point at shanghai in agreement with the velocity vector determined by very long baseline interferometry heki et al 1995 the eastward movement of south china is accommodated by oblique extension along the red river fault at a rate of 10 5 mm yr in the south and by the combination of left lateral strike slip motion on the qinling dabie shan fault and the counterclockwise rotation of the ordos and adjacent blocks in the north the tarim rotates clockwise with respect to dzungaria about a pole located at 44 n 92 e consistent with increasing crustal shortening toward the west throughout the tien shan assuming incompressibility a crustal volume budget over a domain encompassing the tertiary mountain ranges in asia shows that over the last 10 000 years 73 4 of the north south shortening between india and asia has been absorbed by thickening of the lithosphere and 27 has been accommodated by lateral extrusion of continental blocks the present day predominance of thickening in asia results from the relatively slow eastward motion of south china controlled by strike slip faulting in the qinling shan and oblique extension in northeastern china

the role of gravitational potential energy in active deformation in the southwestern united states in orogenic belts intrinsic buoyancy and extrinsic boundary forces act on the continental lithosphere determining the relative importance of these forces is fundamental to understanding the dynamics of continental deformation in the western united states estimates of gravitational potential energy show that given available constraints on rheology buoyancy forces are sufficient to produce the wide variety of rates and styles of observed deformation including absence of deformation this underscores the importance of gravitational potential energy for deformation in orogenic belts

dynamics of the pacific north american plate boundary in the western united states the vertically averaged deviatoric stress tensor field within the western united states was determined with topographic data geoid data recent global positioning system observations and strain rate magnitudes and styles from quaternary faults gravitational potential energy differences control the large fault normal compression on the california coast deformation in the basin and range is driven in part by gravitational potential energy differences but extension directions there are modified by plate interaction stresses the california shear zone has relatively low vertically averaged viscosity of about 1021 pascal seconds whereas the basin and range has a higher vertically averaged viscosity of 1022pascal seconds

lithospheric structure of the qiangtang terrane northern tibetan plateau from complete regional waveform modeling evidence for partial melt we report models of p and s wave velocity and attenuation for the the crust and uppermost mantle of the qiangtang terrane northern tibetan plateau inferred by fitting reflectivity synthetic seismograms to observed complete regional waveforms the data are three component broadband seismograms recorded by the 1991 1992 iris passcal incorporated research institutions for seismology program for array seismic studies of the continental lithosphere tibetan plateau experiment and global seismic network stations in the region the qiangtang terrane has thick crust 65 5 km with p and s wave velocities of 6 1 6 3 and 3 34 3 43 km s respectively yielding an anomalously high crustal poisson s ratio of 0 29 0 02 seismic velocities of the upper mantle of the qiangtang terrane are normal for p waves and slow for s waves 8 10 and 4 35 4 41 km s respectively with a high mantle poisson s ratio of 0 29 0 01 attenuation in the crust and upper mantle is high qp 100 200 and qs 44 89 modeling of the broadband p waveforms suggests that a decrease in mantle velocity occurs at about 160 km depth in the mantle however this is not unambiguously supported by the data and modeling the crust and uppermost mantle of the qiangtang terrane probably contains partial melt based on the high poisson s ratio low shear wave velocities and low q the absence of high frequency sn and the presence of volcanism of mantle lithospheric origin support the presence of partial melt crustal and uppermost mantle structure in the qiangtang terrane is different from that for the lhasa terrane immediately to the south based on results of our previous studies the average crustal p and s wave velocities are 4 faster and 2 slower respectively in the qiangtang terrane relative to the lhasa terrane this yields a significant difference in the crustal poisson s ratio with values of 0 29 for the qiangtang terrane and 0 25 for the lhasa terrane differences in the uppermost mantle p and s wave velocities and poisson s ratios of these two adjacent terranes cannot be explained by temperature differences alone using the mantle temperature estimates of mcnamara et al 1997 we suggest that partial melt of an ultramafic composition beneath the qiangtang terrane fits the velocity and poisson s ratio estimates

preparation and deformation of synthetic aggregates of quartz synthetic quartz aggregates of low porosity have been fabricated from natural quartz powder impure silica gel and high purity silicic acid by hydrothermal isostatic pressing at 300 mpa the resulting specimens have water contents comparable to those of natural quartzites and have been used for deformation tests the flow stress at 1200 1300 k was found to vary in an inverse relationship to the water pressure estimated under the assumption that the measured water content homogeneously filled the measured porosity it is questionable that the natural quartz origin specimens were fully equilibrated in respect to the activity of the water but equilibration has probably been achieved in the amorphous silica origin specimens which crystallized under the experimental conditions in the deformation of the latter materials the experimental activation energy at water pressures approaching the confining pressure was found to be about 150 kj mol 1 for both however the impure gel origin material gave a stress exponent of about 2 3 whereas the high purity silicic acid origin material gave a stress exponent of about 4 in spite of the grain sizes being about 90 m and 20 m respectively it is concluded that while the impure specimens have higher intragranular strengths there is also a significant contribution of grain boundary processes to the strain in them which is absent in the high purity specimens

on the gravitational potential of the earth s lithosphere the mean potential energy of the lithosphere is useful for defining the tectonic reference state trs of the earth and can be used to constrain the ambient state of stress in the plates in the absence of external forces applied at the base or along plate boundaries a lithospheric column with the potential energy of the trs would remain undeformed thus the difference between the potential energy of a lithospheric column and the trs determines whether the column is in an extensional joeutral or compressional state of stress we evaluate and intraplate variations about this mean using a simple first order lithospheric density model this model assumed that the continental geotherm is linear and density variations below a depth of 125 km have negligible influence on and is consistent with observed geoid anomalies across continental margins is estimated to be 2 379 1014 n m 1 which is equivalent to the potential energy of both near sea level continental lithosphere 160 to 220 m for an assumed crustal density c in the range 2800 2700 kg m 3 and cooling oceanic lithosphere at a depth of 4 3 km with the exception of eurasia which has anomalously high mean potential energy 2 383 1014 n m 1 the mean potential energies of the continental plates are nearly identical to the global mean the mean potential of the oceanic plates was found to be a strong function of the mean age of the oceanic lithosphere both the global and plate mean potential energies are relatively insensitive to a wide range in c the potential of the mid ocean ridges 2 391 1014 n m 1 is greater than the global mean which is consistent with the divergent nature of the ridges elevated continental lithosphere with a height of about 70 m has an equivalent potential energy to suggesting that in the absence of external forces continental regions will be in a slightly extensional state of stress the importance of our potential energy formulation is substantiated by the strong correlation between the torque poles associated with the potential energy distributions and the observed plate velocity poles for the south american nazca indo australian and pacific plates

crustal thickening versus lateral expulsion in the indian asian continental collision since the beginning of the continental collision between india and asia there has been about 2500 km of convergence and the northward movement of india has been accommodated by major internal deformation of the asian lithosphere the crustal thickening in and around the tibetan plateau is clearly a direct consequence of this collision but there is considerable debate as to whether a large fraction of the indentation has been accommodated by eastward motion of the lithospheric blocks of southeastern asia and southern china numerical experiments described here test this hypothesis for a range of indentation geometries and rheological models of the lithosphere we employ a thin viscous sheet model of the lithosphere with a depth averaged nonlinear viscous rheology described by a stress strain rate exponent n and including gravitational buoyancy forces scaled by the dimensionless argand number ar the eastern boundary for the collision region is described as a lithostatic boundary the precollision normal stress is determined by static balance and that constant stress is applied throughout the collision the experiments show that during collision the eastern boundary is smoothly displaced to the east at a rate about 1 4 of the indentation rate with only minor variation due to geometry or rheology crustal thickening rates in the region of the plateau are reduced by between 10 and 25 of the corresponding rates determined for similar experiments with a rigid eastern boundary to generalize from the results of these experiments the total north south shortening strain produced by the collision is partitioned between crustal thickening and eastward displacement in the ratio of at least 3 1 and more probably 4 1

the field of crustal velocity in asia calculated from quaternary rates of slip on faults rates of quaternary faulting when treated as measures of strain rates within regions several times larger than the depths of brittle faulting yield a velocity field for asia that matches the average relative velocity between the indian and eurasian plates over the past 2 myr nuvel 1a and the velocity of shanghai with respect to eurasia measured with very long baseline interferometry vlbi by using rates of slip on discrete surfaces faults to define average strain rates we explicitly assume that the large scale deformation can be treated as that of a continuum depending on the jcales over which strain rates are averaged 87 to 93 per cent of the assumed strain ites are fit within assigned 1 sigma uncertainties and 97 to 99 per cent are fit within 2 sigma the internal consistency of these data and the agreement with independent rates nuvel 1a and vlbi test the validity of treating large scale deformation of continents as that of a continuum in the calculated velocity field a large fraction 85 per cent of india s convergence with eurasia is absorbed by crustal or lithospheric thickening among consistently misfit assumed strain rates are those dominated by high rates of strike slip faulting strain associated with 20 30 mm a 1 of slip on the altyn tagh karakax and karakorum faults is inconsistent with the distributions of strain in surrounding regions we conclude that these rates all of which are based on correlations of ages of offset landforms with changes in global climate and not on radiometric dating have been systematically overestimated lateral transport out of india s northward path is relatively minor south china is calculated to move east southeast at less than 10 mm a 1

faulting associated with large earthquakes and the average rate of deformation in central and eastern asia we determine approximate average rates of deformation in portions of eastern and central asia for the last 80 years these rates are based on estimates of source dimensions average displacements seismic moments and orientations of faults of all large earthquakes m 7 5 and as many other major earthquakes m 7 0 as possible in asia since 1900 a major part of this study is the compilation of these data for separate parts of asia we combine seismic moment tensors to estimate the average rotational strain or deformation tensors for these regions for the last 80 years we estimate that shortening across the tien shan has occurred at 11 mm yr between 5 and 19 mm yr left lateral shear on easterly trending planes in mongolia in northwestern tibet gansu and ningxia and on the xianshuihe fault in western sichuan has been very rapid more than 10 mm yr in each region and apparently as rapid as 40 mm yr in northwestern tibet gansu and ningxia slower rates of a few to 10 mm yr appear to characterize east west extension and conjugate north south shortening of tibet the earthquake history of the himalaya is consistent with rates of underthrusting of more than 10 mm yr but is probably less than the convergence rate of india with eurasia of about 50 mm yr the overall strain field is consistent with a large part of india s penetration into eurasia being absorbed by the extrusion of material out of india s way as a result southeast china moved at about 21 mm yr east southeast with respect to eurasia in the last 80 years with an uncertainty of about a factor of 2 in this estimated rate

crustal deformation in the baikal rift from gps measurements three years and four campaigns of global positioning system gps measurements 1994 1997 in the baikal rift zone largest active continental rift system in eurasia show crustal extension at a rate of 4 5 1 2 mm yr in a wnw ese direction a comparison with moment release of large historical earthquakes suggests that elastic strain is currently accumulating in the baikal rift zone along active faults that currently have the potential for a m 7 5 earthquake the gps derived extension rate in the baikal rift zone is at least two times greater than the prediction of most deformation models of asia this result could reflect the dynamic contribution of the pacific eurasia subduction to intracontinental deformation in asia in addition to the effect of the india eurasia collision

partitioning of localized and diffuse deformation in the tibetan plateau from joint inversions of geologic and geodetic observations abstract the spatial complexity of continental deformation in the greater tibetan plateau region can be defined as the extent to which relative motion of the indian and asian plates is partitioned between localized slip on major faults and distributed deformation potency rates provide a quantitative metric for determining the magnitudes of diffuse and on fault crustal deformation which are proportional to strain rates within crustal micro plates and fault slip rates respectively we simultaneously estimate micro plate rotation rates interseismic elastic strain accumulation fault slip rates on major structures and strain rates within 24 tectonic micro plates inferred from active fault maps in the greater tibetan plateau region using quasi static block models constrained by interseismic surface velocities at 731 gps sites and 9 holocene late quaternary geologic fault slip rates the joint geodetic geologic inversion indicates that geologic slip rates are kinematically consistent with differential micro plate motions estimated left lateral slip rates on the altyn tagh west central kunlun and xianshuihe faults are relatively homogeneous along strike 11 5 10 5 and 11 mm yr respectively while segmentation of the eastern kunlun fault by the intersecting elashan and riyueshan faults results in a decreased slip rate consistent with geologic observations the fraction of total potency rate associated with intrablock strain uncorrected for observational noise ranges from 0 27 in the himalayan range block to 0 91 in the aksai chin block monte carlo simulations are used to quantify the likelihood that internal deformation is statistically distinguishable from the uncertainties in geodetic velocities these simulations indicate that internal block deformation is statistically significant only within the himalayan range front where internal deformation accounts for id 0 11 of block potency rate budget west central plateau id 0 68 ganzi yushu xianshuihe 0 58 burma 0 09 jiali 0 38 and aksai chin 0 68 blocks in the other 18 tectonic micro plates within the plateau region estimated internal block potency is not currently distinguishable from the expected contribution of observational noise to residual velocities of the total potency budget within the tibetan plateau 85 is taken up by slip on major faults with the remaining 15 accommodated by internal processes at sub block scale distinguishable from observational noise the localization of the majority of plate boundary activity is also supported by the spatial distribution of modern and historical crustal earthquakes sixty seven percent of the total moment released by earthquakes in the cmt catalog and 90 of historical moment since 1900 have been released within 25 km of the major faults included in the block model representing only 12 of the characteristic half block length scale of 215 km the localization of deformation inferred from geologic geodetic and seismic observations suggests that forces applied to tectonic micro plates drive fault system activity at the india asia collision zone over decadal to quaternary time scales

continental tectonics in the aftermath of plate tectonics the success of plate tectonics required an acceptance of continental drift and thus a reinterpretation of the large scale geological history of most of the earth but the basic tenet of plate tectonics rigid body movements of large plates of lithosphere fails to apply to continental interiors where buoyant continental crust can detach from the underlying mantle to form mountain ranges and broad zones of diffuse tectonic activity

seismic anisotropy beneath the continents probing the depths of geology seismic anisotropy beneath continents is analyzed from shear wave splitting recorded at more than 300 continental seismic stations anisotropy is found to be a ubiquitous property that is due to mantle deformation from past and present orogenic activity the observed coherence with crustal deformation implies that the mantle plays a major if not dominant role in orogenies no evidence is found for a continental asthenospheric decoupling zone suggesting that continents are coupled to general mantle circulation

accommodating large scale intracontinental extension and compression in a single stress field a key example from the baikal rift system the baikal rift system in southern siberia is one of the main intracontinental extensional features on earth the rift system represents the northwestern boundary of the amuria plate and in that respect can be considered as an evolving plate boundary the baikal rift system has been widely studied both in terms of geology and geophysics and many models have been proposed for its formation and evolution however the age of the initiation of deformation and the mechanism driving this deformation are still largely debated while major extension has occurred since the late miocene pliocene the onset of extension seems older than the india asia collision implying that several driving mechanisms may have acted together or in relay through time in this work we review the available data and models for deformation in an area encompassing the baikal rift system the sayan ranges to the west and the transbaikal to the east using a synthesis of this data and our own field and mapping observations we show that the baikal rift system along with transpressional deformation in the sayan ranges and transtension in the transbaikal area can be explained through major left lateral strike slip systems the deformation is strongly controlled by inherited crustal and lithospheric structures and is distributed over a wide area within the western amuria plate that consequently cannot be considered as a rigid block such distributed deformation is likely to have a strong effect on the structure of the future continental margin if extension evolves towards the formation of oceanic crust

rift flank uplift in tibet evidence for a viscous lower crust analysis of high resolution digital topography across 21 plio pleistocene grabens in central and southern tibet shows that high topography occurs along the graben margins profiles across these marginal highlands fit curves describing the flexure of an elastic plate under vertical loading we interpret the features as flexural rift flank uplifts resulting from regional isostatic compensation of low density grabens the short wavelengths of the uplifts 30 40 km imply a low flexural rigidity for the plateau d 2 6 1020 n m similar to values inferred for the basin and range of western north america in addition the grabens themselves are extremely narrow compared to typical continental rifts taken together these observations suggest the existence of a low viscosity lower crust and suggest that isostatic compensation is taking place within this ductile lower crust simple channel flow calculations suggest that crustal viscosities below 1022 pa s are necessary to accommodate this compensation over reasonable timescales