a slab in depth three dimensional geometry and evolution of the indo australian plate geoscientists use plate tectonics to explain many aspects of both continental evolution and evolution of the planet as a whole the subduction of material at convergent plate boundaries forms a fundamental component to the theory of plate tectonics plates continents subduction zones and spreading centers all exhibit motion and geometric evolution so to try and resolve the past geometries of the planet geologists have utilized plate tectonic reconstructions here we present a three dimensional image of the subducted indo australian plate below southeast asia and show that the geometry of the subducted slab at depth is intimately related to the geometric evolution of se asia over the past 50 ma including the collision of india with the asian continent we show how the once semicontinuous subducting indo australian plate has been segmented during collision between india australian and the subduction margin to the north thus we have found that the geometry of the subducted plate should form a key component to the interpretation of the evolution of earth s surface as complexities and evolution of the subducted plate are manifest in the evolution of the overriding plate

shear wave splitting and subcontinental mantle deformation splitting observations are interpreted in terms of the strain induced lattice preferred orientation of mantle minerals especially olivine we consider three hypotheses concerning the origin of the continental anlsotropy 1 strain associated with absolute plate motion as in the oceanic upper mantle 2 crustal stress and 3 the past and present internal deformation of the subcontinental upper mantle by tectonic episodes it is found that the last hypothesis is the most successful namely that the most recent significant episode of internal deformation appears to be the best predictor of q for stable continental regions this is interpreted as fossil anisotropy whereas for presently active regions such as alaska the anisotropy reflects present day tectonic activity in the stable portion of north america there is a good correlation between delay time and lithospheric thickness this is consistent with the anisotropy being localized in the subcontinental lithosphere and suggests that intrinsic anisotropy is approximately constant the acceptance of this hypothesis has several implications for subcontinental mantle deformation first it argues for coherent deformation of the continental lithosphere crust and mantle during orogenies this implies that the anisotropic portion of the lithosphere was present since the deformational episode and rules out the addition of undeformed material to this layer by subsequent underplating or conductive growth of the thermal boundary layer one of the most important issues in the study of orogenies is the need to reconcile the formation of thickened lithosphere with the paradoxically high mantle temperatures often associated with orogenic episodes most efforts to date have focussed on modes of deformation whereby the cold lithospheric mantle is removed by convective instability or delamination and replaced by warm asthenosphere these models however are incompatible with the evidence for preserved coherent lithospheric deformation rather the deformed mantle appears to have been heated in place we suggest that the elevated mantle temperatures may be due to the strain heating accompanying the deformation

a no net rotation model of present day surface motions a significant portion of the earth s surface consists of zones of diffuse deformation the interior regions of these diffuse zones of deformation move at distinctly different velocities from that of adjacent plates and because of their complexities have been ignored in previous no net rotation nnr models we have calculated a new nnr model from a continuous velocity field that incorporates both rigid plate motion and velocity gradients within plate boundary zones when compared with earlier nnr models we find significantly different angular velocities for many plates differences between the nnr model presented here and earlier nnr models can be attributed to the effect of including velocity gradients in diffuse plate boundary zones as well as to the actual differences between geodetically derived surface motions and geologic estimates

contribution of gravitational potential energy differences to the global stress field summary modelling the lithospheric stress field has proved to be an efficient means of determining the role of lithospheric versus sublithospheric buoyancies and also of constraining the driving forces behind plate tectonics both these sources of buoyancies are important in generating the lithospheric stress field however these sources and the contribution that they make are dependent on a number of variables such as the role of lateral strength variation in the lithosphere the reference level for computing the gravitational potential energy per unit area gpe of the lithosphere and even the definition of deviatoric stress for the mantle contribution much depends on the mantle convection model including the role of lateral and radial viscosity variations the spatial distribution of density buoyancies and the resolution of the convection model gpe differences are influenced by both lithosphere density buoyancies and by radial basal tractions that produce dynamic topography the global lithospheric stress field can thus be divided into 1 stresses associated with gpe differences including the contribution from radial basal tractions and 2 stresses associated with the contribution of horizontal basal tractions in this paper we investigate only the contribution of gpe differences both with and without the inferred contribution of radial basal tractions we use the crust 2 0 model to compute gpe values and show that these gpe differences are not sufficient alone to match all the directions and relative magnitudes of principal strain rate axes as inferred from the comparison of our depth integrated deviatoric stress tensor field with the velocity gradient tensor field within the earth s plate boundary zones we argue that gpe differences calibrate the absolute magnitudes of depth integrated deviatoric stresses within the lithosphere shortcomings of this contribution in matching the stress indicators within the plate boundary zones can be corrected by considering the contribution from horizontal tractions associated with density buoyancy driven mantle convection deviatoric stress magnitudes arising from gpe differences are in the range of 1 4 tn m 1 a part of which is contributed by dynamic topography the egm96 geoid data set is also used as a rough proxy for gpe values in the lithosphere however gpe differences from the geoid fail to yield depth integrated deviatoric stresses that can provide a good match to the deformation indicators gpe values inferred from the geoid have significant shortcomings when used on a global scale due to the role of dynamically support of topography another important factor in estimating the depth integrated deviatoric stresses is the use of the correct level of reference in calculating gpe we also elucidate the importance of understanding the reference pressure for calculating deviatoric stress and show that overestimates of deviatoric stress may result from either simplified 2 d approximations of the thin sheet equations or the assumption that the mean stress is equal to the vertical stress

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

