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

extension during continental convergence with application to the tibetan plateau the tibetan plateau is the product of crustal thickening caused by the collision between india and asia and is the largest active example of extensional tectonics in a zone of continental collision throughout most of the tertiary the tectonics of the plateau were dominated by north south shortening a significant proportion of which took place on east west striking thrust faults for the last 5 m y or so the plateau has been thinning by the mechanism of extension on north south trending normal faults previous investigations of the collision have been able to account for the large scale features of the tertiary deformation but have failed to explain the transition in the late tertiary to quaternary strain rate field of the plateau from north south compression to east west extension the transition could in principle be effected by a reduction in the rate of convergence between india and asia or by a uniform reduction in strength of the whole continental lithosphere of asia it could not however be effected by a reduction in strength of the elevated region alone this produces increased compressional strain rates in the weakened zone an alternative explanation for the transition to extension comes from considering the thermal evolution of thickened continental lithosphere the lower part of the lithosphere consists of a thermal boundary layer which when thickened by horizontal shortening is colder and denser than its surroundings convective instability of the thickened thermal boundary layer and its replacement by hot asthenosphere would rapidly raise the surface elevation and gravitational potential energy of the overlying part of the lithosphere the convective instability would happen in a time brief compared with the collision time scale 50 m y in the case of india and asia but would only occur after there had already been substantial thickening of the lithosphere numerical experiments show that for a range of lithospheric parameters the increase in surface height as much as 2 km and of potential energy 5 to 10 1012 n m 1 resulting from convective instability of the lower lithosphere are sufficient for east west extension to replace north south compression as the dominant feature of the stress field 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

neotectonics of the thakkhola graben and implications for recent activity on the south tibetan fault system in the central nepal himalaya the thakkhola graben is one of many north trending rifts that define the neogene structural pattern of the southern tibetan plateau lying at the southern margin of the plateau and extending to the crest of the himalaya the graben provides an opportunity to evaluate the kinematic relationships between east west extensional strain in southern tibet and north south extensional strain in the himalaya neotectonic and structural mapping of the dangardzong fault along the western margin of the graben reveals a southward decreasing component of normal slip coupled with a southward increasing component of right lateral slip that affects pleistocene basin fill sediments we present 14c ages for river terraces in the thakkhola graben that provide a ca 17 2 ka minimum age on the latest stage of dangardzong fault movement near the southern termination of the graben the dangardzong fault apparently offsets the annapurna detachment an early miocene strand of the east striking south tibetan fault system however the dangardzong fault itself terminates against a young i e younger than ca 17 2 ka strand of the south tibetan fault system the dhumpu detachment structural relationships among the dangardzong annapurna and dhumpu faults suggest that the dangardzong structure is a tear fault in the south tibetan allochthon that accommodates differential amounts and rates of displacement along the south tibetan fault system thus although the south tibetan fault system first developed as part of the structural architecture of the himalaya in miocene time at least some strands have been active as recently as the pleistocene in a regional context the south tibetan fault system serves to accommodate the strain gradient between extension in tibet and shortening in the himalaya

coupling and decoupling of crust and mantle in convergent orogens implications for strain partitioning in the crust topographic relief strain partitioning and crust mantle coupling within convergent orogens may be systematically related to the distribution of strength as a function of depth within continental crust this result is obtained by solving analytically for basally driven three dimensional flow within an idealized crust with a laterally invariant depth dependent newtonian viscosity two end member flow states are shown to occur in the first crustal flow is directly coupled to the underlying mantle and is controlled by stresses transmitted vertically through the crust in the second upper crustal to midcrustal flow is decoupled from the underlying mantle and is controlled by stresses transmitted horizontally through the crust this model is applied to deformation of continental crust above plate like mantle that undergoes oblique convergence across a narrow suture if no significant low viscosity zone develops within the crust crustal deformation is coupled to motion of the underlying mantle and a triangular mountain range develops surface shortening occurs at the foot of the mountains normal to the trend of the range while strike slip and extension are partitioned along the axis of the orogen if a low viscosity zone is initially absent but develops during crustal thickening a steep sided flat topped plateau forms crustal deformation is coupled to mantle motions on the plateau flanks but may be strongly decoupled beneath the plateau surface shortening normal to the plateau occurs at the foot of the plateau while surface extension normal to the plateau occurs at the plateau edge if a low viscosity zone is present in the lower crust prior to convergence a wide orogen with low topographic relief develops crustal deformation is decoupled from the mantle except at the edges of the deforming region and surface shortening is partitioned uniformly across the belt this study shows that horizontal and vertical partitioning of strain occur even in the absence of lateral strength variations and that shortening and extension may exist simultaneously in different parts of an orogen similar morphology and strain partitioning can be observed in many real orogenic belts including the eastern alps the tibetan plateau and the lachlan tasman fold belt

