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

numerical modeling of crustal block and fault dynamics earthquakes and slip rates in the tibet himalayan region the tibetan plateau and himalayans have resulted from the continuous indian and eurasian plate convergence following their initial collision about 55 million years ago earthquakes in the region occur mainly in response to the crustal motion and stress localization associated with this convergence to understand the basic features of the crustal motion and seismicity in the tibet himalayan region we develop a numerical model of block and fault dynamics the model structure is composed of six major upper crustal blocks separated by fault planes these blocks are assumed to be perfectly rigid and move as a consequence of the indian plate push and of a flow of the lower crust deformations take place along the fault planes separating the blocks the interaction of the blocks along the fault planes is visco elastic as long as the ratio of the shear stress to the difference between the pore pressure and normal stress remains below a critical strength level when the critical level is exceeded in some part of a fault plane an earthquake stress drop occurs causing also failures in adjacent parts of the fault plane the stress drop affected parts of the fault plane enter in a state of creep immediately after the earthquake and the creep lasts until the stress falls below a certain level we develop several sets of numerical experiments to analyze the earthquake clustering frequency to magnitude relationships earthquake focal mechanisms and fault slip rates in the model large events in the numerical experiments cluster on the fault segments associated with the himalayan frontal thrust as well as at some internal faults of the tibetan plateau the clustering of earthquakes on a given fault is a consequence of the dynamics of the regional fault system rather than that of the fault only we show that variations in the relationship of magnitude to frequency of the events are associated with changes in the motion of the upper crustal blocks and depend on the rheological properties of fault plane zones the focal mechanisms of model events are found to be consistent with that of earthquakes in the region the synthetic moment rate captures the observed rate of seismic moment release we demonstrate in the model that the present crustal motion in the region is indeed governed by the north northeastern movement of india toward eurasia and the movement of the lower crust variations in the rheological properties of fault plane zones and or in the motion of the lower crust influence rates of the crustal block displacements and slips at the faults separating the blocks this can explain the discrepancies in estimates of slip rates over short and long time scales at major faults in the region

new empirical relationships among magnitude rupture length rupture width rupture area and surface displacement source parameters for historical earthquakes worldwide are compiled to develop a series of empirical relationships among moment magnitude m surface rupture length subsurface rupture length downdip rupture width rupture area and maximum and average displacement per event the resulting data base is a significant update of previous compilations and includes the additional source parameters of seismic moment moment magnitude subsurface rupture length downdip rupture width and average surface displacement each source parameter is classified as reliable or unreliable based on our evaluation of the accuracy of individual values only the reliable source parameters are used in the final analyses in comparing source parameters we note the following trends 1 generally the length of rupture at the surface is equal to 75 of the subsurface rupture length however the ratio of surface rupture length to subsurface rupture length increases with magnitude 2 the average surface displacement per event is about one half the maximum surface displacement per event and 3 the average subsurface displacement on the fault plane is less than the maximum surface displacement but more than the average surface displacement thus for most earthquakes in this data base slip on the fault plane at seismogenic depths is manifested by similar displacements at the surface log linear regressions between earthquake magnitude and surface rupture length subsurface rupture length and rupture area are especially well correlated showing standard deviations of 0 25 to 0 35 magnitude units most relationships are not statistically different at a 95 significance level as a function of the style of faulting thus we consider the regressions for all slip types to be appropriate for most applications regressions between magnitude and displacement magnitude and rupture width and between displacement and rupture length are less well correlated and have larger standard deviation than regressions between magnitude and length or area the large number of data points in most of these regressions and their statistical stability suggest that they are unlikely to change significantly in response to additional data separating the data according to extensional and compressional tectonic environments neither provides statistically different results nor improves the statistical significance of the regressions regressions for cases in which earthquake magnitude is either the independent or the dependent parameter can be used to estimate maximum earthquake magnitudes both for surface faults and for subsurface seismic sources such as blind faults and to estimate the expected surface displacement along a fault for a given size earthquake

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 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