constraints on crustal composition beneath a metamorphic core complex results from 3 component wide angle seismic data along the eastern flank of the ruby mountains nevada metamorphic core complexes expose rocks deformed at deep upper to middle crustal levels during extreme crustal extension however mechanisms of crustal extension and exhumation of core complexes remain to be fully understood detailed study of crustal velocity structure and inferences about the composition of the crust beneath core complexes and nearby provide useful constraints on core complex evolution p and s wave velocity structures determined from seismic experiments along the eastern flank of the ruby mountains metamorphic core complex nevada show that the crust can be divided into three main layers corresponding to the upper middle and lower crust we interpreted crustal composition by integrating results of p wave velocities vp s wave velocities vs poisson s ratios seismic anisotropy and reflection character with published geologic maps of the area near surface estimates of vp vs and anisotropy of 1 90 4 8 km s 1 01 2 75 km s 0 25 0 33 and 0 6 2 5 respectively are consistent with surface exposures of unconsolidated to consolidated sedimentary rocks limestone dolomite siltstone sandstone porous sandstone conglomerate and weathered granite results from analysis of reflection responses vp vs and anisotropy also indicate that 1 upper crustal rocks most likely consist of metaquartzite schist granite gneiss and granite granodiorite with vp of 5 80 6 25 km s vs of 3 20 3 72 km s of 0 22 0 25 and anisotropy of 0 6 2 5 2 possible middle crustal rocks are paragranulite felsic granulite felsic amphibolite gneiss granite granodiorite and mica quartz schist with vp of 6 35 6 45 km s vs of 3 70 3 75 km s and of 0 24 and 3 lower crustal rocks most likely consist of granulite rather than amphibolite facies rocks with vp of 6 60 6 80 km s vs of 3 85 3 92 km s of 0 24 0 25 and anisotropy of 3 our principal conclusions are 1 significant addition of gabbroic rocks underplating is unlikely in the lower crust 2 lower crustal rocks were stretched into sub horizontal geometries with aligned minerals during extension creating seismic lamellae in the lower crust 3 present day seismic velocities of highly extended core complex crust and normally extended basin and range crust are similar and 4 orientations of fast shear waves near the surface and in the upper crust are sub parallel to the regional maximum horizontal compressive stress in the nevada part of the basin and range province

the hydraulics of a stratified fluid flowing through a contraction the steady hydraulics of a continuously stratified fluid flowing from a stagnant reservoir through a horizontal contraction was studied experimentally and theoretically as the channel narrows the flow accelerates through a succession of virtual controls at each of which the flow passes from subcritical to supercritical with respect to a particular wave mode when the narrowest section acts as a control the flow is asymmetric about the narrowest section supercritical in the divergent section and selfsimilar throughout the channel with increased flow rate a new enclosed self similar solution was found with level isopycnals and velocity uniform with depth

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

crustal structure of mainland china from deep seismic sounding data abstract since 1958 about ninety seismic refraction wide angle reflection profiles with a cumulative length of more than sixty thousand kilometers have been completed in mainland china we summarize the results in the form of 1 a new contour map of crustal thickness 2 fourteen representative crustal seismic velocity depth columns for various tectonic units and 3 a pn velocity map we found a north south trending belt with a strong lateral gradient in crustal thickness in central china this belt divides china into an eastern region with a crustal thickness of 30 45 km and a western region with a thickness of 45 75 km the crust in these two regions has experienced different evolutionary processes and currently lies within distinct tectonic stress fields our compilation finds that there is a high velocity 7 1 7 4 km s layer in the lower crust of the stable tarim basin and ordos plateau however in young orogenic belts including parts of eastern china the tianshan and the tibetan plateau this layer is often absent one exception is southern tibet where the presence of a high velocity layer is related to the northward injection of the cold indian plate this high velocity layer is absent in northern tibet in orogenic belts there usually is a low velocity layer lvl in the crust but in stable regions this layer seldom exists the pn velocities in eastern china generally range from 7 9 to 8 1 km s and tend to be isotropic pn velocities in western china are more variable ranging from 7 7 to 8 2 km s and may display azimuthal anisotropy

mantle flow deflected by interactions between subducted slabs and cratonic keels subducting slabs can influence mantle flow but the importance of neighbouring continental cratons is little understood geodynamical modelling constrained by seismic data that identify regions of mantle flow beneath the caribbean south american plate margin shows that the deep rooted south american craton acts to deflect and enhance mantle flow into a narrow channel

teleseismic shear wave splitting in se tibet insight into complex crust and upper mantle deformation we measured shear wave splitting of teleseismic xks phases i e sks skks and pks recorded by more than 300 temporary chinarray stations in yunnan of se tibet the first order pattern of xks splitting measurements shows that the fast polarization directions change at 26 27 n from dominant n s in the north to e w in the south while splitting observations around the eastern himalayan syntax well reflect anisotropy in the lithosphere under left lateral shear deformation the dominant e w to the south of 26 n is consistent with the maximum extension in the crust and suggest vertically coherent pure shear deformation throughout the lithosphere in yunnan however the thin lithosphere 80 km could account for only part 0 7 s of the observed splitting delay times t 0 9 1 5 s anisotropy in the asthenosphere is necessary to explain the nw se and nearly e w in these regions the ne sw can be explained by the counter flow caused by the subduction and subsequent retreat of the burma slab the e w is consistent with anisotropy due to the absolute plate motion in se tibet and the eastward asthenospheric flow from tibet to eastern china accompanying the tectonic evolution of the plateau our results provide new information on different deformation fields in different layers under se tibet which improves our understanding on the complex geodynamics related to the tectonic uplift and southeastward expansion of tibetan material under the plateau