evidence for tibetan plateau uplift before 14 myr ago from a new minimum age for east west extension important changes in south asian climate occurred in the late miocene epoch 8 myr ago 1 2 and these have been attributed by some researchers to uplift of the tibetan plateau at about the same time3 5 unfortunately this link has been difficult to test because the timing of plateau uplift remains poorly constrained by independent evidence one way to determine the minimum age of uplift is to establish the initiation age of the north striking normal fault systems in southern tibet that are widely regarded6 10 as being related to gravitational collapse of the tibetan plateau here we report an40ar 39ar age of 14 myr for hydrothermal mica from an extensional fracture belonging to such a fault system in north central nepal this age implies that east west extension began before 14 myr ago in at least some parts of the tibetan plateau suggesting that the plateau attained its high mean elevation well before late miocene time

the role of a heterogeneous inclusion during continental collision abstract we study a simple model of the deformation due to collision between two continental plates the model is patterned after the geometry of the india asia collision but it may be applicable to more general situations the collision is modelled by the indentation of a viscoplastic plate asia by a rigid punch india a finite element approach capable of dealing with finite deformation heterogeneity and non linear rheologies is used to study the role of an inclusion of increased rigidity on the viscoplastic flow of the plate the heterogeneous zone was chosen to approximate the shape and position of the tarim basin behind tibet the lateral boundaries of the plate were defined so as to approximate those that are believed to apply for the asian plate i e fixed boundaries to the west and north free boundaries to the west and south several geometries and contrasts of rheology were analyzed to verify the general validity of our results the purpose was to determine the main kinematic features of the flow not to model the particular case of tarim it is found that the heterogeneity affects the flow in a drastic way even for small contrasts of rigidity and for small size heterogeneities the main effect of the heterogeneity is to create two important shear zones or internal boundary layers the most important is in front of the inclusion separating the main flow from the inclusion as the flow proceeds the inclusion rotates and the shear zone progressively bends in qualitative agreement with the bending of the large strike slip faults like the altyn tagh that separates the tarim from tiber our model of continuous deformation does not allow for the formation of faults or fractures but the shear zones may be understood as the initial stages for the generation of fractures or velocity discontinuities another shear zone develops to the left of the inhomogeneity near the fixed boundary and in the wake of the punch the role of these shear zones is to isolate the heterogeneity from the main flow to facilitate its rigid rotation inside the flow of the softer plate a final observation is that in the presence of the heterogeneity the lower left hand corner rotates almost like a rigid block producing a sharp bend on the left hand side of the plate as the inclusion rotates pushed by the punch this rigid block rotates clockwise to open up the space necessary for the homogeneity

lithospheric scale folding numerical modelling and application to the himalayan syntaxes we describe the eastern and western himalayan syntaxes which are large scale antiforms situated at geodynamically similar locations and the metamorphic evolution of which is coeval in the india asia collisional history to understand the mechanical plausibility of the structural interpretation we present two dimensional finite element modelling of lithospheric folding the models reveal the coeval development of adjacent synformal basins analogous to the peshawar and kashmir basins on both sides of the western syntaxis similarities between geological data and calculated models indicate that lithospheric buckling is a basic response to large scale continental shortening and an efficient mountain building process