crustal deformation of the lhasa terrane tibet plateau from project indepth deep seismic reflection profiles international deep profiling of tibet and the himalaya indepth deep reflection data in the yangbajain damxung graben of southern tibet yield evidence of magmatism and deformation beneath the southern lhasa terrane and yarlung zangbo suture shallow reflections and low velocity first arrivals indicate a thin quaternary graben fill of generally less than a few hundred meters thickness underlying stratified reflections extending to a maximum of about 11 5 km depth likely originate from deformed paleozoic cenozoic supracrustal strata of the lhasa terrane a prominent undulatory band of reflections within the crystalline basement extends beneath the length of the yangbajain damxung graben depth ranges from 12 km to 18 km this horizon has been interpreted to mark the top of a midcrustal partial melt zone underlying southern tibet the undulatory character of this horizon additionally suggests that it may be tectonically deformed a somewhat deeper subhorizontal wide angle reflection extends southward beneath the outcrop of the yarlung zangbo suture indicating that the suture is cut off or superposed by a younger structure at depth three gently north dipping reflections at 19 km 24 km and 27 km depth beneath the gangdese batholith are suggestive of a midcrustal duplex and may mark the northward downdip extension of the late oligocene early miocene gangdese thrust system a prominent 40 north dipping reflection imaged in the deep crust beneath the gangdese batholith between 40 and 60 km depth might mark the downdip expression of the yarlung zangbo suture or alternatively a younger reverse fault in the lower crust or both viewed in aggregate the reflection data are suggestive of moderate postcollisional shortening of the upper crust of the lhasa terrane accompanied by melting of the middle crust although the data are permissive of wholesale underthrusting or fluid injection of indian continental crust beneath the lhasa terrane they show no direct evidence for this having occurred

numerical models of the geodynamo and observational constraints 1 the past few years have seen the emergence of a large number of numerical simulations of the geodynamo in parallel both new and old geomagnetic archeomagnetic and paleomagnetic observations have been interpreted as actual geomagnetic features and used as constraints for dynamo models naturally model predictions should be tested against actual characteristics of the geomagnetic field despite huge differences sometimes in excess of a billion between the values of parameters used in the simulations and those estimated for the earth it is intriguing that many available simulations succeed in producing largely axial dipolar magnetic fields with weaker nondipolar structures in agreement with the first order characteristics of the geomagnetic field yet when considering finer characteristics there are significant differences and failures to actually produce a number of fundamental characteristic features in this presentation we first review numerical results obtained to date then we attempt to summarize which field characteristics derived from observational data sets can be considered robust on the basis of simple criteria used to evaluate the degree of confidence that can be placed in each datum we sort presumably characteristic geomagnetic features into three categories robust controversial and unlikely we conclude that numerical models should be illustrated with a number of key predictions averaged over at least 10 dipole diffusion times these predictions should be tested against the subset of robust observations only controversial observations should await additional confirmation

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

numerical modelling of earthquake flow in the southeastern carpathians vrancea effect of a sinking slab abstract the earthquake prone vrancea region is modelled as a system of perfectly rigid blocks separated by infinitely thin plane faults the interaction of the blocks along the fault planes and with the surrounding medium is assumed to be a viscous elastic the displacements of the block system are caused by motions of boundary blocks the velocities of the motions are found from a model of mantle flows induced by a sinking slab beneath the vrancea region when a ratio of stress to pressure for some portion of a fault plane exceeds a certain strength level a stress drop earthquake occurs as a result of the numerical simulation a catalog of synthetic earthquakes is produced several numerical experiments for various model parameters show that the spatial distribution of synthetic events is significantly sensitive to the directions of the block movements small variations in a slab rotation control the pattern of the synthetic seismicity the results of the analysis show that the catalogs obtained by the simulation of the block structure dynamics have certain features similar to those of the real earthquake catalog of the vrancea region

