dynamic topography produced by lower crustal flow against rheological strength heterogeneities bordering the tibetan plateau summary dynamic stresses developed in the deep crust as a consequence of flow of weak lower crust may explain anomalously high topography and extensional structures localized along orogenic plateau margins with lubrication equations commonly used to describe viscous flow in a thin gap geometry we model dynamic stresses associated with the obstruction of lower crustal channel flow due to rheological heterogeneity dynamic stresses depend on the mean velocity u viscosity and channel thickness h uniquely through the term u h 2 these stresses are then applied to the base of an elastic upper crust and the deflection of the elastic layer is computed to yield the predicted dynamic topography we compare model calculations with observed topography of the eastern tibetan plateau margin where we interpret channel flow of the deep crust to be inhibited by the rigid sichuan basin model results suggest that as much 1500 m of dynamic topography across a region of several tens to a hundred kilometres wide may be produced for lower crustal material with a viscosity of 2 10 18 p asfl owing in a 15 km thick channel around a rigid cylindrical block at an average rate of 80 mm yr 1

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

plate tectonic history basin development and petroleum source rock deposition onshore china abstract china comprises a mosaic of distinct continental fragments separated by fold belts these fold belts are suture zones resulting from the accretion of various fragments formerly separated by intervening areas of oceanic crust the major sedimentary basins onshore china can be classified into four groups those in western china are flexural developing as a result of north south compression in contrast those in the east are extensional and related to development of the pacific oceanic margin in central china basins have a more problematic origin those of north central china ordos sichuan are flexural basins controlled by eastward directed thrusting along their western margin in contrast basins further south chuxiong shiwandashan are predominantly extensional and related to major strike slip movements by synthesizing basin stratigraphies across china in tectonostratigraphic terms and in particular comparing the nature and timing of unconformities it is possible to formulate a coherent model for the palaeoreconstruction of china we identify five major tectonostratigraphic breaks which equate with the collision of the following continental fragments tarim north china carboniferous permian south china block permian triassic qiantang late triassic early jurassic lhasa block late jurassic early cretaceous and india early tertiary prior to permian times the southern margin of eurasia ran approximately along the northern border of modern china the late carboniferous collision of tarim north china with eurasia resulted in the development of a flexural basin junggar and deposition of non marine clastics to the south of the suture shallow marine deposition continued in the late permian early triassic the progressive collision from east to west of the south china block with the north china block resulted in a change to fluvial lacutrine sedimentation across the entire north china tarim block open marine carbonate deposition in the north of the south china block passed southward into a deeper marine clastic sequence deposited in a backarc basin further south a subduction zone existed along the southeastern margin of the south china block in western china northward subduction throughout the triassic resulted in the development of the songban ganzi accretionary prism with retroarc thrusting resulting in flexure and the first development of the tarim basin oblique collision of the qiantang block in the late triassic along the east of the south china block resulted in east west directed thrusting which initiated the suchuan and ordos basins continued strike slip deformation along the south western margin of the south china block resulted in the development of basins with a significant extensional component such as chuxiong the collision of the qiantang block with the southern edge of the tarim basin early jurassic resulted in a renewed clastic influx in both the tarim and junggar basins along the eastern pacific margin a compressional arc and retroarc basin in the south passed northwards into an extensional arc system subduction rollback of the extensional arc initiated rifting in the late jurassic in the eren and songliao basins the late jurassic early cretaceous collision of the lhasa block in the west rejuvenated the thrust systems bordering the western basin and resulted in a renewed clastic influx in the southeast the compressional arc phase culminated in widespread thrusting and folding of early cretaceous age in the northeast extension continued with the progressive migration of the rift system southward with time the arrival of the indian block in the early tertiary rejuvenated the bounding thrust belts of all the western basins in the east the change in convergence of the pacific plate to a more westerly direction is marked by extension and widespread rifting along the entire length of eastern china throughout most of china mesozoic and cenozoic deposition occurred in predominantly non marine environments source rocks in such settings comprise principally mudstones deposited in lakes organic rich mudstones these can accumulate in both deep and shallow lakes in order to accumulate substantial volumes the lake must be significant in space and time in china lacustrine orms occur in both rift and flexural basins lacustrine orms deposited under humid climatic conditions are restricted to the period of maximum tectonic subsidence in the flexural basins of western china source rock deposition follows basin initiation by 20 30 ma in the extensional basins of eastern china source rock deposition takes place 5 15 ma after basin initiation by contrast semi arid and arid climate lacustrine orms whilst being best developed during the period of maximum tectonic subsidence occur at all stages in the basin history

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

evaluation of permeability and non darcy flow in vuggy macroporous limestone aquifer samples with lattice boltzmann methods 1 lattice boltzmann flow simulations provide a physics based means of estimating intrinsic permeability from pore structure and accounting for inertial flow that leads to departures from darcy s law simulations were used to compute intrinsic permeability where standard measurement methods may fail and to provide better understanding of departures from darcy s law under field conditions simulations also investigated resolution issues computed tomography ct images were acquired at 0 8 mm interscan spacing for seven samples characterized by centimeter scale biogenic vuggy macroporosity from the extremely transmissive sole source carbonate karst biscayne aquifer in southeastern florida samples were as large as 0 3 m in length 7 9 cm scale length subsamples were used for lattice boltzmann computations macroporosity of the subsamples was as high as 81 matrix porosity was ignored in the simulations non darcy behavior led to a twofold reduction in apparent hydraulic conductivity as an applied hydraulic gradient increased to levels observed at regional scale within the biscayne aquifer larger reductions are expected under higher gradients near wells and canals thus inertial flows and departures from darcy s law may occur under field conditions changes in apparent hydraulic conductivity with changes in head gradient computed with the lattice boltzmann model closely fit the darcy forchheimer equation allowing estimation of the forchheimer parameter ct scan resolution appeared adequate to capture intrinsic permeability however departures from darcy behavior were less detectable as resolution coarsened

density structure and isostatic state of the crust in the longmenshan and adjacent areas abstract in this study we ve analyzed some new gravity data obtained from two profiles across the longmenshan lms fault system in order to determine the density structure and isostatic state of the crust beneath the songpan ganzi block sg and the sichuan basin sc according to our research bouguer gravity anomalies along the two profiles range from 450 mgal in the sg to 80 mgal in the sc and density structural models show obvious differences between the sg and the sc compared with the sc the sg has lower average density in the crust and it has a low density layer with thickness of 5 km in the middle crust in the sc area however the crustal density increases gradually with depth and there exists no low density layer in the crust the moho interface decreases from 60 km beneath the sg in east tibet to 40 km beneath the sc in the lms region the topography bouguer gravity and moho depth change abruptly in a narrow tectonic belt the crustal isostatic state which is derived from the airy compensation model based on the topography gravity and deep seismic sounding data indicates that the sg on the northwest side of the lms is not in an isostatic state compared with the isostatic crust beneath the sc the isostatic equilibrium of the crust under the northeast section of the lms is not achieved as that in the southwest section of the lms based on the research we infer that the non homogenous collision between the tibetan plateau and the yangtze craton has caused the strong variation of the moho depth and the special isostatic state of the crust in this area