geodetic measurement of crustal motion in southwest china global positioning system measurements performed over a 2 4 yr period confirm the xianshuihe xiaojiang fault system as one of the primary active structures in southwest china stations southwest of these faults show southerly motions of 5 15 mm yr relative to the western sichuan basin and stations in a 200 km geodetic network located northwest of the basin move at only 0 5 mm yr these results are consistent with clockwise rotation of southwestern sichuan and western yunnan about the eastern tibetan syntaxis accommodated by left lateral slip of 12 4 mm yr on the xianshuihe xiaojiang fault system the results imply that if eastward extrusion of crustal material from the plateau occurs at present it is not accommodated by east west shortening along the margin of the plateau and must involve wholesale eastward motion of low lying regions to the east

the process and mechanism of the rise of the qinghai tibet plateau abstract during the period from the late cretaceous to eocene the qinghai tibet plateau started a new stage of continental collision crustal shortening and thickening and continuous rise the rising process of the plateau progressed through three stages rise due to subduction collision rise due to convergence compression and rise due to isostatic adjustment the crustal thickening and shortening of the plateau occurred through different modes of motion of the layers of the earth s crust the model of the plateau rise may be generalized as intracontinental convergence crustal thickening isostatic adjustment

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

global rates of geologic cycling tectonic diffusion of upper crustal lithosomes abstractage frequency distributions of various rock types now exposed at earth s surface primarily reflect characteristic crustal depths at which they originate and the cumulative influence of vertical tectonic movements that displace them abundances of mid to upper crustal lithosomes increase rapidly to some modal age and then gradually decrease with increasing age paucity of exposures younger than the modal age reflects an intrinsic interval of geologic time needed for tectonic movements to exhume rock units while increasing rarity of exposures older than the modal age reflects the overall dispersion uplift and erosion subsidence and burial from initial lithosome positions relative to earth s surface we develop a random walk model that describes the geologic cycling of crustal rock as a tectonically driven diffusive system because at a globally averaged spatial and temporal scale tectonic subsidence burial and uplift erosion are essentially random processes from this we produce solutions to

northeastward growth of the tibet plateau deduced from balanced reconstruction of two depositional areas the qaidam and hexi corridor basins china we address the problem of late cenozoic uplift erosion and growth of northeastern tibet by reconstructing from isopach maps and drill holes the volumes and masses of sediments deposited in the qaidam and hexi corridor basins since 35 myr ago the mass budget is based on simple geometrical assumptions such as regional similarity of the thickness ratios between strata of different ages in the qaidam where our record extends back to the oligocene the budget shows a huge rise of the accumulation rates after the beginning of the pliocene 5 3 ma the early pliocene seems to be the period of maximum deposition with accumulation rates in excess of 1 mm yr 1 2 7 kg m 2 yr 1 of compacted rocks throughout the basin there also seems to be a southeastward shift of the largest depocenters between the upper pliocene 3 4 1 6 ma and the quaternary in the hexi corridor sedimentation is confined to small foreland flexural depressions associated with the frontal thrusts of the qilian shan and occurs at an average rate one order of magnitude smaller than in the qaidam basin the accumulation rate is maximum in the quaternary the sedimentation history appears to support a plateau building mechanism resulting from the combination of two geologicallly common processes crustal scale thrusting and sedimentary basin infilling the time needed to completely fill the qaidam basin and make its catchment a plateau closely resembling that of the highest part of tibet qangtang is of the order of 9 myr the mechanism now at work north of the kunlun which involves rapid infilling of broad flat areas separated by relatively narrow mountain ranges has thus probably been important in producing the high smooth topography that characterizes much of central tibet

a propagation model for the internal solitary waves in the northern south china sea a two dimensional regularized long wave equation model is developed to study the dynamic mechanisms of the propagation and evolution of the internal solitary waves isws in the northern south china sea scs it is shown that the bottom topography would cause the polarity reversal of isws the change of the local wave crestline shape and some diminution in wave amplitude even if the isws are induced at the small sill channel along the luzon strait they could propagate westward with their crestlines covering a large area in the latitudinal direction in the northern scs when there are two trains of isws propagating from the same source site with a time lag but different amplitudes of initial solitons the latter train of isws with a larger amplitude may catch then swallow the former one with a smaller amplitude and the wave amplitude of the merged isw train decreases while the wave number increases when there are two trains of isws propagating from the different source sites at the same time with the same amplitude of initial solitons the crestlines of the two isw trains may meet and a new leading soliton is induced at the connection point once the isw trains collide with the island before the island a weak isw train is reflected behind the island the former crestlines of the isw train are torn by the island into two new trains which may reconnect after passing around the island the propagation direction the wave amplitude and the reconnection point of the new merged isw train behind the island depend on the relative orientation of the original soliton source site to the island

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