gap winds in a fjord part i observations and numerical simulation abstract gap winds in howe sound british columbia are described and placed in context by reviewing studies of similar phenomena in other locations an observational program consisting of a surface mesonetwork and vertical soundings shows that gap winds vary considerably along and across the channel as well as vertically wind strength generally increases down channel and strongest winds are found below 1000 m depth results from application of a 3d mesoscale numerical model to a gap wind case compare reasonably well with observations model output reveals more details of horizontal and especially vertical flow structure than is possible from observations model vertical cross sections and froude number output indicate similarity with hydraulic flow this is further substantiated by a force balance analysis of model output

seismic wide angle constraints on the crust of the southern urals a wide angle seismic reflection refraction data set was acquired during spring 1995 across the southern urals to characterize the lithosphere beneath this paleozoic orogen the wide angle reflectivity features a strong frequency dependence while the lower crustal reflectivity is in the range of 6 15 hz the pmp is characterized by frequencies below 6 hz after detailed frequency filtering the seismic phases constrain a new average p wave velocity crustal model that consists of an upper layer of 5 0 6 0 km s which correlates with the surface geology 5 7 km depths at which the velocities increase to 6 2 6 3 km s 10 30 km depths at which on average the crust is characterized by velocities of 6 6 km s and finally the lower crust from 30 35 km down to the moho which has velocities ranging from 6 8 to 7 4 km s two different s wave velocity models one for the n s and one for the e w were derived from the analysis of the horizontal component recordings crustal sections of poisson s ratio and anisotropy were calculated from the velocity models the poisson s ratio increases in the lower crust at both sides of the root zone a localized 2 3 anisotropy zone is imaged within the lower crust beneath the terranes east of the root this feature is supported by time differences in the sms phase and by the particle motion diagrams which reveal two polarized directions of motion velocities are higher in the central part of the orogen than for the siberian and eastern plates these seismic recordings support a 50 56 km crustal thickness beneath the central part of the orogen in contrast to moho depths of 45 km documented at the edges of the transect the lateral variation of the pmp phase in frequency content and in waveform can be taken as evidence of different genetic origins of the moho in the southern urals

azimuthal anisotropy and mantle flow underneath the southeastern tibetan plateau and northern indochina peninsula revealed by shear wave splitting analyses abstract seismic azimuthal anisotropy beneath the transitional region between the southeastern tibetan plateau and the indochina peninsula an area in which the fast orientation of mantle anisotropy changes to dominantly e w from mostly n s on the plateau at the north is quantified using splitting of the sks skks and pks phases among the 50 stations with one or more splitting measurements 22 possess an azimuthal coverage that is adequate for the identification and characterization of complex anisotropy the resulting splitting parameters from 15 such stations exhibit systematic back azimuthal variations with a 90 periodicity which is consistent with a two layered anisotropy model the upper layer parameters are consistent with crustal anisotropy measurements obtained independently based on the sinusoidal moveout of p to s conversions from the moho with the fast orientations being mostly parallel to major shear zones the lower layer fast orientations and the fast orientations at stations with azimuthally invariant splitting parameters are mostly e w which is significantly different from the dominantly n s trend of the surface expression of major structural fabrics in the area they are also inconsistent with the absolute plate motion directions when combined with results from seismic tomography and focal mechanism solutions the observed azimuthal anisotropy can be adequately explained by the movement of the lithosphere relative to the underlain asthenosphere most likely associated with the westward rollback of the subducted indian plate

correlated crust and mantle strain fields in tibet vertical to subvertical planes of shear within the active crustal deformation field in tibet align with fast directions of shear wave polarization this observation suggests that the present day velocity gradient tensor field within the crust of tibet correlates with the lithospheric mantle velocity gradient tensor field beneath tibet this inference requires the following 1 the 100 axes of olivine are aligned within established zones of lithospheric mantle fabric or anisotropy built up through finite strain 2 present day shear in the mantle lithosphere is occurring within these established zones of fabric parallel to the direction of olivine 100 axis alignment 3 the zones of mantle lithospheric fabric or zones of mantle shear are aligned with zones of crustal shear faults the correlation of crustal and mantle strain fields most simply results from the fact that both crust and mantle lithosphere are under the influence of similar velocity boundary conditions furthermore the observations confirm the distributed nature of lithospheric mantle deformation beneath east central tibet and suggest that left lateral shear has been the dominant component of finite shear there

correction to seismic anisotropy of oceanic upper mantle shear wave splitting methodologies and observations by cecily j wolfe and paul g silver

shear wave splitting continental keels and patterns of mantle flow in this study we investigated the origin of seismic anisotropy in the mantle beneath north america in particular we evaluated whether shear wave splitting patterns in eastern north america are better explained by anisotropy caused by lithospheric deformation anisotropy due to mantle flow beneath the lithosphere or a combination of both we examined new measurements of shear wave splitting from the missouri to massachusetts broadband seismometer array moma the north american mantle anisotropy and discontinuity experiment nomad as well as splitting parameters from several previous studies we developed a simple finite difference model that approximates mantle flow around a complex three dimensional continental lithospheric keel to evaluate potential anisotropy from mantle flow beneath the lithosphere in eastern north america we compared shear wave splitting observations to predicted splitting parameters calculated using this mantle flow model our results indicate that a significant portion of observed shear wave splitting in eastern north america can be explained by mantle flow around the continental keel however shear wave splitting patterns in a few regions of eastern north america indicate that a component of lithospheric anisotropy must exist particularly in regions containing the largest keel thicknesses for eastern north america as well as for splitting observations in australia europe and south america we favor a model in which anisotropy is controlled by a combination of both lithospheric deformation and subcontinental mantle flow