shear wave splitting in the upper mantle wedge above the tonga subduction zone summary intermediate period seismograms of shear waves from deep earthquakes are examined for shear wave splitting diagnostic of seismic velocity anisotropy in the upper mantle wedge between the subducting tonga slab and stations in the fiji islands earthquakes with near vertical incidence angles are selected to minimize complications of converted phases and the free surface horizontal component seismograms for 25 deep earthquakes recorded at station lak on the lau ridge show splitting of up to 0 9 s with an average of 0 45 0 25 s between fast and slow s wave components anisotropy in the upper mantle is the most likely cause of the splitting the ese orientation of the fast velocity direction at station lak corresponds to an axis of symmetry of the anisotropic mantle and is approximately perpendicular to magnetic lineations in the lau basin and the northeastern portion of the south iiji basin the observations are consistent with a model of flow alignment of olivine 100 and orthopyroxene 001 axes parallel to the direction of basin extension seismograms for eight deep earthquakes at station sva on viti levu show up to 1 0 s splitting with an average of 0 65 0 29 s and a weak nne preferred orientation of fast velocity polarizations of shear wave particle motions corrected for anisotropy are within 36 5 26 of polarixations predicted for nine moment tensor solutions the magnitude of shearwave splitting and estimates of the thickness of the lithosphere in the fiji region imply that the anisotropy cannot be contained in the lithosphere alone either the anisotropy is distributed in both the lithosphere and asthenosphere or the asthenosphere alone is anisotropic

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

structure of the crust beneath the southeastern tibetan plateau from teleseismic receiver functions abstract southeastern tibet marks the site of presumed clockwise rotation of the crust due to the india eurasian collision and abutment against the stable sichuan basin and south china block knowing the structure of the crust is a key to better understanding of crustal deformation and seismicity in this region here we analyze recordings of teleseismic earthquakes from 25 temporary broadband seismic stations and one permanent station using the receiver function method we find that the crustal thickness decreases gradually from the tibetan plateau proper to the sichuan basin and yangtze platform but that significant intra crustal heterogeneity exists on shorter lateral scales 0 30 are detected beneath the southeastern margin of tibet but in the sichuan basin and southeastern yangtze platform the values are close to the global average combined with high regional heat flow and independent evidence for mid crustal layers of high electric conductivity the large intra crustal s wave velocity reduction 12 19 and the intermediate to high average crustal poisson s ratios are consistent with partial melt in the crust beneath parts of southeastern tibet these results could be used in support of deformation models involving intra crustal flow with the caveat that significant lateral variation in location and strength of this flow may occur

crustal structure and deformation of the se tibetan plateau revealed by receiver function data we analyze a large amount of receiver function data recorded by regional seismic networks of the china earthquake administration to estimate crustal structure and deformation beneath the southeast margin of the tibetan plateau and its surrounding areas we develop a comprehensive analysis method that facilitates robust extraction of azimuthal seismic anisotropy from receiver function data the method includes an estimate of fast polarization direction and splitting time by a joint analysis of radial and transverse receiver function data and an evaluation of measurement reliability by statistical and harmonic analysis we find significant seismic anisotropy with a splitting time of 0 5 0 9 s beneath the se margin of the tibetan plateau both the splitting time and fast polarization direction are comparable to those estimated from sks skks data suggesting that crustal anisotropy is the main cause of shear wave splitting of the sks skks wave this also suggests that deformation in the upper mantle is either weak or predominantly vertical and is obviously different from the one in the crust a vertical flow in the upper mantle combined with the observation of a thin lithosphere beneath the study area leads to the inference that part of the mantle lithosphere may have been removed and is descending into deep mantle stations located in the surrounding areas on the other hand exhibit very little to no crustal anisotropy the estimated moho depth and vp vs ratio also show a distinct difference between the se tibetan plateau and the surrounding regions stations on the tibetan plateau have a vp vs ratio of 1 79 which is substantially higher than those measured at the yunnan guizhou yungui plateau 1 69 our observations here are consistent with the scenario that the se tibet has been built by lower crustal flow they also suggest that the mantle lithosphere beneath the margin may have been mechanically decoupled from the upper crust

crustal structure across longmenshan fault belt from passive source seismic profiling 1 we analyse receiver functions from 29 broad band seismographs along a 380 km profile across the longmenshan lms fault belt to determine crustal structure beneath the east tibetan margin and sichuan basin the moho deepens from about 50 km under songpan ganzi in east tibet to about 60 km beneath the lms and then shallows to about 35 km under the western sichuan basin the average crustal vp vs ratios vary in the range 1 75 1 88 under songpan ganzi in east tibet 1 8 2 0 under the lms and decrease systematically across the nw part of the sichuan basin to less than 1 70 a negative phase arrival above the moho under songpan ganzi and sichuan basin is interpreted as a ps conversion from the top of a low velocity layer in the lower crust the very high crustal vp vs ratio and negative polarity ps conversion at the top of lower crust in east tibet are inferred to be seismic signatures of a low viscosity channel in the eastern margin of the tibetan plateau the lateral variation of moho topography crustal vp vs ratio and negative polarity ps conversion at the top of the lower crust along the profile seem consistent with a model of lower crust flow or tectonic escape

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

estimating crustal seismic anisotropy with a joint analysis of radial and transverse receiver function data summary we developed an integrated technique for estimating crustal anisotropy with a horizontal axis using radial r and transverse t receiver functions the technique includes computing three individual and one joint objective function jof and a reliability analysis of the estimated anisotropy the individual objective functions iofs are designed to 1 maximize the peak energy of the stacked r receiver function after a cosine moveout correction in the ps arrival time 2 to maximize the correlation of the radial receiver functions after a full correction of anisotropy or 3 to minimize the total energy of transverse receiver functions stacked after a removal of crustal anisotropy the jof was computed by a weighted average of the three iofs and the reliability analysis uses the principle that stacking coherent signals can lead to an increase of signal to noise ratio we applied the technique to synthetic receiver functions generated with 30 60 per cent white noise from a variety of anisotropic and heterogeneous models the synthetic tests indicate that the proposed technique has good capability to recover the input models despite the presence of random and other coherent noises such as those caused by inhomogeneous structures in the data the technique can always provide accurate estimates of crustal anisotropy we applied the technique to two permanent seismic stations in western china and found significant crustal anisotropy beneath one station located at the northern edge of the tibetan plateau the observed fast polarization direction at this station follows the direction of the maximum horizontal tensile stress suggesting that the observed seismic anisotropy is likely caused by mineral alignment in the lower crust the station situated in the sichuan basin on the other hand shows little to no seismic anisotropy which may suggest that the crust beneath the basin is nearly rigid with very little deformation the developed technique can be applied to any broadband seismic stations that have a good backazimuthal coverage of teleseismic events