three dimensional elastic wave velocity structure of the western and central tien shan arrival times of compressional p and shear s waves generated by earthquakes at local and teleseismic distances and recorded by seismographs located in the western and central tien shan are used to determine one and three dimensional elastic wave velocity structures of the crust and upper mantle beneath the mountain belt the best fit one dimensional structures suggest that the average depth of the mohorovicic discontinuity in this area is 50 km the three dimensional structure of the upper crust reveals thick sediments within each of the major depressions in the region a 7 km thick wedge of sediment beneath the chu depression is outlined at depth by a south dipping plane of seismic activity suggesting the presence of an active decollemont these low velocities extend continuously to the southeast toward issyk kul suggesting a structural relationship between the two however rather than being consumed it appears that issyk kul is overthrusting the surrounding ranges the low velocity sediments in the fergana basin reach depths of 10 km and are bounded on three sides by amorphous bands of seismicity velocities at midcrustal depths generally are lower beneath the central tien shan than beneath the western tien shan this pattern becomes more evident in the uppermost mantle with p velocity contrasts of as much as 10 across a boundary that corresponds roughly to the geographical position of the talasso fergana fault the low velocities beneath the central tien shan exceed 150 km depth but do not appear to be deeper than 300 km depth there is no evidence for a lithospheric root beneath this part of the range rather the low velocities imply the presence of a positive buoyancy force uplifting the mountains evidence that this low velocity region existed before the collision suggests that the tien shan may not owe its rejuvenation simply to its location at the northern edge of a strong tarim basin but rather to an anomalous upper mantle that was easier to deform than the surrounding lithosphere

paleogeographic and tectonic evolution of southern tibet since middle cretaceous time new paleomagnetic data and synthesis during the second french chinese paleomagnetic campaign in tibet in the summer of 1981 we sampled the middle cretaceous red beds of the takena formation and the cenozoic volcanics of the lingzizong formation a detailed analysis of about 400 specimens of red beds reveals a stable component of normal magnetization n in the north of the lhasa block north of the nyainqentanglha range and two stable components in the south lt ht where the red beds are more strongly folded and are unconformably overlain by the lingzizong volcanics all three components are carried by hematite blocking temperatures range from 400 to 600 c for lt and up to the neel temperature of hematite for both n and ht lt is shown to postdate part of the folding of the takena whereas n and ht both indicate a positive fold test at the 99 probability level in addition a small less than 10 of the natural remanent magnetization antiparallel component is superimposed in most samples and may be consistent with a model of long term acquisition of a chemical remanent magnetization for the red beds analysis of the volcanics reveals a single stable component with mixed polarities carried either by hematite or magnetite the paleomagnetic directions indicate that 1 the lhasa block stood at 12 5 3 n by upper cretaceous and at 13 5 6 5 n in the paleocene eocene thus confirming the limited displacement of the margin of eurasia between albian aptian deposition of the red beds and eocene onset of the india asia collision these values imply crustal shortening by as much as 1900 850 km north of the block since the upper cretaceous and by 1400 1200 km south of the block since the onset of the collision 2 southern tibet formed a single tectonic unit by middle cretaceous time and 3 differential rotations are observed within the block and are related to the collision these rotations are discussed in the frame of a recent propagating extrusion tectonics model this study confirms the close paleogeographic relationship between indochina khorat plateau and southern tibet prior to the india asia collision