interpretation of sks waves using samples from the subcontinental lithosphere abstract the seismic properties of five south african kimberlite nodules have been calculated from petrofabric measurements these garnet lherzolite and harzburgite composition nodules 56 80 olivine are known from previous studies to have originated at depths of 120 170 km in the subcontinental lithosphere their in situ seismic properties have been calculated by extrapolating the elastic constants to the appropriate temperature 900 1050 c and pressure 3 0 3 5 gpa conditions the average of the five samples has a maximum s wave anisotropy v s of 3 7 which is sufficiently high to explain previously reported teleseismic sks delay times t in continental shield areas of up to 1 7 s the biggest delay time is observed when the foliation is vertical the lineation horizontal and the fastest s wave polarized parallel to the lineation i e an orientation for transcurrent motion delay times of over 1 s can only be generated by this orientation a strong mantle crust mechanical coupling is suggested for such situations the potential use of seismic anisotropy as an indicator of strain in the lithosphere has been investigated using experimental and simulated fabric data for olivine the most abundant phase 70 in the upper mantle the v p and v s seismic anisotropy increases with fabric strength and finite strain recrystallization tends to reduce fabric strength and seismic anisotropy resulting in saturation values for experimental and simulated fabrics which correspond to approximately 50 60 strain or 8 9 v s for pure olivine aggregates a survey of 25 naturally deformed peridotites of oceanic origin suggests an average maximum s wave anisotropy for the olivine component of 9 or about 8 for a lherzolite or harzburgite rock composition when the orthopyroxene component is taken into account the ophiolite samples are twice as anisotropic as the kimberlite nodules if an average anisotropy value is representative of the lithosphere then the sks delay times represent variations in anisotropic layer thickness for delay times over 1 s as the orientation is constrained to be foliation vertical and lineation horizontal the magnitude of s wave anisotropy is less sensitive than v p to variations in olivine volume fraction in range 50 100 and to deformation or fabric intensity further suggesting that s wave anisotropy is particularly apt for the determination of the anisotropic lithosphere thickness for smaller delay times there is some trade off between structural orientation and thickness

mantle conveyor beneath the tethyan collisional belt abstract collisional belts are generated by the arrival of continental lithosphere into a subduction zone the tethyan suture from the bitlis to the himalayas is a prime example where the arabian and indian plates collided with eurasia during the cenozoic while the kinematics of this process are well established its dynamics are more uncertain india and arabia intriguingly keep advancing in spite of large collisional resisting forces and in the absence of a substantial upper mantle slab driving force at present day we perform global mantle circulation computations to test the role of deep mantle flow as a driving force for the kinematics of the tethyan collisional belt evaluating different boundary conditions and mantle density distributions as inferred from seismic tomography or slab models our results show that mantle drag exerted on the base of the lithosphere by a large scale convective conveyor belt with an active upwelling component is likely the main cause for the ongoing indentation of the indian and arabian plates into eurasia

the kaapvaal craton seismic anisotropy petrophysical analyses of upper mantle kimberlite nodules a dense network of seismic stations has been deployed on the kaapvaal craton south africa to investigate the upper mantle seismic structures in order to bring independent petrophysical constraints we analyze a direct sampling of the cratonic upper mantle and determine the seismic properties of 48 mantle nodules brought up to the earth s surface by kimberlite eruptions seismic properties of these nodules are calculated from the olivine and pyroxene crystal preferred orientations and the single crystal elastic constants despite variations in the nodules compositions microstructures and crystallographic preferred orientations seismic anisotropy is rather homogeneous throughout the craton mean s wave anisotropy is weak 2 64 which is compatible with the small measured sks wave splitting mean delay time of 0 62 s

lithospheric thickness of the chinese continent abstract we invert for the upper mantle temperatures of the chinese continent in the depth range of 70 240 km from a recent s velocity model the depth where temperatures intersect a mantle adiabat with a potential temperature of 1300 c is in close correspondence with the top of the seismic low velocity zone for most regions this correspondence implies that seismic lithosphere estimated from short time scale seismic information may be equivalent to the long time scale geodynamical lithosphere defining the 1300 c adiabat as coinciding with the lithospheric base we estimate the seismic thermal lithosphere thickness the estimated thickness shows obvious dependence on the tectonic settings beneath eastern china which mainly belongs to the circum pacific tectonic domain it has a thickness of 100 km and beneath the qinghai tibet plateau and south to the tarim craton which belong to the tethyan tectonic domain it has a thickness of 160 220 km the lithospheric thicknesses of the three large para platforms cratons range from 170 km for the western yangtze 140 km for tarim and 100 km for sino korean the three cratons may have been reshaped by phanerozoic tectonic activities and are thinner than most cratons in other continents