fine structure of the lowermost crust beneath the kaapvaal craton and its implications for crustal formation and evolution high quality data from a dense seismic array covering an area of approximately 60 40 km2 are used to obtain tight quantitative estimates of the fine scale velocity and density structure of the lowermost crust and the crust mantle boundary moho beneath the kaapvaal craton in the vicinity of kimberley south africa results based on a modified receiver function waveform analysis of moho conversions and crustal reverberation phases show that the crust beneath the array is thin 35 4 km with an average poisson s ratio of 0 254 the minimum s wave velocity contrast across the moho is 17 3 while the contrast in density is 15 4 the density contrast across the moho is particularly diagnostic for an assumed uppermost mantle density beneath kimberley of 3 3 gm cc as determined from mantle xenoliths the density of the lowermost crust is 2 86 gm cc indicating rocks of felsic to intermediate composition analysis of waveform broadening of the crustal reverberation phases relative to that of the direct p wave shows the thickness of the moho transition zone to be less than 0 5 km and the maximum variation in crustal thickness over the region of the array to be less than 1 km the flat and almost perfectly sharp moho together with the absence of a mafic lower crust suggests large scale crustal reworking in the period between crustal formation and the time of cratonic stabilization

structural heterogeneity of the longmenshan fault zone and the mechanism of the 2008 wenchuan earthquake ms 8 0 we determined detailed 3 d images of p and s wave velocity vp vs and poisson s ratio in and around the longmenshan lms fault zone by using a large number of p and s wave arrival times from the aftershocks of the 2008 wenchuan earthquake ms 8 0 and other local events our results show that the structure of the lms fault zone north of the wenchuan main shock is very different from that south of the main shock the southern section of the lms fault zone contains a broad zone with low vp low vs and high anomalies while the northern segment exhibits more scattered heterogeneities corresponding to most of the aftershocks which occurred there a prominent low vp low vs and high anomaly exists directly beneath the wenchuan main shock hypocenter suggesting that in addition to compositional variations fluid filled fractured rock matrices may exist in the lms fault zone which may have influenced the generation of the large wenchuan earthquake our tomographic results provide sound seismic evidence for the hypothesis that an upward intrusion of lower crustal flow occurred along the lms fault zone in addition most small earthquakes before the 2008 wenchuan main shock occurred around the guanxian jiangyou fault while the wenchuan aftershocks were mainly concentrated on the yingxiu beichuan fault suggesting that the rupture process of the wenchuan earthquake may belong to an out of sequence thrusting event a suggestion which is in good agreement with the results from geological surveys and also quite similar to the rupture processes of the 1999 chi chi earthquake m 7 5 and the 2005 kashmir earthquake m 7 6 a few aftershocks occurred close to the blind guangyuan dayi fault in the sichuan basin suggesting that this blind fault was also ruptured by the wenchuan earthquake consistent with geological surveys

contemporary crustal deformation around the southeast borderland of the tibetan plateau 1 we derive a detailed horizontal velocity field for the southeast borderland of the tibetan plateau using gps data collected from the crustal motion observation network of china between 1998 and 2004 our results reveal a complex deformation field that indicates that the crust is fragmented into tectonic blocks of various sizes separated by strike slip and transtensional faults most notably the regional deformation includes 10 11 mm yr left slip across the xianshuihe fault 7 mm yr left slip across the anninghe zemuhe xiaojiang fault zone 2 mm yr right slip across a shear zone trending northwest near the southern segment of the lancang river fault and 3 mm yr left slip across the lijiang fault deformation along the southern segment of the red river fault appears not significant at present time the region south and west of the xianshuihexiaojiang fault system whose eastward motion is resisted by the stable south china block to the east turns from eastward to southward motion with respect to south china resulting in clockwise rotation of its internal subblocks active deformation is detected across two previously unknown deformation zones one is located 150 km northwest of and in parallel with the longmenshan fault with 4 6 mm yr right slip and another is continued south southwestward from the xiaojiang fault abutting the red river fault with 7 mm yr left slip while both of these zones are seismically active the exact locations of faults responsible for such deformation are yet to be mapped by field geology comparing our gps results with predictions of various models proposed for tibetan plateau deformation we find that the relatively small sizes of the inferred microblocks and their rotation pattern lend support to a model with a mechanically weak lower crust experiencing distributed deformation underlying a stronger highly fragmented upper crust

evaluation of anisotropy by shear wave splitting the polarizations of three component shear wavetrains carry unique information about the internal structure of the rock through which they pass specifically commonly observed shear wave splitting may contain information about the orientation of crack distributions this information cannot usually be recovered from shear waves recorded at the free surface because of interference with the interaction of the shear wave with the surface even for nearly vertical incidence however shear wave splitting in synthetic three component vertical seismic profiles in some cases may be interpreted directly in terms of the direction of strike of vertical cracks because the human eye is not skilled at identifying the phase relationships between three component seismograms played out conventionally as parallel time series the polarizations are displayed in orthogonal sections of the particle displacements to facilitate recognition and evaluation of the shear wave splitting estimating the orientations of cracks an

co seismic ruptures of the 12 may 2008 ms 8 0 wenchuan earthquake sichuan east west crustal shortening on oblique parallel thrusts along the eastern edge of tibet abstract the m s 8 0 wenchuan earthquake which devastated the mountainous western rim of the sichuan basin in central china produced a surface rupture over 200 km long with oblique thrust dextral slip and maximum scarp heights of 10 m it thus ranks as one of the world s largest continental mega thrust events in the last 150 yrs field investigation shows clear surface breaks along two of the main branches of the ne trending longmen shan thrust fault system the principal rupture on the nw dipping beichuan fault displays nearly equal amounts of thrust and right lateral slip basin ward of this rupture another continuous surface break is observed for over 70 km on the parallel more shallowly nw dipping pengguan fault slip on this latter fault was pure thrusting with a maximum scarp height of 3 5 m this is one of the very few reported instances of crustal scale co seismic slip partitioning on parallel thrusts this out of sequence event with distributed surface breaks on crustal mega thrusts highlights regional ew directed present day crustal shortening oblique to the longmen shan margin of tibet the long rupture and large offsets with strong horizontal shortening that characterize the wenchuan earthquake herald a re evaluation of tectonic models anticipating little or no active shortening of the upper crust along this edge of the plateau and require a re assessment of seismic hazard along potentially under rated active faults across the densely populated western sichuan basin and mountains