the effects of accretion erosion and radiogenic heat on the metamorphic evolution of collisional orogens petrological and thermochronological data provide our best record of the thermal structure of deeply eroded orogens and in principle might be used to relate the metamorphic structure of an orogen to its deformational history in this paper we present a two dimensional thermal model of collisional orogens that includes the processes of accretion and erosion to examine the p t evolution of rocks advected through the orogen calculated metamorphic patterns are similar to those observed in the field metamor phic temperatures depths and ages generally increase with distance from the toe of the orogen p t paths are anti clockwise with rocks heating during burial and early stages of unroofing followed by cooling during late stage unroofing the results indicate that peak metamorphic temperatures within the core of a collisional orogen and the distance from the toe of an orogen to the metamorphic core can be related to the relative rates of accretion erosion and plate convergence model orogens displaying high metamorphic temperatures 600 c are associated with low ratios of accretion rate to plate convergence velocity and with high heat flow through the foreland model orogens with metamorphic cores far from the toe of the orogen are associated with high ratios of accretion rate to erosion rate calculated metamor phic gradients mimic steady state geotherms and inverted thermal gradients can be preserved in the metamorphic record suggesting reconsideration of the concept that the metamorphic record does not closely reflect geothermal gradients within an orogen

crustal structure of the himalayan orogen at 90 east longitude from project indepth deep reflection profiles international deep profiling of tibet and the himalaya indepth deep seismic reflection profiles show that the indian lithosphere is underthrusting the central himalaya along a gently north dipping decollement that is traceable northward beneath the tethyan belt to 28 6 n and to a depth of about 45 km the decollement carries in its hanging wall a near crustal thickness slice of internally deformed indian continental basement and cover represented in outcrop by the greater himalayan belt and structurally overlying tethyan belt geometric relationships suggest that the decollement probably ramps downward beneath the northern tethyan belt to near the base of the crust and that the hanging wall crustal slice was detached from indian mantle lithosphere that presently underlies southern tibet to the north the north himalayan anticlinorium kangmar dome appears to be a large duplex ramp anticline in the hanging wall of the decollement the south tibetan detachment appears to have been imaged in two areas in the wagye la area it appears to follow the tethyan belt greater himalayan belt contact northward in the subsurface and to have been folded over the kangmar dome in the zherger la area it appears to have been active relatively recently as a 30 north dipping normal fault that cuts deeply into the greater himalayan belt allochthon palinspastic reconstruction indicates a minimum of about 326 km of shortening of the indian basement across the himalaya since the initiation of the main central thrust hence indian continental lithosphere largely stripped of its overlying crust can extend northward beneath the tibetan plateau to at least 32 n where earthquake seismological observations show that the properties of the upper mantle change markedly e g mcnamara et al 1995 taken together the seismic reflection and earthquake seismological observations lend support for the view that the indian mantle lid mechanically underplates roughly the southern half of the tibetan plateau the indepth reflection data do not show that indian crust attached to this lithosphere extends north of southernmost tibet kangmar dome nor do they negate this possibility similarly the reflection data do not yield direct evidence for or against lower crustal subduction beneath the himalaya

finite strain calculations of continental deformation 1 method and general results for convergent zones a finite element method is used to calculate stresses and strain rates in a thin viscous sheet representing the continental lithosphere that is bordered in one part by an indenting boundary we calculate solutions for velocity fields crustal thickness distributions strain rates rotation time integrated vorticity and finite strain ellipsoids to show the dependence of these quantities on the controlling parameters the stress strain exponent n and the strength of the sheet relative to the gravitational forces which may be expressed as a dimensionless number the argand number ar for n 3 and ar between 1 and 10 a plateau of thickened crust is formed in front of the indenter the strain histories of individual elements within the viscous sheet may be quite complex and the finite strain at the end of deformation may bear no obvious relation to the principal stress orientations during the deformation the strain rate fields are interpreted in terms of the style of faulting that would be expected in the brittle upper crust if it were coupled to the ductile but stronger upper mantle in general the length scale of deformation increases with ar and time since the start of the collision and decreases with n for ar 3 the increase in length scale is accompanied by a change in the style of deformation from predominantly crustal thickening with thrust faulting in the brittle layer to predominantly transcurrent deformation strike slip faulting in the brittle layer