gap winds and wakes sar observations and numerical simulations the nature of terrain induced gap winds and wakes in the atmosphere is examined using surface wind data from synthetic aperture radar sar and the shallow water equations the shallow water model is used to predict the types of wake jet wind patterns that might occur behind an idealized pair of bell shaped hills with a gap between them a regime diagram is constructed based on the width of the gap and the upstream froude number specific predictions of the model are found to compare moderately well with sar data from four examples of airflow near unimak island in the aleutian chain the model predicts the observed wakes and jets including jets that exceed the upstream speed theoretical analysis considers the relative importance of rising terrain and narrowing valley walls in the acceleration of gap winds wind speeds in the wake region are controlled by the bernoulli function and regional pressure gap winds therefore are streams of air that have avoided bernoulli loss over the terrain by passing through gaps the speed of gap winds can exceed the upstream speed only in ridgelike situations when the regional leeside pressure is lower than the upstream pressure

subduction in the indo burma region is it still active the indo burma region which includes longitude 92 96 e latitude 20 26 n is a subduction zone where the indian plate underthrusts the south eastern asian plate but the nature of subduction is complex i examine here the distribution of the p t and b axes of 37 earthquakes the time period 1977 1995 depths 0 153 km and magnitudes mb 4 8 6 8 with respect to the geometry of the wadati benioif zone the analysis shows that the t axes are clustered close to the down dip direction of the subducting slab implying a predominant down dip tensional stress regime within the slab which is typical of intermediate depth seismicity of active subduction zones these results suggest that the subduction in the indo burma region is possibly active

diurnal and semidiurnal tides in global surface pressure fields global surface pressure data from 1976 to 1997 from over 7500 land stations and the comprehensive ocean atmosphere data set have been analyzed using harmonic and zonal harmonic methods it is found that the diurnal pressure oscillation s1 is comparable to the semidiurnal pressure oscillation s2 in magnitude over much of the globe except for the low latitude open oceans where s2 is about twice as strong as s1 over many land areas such as the western united states the tibetan plateau and eastern africa s1 is even stronger than s2 this is in contrast to the conventional notion that s2 predominates over much of the globe the highest amplitudes 1 3 mb of s1 are found over northern south america and eastern africa close to the equator here s1 is also strong 1 1 mb over high terrain such as the rockies and the tibetan plateau the largest amplitudes of s2 1 0 1 3 mb are in the tropics over south america the eastern and western pacific and the indian ocean here s1 peaks around 0600 0800 lst at low latitudes and around 1000 1200 lst over most of midlatitudes while s2 peaks around 1000 and 2200 lst over low and midlatitudes here s1 is much stronger over the land than over the ocean and its amplitude distribution is strongly influenced by landmasses while the land sea differences of s2 are small the spatial variations of s1 correlate significantly with spatial variations in the diurnal temperature range at the surface suggesting that sensible heating from the ground is a major forcing for s1 although s2 is much more homogeneous zonally than s1 there are considerable zonal variations in the amplitude of s2 which cannot be explained by zonal variations in ozone and water vapor other forcings such as those through clouds reflection and absorption of solar radiation and latent heating in convective precipitation are needed to explain the observed regional and zonal variations in s2 the migrating tides and 1 s 1 predominate over other zonal wave components however the nonmigrating tides are substantially stronger 2 s 2 than previously reported the amplitudes of both the migrating and nonmigrating tides decrease rapidly poleward with a slower pace at middle and high latitudes

seismic anisotropy of oceanic upper mantle shear wave splitting methodologies and observations we develop methodologies to obtain accurate measurements of shear wave splitting and apply these techniques to examine the pattern of oceanic upper mantle anisotropy to obtain high quality estimates of receiver splitting at island stations we devise a stacking method that finds the optimum splitting parameters and 95 error region for teleseismic shear wave phases from a suite of earthquake events to obtain additional measurements in oceanic regions we develop techniques to measure splitting parameters and errors for ss phases that sample upper mantle anisotropy at their bounce points however we find the data are often of low resolution and anomalous characteristics are sometimes found that make splitting difficult to interpret ten second splitting of ss is observed across the banjo seismic array but we cannot unambiguously attribute this signal to mantle anisotropy at the bounce the receiver splitting methods are used to assess the adequacy of a two anisotropic layer model for the pacific region with fast polarization azimuth in the lithosphere oriented in the fossil spreading direction and in the asthenosphere oriented in the absolute plate motion direction this model has been proposed to explain surface wave data in the pacific ocean but our splitting results demonstrate that oceanic anisotropy patterns are more heterogeneous than would be predicted while the island splitting measurements could reflect the influence of individual hotspot upwellings hotspot effects do not appear to be universally dominant we propose that splitting observations alternatively indicate broad scale differences in the underlying character of oceanic upper mantle anisotropy associated with coherant patters of lithospheric structure and asthenospheric flow in particular splitting surface wave models and regional studies all support a model where lithospheric anisotropy throughout the south pacific has been erased or reoriented toward the absolute plate motion direction whereas more limited observations in the north pacific indicate that the fossil lithospheric signature appears to be preserved

shear wave splitting in the southern margin of the ordos block north china 1 the southern margin of the ordos block in north china marks the transition from a stable cratonic block to active transextensional deformation in the weihe graben driven by the northeastward expansion of the tibet plateau we deployed 29 portable digital stations across this boundary and studied the upper mantle anisotropy using sks splitting the anisotropy is weak t 0 5 s at stations within the ordos block and in the north weihe graben but strong t 1 23 0 39 s in the south weihe graben and qinling orogen further south the directions of the fast s wave are 98 7 subparallel to the orientations of surface structures and the gps site velocities relative to stable asian continent we suggest that this e w oriented anisotropy is the result of both the collision between the north and south china blocks in the mesozoic and the eastward upper mantle flow driven by the indo asian collision