s receiver function analysis of the crustal and lithospheric structures beneath eastern tibet abstract s receiver functions from 51 permanent broad band stations are used to investigate the structures of the crust and lithosphere beneath eastern tibet and sichuan basin by stacking the s receiver functions the sp converted phases from the moho and the lithosphere asthenosphere boundary lab are clearly visible the moho depth ranges from 40 to 50 km beneath sichuan basin and from 55 to 60 km on the west side of the longmenshan lms fault the largest depth is 70 km occurring at the northern sichuan yunnan terrane the lab depth beneath sichuan basin ranges from 130 to 170 km while that on the west side of lms ranges from 100 to 120 km from the east to the west across lms there is a sudden rise of the lab followed by a gradual dip western sichuan basin is found to be surrounded by thinner lithosphere extending northeastward along the lms and southward to indochina the thinner lithosphere may indicate that the lower part of eastern tibet lithosphere is heated and delaminated by the underlying hot asthenospheric flow thus reflecting the size and direction of the flow it appears that asthenospheric flow from tibet is resisted by the cold rigid sichuan basin thus branching into northeastward and southeastward parts the complicated topography of the lithosphere likely results from this dynamic process

shear wave splitting and implications for mantle flow beneath the melt region of the east pacific rise shear wave splitting across the fast spreading east pacific rise has been measured from records of sks and skks phases on the ocean bottom seismometers of the mantle electromagnetic and tomography melt experiment the direction of fast shear wave polarization is aligned parallel to the spreading direction delay times between fast and slow shear waves are asymmetric across the rise and off axis values on the pacific plate are twice those on the nazca plate splitting on the pacific plate may reflect anisotropy associated with spreading induced flow above a depth of about 100 km as well as a deeper contribution from warm asthenospheric return flow from the pacific superswell region

receiver functions of seismic waves in layered anisotropic media application to the estimate of seismic anisotropy abstract we investigate the effect of seismic anisotropy on p wave receiver functions calculating synthetic seismograms for p wave incidence on multilayered anisotropic structure with hexagonal symmetry the main characteristics of the receiver functions affected by the anisotropy are summarized as 1 appearance of seismic energy on radial and transverse receiver functions 2 systematic change of p to s ps converted waveforms on receiver functions as ray back azimuth increases and 3 reversal of the ps phase polarity on the radial receiver function in a range of the back azimuth another important influence is shear wave splitting of the ps converted waves and other later phases reverberated as s wave by numerical experiments using synthetic receiver functions we demonstrate that the waveform cross correlation analysis is applicable to splitting ps phases on receiver functions to estimate the seismic anisotropy of layer structure advantages to utilizing the ps phases are 1 they appear more clearly on receiver functions than on seismograms and 2 they inform us about what place along the seismic ray path is anisotropic real analysis of shear wave splitting is executed to the moho generated ps phases that are identified on receiver functions at six seismic stations in the chugoku district southwest japan the time lags between the two arrivals of the split ps phases are estimated at 0 2 0 7 sec and the polarization directions of the fast arrival components are from north south to northeast southwest this result is consistent with recent results of shear wave splitting measurements and the trend of linear epicenter distributions of crustal earthquakes and active fault strikes in the chugoku district

crustal anisotropy from moho converted ps wave splitting analysis and geodynamic implications beneath the eastern margin of tibet and surrounding regions abstract in addition to crustal thickening distinctly different mechanisms have been suggested to accommodate the huge convergences caused by the continental collision between india and eurasia as the transition zone between the two grand tectonic domains of asia the tethys and the pacific east tibet and its surrounding regions are the ideal places to study continental deformation pervasive rock deformation may produce anisotropy on the scale of seismic wavelengths thus seismic anisotropy provides insight into the deformation of the crust and mantle beneath tectonically active domains in this study we calculated receiver function pairs of radial and transverse components at 98 stations located in sichuan and yunnan provinces china we selected 7423 pairs with high signal to noise ratio snr and unambiguous moho converted ps phases pms to measure the pms splitting owing to the crustal anisotropy both the crustal thickness and the average crustal vp vs ratio were calculated simultaneously by the h k stacking method the geodynamic implications were also investigated in relation to surface geological features gps velocities absolute plate motion apm sks skks splitting and other seismological observations in addition to the fast polarization directions fpds of the crustal anisotropy we observed a conspicuous sharper clockwise rotation around the eastern himalayan syntaxis than was revealed by gps velocities the distributed fpds within and near the main active fault zones also favored the directions parallel to the faults this implied that the deformation of a continuous medium revealed by gps motions is a proxy for the deformation of the brittle shallow crust only while the main active faults and the deep crustal interiors both play important roles in the deep deformation our results suggest that the deformation between the crust and upper mantle within the northernmost section of the indochina block is decoupled due to the large difference in the directions between the observations related to the crust gps and crustal anisotropy and mantle apm and mantle anisotropy focusing on the transition zone between the plateau and the south china and indochina blocks we suggest that the motion of the central yunnan sub block is a southeastward extrusion by way of tectonic escape there is less deformation in the deep crust and the motion is controlled by the active boundary faults of the ailaoshan red river shear zone to the west and the xianshuihe xiaojiang fault to the east the lower crustal flow within the plateau southeastward reached the lijiang xiaojinhe fault but further south it was obstructed by the central yunnan sub block

reconciling lithospheric deformation and lower crustal flow beneath central tibet a viscous region with deformable boundaries is used to model simultaneous crustal flow and lithospheric coupling in northern tibet this model suggests that 1 the deformation of northern tibet is different from that of southern tibet because of structural differences between the regions 2 the viscosity contrast between the crust and the mantle lithosphere is relatively small beneath northern tibet and 3 crustal flow is compatible with crust mantle coupling under these conditions