mechanical model for the tectonics of doubly vergent compressional orogens a mechanical model of crustal shortening and deformation driven by the relative convergence of rigid underlying mantle plates explains many features of convergent orogens results based on numerical models and supported by sandbox models show that a coulomb crustal layer subject to basal velocity boundary conditions corresponding to asymmetric detachment and subduction of the underlying mantle passes through three stages of orogenic growth 1 block uplift bounded by step up shear zones 2 development of a low taper wedge over the underthrusting mantle plate and 3 development of a low taperwedge overlying the overthrusting mantle plate and verging in the opposite direction when modified by isostasy basal viscous flow surface erosion and denudation and sedimentation the resultant model orogens exhibit a variety of styles with characteristics in common with small rapidly denuded orogens large orogens with plateaus and extensional characteristics and active subduction margins with doubly vergent accretionary wedges and deformed fore arc basins

active thrusting and folding in the qilian shan and decoupling between upper crust and mantle in northeastern tibet fieldwork south of the city of gaotai gansu province china shows that active shortening of surface sediments in the foothills of the yumu shan a large fore mountain of the qilian shan at the northeastern edge of tibet involves both overthrusting and flexural slip folding north of this mountain we found and mapped a prominent north facing thrust scarp that offsets a holocene fan sloping gently 3 4 to the north part of this scarp appears to be related to the m 7 5 180 a d earthquake that may have led to the demise of the han dynasty city of luo tuo chen in the hexi corridor a set of 10 100 150 m long profiles measured across this scarp 3 2 m high on the average can be made to fit the diffusion degraded morphology of a surface break related to the 180 a d event using a value of about 3 3 m 2 10 3 yr for the mass diffusivity of fanglomerates in this part of gansu province smaller mountain facing scarps on a terrace capped foothill result from bedding slip concurrent with active folding of underlying steeply northdipping plioquaternary sandstone and conglomerate beds holocene uplift rates along the yumu shan which is only one of the qilian shan ranges are estimated to be between 0 4 and 1 9 mm yr which implies that much of the mountain formed in the quaternary the periclinal structure of the plioquaternary envelope under which the paleozoic core of the yumu shan plunges towards the west suggests that the whole 3200 m high mountain is a basement ramp anticline mountains striking parallel to the yumu shan with similar structure and comparable or greater sizes north and south of the hexi corridor probably also correspond to recent crustal ramp anticlines this implies that the wide mountainous upper crustal wedge making the northeastern edge of the tibet qinghai plateau is detached from the underlying lower crust and 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

large scale crustal deformation of the tibetan plateau the topography velocity and strain fields calculated from a three dimensional newtonian viscous model for large scale crustal deformation are generally in good agreement with results from geological geodetic and earthquake studies in and around the tibetan plateau provided that the model rheology incorporates a weak zone within the deep crust beneath the plateau equivalent to a viscosity of 1012 pa s within a 250 m thick channel or 1018 pa s within a 15 km thick channel model studies and observations show a plateau at 5 km elevation with steep topographic gradients across the southern and northern plateau margins and more gentle gradients across the southeastern and northeastern margins rapid shortening strain is concentrated along the lower portions of the northern and southern plateau margins at rates 20 mm yr model results show north south shortening 10 mm yr in reasonable agreement with gps data 5 8 mm yr of north south shortening across the northern two thirds of the plateau and east west stretching 10 15 mm yr across the eastern half of the high plateau in reasonable agreement with seismic geologic and gps data upper crustal material moves eastward from the plateau proper into a lobe of elevated topography that extends to the south and east clockwise rotation of material around the east himalayan syntaxis at rates up to 10 mm yr occurs partly as a result of dextral shear between indian and asian mantle at depth and partly as a result of gravitational spreading from the high plateau to the south and east there is little difference in model surface deformation for assumptions of moderately weak or extremely weak lower crust except in southern and northern tibet where margin perpendicular extension occurs only for the case of an extremely weak lower crust our results suggest that the tibetan plateau is likely to have gone through a two stage development the first stage produced a long narrow high orogen whose height may have been comparable to the modern plateau the second stage produced a plateau that grew progressively to the north and east east west stretching eastward plateau growth and dextral rotation around the east himalayan syntaxis probably did not begin until well into the second stage of plateau growth perhaps becoming significant after 20 m y of convergence