subduction of the indian lithosphere beneath the tibetan plateau and burma abstract subduction of the indian lithosphere under eurasia plays an important role in the tectonic evolution of the tibetan plateau and surrounding regions to improve our knowledge of pertinent mantle structures through tomographic imaging we combine p wave arrival time data from temporary arrays in tibet and stations of the chinese seismograph network with reprocessed data from the international seismological centre the new images reveal considerable spatial variations in mantle structure along the collision zone and the horizontal distance over which presumed continental indian lithosphere slides northward beneath the plateau decreases from west where it underlies the himalayas and the entire plateau to east where no indication is found for present day underthrusting beyond the himalayan block and indus tsangpo suture p wavespeed appears low in the shallow mantle beneath much of central and eastern tibet these observations suggest that indian lithosphere underlies only the southwestern part of the plateau and that the central and northeastern part is underlain by lithosphere of asian origin the continental parts of the indian plate that currently underthrust or subduct below the plateau appear generally detached from the oceanic slab fragments that subducted longer ago and that have been detected deeper in the mantle except perhaps in the central part of the collision zone 80 90 e this suggests that ongoing india eurasia collision must increasingly be driven by other forces such as subduction of the indian plate further west and east e g beneath indochina our images reveal structures associated with eastward subduction along the burmese arc southeast of the eastern himalayan syntaxis the tengchong volcanic complex in southwest china is marked by slow wave propagation to 150 km depth suggesting a causal link to subduction along the burmese arc the low velocities beneath the red river fault region extend to greater depth and may be related to upper mantle processes further southeast

diurnal winds in the himalayan kali gandaki valley part i observations abstract the diurnal wind system of the kali gandaki valley in nepal was explored in september and october 1998 in a field campaign using pilot balloons as the main observational tool this valley connects the plateau of tibet with the indian plains the river crosses the himalayas forming the deepest valley on earth intense upvalley winds blow up this valley during the day observations were made along the river at various spots selected between the exit point from the himalayas and the source close to the plateau of tibet the strongest upvalley winds were found between marpha and chuksang with typical speeds of 15 20 m s 1 the upvalley wind sets in first at the ground but an upvalley wind layer of 1000 2000 m depth forms quickly after the onset this deep inflow layer persists up to lo manthang a town located a few kilometers south of the plateau of tibet deceleration in the late afternoon and evening also appears to commence near the ground weak drainage flow forms late in the night the causes o

three dimensional velocity structure of crust and upper mantle in southwestern china and its tectonic implications using p and s arrival times from 4625 local and regional earthquakes recorded at 174 seismic stations and associated geophysical investigations this paper presents a three dimensional crustal and upper mantle velocity structure of southwestern china 21 34 n 97 105 e southwestern china lies in the transition zone between the uplifted tibetan plateau to the west and the yangtze continental platform to the east in the upper crust a positive velocity anomaly exists in the sichuan basin whereas a large scale negative velocity anomaly exists in the western sichuan plateau consistent with the upper crustal structure under the southern tibetan plateau the boundary between these two anomaly zones is the longmen shan fault the negative velocity anomalies at 50 km depth in the tengchong volcanic area and the panxi tectonic zone appear to be associated with temperature and composition variations in the upper mantle the red river fault is the boundary between the positive and negative velocity anomalies at 50 km depth the overall features of the crustal and the upper mantle structures in southwestern china are a low average velocity large crustal thickness variations the existence of a high conductivity layer in the crust or and upper mantle and a high heat flow value all these features are closely related to the collision between the indian and the asian plates

pn anisotropic tomography and dynamics under eastern tibetan plateau we present a new anisotropic tomographic model of the uppermost mantle around eastern tibet using pn traveltime data from a newly deployed temporary seismic array and recent observation bulletins of chinese provincial networks our results are generally consistent with previous results but provide new insights into the dynamics of tibetan plateau prominent high velocity high v anomalies are visible under alashan block and qaidam and sichuan basins which clearly outline their tectonic margins a distinct high v zone representing the double sided subduction of indo eurasian plates is imaged from lhasa block to the south of qaidam basin a pronounced low velocity low v zone is observed from songpan ganzi block to southern chuan dian diamond block suggesting the existence of hot material upwelling there crustal strong earthquakes frequently occurred around high v anomalies or transition zones from high v to low v anomalies suggesting that these earthquakes could be related to lateral heterogeneities in the mantle the pn fast direction approximately rotates around eastern himalayan syntaxis and it is tangential to the margins of sichuan basin suggesting that the mantle material flow of tibetan plateau may have affected east china in the yunnan region to the south of 26 n the pn fast direction is different from sks splitting results indicating that the mantle lithosphere could be mechanically decoupled at certain depth below the uppermost mantle which might be attributable to the subduction of indian or burma slab although the correlation between anisotropy and velocity is complicated anisotropy strength could be associated with the pattern of velocity anomalies in the region

evidence of active mantle flow beneath south china 1 the india eurasia collision is responsible for producing the himalayan mountains and tibetan plateau and has been hypothesized to have significant far field influences including driving the baikal rift and the eastward extrusion of south china however quantification of lithospheric buoyancy forces and integrated effect of tractions acting at base of the lithosphere are unable to explain the observed surface motions within south china we present 198 new sks shear wave splitting observations beneath south china and invert these data along with published gps data to solve for the subasthenospheric flow field beneath south china to assess the role of small scale convection here we find a 15 20 mm yr southwestward directed mantle flow toward the burma slab this flow is consistent with the mantle response of slab retreat over the past 25 ma and counter flow due to subduction of burma sunda slabs demonstrating the importance of localized mantle convection on present day plate motions