lithospheric thickness and upper mantle deformation beneath the ne tibetan plateau inferred from s receiver functions and sks splitting measurements summary we present s receiver functions and sks splitting measurements from the china seismograph network located in the qinghai and gansu provinces teleseismic data are used to interpret the structure of the lithosphere asthenosphere boundary lab and upper mantle deformation beneath the northeastern tibetan plateau netp in regions north of the east kunlun fault based on our observations the lab lies at a depth of 125 135 km beneath the northeastern songpan ganzi block and the west qinling orogen between 145 and 175 km beneath the kunlun and qilian orogen and deepens below the qaidam basin 175 190 km ordos craton 170 km and alashan platform 200 km the netp is characterized by a nearly uniform fast nw se s wave direction these observations are different from those to the south of the kunlun fault where fast s directions are rotating clockwise from the inner plateau the change in fast directions across the kunlun fault implies a sudden variation of upper mantle deformation shear wave splitting delay times vary from 0 8 to 1 9 s data from beneath regions north of the kunlun ayimaqin suture showed that delay time was positively correlated with lithospheric thickness with an increase of 0 7 s per 100 km this indicates that the anisotropy may develop in the uppermost mantle such as the lithosphere beneath the netp

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

seismic anisotropy in eastern tibet from shear wave splitting reveals changes in lithospheric deformation abstract knowledge about seismic anisotropy can provide important insight into the deformation of the crust and upper mantle beneath tectonically active regions here we focus on the southeastern part of the tibetan plateau in sichuan and yunnan provinces sw china we measured shear wave splitting of core refracted phases sks and skks at a temporary array of 25 iris passcal stations we calculated splitting parameters using a multi channel and a single record cross correlation method multiple layers of anisotropy cannot be ruled out but are not required by the data a fresnel zone analysis suggests that the shallow mantle between 60 and 160 km depth is the most likely source of anisotropy the polarization directions reveal a pronounced transition from primarily north south in the north sichuan to mostly east west orientations in the south yunnan in the southern part of the study region that is south of 26 n the fast polarization directions do not correlate well with known surface features and geodetic estimates of the crustal displacement fields whereas gps campaigns provide evidence suggesting north south crustal flow across the red river fault the pattern of anisotropy argues against such flow in the upper mantle these observations support models that allow differential movement of upper crust relative to lithospheric mantle in the northern part of the study region the relationships are more ambiguous and coherent deformation of the crust and mantle lithosphere cannot be excluded the interpretation of the shear wave splitting results is non unique but we suggest that the observed n s transition reflects a fundamental change in deformation regime across our study region it may be related to lateral variations in lithospheric rheology or may mark a transition from the direct impact of the continental collision to dominance of the far field strain field associated with regional subduction processes understanding the nature of the lateral change in deformation regime may prove critical for our understanding the geotectonic evolution of eastern tibet

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

high resolution mantle tomography of china and surrounding regions 1 a high resolution p wave tomographic model of the crust and mantle down to 1100 km depth under china and surrounding regions is determined by using about one million arrival times of p pp pp and pcp waves from 19 361 earthquakes recorded by 1012 seismic stations the subducting pacific slab is imaged clearly as a high velocity zone from the oceanic trenches down to about 600 km depth and intermediate depth and deep earthquakes are located within the slab the pacific slab becomes stagnant in the mantle transition zone under east china the western edge of the stagnant slab is roughly coincident with a surface topographic boundary in east china the active changbai and wudalianchi intraplate volcanoes in northeast china are underlain by significant slow anomalies in the upper mantle above the stagnant pacific slab these results suggest that the active intraplate volcanoes in ne china are not hot spots but a kind of back arc volcano associated with the deep subduction of the pacific slab and its stagnancy in the transition zone under the mariana arc however the pacific slab penetrates directly down to the lower mantle the active tengchong volcano in southwest china is related to the eastward subduction of the burma microplate the subducting indian and philippine sea plates are also imaged clearly the indian plate has subducted down to 200 300 km depth under the tibetan plateau with a horizontal moving distance of about 500 km high velocity anomalies are revealed in the upper mantle under the tarim basin ordos and sichuan basin which are three stable blocks in china

an integrated global model of present day plate motions and plate boundary deformation summary in this paper we present a global model gsrm 1 of both horizontal velocities on the earth s surface and horizontal strain rates for almost all deforming plate boundary zones a model strain rate field is obtained jointly with a global velocity field in the process of solving for a global velocity gradient tensor field in our model we perform a least squares fit between model velocities and observed geodetic velocities as well as between model strain rates and observed geological strain rates model velocities and strain rates are interpolated over a spherical earth using bi cubic bessel splines we include 3000 geodetic velocities from 50 different mostly published studies geological strain rates are obtained for central asia only and they are inferred from quaternary fault slip rates for all areas where no geological information is included a priori constraints are placed on the style and direction but not magnitude of the model strain rate field these constraints are taken from a seismic strain rate field inferred from normalized focal mechanisms of shallow earthquakes we present a global solution of the second invariant of the model strain rate field as well as strain rate solutions for a few selected plate boundary zones generally the strain rate tensor field is consistent with geological and seismological data with the assumption of plate rigidity for all areas other than the plate boundary zones we also present relative angular velocities for most plate pairs we find that in general there is a good agreement between the present day plate motions we obtain and longterm plate motions but a few significant differences exist the rotation rates for the indian arabian and nubian plates relative to eurasia are 30 13 and 50 per cent slower than the nuvel1a estimate respectively and the rotation rate for the nazca plate relative to south america is 17 per cent slower on the other hand caribbean north america motion is 76 per cent faster than the nuvel 1a estimate while crustal blocks in the india eurasia collision zone move significantly and self consistently with respect to bounding plates only a very small motion is predicted between the nubian and somalian plates by integrating plate boundary zone deformation with the traditional modelling of angular velocities of rigid plates we have obtained a model that has already been proven valuable in for instance redefining a no net rotation model of surface motions and by confirming a global correlation between seismicity rates and tectonic moment rates along subduction zones and within zones of continental deformation

elastic thickness mechanical anisotropy and deformation of the southeastern tibetan plateau abstract the southeastern tibetan plateau se tibet located between the eastern himalayan syntaxis and sichuan basin is a heterogeneous structure combining different tectonic regimes it encompasses a gently sloping long wavelength margin and is characterized by intense deformation and active seismicity mechanical strength of the lithosphere is one of the parameters that controls deformation and evolution of any region in particular anisotropy of the strength is important for understanding the relation between tectonic processes and lithosphere deformation as a proxy the effective elastic thickness of the lithosphere t e can be used to study both spatial distribution and anisotropy of the lithospheric strength here we use the fan wavelet coherence method to determine variations of t e and its anisotropy based on the topography and bouguer gravity anomaly data the coherence analysis shows that low t e 5 30 km and high t e anisotropy are widespread in se tibet indicating a weak and highly mechanical anisotropic lithosphere this also implies that the style of lithospheric deformation is continuous and may be dominated by lower crustal flow compared with the crustal stress and lithospheric strain indicators we find that the t e anisotropy is predominantly governed by the present day tectonic stress in se tibet this may indicate a post collision compression tectonics within the plateau due to the india eurasia plate convergence and a resistance by the rigid sichuan basin however the anisotropy reflects extension stresses in approximately wnw ese direction in yunnan and northern indochina off the plateau which may relate to the burma plate subduction processes we infer that the present day tectonic stresses and mantle processes associated with the india eurasia collision continuously modify and rework the mechanical structure and both govern the lithospheric deformation in se tibet