did the indo asian collision alone create the tibetan plateau it is widely believed that the tibetan plateau is a late cenozoic feature produced by the indo asian collision however because tibet was the locus of continental accretion and subduction throughout the mesozoic crustal thickening during that time may also have contributed to growth of the plateau this portion of the geologic history was investigated in a traverse through the central lhasa block southern tibet together with earlier studies our mapping and geochronological results show that the lhasa block underwent little north south shortening during the cenozoic rather our mapping shows that 60 crustal shortening perhaps due to the collision between the lhasa and qiangtang blocks occurred during the early cretaceous this observation implies that a significant portion of southern tibet was raised to perhaps 3 4 km elevation prior to the indo asian collision

the namche barwa syntaxis evidence for exhumation related to compressional crustal folding abstract geological observation in the eastern end of the himalayas shows that the asia india suture is folded and continues southward into india and burma where the continuation of the transhimalayan plutonic complex can be identified metamorphic rocks derived from india occur structurally below the suture in the core of a regional antiform isotopic and fission track dating are consistent with the geological interpretation in addition they establish cooling exhumation of rocks from c 30 km depth within the last 4 ma rapid exhumation of deep crustal levels in mountain systems is generally related to extension subsequent to thermal softening of a thickened continental lithosphere however in the eastern himalayas we argue that exhumation is caused by c 10 mm yr 1 erosion coeval with crustal scale folding the general history of this syntaxis resembles the evolution of the western himalayas syntaxis in pakistan

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

gps measurements of present day convergence across the nepal himalaya the high elevations of the himalaya and tibet result from the continuing collision between india and asia which started more than 60 million years ago1 4 from geological and seismic studies of the slip rate of faults in asia5 it is believed that approximately one third of the present day convergence rate between india and asia 58 4mmyr 1 is responsible for the shortening uplift and moderate seismicity of the himalaya great earthquakes also occur infrequently in this region releasing in minutes the elastic strain accumulated near the boundary zone over several centuries and accounting for most of the advance of the himalaya over the plains of india the recurrence time for these great earthquakes is determined by the rate of slip of india beneath tibet which has hitherto been estimated indirectly from global plate motions6 from the slip rates of faults in asia7 8 from seismic productivity9 and from the advance of sediments on the northern ganges plain10 here we report geodetic measurements using the global positioning system gps of the rate of contraction across the himalaya which we find to be 17 52 2 mm yr 1 from the form of the deformation field we estimate the rate of slip of india beneath tibet to be 20 5 2 mmyr 1 strain sufficient to drive one or more great himalayan earthquakes with slip similar to that accompanying the magnitude 8 1 bihar nepal 1934 earthquake may currently be available in western nepal

global positioning system measurements of indian plate motion and convergence across the lesser himalaya we use global positioning system gps measurements acquired from 1991 to 1995 to constrain the motion of sites in bangalore in southern india and kathmandu nepal relative to a global gps network these measurements permit estimates of the northward motion of the indian plate and convergence between the southern himalaya and the indian subcontinent the velocities of bangalore and kathmandu in the itrf92 reference frame agrees with that predicted by the nnr nuvel1a plate motion model for indian plate motion and differ from that predicted for the australian plate confirming the independent motion of the indian and australian plate fragments no significant motion was detected between bangalore and kathmandu during the three years from 1991 1994 even though kathmandu is located in the hanging wall of the active himalayan thrust system the himalayan thrust system is thought to accommodate 18 7 mm yr of convergence and has been the source of several historic m 8 earthquakes the absence of motion of kathmandu relative to the indian plate can be explained if the thrust system is presently locked south of the greater himalaya our preferred model has no steady slip on the detachment south of the greater himalaya and steady slip at a rate greater than 6 mm yr 1 3 of the long term convergence rate can be ruled out at 95 confidence level

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