comparison of shear wave splitting and finite strain from the india asia collision zone we investigate whether observations of shear wave splitting in asia may be related to distributed strain of the lithosphere or to strain caused by motion of the asian lithosphere over the deep mantle the hs2 nuvel angular velocity for eurasia and that of minster and jordan 1978 predict absolute plate motions that are a factor of at least 4 slower than the nne relative motion of india with respect to stable eurasia shear of the upper mantle beneath the south asian lithosphere is therefore likely to be dominated by the velocity field caused by the internal deformation of asia rather than by plate motion if that shear causes sks splitting fast polarization directions in tectonic asia should be aligned approximately in the direction of the motion of asian lithosphere relative to the deep mantle which for southern asia is nne in general however fast polarization directions are not aligned with this direction and in some cases are nearly orthogonal to it we calculate finite deformation of the asian lithosphere in response to its indentation by india and from the deformation predict the orientations of fast polarization directions aligned with the principal axes of elongation these predictions agree with the gross features of the anisotropy in tibet and the tien shan orientations of anisotropy in the baikal region of siberia are consistent with the orientations of the local principal axes of active strain orientations of anisotropy in mongolia however are inconsistent both with the azimuth of the velocity of this region with respect to the deeper mantle and with the orientation of the calculated maximum elongation of the region we conclude that sks splitting observations from asia agree better with those predicted by calculations of internal lithospheric deformation due to indentation by india than with strain predicted from motion of the lithosphere over the deeper mantle these results place constraints on the depth range where shear wave splitting occurs which we propose lies predominantly within the lithospheric mantle

the african landscape through space and time it is generally accepted that cenozoic epeirogeny of the african continent is moderated by convective circulation of the mantle nevertheless the spatial and temporal evolution of africa s basin and swell physiography is not well known here we show how continental drainage networks can be used to place broad constraints on the pattern of uplift through space and time first we assemble an inventory of 710 longitudinal river profiles that includes major tributaries of the 10 largest catchments river profiles have been jointly inverted to determine the pattern of uplift rate as a function of space and time our inverse model assumes that shapes of river profiles are controlled by uplift rate history and modulated by erosional processes which can be calibrated using independent geologic evidence e g marine terraces volcanism and thermochronologic data our results suggest that modern african topography started to develop 30 myr ago when volcanic swells appeared in north and east africa during the last 15 20 myr subequatorial africa was rapidly elevated culminating in the appearance of three large swells that straddle southern and western coasts our results enable patterns of sedimentary flux at major deltas to be predicted and tested we suggest that the evolution of drainage networks is dominated by rapid upstream advection of signals produced by a changing pattern of regional uplift an important corollary is that with careful independent calibration these networks might act as useful tape recorders of otherwise inaccessible mantle processes finally we note that there are substantial discrepancies between our results and published dynamic topographic predictions

diurnal winds in the himalayan kali gandaki valley part ii modeling abstract the penn state ncar mesoscale model mm5 is used to simulate and better understand the wind observations in the kali gandaki valley reported in the first part of this paper the kali gandaki river originates in nepal near tibet flows southward through the mustang basin crosses the himalayas in a gorge and descends to the lowlands of nepal extremely strong diurnal upvalley flow in the gorge and the basin alternates with rather weak drainage flow in the night as proposed in part i the mustang basin and the tibetan plateau can be considered as an elevated heat source driving the upvalley flow during the day however the extreme strength of the diurnal upvalley winds and the order of magnitude asymmetry between day and night cannot be explained with a simple plateau circulation theory the model is successful in simulating almost all aspects of the observations the simulations strongly suggest that the observed acceleration of the upvalley winds near the entrance to the mustang basin is linked

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

accretionary tectonics of burma and the three dimensional geometry of the burma subduction zone the geometry of the burma wadati benioff zone wbz has been determined by fitting a trend surface parameterized with eight effective degrees of freedom to 184 well located hypocenters the dip of this surface which passes through the middle of the wbz varies from about 50 in the north near the eastern himalayan syntaxis to about 30 in the bay of bengal area the eastern edge of the indo burman ranges closely follows the map projection of the 60 km depth contour of the wbz the curvature of the indo burman ranges is controlled by the geometry of the interface between the more steeply dipping part of the indian plate and the leading edge of the overriding burma platelet shallow earthquakes beneath the indo burman ranges are primarily confined to the underthrusting indian plate their focal mechanisms indicate strike slip faulting and north south shortening parallel to the eastern margin of the indian plate

gap winds in a fjord part ii hydraulic analog abstract a simple shallow water model of gap wind in a channel that is based upon hydraulic theory is presented and compared with observations and output from a 3d mesoscale numerical model the model is found to be successful in simulating gap winds the speed and depth of gap wind flow is strongly controlled by topography horizontal or vertical channel contractions can act to force strong shallow supercritical flow downwind and light deep subcritical flow upwind force balance analysis of the hydraulic model output confirms mesoscale model results and indicates that the prime force balance in gap wind is between external pressure gradient and friction for supercritical flow and between external pressure gradient and height pressure gradient for subcritical flow this force balance changes near channel controls when the balance is between advection and height pressure gradient sensitivity analyses show positive sensitivity of gap wind speed to changes in discharge and external pressure gradient nega

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