oblique stepwise rise and growth of the tibet plateau two end member models of how the high elevations in tibet formed are i continuous thickening and widespread viscous flow of the crust and mantle of the entire plateau and ii time dependent localized shear between coherent lithospheric blocks recent studies of cenozoic deformation magmatism and seismic structure lend support to the latter since india collided with asia 55 million years ago the rise of the high tibetan plateau likely occurred in three main steps by successive growth and uplift of 300 to 500 kilometer wide crustal thrust wedges the crust thickened while the mantle decoupled beneath gently dipping shear zones did not sediment infilling bathtub like of dammed intermontane basins formed flat high plains at each step the existence of magmatic belts younging northward implies that slabs of asian mantle subducted one after another under ranges north of the himalayas subduction was oblique and accompanied by extrusion along the left lateral strike slip faults that slice tibet s east side these mechanisms akin to plate tectonics hidden by thickening crust with slip partitioning account for the dominant growth of the tibet plateau toward the east and northeast

azimuthal shear wave velocity anisotropy in the basin and range province using moho ps converted phases an eight station large aperture seismometer array was deployed in the basin and range province of west central nevada for the 1988 1989 program for array seismic studies of the continental lithosphere passcal passive source seismic experiment during the 10 months of data collection a total of 100 teleseismic events were recorded 30 source equalized p wave receiver functions show waveform variations in timing amplitude and polarity in the moho ps phase that is diagnostic of shear wave velocity anisotropy information regarding the crustal anisotropic component can be recovered through the ps phase since it is generated at the moho and is sensitive to shear wave splitting in the crust only we utilize single event stacked and array beam formed receiver functions to minimize the influence of structural variations and scattering on the observed converted phases analysis of shear wave splitting of the ps conversion from the crust mantle boundary indicates a fast azimuth of anisotropy oriented approximately nw se with observed time delays of 0 20 s this strongly suggests coherent crustal anisotropy with a fast direction perpendicular to the strike of fault block mountain ranges this fabric is most likely due to the preferred crystallographic orientation of seismically anisotropic minerals in the middle to lower crust rather than the distribution of basin and range upper crustal fracture systems

receiver function study of the crustal structure of the southeastern caribbean plate boundary and venezuela 1 we have investigated crustal thickness and composition across the southeastern caribbean plate boundary with the receiver function technique we used teleseismic data recorded by a temporary broadband array deployed under the bolivar project and the permanent national seismic network of venezuela we used the primary p to s conversion and crustal reverberations to estimate crustal thickness and average crustal vp vs ratio over the region we observe large variations in crustal thickness and poisson s ratio estimated moho depth ranges from 16 km beneath the southeastern caribbean sea to 52 km beneath northeastern venezuela and the venezuelan andes there is a good correlation between crustal structure and tectonic terranes data from the precambrian guayana shield suggest that the underlying crustal structure is relatively uniform with a moderate thickness 37 km and an intermediate composition a thick crust is found below the foreland basins the two mountain systems in northern venezuela the serrania del interior and the serrania del falcon have a thin crust with arc composition and are likely dynamically supported by elastic rebound or underthrusting of the oceanic plateau that characterizes the southern caribbean on the other hand the venezuelan andes and perija range on the western side of the country are probably isostatically balanced by thick crustal roots

seismic evidence for a 920 km discontinuity in the mantle analysis of hundreds of seismograms from a seismic array in the japanese islands reveals a seismic discontinuity at a depth of about 920 km in the mantle beneath the tonga subduction zone this discontinuity is also evident beneath subduction zones in the japan and flores seas the discontinuity probably represents either a change in mantle chemical composition or the signature of the subducted slab s garnet layer

shear wave anisotropy beneath the tibetan plateau eleven broadband digital seismic stations were deployed across the central tibetan plateau in the first extensive passive source experiment attempted within the tibetan plateau one year of recording resulted in 186 event station pairs which we analyze to determine the characteristics of shear wave splitting in the upper mantle beneath the array measurements of the fast polarization direction and delay time t for sks and direct s arrivals reveal systematic variations along the north south oriented array in the north central region of the plateau very large delay times are observed at three stations the largest of which is budo with of t 2 4 s however at tunl which is off the northern edge of the plateau and 110 km from budo and at sites in the south central plateau t decreases by nearly a factor of 3 we also observe systematic rotation of from about 45 ne to 90 e w from south to north along the array a previously identified zone of inefficient sn propagation correlates well with our region of large t observations the large delay times suggest that a relateively high number of anisotropic crystals are preferentially alligned within the mantle lid beneath the north central portion of the tibetan plateau in most cases fast polarization directions appear to be parallel to surface geologic features suggesting as much as 200 km of the upper mantle has been involved in the collisional deformation that has produced the tibetan plateau

indepth wide angle reflection observation of p wave to s wave conversion from crustal bright spots in tibet three component wide angle seismic data acquired in southern tibet during project indepth show strong p to s converted reflections from reflectors that are aligned at a depth of 15 kilometers beneath the northern yadong gulu rift these converted reflections are locally higher in amplitude than the corresponding p wave reflections modeling of reflection mode conversion as a function of incidence angle indicates that this condition obtains for a reflector that is a solid over fluid interface it is not typical of a solid solid interface the likely candidates for a fluid trapped within the crystalline crust of southern tibet are granitic magma and water brine