seismic evidence for stepwise thickening of the crust across the ne tibetan plateau we performed a receiver function analysis on teleseismic data recorded along two 550 km long profiles crossing the northeastern tibetan plateau results from time to depth migration grid search vp vs determination and simulated annealing inversion of waveforms reveal that the crust thickens from 50 km near the northern edge of the plateau to 80 km south of the jinsha suture in the qiang tang block crustal thickening occurs in staircase fashion with steps located beneath the main reactivated sutures the vp vs ratio close to the global continental average does not suggest widespread partial melting but rather a more usual separation between an upper felsic and a lower mafic part within the northeastern tibetan crust

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

geomorphic and stratigraphic evidence of quaternary diapiric activity enhanced by fluvial incision navarres salt wall and graben system se spain abstract the structural and stratigraphic geological features associated with salt structures have received extensive attention from geologists however there is limited work on the geomorphic and quaternary record of diapiric activity despite its practical implications e g salt and hydrocarbon extraction geostorage activities this work analyses the quaternary geomorphic and stratigraphic evidence of diapiric activity in the navarres salt wall se spain developed along the axis of a neogene graben this salt system is located in a region characterized by a peculiar network of orthogonal grabens that control the configuration and evolution of the drainage network the protruding salt walls are spatially associated with the erosionally unloaded and deeply entrenched graben sections situated close to the regional base level evidence of recent current activity in the navarres salt wall include 1 internally drained areas in a marginal withdrawal basin with long sustained aggradation 2 distorted drainage network defeated and deflected streams wind gaps knickpoints changes in fluvial style which changed from an initial axial longitudinal pattern to a longitudinal marginal distribution at the flanks of the salt wall 3 development of diapiric fault scarps at the edge of the salt wall 4 tilted fluvial terraces dipping away from the salt wall and with localized sediment thickening the available data indicate along strike variability in deformation style and long term vertical deformation rates these range from 0 09 mm yr to significantly higher values 0 2 0 4 mm yr in the section where the salt wall displays a more prominent relief and is affected by deep fluvial entrenchment the regional analysis and the characterization of the navarres salt wall suggest that quaternary diapirism in the region is enhanced by erosional unloading related to fluvial entrenchment which expands from the regional base level through the upstream propagation of an incision wave

rapid slip along the central altyn tagh fault morphochronologic evidence from cherchen he and sulamu tagh consistent with the d 18 o record from the guliya ice cap in the west kunlun the features of interest were all formed by glacial and fluvial processes subsequent to marine isotope stage 5e with the youngest features having formed during the early holocene optimum this near field morphochronological slip rate is averaged over many earthquake cycles and is hence little affected by interseismic strain it is kinematically consistent with other somewhat lower geomorphic slip rate measurements to the east the average rate and lower bounds obtained from alternate interpretational models 18 4 mm yr cannot be reconciled with the most rece geodetic measurements 7 mm yr suggesting that interseismic strain and interactions with adjacent faults may lead to disparate geologic and geodetic rate estimates this late pleistocene holocene morphochronologic rate would imply that at this longitude the altyn tagh fault on the north edge of tibet might absorb almost as much of india s convergence relative to siberia as the himalayan main frontal thrust does on the southern edge of the plateau index terms 1035 geochemistry geochronology 1208 geodesy and gravity crustal movements intraplate 8110 8102 tectonophysics continental contractional orogenic belts 8107 tectonophysics continental neotectonics keywords slip rates cosmogenic dating indo asian collision

holocene left slip rate determined by cosmogenic surface dating on the xidatan segment of the kunlun fault qinghai china cosmogenic dating using in situ sup 26 al and sup 10 be in quartz pebbles from alluvial terrace surfaces constrains the late holocene slip rate on the xidatan segment of the kunlun fault in northeastern tibet two terrace risers offset by 24 3 and 33 4 m having respective ages of 1799 388 and 2914 471 yr imply a slip rate of 12 1 2 6 mm yr the full range of ages obtained le 22 8 k y most of them between 6 7 and 1 4 k y confirm that terrace deposition and incision hence landform evolution are modulated by post glacial climate change coupled with minimum offsets of 9 12 m this slip rate implies that great earthquakes m approximately 8 with a recurrence time of 800 1000 yr rupture the kunlun fault near 94 e