synchronous anatexis metamorphism and rapid denudation at nanga parbat pakistan himalaya the nanga parbat haramosh massif of the western himalaya is a north trending half window of indian crust that provides spectacular exposures of precambrian basement gneisses that have been overprinted by himalayan metamorphism we report here petrologic data and u pb dates on zircon and monazite which document that nanga parbat gneisses underwent a pliocene pleistocene episode of high grade metamorphism and anatexis during an interval in which the nanga parbat massif was undergoing rapid denudation at mean rates of 5 mm yr we speculate that by initiating decompression melting this denudation may be at least partly responsible for the anatexis and high grade metamorphism

tectonic exhumation of the nanga parbat massif northern pakistan abstract structural analysis in the nanga parbat region of northern pakistan indicates that extensional deformation has been in part responsible for the exhumation of gneissic rocks of the indian plate basement complex that dominates the nanga parbat massif this massif has been mapped as a syntaxis at the boundary between the indian and asian plates in the northwest himalaya and has received recent attention for evidence of rapid exhumation within the last 10 ma field work along the main mantle thrust mmt and within the indian plate rocks in the area between babusar pass and toshe gali southwest of the peak of nanga parbat provided us with evidence from shear fabrics and a shallowly wsw plunging stretching lineation for a dominant phase of extensional deformation during this deformation rock of the kohistan sequence and the indian plate cover sequence moved to the wsw relative to rocks of the indian plate basement including the nanga parbat gneiss a post metamorphic ductile shear fabric pervades the mmt contact the indian plate cover sequence below the mmt but decreases in intensity in the indian plate basement and the nanga parbat gneiss using 39 ar 40 ar geochronology of hornblendes we have constrained the age of pre deformational amphibolite facies metamorphism to be 40 56 ma published cooling ages from the babusar pass region further document that the extensional deformation was completed by 20 ma we interpret this ductile extensional shear as the mechanism responsible for significant tectonic exhumation of the gneissic and granitic rocks of the nanga parbat massif

age of initiation of collision between india and asia a review of stratigraphic data abstract the collision of india with asia is perhaps the most profound tectonic event to have occurred in past 100 ma it is responsible for the uplift of the himalayas and tibetan plateau and has been argued to have been responsible for geological geochemical and climatological consequences of global extent yet the age of initiation of this collision remains poorly constrained the literature is replete with estimates that range from the late cretaceous 65 ma to latest eocene

velocities and propagation characteristics of pn and sn beneath the himalayan arc and tibetan plateau possible evidence for underthrusting of indian continental lithosphere beneath tibet new seismological observations on velocities and propagation characteristics of pn and sn waves beneath tibet can be interpreted although not uniquely to indicate the shallow angle underthrusting of the greater indian continental lithosphere beneath the tibetan plateau this is inferred from available geological and geophysical data and is now supported by new seismological observations the most significant observation is that high frequency sn waves shear waves that travel in the mantle lithosphere propagate efficiently in the uppermost mantle beneath the tibetan plateau except beneath the north central part of tibet the chang thang terrane furthermore sn propagates efficiently across the stable block of the tarim basin along the tien shan and himalayan mountains and across the indian shield apparent velocities of pn and sn waves that traverse the uppermost mantle beneath the tibetan plateau are 8 42 and 4 73 km s respectively these velocities are very similar to those beneath the himalayan mountains and the indian shield the efficient propagation and the velocities of pn and sn waves beneath tibet are also similar to those commonly observed beneath shield and stable continental regions these results are consistent with the underthrusting model and do not favor the class of alternative models in which the tibetan plateau is formed by shortening and thickening of hot crustal and uppermost mantle material in response to the convergence of india and asia

uplift of tibetan plateau a history of the elevation and crustal thickness of the tibetan plateau since the continental collision at 45 ma has been computed using a kinematic model based on plate tectonic reconstructions and conservation of crustal volume in our reconstructions the indochina block was attached to the southern margin of eurasia before being extruded to the southeast as evidenced by the opening of the south china sea this opening affected the mass input into tibet so that tibet did not uplift to an elevation higher than 1 km until early miocene when spreading of the south china sea ceased the model then gives a steady rise of tibet from middle miocene to its present elevation