seismic anisotropy beneath the eastern margin of the tibetan plateau from sks splitting observations abstract the eastern margin of the tibetan plateau is one of the areas of the entire region with the most complex deformation seismic anisotropy reveals the imprints of past and present deformation in the lithosphere and asthenosphere which could help elucidate the deformation mechanism this study investigated the characteristics of seismic anisotropy beneath 512 seismic stations based on sks phases from 200 earthquake the measured splitting delay time varied from 0 2 1 4 s average 0 7 s in the central sichuan basin scb the average delay time is only 0 35 s consistent with the characteristics of the relatively thick and rigid lithosphere of the stable scb however the northeastern scb has relatively stronger anisotropy with a fast polarization direction fpd with ne sw orientation which is markedly different from the weak anisotropy with nw se orientation in the central scb we speculate that owing to e w compressive stress caused by eastward escape of material from the tibetan plateau the lithosphere in this area has undergone considerable deformation along the right lateral maximum shear direction in the chuandian diamond block the fpd of several adjacent stations was approximately e w i e almost perpendicular to the fpd in the northern and southern regions which might be related to partial melting in the lithosphere and to compressive stress the fpd in the west qinling orogenic belt is approximately parallel to the strike of the surface faults and the average delay time is approximately 1 3 s as the lithosphere is insufficient to produce such strong anisotropy in this region we infer that the lithosphere and asthenosphere both contribute to the seismic anisotropy delay time in the southern songpan ganzi block is much smaller than in the north which might be attributable to low velocity mantle upwelling and complex tectonic stress in this area

crustal thickening and lateral extrusion during the indo asian collision a 3d viscous flow model abstract the 2000 km crustal shortening from the indo asian collision in the past 40 70 million years has been accommodated mainly by crustal thickening in the himalayan tibetan plateau and lateral extrusion of blocks of asian continents however the spatiotemporal evolution of crustal thickening and lateral extrusion hence the far field impact of the indo asian collision on asian tectonics remains controversial here we present a 3d viscous flow model that simulates the partitioning of crustal material between thickening and lateral extrusion during the indo asian collision the model assumes conservation of crustal mass and is constrained by the history of the indo asian plate convergence and by the present day topography of the himalayan tibetan plateau the results show that much of the early collision was absorbed by crustal thickening within the himalayan tibetan plateau however lateral extrusion of crustal material has gradually become dominant in the past 10 20 myr indicating increasing influence of the indo asian collision on asian tectonics

heterogeneity and anisotropy of the lithosphere of se tibet from surface wave array tomography 1 understanding the geotectonic evolution of the southeastern tibetan plateau requires knowledge about the structure of the lithosphere using data from 77 broadband stations in sw china we invert rayleigh wave phase velocity dispersion curves from ambient noise interferometry t 1 0 40 s and teleseismic surface waves t 2 0 150 s for 3 d heterogeneity and azimuthal anisotropy in the lithosphere to 150 km depth our surface wave array tomography reveals 1 deep crustal zones of anomalously low shear wave speed and 2 substantial variations with depth of the pattern of azimuthal anisotropy upper crustal azimuthal anisotropy reveals a curvilinear pattern around the eastern himalayan syntaxis with fast directions generally parallel to the main strike slip faults the mantle pattern of azimuthal anisotropy is different from that in the crust and varies from north to south the tomographically inferred 3 d variation in azimuthal anisotropy helps constrain the source region of shear wave splitting south of 26 n off the high plateau most of the observed splitting can be accounted for by upper mantle anisotropy but for stations on the plateau proper with thick crust crustal anisotropy cannot be ignored on long wavelengths the pattern of azimuthal anisotropy in the crust differs from that in the mantle this is easiest explained if deformation varies with depth the deep crustal zones of low shear wave speed and presumably mechanical strength may represent loci of ductile deformation but their lateral variation suggests that in se tibet localized crustal channel flow and motion along the major strike slip faults are both important

crustal structure and deformation under the longmenshan and its surroundings revealed by receiver function data abstract the convergence of india and eurasia and the obstruction from the rigid sichuan basin cause the longmenshan lms to have the steepest topographic gradient at the eastern margin of the tibetan plateau however the mechanisms of surface uplift are still controversial in this paper we estimate the crustal structure and deformation under the lms and its surroundings by analyzing a large amount of receiver function data recorded by regional seismic networks of the china earthquake administration we apply a comprehensive splitting measurement technique on ps conversion phase at the moho moho ps splitting to calculate crustal anisotropy from azimuthal variations of receiver functions our results show that most of the seismic stations beneath the lms area exhibit significant seismic anisotropy with the splitting time of 0 22 0 94 s and a fast polarization direction of nw se while less or even no crustal anisotropy has been observed under the sichuan basin comparing the fast polarization directions of moho ps splitting with the indicators of lithospheric deformation such as shear wave splitting absolute plate motion and global positioning system imply a consistent tendency of deformation between the lower crust and upper mantle but decoupling deformation in the crust beneath the lms area we further compare moho ps splitting time to that estimated from previous sks splitting indicating that crustal anisotropy is an important source of the sks splitting time in this study area in addition a thick crust 50 km with high vp vs values 1 74 1 86 is also observed using the h stacking method these seismic observations are consistent with the scenario that the lms area has been built by the lower crustal flow combined with the seismic reflection refraction profile and geology studies we further suggest that the lower crustal flow may extrude upward into the upper crust along the steeply dipping strike faults under the lms area resulting in the surface uplift of the lms

constraining p wave velocity variations in the upper mantle beneath southeast asia we have produced a p wave model of the upper mantle beneath southeast se asia from reprocessed short period international seismological centre isc p and pp data short period p data of the annual bulletin of chinese earthquakes abce and long period pp p data we used 3d sensitivity kernels to combine the datasets and mantle structure was parameterized with an irregular grid in the best sampled region our data resolve structure on scale lengths less than 150 km the smearing of crustal anomalies to larger depths is reduced by a crustal correction using an a priori 3d model our tomographic inversions reveal high velocity roots beneath the archean ordos plateau the sichuan basin and other continental blocks in se asia beneath the himalayan block we detect high seismic velocities which we associate with subduction of indian lithospheric mantle this structure is visible above the 410 km discontinuity and may not connect to the remnant of the neo tethys oceanic slab in the lower mantle our images suggest that only the southwestern part of the tibetan plateau is underlain by indian lithosphere and thus that the upper mantle beneath northeastern tibet is primarily of asian origin our imaging also reveals a large scale high velocity structure in the transition zone beneath the yangtze craton which could have been produced in multiple subduction episodes the low p wave velocities beneath the hainan island are most prominent in the upper mantle and transition zone they may represent counter flow from the surrounding subduction zones and may not be unrelated to processes beneath eastern tibet 2005 elsevier b v all rights reserved