active faulting and tectonics of burma and surrounding regions this study of the active tectonics of burma and surrounding regions is based mostly on an interpretation of landsat imagery and on fault plane solutions for shallow and intermediate earthquakes fold axes in the indoburman ranges and an east dipping inclined zone of intermediate depth earthquakes suggest that a slab of oceanic lithosphere was recently subducted to the east under the indoburman ranges fault plane solutions of shallow earthquakes however do not show underthrusting at the present time instead p axes are oriented roughly north south parallel not perpendicular to the fold axes in the ranges p axes for fault plane solutions of intermediate shocks also trend north south parallel to the strike of the subducted slab and may indicate that the hanging slab is being dragged north by india through the surrounding asthenosphere the sagaing fault probably accommodates most of the right lateral slip of india past indochina but india s northward movement also causes internal deformation within burma thailand and yunnan differences between the tectonic fabric seen on the landsat images and the tectonic regime compatible with the fault plane solutions throughout the region imply major changes in the directions of relative motion in the orientations of strain and in the style of deformation in the late cenozoic presumably these changes reflect the transition from subduction to collision and strike slip motion as india penetrated into eurasia past burma large parts of burma and thailand may have rotated several tens of degrees clockwise in this process

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

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

the role of water in the deformation of dunite specimens of two natural dunites anita bay of 100 m grain size and aheirn of 900 m grain size have been deformed in experiments at temperatures between 1200 and 1400 c and constant strain rates between 10 4 and 10 6 s 1 in a gas medium deformation apparatus at 300 mpa confining pressure prior to these experiments the specimens were given a heat treatment in a controlled oxygen fugacity furnace at 1200 c for 60 hours to dehydrate the trace layer silicate minerals present and drive off adsorbed water the predried specimens are much stronger than specimens tested under wet conditions in an earlier study the predried specimens also do not show the apparent grain size dependence seen under wet conditions infrared absorption measurements on the present specimens and those of the earlier study indicate that for a given rock flow stress can be correlated with the intensity of a residual water related broad absorption the addition of water to a predried specimen results in both a mechanical reweakening and an intensifying of the broadband infrared component the weakening of the dunites under wet conditions is probably strongly influenced by the presence of small amounts of a water related species in the grain boundaries but intragranular weakening may also be contributing it is suggested that the concentration of the water related species in the grain boundaries is governed by its concentration in the small melt fraction present the amounts of water that have been found to be of importance for the strength of these dunites are of the order of 0 01 wt amounts much less than most estimates of the water content of the earth s undepleted upper mantle such water may have an important role in determining upper mantle flow behavior stress strain curves for all experiments are available with entire manuscript on microfiche order from american geophysical union 2000 florida avenue n w washington dc 20009 document b84 008 2 50 payment must accompany order

insar observations of low slip rates on the major faults of western tibet two contrasting views of the active deformation of asia dominate the debate about how continents deform i the deformation is primarily localized on major faults separating crustal blocks or ii deformation is distributed throughout the continental lithosphere in the first model western tibet is being extruded eastward between the major faults bounding the region surface displacement measurements across the western tibetan plateau using satellite radar interferometry insar indicate that slip rates on the karakoram and altyn tagh faults are lower than would be expected for the extrusion model and suggest a significant amount of internal deformation in tibet

lhasa block and bordering sutures a continuation of a 500 km moho traverse through tibet a north south wide angle fan profile through the yarlung zangbo suture across the lhasa block and through the bangong nujiang suture shows several significant and sharp changes in crustal thickness both sutures appear as vertical 20 km steps in the moho and may have been the loci for eastward strike slip motion of the tibetan lithosphere