fault slip and gps velocities across the shan plateau define a curved southwestward crustal motion around the eastern himalayan syntaxis characterizing the 700 km wide system of active faults on the shan plateau southeast of the eastern himalayan syntaxis ehs is critical to understanding the geodynamics and seismic hazard of the large region that straddles neighboring china myanmar thailand laos and vietnam here we evaluate the fault styles and slip rates over multi timescales reanalyze previously published short term gps velocities and evaluate slip rate gradients to interpret the regional kinematics and geodynamics that drive the crustal motion relative to the sunda plate gps velocities across the shan plateau define a broad arcuate tongue like crustal motion with a progressively northwestward increase in sinistral shear over a distance of 700 km followed by a decrease over the final 100 km to the syntaxis the cumulative gps slip rate across the entire sinistral slip fault system on the shan plateau is 12 mm yr our observations of the fault geometry slip rates and arcuate southwesterly directed tongue like patterns of gps velocities across the region suggest that the fault kinematics is characterized by a regional southwestward distributed shear across the shan plateau compared to more block like rotation and indentation north of the red river fault the fault geometry kinematics and regional gps velocities are difficulty to reconcile with regional bookshelf faulting between the red river and sagaing faults or localized lower crustal channel flows beneath this region the crustal motion and fault kinematics can be driven by a combination of basal traction of a clockwise southwestward asthenospheric flow around the ehs and gravitation or shear driven indentation from north of the shan plateau

on the problem of violent valley winds abstract observational results of a one month mesoscale experiment in a valley are used to emphasize the prominent part played by an inversion layer in air flow dynamics a model based on the analogy between shallow water flow and air flow beneath an inversion is applied using the actual topography of the valley theoretical arguments for the appearance of violent winds are presented and the model results are compared with the in situ measurements the qualitatively good agreement for the wind speed and the height of the inversion layer allows for further confidence in the model predictions for example the existence of a sharp transition zone or jump with a strong horizontal gradient of wind speed and temperature present in the model is consistent with aircraft measurements

diurnal winds in the himalayan kali gandaki valley part iii remotely piloted aircraft soundings abstract in 1998 a field campaign has been conducted in the north south oriented kali gandaki valley in nepal to explore the structure of its extreme valley wind system piloted ballon pibal observations were made to map the strong upvalley winds as well as the weak nocturnal flows part i the stratification of the valley atmosphere was not explored in part ii of this multipart paper numerical simulations are presented that successfully simulate most of the wind observations moreover the model results suggest that the vigorous upvalley winds can be seen as supercritical flow induced by contractions of the valley here the results of a further campaign are reported where remotely piloted airplanes were used to obtain vertical profiles of temperature and humidity up to heights of 2000 m above the ground such profiles are needed for an understanding of the flow dynamics in the valley and for a validation of the model results this technique is novel in some respects and turned out to be highly rel

vertical coherence of deformation in lithosphere in the eastern himalayan syntaxis using gps quaternary fault slip rates and shear wave splitting data we present 59 new sks skks and combine them with 69 previously published data to infer the mantle deformation field in se tibet the dense set of anisotropy measurements in the eastern himalayan syntaxis ehs are oriented along a ne sw azimuth and rotate clockwise in the surround regions we use gps measurements and geologic data to determine a continuous surface deformation field that is then used to predict shear wave spitting directions at each station comparison of splitting observations with predictions yields an average misfit of 11 7 illustrating that deformation is vertically coherent consistent with previous studies within the central ehs in areas directly surrounding the namche barwa metamorphic massif the average misfit of 11 stations increases to 60 8 and vertical coherence is no longer present the complexity of the mantle anisotropy and surface observations argues for local alteration of the strain fields here associated with recent rapid exhumation of the indian crust

seismic anisotropy beneath tibet evidence for eastward extrusion of the tibetan lithosphere strong seismic anisotropy beneath tibet has recently been reported from the study of sks shear wave splitting the fast split waves are generally polarized in an easterly direction close to the present day direction of motion of the tibetan crust relative to stable eurasia as deduced from holocene slip rates on the major active faults in and around tibet this correlation may be taken to suggest that the whole tibetan lithosphere is being extruded in front of indenting india and that the anisotropic layer is the deforming asthenosphere that accommodates the motion of the tibetan lithosphere relative to the fixed mantle at depth uncertainties about this motion are at present too large to bring unambiguous support to that view assuming that this view is correct however a simple forward model is used to compute theoretical delay times as a function of the thickness of the anisotropic layer the observed delay times would require a 50 100 km thick anisotropic layer beneath south central tibet and an over 200 km thick layer beneath north central tibet where particularly hot asthenosphere has been inferred this study suggests that the asthenospheric anisotropy due to present absolute block motion might be dominant under actively deforming continents

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

the mantle flow field beneath western north america although motions at the surface of tectonic plates are well determined the accompanying horizontal mantle flow is not we have combined observations of surface deformation and upper mantle seismic anisotropy to estimate this flow field for western north america we find that the mantle velocity is 5 5 1 5 centimeters per year due east in a hot spot reference frame nearly opposite to the direction of north american plate motion west southwest the flow is only weakly coupled to the motion of the surface plate producing a small drag force this flow field is probably due to heterogeneity in mantle density associated with the former farallon oceanic plate beneath north america

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