lithospheric flexure in the sichuan basin and longmen shan at the eastern edge of tibet 1 the mountain range at the steep eastern edge of the tibetan plateau the longmen shan was deformed by a mw7 9 earthquake with oblique thrust and strike slip motion in 2008 the tectonic processes and structure of the lithosphere beneath this range have been controversial gravity measurements reflect the distribution of mass within the earth including the large load of rock above the geoid in mountain ranges we investigate the response of the lithosphere to the load of the longmen shan and estimate the flexural rigidity or effective elastic thicknessteusing new gravity data acquired by recent satellites and the combined goco2s and eigen 6c datasets the free air gravity anomalies over the longmen shan show that its mass is supported by flexure of the adjacent sichuan basin lithosphere similar to the flexural support of the himalayas analysis of a stacked profile of the free air anomalies shows that the effective elastic thickness of the basin lithosphere is greater than 10 km but has a broad minimum misfit function with no upper limit on the thickness two dimensional admittance analysis shows thete of easternmost tibet is very low approximately 7 km

geodynamics of the southeastern tibetan plateau from seismic anisotropy and geodesy ongoing plate convergence between india and eurasia provides a natural laboratory for studying the dynamics of continental collision a fi rst order process in the evolution of continents regional climate and natural hazards in southeastern tibet the fast directions of seismic anisotropy determined using shear wave splitting analysis correlate with the surfi cial geology including major sutures and shear zones and with the surface strain derived from the global positioning system velocity fi eld these observations are consistent with a clockwise rotation of material around the eastern himalayan syntaxis and suggest coherent distributed lithospheric deformation beneath much of southeastern tibet at the southeastern edge of the tibetan plateau we observe a sharp transition in mantle anisotropy with a change in fast directions to a consistent e w direction and a clockwise rotation of the surface velocity surface strain fi eld and fault network toward burma around the eastern himalayan syntaxis the coincidence between structural crustal features surface strain and mantle anisotropy suggests that the deformation in the lithosphere is mechanically coupled across the crust mantle interface and that the lower crust is suffi ciently strong to transmit stress at the southeastern margin of the plateau in yunnan province a change in orientation between mantle anisotropy and surface strain suggests a change in the relationship between crustal and mantle deformation lateral variations in boundary conditions and rheological properties of the lithosphere play an important role in the geodynamic evolution of the himalayan orogen and tibetan plateau and require the development of three dimensional models that incorporate lateral heterogeneity

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

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

development of shape and lattice preferred orientations application to the seismic anisotropy of the lower crust abstract we review the physical basis of the development of fabrics in plastic and viscous flow and illustrate the typical fabrics formed by these processes in the main rock forming silicates of the lower crust feldspar quartz pyroxene and amphibole the orientation process in plastic deformation where a single slip is dominant is recalled and the role of the constraint of neighbouring grains is emphasized the fabric development of anisometric crystals in viscous flow is discussed as a function of the main controlling parameters shear strain aspect ratio and interference between crystals the same sense of fabric asymmetry is introduced by plastic and viscous flow between the flow plane and the shape preferred orientation and hence coherent kinematic analysis can be undertaken in both modes of flow in order to assess the role of such fabrics in the seismic laminations of the lower continental crust we have calculated the seismic p wave properties of typical fabrics for hypothetical monomineralic and polymineralic rocks the calculations show that the strongest anisotropies develop in monomineralic rock with values between 5 and 16 compared with 5 and 8 for typical rock compositions the strongest anisotropies for layered monomineralic rocks generated by fabrics is only 6 compared to the 14 suggested by model studies of the observed seismic laminations we suggest that other effects such as compositional layering and or constructive interference of seismic waves are responsible for augmenting the apparent anisotropy

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

structural features of the subducting slab beneath the kii peninsula central japan seismic evidence of slab segmentation dehydration and anisotropy 1 seismic evidence of lithospheric slab segmentation dehydration and anisotropic rocks beneath the kii peninsula central japan is detected with receiver function rf analyses by stacking radial and transverse rfs with moveout perturbations we estimate the downdip azimuth dip angle and depth of the slab moho beneath each seismic station comparing the downdip azimuth to the slab convergence azimuth we can classify the slab into three segments the segments correspond to both the fault zones of past megathrust earthquakes and to the spatial distribution of intraslab seismicity the slab moho lies just above the intraslab seismicity in the southwestern kii peninsula slab moho lies deeper at the eastern part of the peninsula at the bottom of seismicity this suggests that earthquakes occur only within the slab mantle under the southern kii peninsula the dip angle of slab moho becomes steeper at 40 to 50 km depth perhaps caused by dehydration transformation of basalt to eclogite at some stations located in the western part of the kii peninsula rf estimated dip angle of the interface is much steeper than that estimated by spatial trends in seismicity since the local seismicity distribution reflects well the larger scale geometry of the slab the rf estimated dip angle may imply anisotropic wave speeds within strongly sheared and or hydrated media around the slab moho these detailed features of the slab may be one key in distinguishing the source regions of megathrust earthquakes in the nankai region

amphibole and lower crustal seismic properties abstract measurements of seismic anisotropy frequently offer insights of the deformation kinematics in situ within the earth in the continental crust seismic anisotropy is commonly ascribed to regionally aligned lattice preferred orientation lpo of mica crystals however determinations of composition suggest that the deep continental crust contains rather little mica in this contribution the role of amphibole a far more plausible constituent especially for deep continental crust of basic composition is assessed a field analogue approach is adopted where the lpo are measured correlated with the strain state of the rocks and used to calculate the seismic properties the seismic properties may be up scaled through the specimen using the modal proportions of the constituent mineral phases to examine the role of composition in controlling the intensity of seismic anisotropy seismic properties are computed using an array of modal compositions so called rock recipes the case history comes from the central block of the lewisian of nw scotland where a regional basic dyke swarm the scourie dykes is locally deformed within amphibolite facies shear zones conditions that are representative of the deep crust shear strain profiles across a 100 m wide shear zone using deflected banding in the host gneisses yield values of 0

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

molybdenum isotope evidence for global ocean anoxia coupled with perturbations to the carbon cycle during the early jurassic relatively brief periods of severe paleoenvironmental change during the jurassic and cretaceous were associated with the widespread accumulation of organic rich marine deposits termed oceanic anoxic events oaes these intervals involved abrupt global warming of 5 10 c higher rates of continental weathering elevated extinction rates and large scale perturbations to the global carbon cycle the major oaes also overlapped temporally the emplacement of large igneous provinces however despite being known as oaes the extent of seawater anoxia at those times is undefined and the causative processes remain unclear here we show how changes in seawater molybdenum isotope ratios a proxy for seawater anoxia during the toarcian early jurassic oae define the onset and expansion of oxygen deficient conditions our data also place constraints on the areal extent of marine anoxia during the event and demonstrate that anoxia expanded and contracted periodically broadly in line with precession driven changes in 13 corg despite their intermittent occurrence over geological history oaes have an important contemporary relevance because the magnitude and high rates of environmental change then were broadly similar to those occurring at the present day