palaeomagnetism and 40ar 39ar geochronology of upper palaeogene volcanic rocks from central tibet implications for the central asia inclination anomaly the palaeolatitude of tibet and post 50 ma shortening within asia summary palaeomagnetic results from upper palaeogene volcanics from three localities in the qiangtang terrane of central tibet record a maximum likelihood palaeolatitude of 28 7 3 7 n these are the first volcanic based palaeomagnetic results from the plateau interior that appear to average secular variation well enough to yield a reliable palaeolatitude estimate this palaeolatitude estimate is significantly lower than the 36 3 n predicted for central tibet by late palaeogene eurasia reference poles consistent with the so called central asia inclination anomaly but it is consistent with values predicted from upper palaeogene reference poles for mongolia the mean inclination of coeval sedimentary rocks from central tibet is 15 5 shallower than our volcanic results similar to the magnitude of sedimentary inclination shallowing estimated from elongation inclination analysis of remanent directions in clastic rocks from the same region we conclude that most of the inclination anomaly observed in volcanic rocks and that which remains after removal of sedimentary flattening in clastic rocks probably is due to the combined effect of reference poles that are not representative of stable asia and by cenozoic intracontinental shortening within central asia on the order of several hundred not more than 1000 kilometers our compilation of palaeomagnetic data from the qiangtang terrane further suggest that the southern margin of the lhasa terrane the southern margin of pre india asia was located as far south as 20 n throughout eocene time these palaeolatitudes in conjunction with the upper crustal shortening record for asia indicate that only approximately one third of the post 50 ma india asia convergence is partitioned in asian lithosphere low latitudes also imply that much of the proto tibetan plateau occupied southerly tropical to subtropical latitudes which has important implications for the weathering and erosion of the tibetan plateau palaeoaltimetry proxies and early cenozoic climate dynamics

multiple episodes of continental subduction during india asia convergence insight from seismic tomography and tectonic reconstruction abstract high wavespeed tomographic anomalies shallower than 1100 km beneath the india asia collision zone are interpreted as continental slabs subducted during collision combining anomaly positions with paleogeographic reconstructions of india we constrain the spatio temporal evolution of multiple episodes of continental subduction likely related to these anomalies this study highlights the different evolution at lithospheric scales of the western and eastern parts of the collision zone the evolution of the western part is characterized by two episodes of steep subduction of the northern margin of india the first episode involving an area with a lateral extent as large as 1500 km started at about 40 30 and ended by a slab break off process at 15 ma the second episode consists on subduction beneath the hindu kush mountains since 8 ma to the east of the collision zone no anomaly related to steep subduction along the northern edge of india is found we interpret two tomographic anomalies beneath the eastern border of the indian plate beneath burma and beneath the andaman sea as the result of two successive episodes of southeastward extrusion followed by subduction we suggest that both extruded portions were initially located along the northern margin of india and that they slid around the eastern syntaxis then southward along the eastern boundary of indian plate both portions subducted along the eastern border of india south of the eastern syntaxis we provide a rough estimate of the amount of indian lithosphere consumed during these subduction and extrusion episodes by comparing this amount with the total amount of indian lithosphere at the onset of collision we conclude that these processes accommodated most of the india asia convergence during collision

simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection through the assimilation of present day mantle seismic structure adjoint methods can be used to constrain the structure of the mantle at earlier times i e mantle initial conditions however the application to geophysical problems is restricted through both the high computational expense from repeated iteration between forward and adjoint models and the need to know mantle properties such as viscosity and the absolute magnitude of temperature or density a priori we propose that an optimal first guess to the initial condition can be obtained through a simple backward integration sbi of the governing equations thus lessening the computational expense given a model with known mantle properties we show that a solution based on an sbi generated first guess has smaller residuals than arbitrary guesses mantle viscosity and the effective rayleigh number are crucial for mantle convection models neither of which is exactly known we place additional constraints on these basic mantle properties when the convection induced dynamic topography on earth s surface is considered within an adjoint inverse method besides assimilating present day seismic structure as a constraint we use dynamic topography and its rate of change in an inverse method that allows simultaneous inversion of the absolute upper and lower mantle viscosities scaling between seismic velocity and thermal anomalies and initial condition the theory is derived from the governing equations of mantle convection and validated by synthetic experiments for both one layer viscosity and two layer viscosity regionally bounded spherical shells for the one layer model at any instant of time the magnitude of dynamic topography is controlled by the temperature scaling while the rate of change of topography is controlled by the absolute value of viscosity for the two layer case the rate of change of topography constrains upper mantle viscosity while the magnitude of dynamic topography determines the temperature scaling lower mantle viscosity when upper mantle lower mantle density anomaly dominates the flow field this two stage scheme minimizes the tradeoff between temperature and lower mantle viscosity for both cases we show that the theory can constrain mantle properties with errors arising through the adjoint recovery of the initial condition for the two layer model this error is manifest as a tradeoff between the temperature scaling and lower mantle viscosity

nonlinear dynamics granular media and dynamic earthquake triggering the landers earthquake of 1992 triggered a series of earthquakes in other parts of california focusing attention on the mechanism of remote triggering laboratory experiments on granular material under pressure suggest that triggering is a result of seismic waves impinging on a fault and inducing elastic nonlinear behaviour of the fault core accompanied by instantaneous weakening and failure the resulting softening to weakening model fits in well with field observations of real earthquakes the 1992 magnitude 7 3 landers earthquake triggered an exceptional number of additional earthquakes within california and as far north as yellowstone and montana1 2 3 since this observation other large earthquakes have been shown to induce dynamic triggering at remote distances for example after the 1999 magnitude 7 1 hector mine1 and the 2002 magnitude 7 9 denali4 earthquakes and in the near field as aftershocks5 the physical origin of dynamic triggering however remains one of the least understood aspects of earthquake nucleation1 2 3 4 5 the dynamic strain amplitudes from a large earthquake are exceedingly small once the waves have propagated more than several fault radii for example a strain wave amplitude of 10 6 and wavelength 1 m corresponds to a displacement amplitude of about 10 7 m here we show that the dynamic elastic nonlinear behaviour of fault gouge perturbed by a seismic wave may trigger earthquakes even with such small strains we base our hypothesis on recent laboratory dynamic experiments conducted in granular media a fault gouge surrogate6 7 from these we infer that if the fault is weak8 9 10 seismic waves cause the fault core modulus to decrease abruptly and weaken further if the fault is already near failure this process could therefore induce fault slip

global slab deformation and centroid moment tensor constraints on viscosity we analyze moment tensor solutions from deep subduction zone earthquakes to determine global slab deformation patterns inferred strain rates are compared to predicted deformation patterns from fluid models to help constrain the first order radial and lateral viscosity structure of the earth while all slabs that reach the lower mantle are compressed at their tip intermediate depth patterns are more complex we compute 3 d spherical flow with various slab rheologies and compare the angular misfit between the compressive eigenvectors of the resultant stress field and global centroid moment tensor gcmt solutions we find that upper mantle slab viscosities of 10 100 and lower mantle viscosities of 30 100 times the upper mantle produce the best match to gcmts a 0 1 viscosity reduction in the asthenosphere seems preferred slab geometry and lower mantle viscosity exert significant control on deformation inclusion of the phase changes at 410 km and 660 km increases extensional deformation at intermediate depth and compressional deformation at the lower mantle improving the match to gcmts for strong slabs our conclusions are fairly insensitive to surface boundary conditions however models which include net rotations of the surface with respect to the lower mantle produce compression at intermediate depths for west directed slabs and extension for east directed slabs without allowing for regional variations these models yield the best match to gcmts while significant deviations between model and seismicity remain our results show that seismicity provides an underutilized constraint for slab dynamics

india eurasia collision chronology has implications for crustal shortening and driving mechanism of plates the motion of the indian plate is determined in an absolute frame of reference and compared with the position of the southern margin of eurasia deduced from palaeomagnetic data in tibet the 2 600 900 km of continental crust shortening observed is shown to have occurred in three different episodes subduction of continental crust intracontinental thrusting and internal deformation and lateral extrusion the detailed chronology of the collision and plate reorganizations in the indian and pacific oceans supports the hypothesis that slab pull is a dominant driving mechanism of plate tectonics

cenozoic plate reconstruction of southeast asia abstract the india eurasia collision alone set up a series of chain reactions and caused the formation and destruction of sedimentary basins within the domain of the collision belt changes in the rate and angle of convergence between the india and eurasia plates reflect different stages of tectonic development in southeast asia for example extrusion of the indochina block induced the consumption of the pre existing proto south china sea along northeastern kalimantan and led to the eventual opening of the south china sea along the south china margin subsequent motions of the sino burma thailand malay peninsula sumatra and kalimantan blocks have produced a succession of basins stretching from north sumatra to central thailand and on to the natuna area we present reconstructions of the southeast asia region at 60 ma 50 ma 40 ma 30 ma 20 ma 15 ma 10 ma and 5 ma it is clear from the reconstructions that the impact between greater india and southeast asia took place in the northwestern part of southeast asia the duration for the impact was probably from the middle eocene to early miocene this timing is in good agreement with the dating of the main red river fault motion wu et al 1989 scharer et al 1990 since the impact between greater india and southeast asia was basically west of the burma block there is no reason to assume that the sumatra malay peninsula and kalimantan should extrude to the southeast first along the left laterally displaced mae ping and three pagodas fault zones as suggested by peltzer and tapponnier 1988 if the opening of the central thailand basins the gulf of thailand and the malay basin are taken into consideration a dextral megashear zone is required to compensate the relative motion between indochina and the malay peninsula this dextral megashear zone might even extend into western kalimantan and serve as a boundary between the indochina block and kalimantan

remnants of a cretaceous intra oceanic subduction system within the yarlung zangbo suture southern tibet extensive field investigations along the yarlung zangbo suture zone in southern tibet reveal the presence of now fragmented remnants of a south facing intra oceanic subduction system this system developed within tethys during the cretaceous the associated arc forearc ophiolite and subduction complex were emplaced onto the leading edge of india at the end of the cretaceous rapid sedimentation in oblique slip basins and disruption of water saturated sediments into melange was widespread and concomitant with ophiolite emplacement we describe the tectonic entities that developed during this previously unrecognized phase of tethys tibet evolution and present a new model for the evolution of this portion of tibet

the himalayan arc large scale continental subduction oroclinal bending and back arc spreading abstract palaeomagnetic results from the himalayan arc and southern tibet compared with simulated apparent polar wander paths for the indian plate show a consistent pattern of rotations of the himalayan arc relative to the indian shield varying gradually from 45 clockwise in the northwestern himalaya to slightly counterclockwise in the lhasa region this pattern is consistent with continental underthrusting of greater india beneath the tibetan plateau since the early miocene over at least 650 km at the longitude of western nepal and oroclinal bending since the latest miocene available palaeomagnetic observations support the steady state model for the formation of the himalayan arc with refinements as follows 1 collision between greater india s northern boundary and southern asia occurred at equatorial latitudes with progressive suturing from palaeocene in the northwestern himalaya until early eocene in the eastern himalaya 2 continuing convergence and indentation of greater india into southern asia over about 2000 km up to the early miocene resulted in southeastward extrusion of indochina and 3 neogene counterclockwise rotational underthrusting of greater india along the main central thrust with pliocene quaternary oroclinal bending of the himalayan arc

the zedong terrane a late jurassic intra oceanic magmatic arc within the yarlung tsangpo suture zone southeastern tibet an overturned sequence of igneous and volcaniclastic rocks crops out along the yarlung tsangpo suture zone ytsz in southeastern tibet these rocks are remnants of an intra oceanic island arc the zedong terrane that once lay between india and asia 40ar 39ar dating and u pb ion microprobe analyses reported here reveal that this arc was active during the jurassic u pb dating of zircon from a dacite breccia from the middle portion of the zedong terrane yields an age of 161 2 3 ma 1 40ar 39ar dating of hornblende from a cross cutting andesite dyke yields the youngest age of 152 2 3 3 ma 1 40ar 39ar results from hornblende from andesite dykes and an andesite breccia from the upper portion yield a mean age of 156 0 3 ma 1 additional u pb and 40ar 39ar dating of zircon and hornblende from quartz diorite yields a mean age of 156 8 0 8 ma 1 geochronological data reported here and other published work indicate that the intra oceanic subduction system to which the zedong terrane belonged was active from at least late jurassic to early to mid cretaceous

application of three dimensional lagrangian residual transport three dimensional lagrangian residual circulation computed from an intratidal hydrodynamic model is used to drive an intertidal water quality transport model the sum of the eulerian residual and stokes drift is used as a first order approximation for the lagrangian residual current the stokes drift approximates residual currents induced by the nonlinear interactions of the tidal currents computed residual velocities show close agreement with a two dimensional vertical longitudinal analytical result the methodology is applied to chesapeake bay for the entire year of 1985 salinity computed with the intertidal transport model i e with lagrangian residuals shows good agreement with observed salinity data and that computed with the intratidal hydrodynamic model the basic effect of the stokes drift is to transport water mass e g salinity up estuary the procedures provide a practical means of interfacing intratidal hydrodynamic model information to water quality models for making long term in

coseismic and postseismic velocity changes measured by repeating earthquakes 1 repeating earthquakes that rupture approximately the same fault patch and have nearly identical waveforms are a useful tool for measuring temporal changes in wave propagation in the earth s crust since source and path effects are common to all earthquakes in a repeating earthquake sequence multiplet differences in their waveforms can be attributed to changes in the characteristics of the medium we have identified over 20 multiplets containing between 5 and 40 repeating events in the aftershock zones of the 1989 loma prieta and 1984 morgan hill california earthquakes postmain shock events reveal delays of phases in the early s wave coda of as much as 0 2 s relative to premain shock events the delay amounts to a path averaged coseismic velocity decrease of about 1 5 for p waves and 3 5 for s waves since most of the multiplets are aftershocks and follow omori s law we have excellent temporal sampling in the immediate postmain shock period we find that the amplitude of the velocity decrease decays logarithmically in time following the main shock in some cases it returns to the premain shock values while in others it does not similar results are obtained for the morgan hill main shock because the fractional change in s wave velocity is greater than the fractional change in p wave velocity it suggests that the opening or connection of fluid filled fractures is the underlying cause the magnitude of the velocity change implies that low effective pressures are present in the source region of the velocity change our results suggest that the changes are predominantly near the stations and shallow but we cannot exclude the possibility that changes occur at greater depth as well if the variations are shallow we may be detecting the lingering effects of nonlinearity during main shock strong ground motion if the variations are deep it suggests that pore pressures at seismogenic depths are high which would likely play a key role in the earthquake process

age spreading rates and spreading asymmetry of the world s ocean crust we present four companion digital models of the age age uncertainty spreading rates and spreading asymmetries of the world s ocean basins as geographic and mercator grids with 2 arc min resolution the grids include data from all the major ocean basins as well as detailed reconstructions of back arc basins the age spreading rate and asymmetry at each grid node are determined by linear interpolation between adjacent seafloor isochrons in the direction of spreading ages for ocean floor between the oldest identified magnetic anomalies and continental crust are interpolated by geological estimates of the ages of passive continental margin segments the age uncertainties for grid cells coinciding with marine magnetic anomaly identifications observed or rotated to their conjugate ridge flanks are based on the difference between gridded age and observed age the uncertainties are also a function of the distance of a given grid cell to the nearest age observation and the proximity to fracture zones or other age discontinuities asymmetries in crustal accretion appear to be frequently related to asthenospheric flow from mantle plumes to spreading ridges resulting in ridge jumps toward hot spots we also use the new age grid to compute global residual basement depth grids from the difference between observed oceanic basement depth and predicted depth using three alternative age depth relationships the new set of grids helps to investigate prominent negative depth anomalies which may be alternatively related to subducted slab material descending in the mantle or to asthenospheric flow a combination of our digital grids and the associated relative and absolute plate motion model with seismic tomography and mantle convection model outputs represents a valuable set of tools to investigate geodynamic problems

subduction of diverging plates and the principles of slab window formation consumption of an ocean basin by subduction commonly brings a sea floor spreading ridge toward a deep sea trench if plate divergence and convergence continue after the ridge intersects the subduction zone a slab window forms between the subducted parts of the diverging oceanic plates producing anomalous thermal physical and chemical effects in the surrounding asthenospheric mantle in turn these conditions alter the tectonic and magmatic evolution of the overriding plate usually disturbing ordinary fore arc and arc regimes differential lithospheric stresses on opposite sides of the triple junction contribute to disturbances in the overriding plate anomalous magmatism from fore arc to back arc may be accompanied by fore arc metamorphism strike slip faulting uplift extension and in extreme cases rifting the shape and size of the window are controlled mainly by the pre subduction ridge transform trench configuration slab dip angles and vectors of plate convergence subducted ridge segments expand into windows whose margins approximately parallel the motion vectors between the triple junction and the subducting plates subducted transform faults continue to be active usually as oblique slip faults until the plates separate as transform faults subduct they become longer on the plate which occupies the acute angle between ridge and trench and shorter on the other plate trains of isolated windows produced by subduction of a segmented ridge transform system progressively expand during descent commonly merging together to form a composite slab window oblique subduction of a highly segmented ridge is likely to produce two or more fraternal slab windows one at each site of ridge trench intersection above a slab window arc volcanism diminishes and may be replaced by volcanism of mid ocean ridge or rift affinity the change in chemical character reflects various processes including elevated heat flow decreasing hydration of the upper mantle juxtaposition of supra and sub slab mantle reservoris asthenospheric upwelling and melting of the trailing plate edges if the slab window migrates the anomalous magmatic regime may be replaced by renewed arc volcanism identifying the effects of slab windows in ancient convergent margin assemblages requires an understanding of slab window principles and implications

subducted slabs and the geoid constraints on mantle rheology and flow geoid anomalies are primarily the result of the density contrasts driving mantle convection and plate motion the total geoid anomaly resulting from a given density contrast in a convecting earth is affected by the mass anomalies associated with the flow induced deformation of the upper surface and internal compositional boundaries as well as by the density contrast itself these boundary deformations and hence the total gravity field depend on the radial distribution of effective viscosity if the internal density contrasts can be estimated as is the case for subducted slabs useful constraints can be placed on the depth and on the variation of viscosity with depth of the convecting system the degree 4 9 components of the observed long wavelength geoid are highly correlated with those predicted by a density model for seismically active subducted slabs the positive sign of the correlation requires that the effective viscosity increases with depth by a factor of 30 or more the amplitude of the correlation cannot be explained by the density contrasts associated with just the seismically active parts of subducted slabs however the amplitude can be explained if the density contrasts associated with subduction extend into the lower mantle or if old lithosphere is piled up at the base of the upper mantle beneath subduction zones to a thickness in excess of 350 km over horizontal distances of thousands of kilometers mantlewide convection in a mantle that has a viscosity increasing with depth provides a simple explanation of the long wavelength geoid anomalies over subduction zones

tethyan subducted slabs under india abstract tomographic imaging of the mantle under tibet india and the adjacent indian ocean reveals several zones of relatively high p wave velocities at various depths under the hindu kush region in northeastern afghanistan and southern tajikistan a regional northward dipping slab is seen in the entire upper 600 km of the mantle and is apparently still attached to the lithosphere of the indian plate under northern pakistan this same slab shows a roll over structure with the deeper portion overturned and dipping southward as can also be seen in the distribution of earthquake hypocenters farther east southeast e g in the vicinity of nepal a well resolved anomaly below 450 km depth is connected to the slab under the hindu kush but seems to be separated from the lithosphere above 350 km these upper mantle anomalies are interpreted as the remnants of delaminated sub continental lithosphere that went down when greater india continued to converge northward with asia after 45 ma the deeper high velocity anomalies under the indian sub continent appear clearly separated from the shallower ones as well as from each other and are inferred to represent remnants of oceanic lithospheric slabs that have sunk into the lower mantle and were subsequently overridden by the indian plate they occur at depths between 1000 and 2300 km and occasionally descend down to the core mantle boundary the anomalies form three parallel wnw ese striking zones we interpret the two southern zones as remnants of oceanic lithosphere that was subducted when the neo tethys ocean closed between india and tibet in the cretaceous and earliest tertiary the northern deep mantle zone under northern afghanistan the himalayas and the lhasa block in southern tibet may represent the last subducted remnant of the paleo tethys ocean which is thought to have closed before the hauterivian stage of the early cretaceous the middle zone continues southeastward as a rather straight high velocity zone towards sumatra where it becomes convex southward and parallel to the subduction zone under the sunda arc comparison of this straight middle zone near india with the shallower upper 600 1000 km northern zone which displays a cusp like shape near the yunnan sw china syntaxis of the eastern himalayas supports the notion that the shallow northern zone represents later subduction than the deeper middle zone the suggestion of a counterclockwise rotation 20 of the indian plate during tertiary indentation of asia is supported by these features the present day latitudes of 5 35 n of the deep slabs under india and adjacent areas correspond to the approximate paleolatitudes of the cretaceous subduction zones the slab remnants in the middle mantle occur therefore near the ancient locations where they started their downward journey which implies that lateral movements in the deeper mantle were not large

discussion on the role of the red river shear zone yunnan and vietnam in the continental extrusion of se asia journal vol 163 2006 1025 1036 p h leloup p tapponnier r lacassin write in his recent paper searle 2006 acknowledges that the 1000 km long ailao shan red river shear zone is a large miocene left lateral shear zone but speculates that left lateral slip started after 21 ma and claims that the total finite offset remains unknown from this he concludes that continental extrusion was only a relatively minor tectonic factor during the india asia collision as long argued by other workers e g england houseman 1986 cobbold davy 1988 dewey et al 1989 houseman england 1993 we summarize below the field and geochronological evidence that makes us maintain a viewpoint in better accordance with facts timing of left lateral shear along the ailao shan red river shear zone the ailao shan red river shear zone is composed mostly of high grade metamorphic rocks and deformed granitoids with ubiquitous evidence for left lateral shear parallel to the belt e g tapponnier et al 1986 1990 leloup et al 1993 1995 2001 the crystallization of the granitoids has been dated between 22 and 35 ma fig 1d f for example by scharer et al 1990 1994 and zhang scharer 1999 leading those workers to propose that left lateral shear started at least c 35 ma ago in contrast searle 2006 claims that all the deformed granitoids found within the shear zone predate left lateral shear and that their crystallization ages should thus be interpreted to provide an upper limit for the onset of deformation rather than a lower limit thus suggesting a maximum age of 21 ma for this deformation a clear understanding of the p t history and in situ deformation history of the intrusions as well as of their relationships with surrounding paragneisses is fundamental for interpreting correctly the geochronological

slow dynamics in the nonlinear elastic response of berea sandstone a typical resonance curve measured accel eration versus drive frequency made on a thin bar of rock shows peak bending with a softening nonlinear modulus as drive levels are increased the shapes of these nonlinear resonance curves were found in earlier work to depend somewhat on sweep rate these slow dynamics are now examined and quantified we have measured slow dynamics in a 0 3 m long 50 mm diam eter bar of berea sandstone under ambient conditions peak strain levels during the experiments ranged from 10 to 10 at driving frequencies near 4 khz the fun damental longitudinal resonance frequency of the bar slow dynamics begin to manifest themselves at strain amplitudes above 10 a at ambient conditions and at the onset of nonlinear peak bending strains above this value condition the rock altering its response for min utes to hours after the drive has been turned off ter the pressure on the rock was changed the work presented here is the first to specifically document the slowly varying time response to periodic excitation i e slow dynamics of a berea sandstone sample moreover we demonstrate that the slow dynamics produce promi nent effects even at room conditions slow dynamics manifest themselves in a sandstone sample during high intensity longitudinal elastic reso nance we first report the effects of high strain on a sample by comparing conditioned resonance curves e g after the sample is driven at high intensity with a set of fully recovered curves e g after the sample is left untouched for a long period we illustrate that the first resonance curve made on a recovered sample differs remarkably from subsequent curves we present several results that illustrate aspects of the sample s memory and concludiwith a discussion of measurements to de termine how the rock recovers and to understand what

fluid effects on velocity and attenuation in sandstones measurements of ultrasonic velocity and attenuation in sandstones with a variety of saturating fluids are compared with the predictions of biot s theory of porous media acoustics the data show systematic deviations from biot theory as a function of pore fluid viscosity two non biot attenuation mechanisms are needed to explain the data local flow in microcracks and enhanced attenuation from pore wall roughness a model calculation of attenuation with pore wall roughness shows that attenuation in rocks can be many times larger than biot predictions this attenuation enhancement explains the inability to detect a biot slow wave in fluid saturated rocks and the observed correlation between ultrasonic attenuation and clay content in sandstones the non biot attenuation mechanisms are negligible when extrapolated to the lower frequencies of geophysical field measurements so ultrasonic attenuation in rocks provides no information about attenuation at the field scale

in situ measurements of nonlinear and nonequilibrium dynamics in shallow unconsolidated sediments abstract we present results and in situ measurements from a field experiment in which a large seismic shaker truck is used to induce nonlinear and nonequilibrium dynamics in shallow unconsolidated sediments an array of accelerometers was deployed adjacent to the shaker truck to record strong ground motions exceeding 1 g we determined high strain rayleigh wave dispersion across the array in the band from 5 to 60 hz rayleigh wave phase velocities at frequencies above 10 hz are wave amplitude dependent and a function of the driving force amplitude as driving force was increased phase velocity decreased consistent with nonlinear dynamics we demonstrate the existence of a temporary nonequilibrium state occurring in the near surface soils during and after the induced nonlinear behavior nonlinear conditioning is demonstrated by measuring changes in rayleigh wave phase velocity for input signals with the same applied driving force amplitude a logarithmic slow dynamic recovery process is observed by analyzing the temporal variation in velocity of the noise field produced by the shaker truck when sitting idle measurements from weak motion seismic surveys taken before and after shaking show compelling evidence that induced nonlinear behavior in the shallow sediments is reversible we demonstrate that an active source field based approach has the potential to expand our knowledge of how sediments respond to strong ground motions provide additional insight into the poorly understood slow dynamic recovery process and possibly even lead to a new site specific and noninvasive technique for characterizing the nonlinear properties of sedimentary deposits

harmonic analysis of tides and tidal currents in south san francisco bay california abstract water level observations from tide stations and current observations from current meter moorings in south san francisco bay south bay california have been harmonically analysed at each tide station 13 harmonic constituents have been computed by a least squares regression without inference tides in south bay are typically mixed there is a phase lag of approximately 1 h and an amplification of 1 5 from north to south for a mean semi diurnal tide because most of the current meter records are between 14 and 29 days only the five most important harmonics have been solved for east west and north south velocity components the eccentricity of tidal current ellipse is generally very small which indicates that the tidal current in south bay is strongly bidirectional the analyses further show that the principal direction and the magnitude of tidal current are well correlated with the basin bathymetry patterns of eulerian residual circulation deduced from the current meter data show an anticlockwise gyre to the west and a clockwise gyre to the east of the main channel in the summer months due to the prevailing westerly wind opposite trends have been observed during winter when the wind was variable

re examination of the magnetic lineations of the gascoyne and cuvier abyssal plains off nw australia this thesis addresses causes of two forms of secondary hawaiian volcanism rejuvenated onshore eruptions and offshore hawaiian arch flows it is proposed that secondary volcanism is generated as a direct consequence of lithospheric flexural uplift that surrounds new volcanic shields as they grow this uplift causes decompression of the underlying asthenosphere which is assumed to be chemically and isotopically heterogeneous near its solidus and derived from the hawaiian mantle plume uplift is modeled as the axisymmetric response of an elastic plate to a volcanic point load that grows linearly in time to model flow in the asthenosphere the rate of flexure of the lithosphere is taken as the upper boundary condition on an isoviscous incompressible fluid half space the first feature of secondary volcanism this model explains is the observed spatial gap between secondary volcanism and active shields best agreement is found with the majority of the observed spatial gaps with a lithosphere of effective elastic thickness te 25 35 km secondly this work demonstrates that the flexural model can produce observed crustal production if some magma focusing toward individual eruption sites occurs from the mantle over an area two to ten times the eruption area the third feature this model addresses is that secondary lavas are isotopically distinct from shield lavas in this model melting the same two component mantle forms the secondary and shield lavas but the components are sampled by melting at rates that differ between the locations as predictable functions of depth flexural decompression produces melts that preferentially sample the mantle component that begins melting shallowest and which is associated with sr and nd isotope ratios most like those of secondary lavas melting in the center of a mantle plume is assumed to generate shield volcanism and is predicted to preferentially sample the component that begins melting deepest which is associated with sr and nd isotope ratios more similar to shield lavas models therefore successfully predict the observed mean difference in 87sr 86sr and 143nd 144nd compositions between the secondary and shield lavas the fourth feature addressed is that secondary lavas are alkalic and shield lavas are dominantly tholeiitic to explain this difference the mean extent of partial melting is computed and it is found that a model plume composed mostly of depleted peridotite 90 and some pyroxenite 10 will yield a lower extent of melting for secondary lavas than shield lavas this particular model assumes lithospheric thicknesses 90 100 km and plume potential temperatures mean of 1550 c that are consistent with independent studies of the hawaiian hotspot thus asthenospheric melting by flexural decompression is a viable mechanism of intraplate volcanism which can explain many general characteristics of secondary hawaiian volcanism

paleogene island arc collision related conglomerates yarlung tsangpo suture zone tibet abstract coarse clastic rocks in the liuqu conglomerate formed in both terrestrial and subaqueous settings record a paleogene phase in the tectonic evolution of tibet facies changes are commonly abrupt with rapid changes in clast types grain size and stratal patterns sediments were derived from the leading northern edge of the indian margin and a late jurassic cretaceous intraoceanic island arc that lay within tethys the coarse clastic sedimentary rocks of the liuqu conglomerates are extremely proximal but are locally offset relative to their original source terranes they record aspects of the history of collision between these terranes and are interpreted to have been deposited in oblique slip basins that developed along the zone of collision the absence of clasts derived from terranes to the north of the yarlung tsangpo suture suggests that basins associated with deposition of the liuqu conglomerate developed prior to the final collision between india and asia

birth death and resurrection the life cycle of suprasubduction zone ophiolites 1 suprasubduction zone ssz ophiolites display a consistent sequence of events during their formation and evolution that suggests that they form in response to processes that are common to all such ophiolites this sequence includes the following 1 birth which entails the formation of the ophiolite above a nascent or reconfigured subduction zone this stage is typically characterized by the eruption of arc tholeiite lavas and the formation of layered gabbros and sheeted dike complex 2 youth during which is continued melting of refractory asthenosphere depleted during birth occurs in response to fluid flux from the subducting slab with extensional deformation of the older plutonic suite eruption of refractory lavas and the intrusion of wehrlite pyroxenite 3 maturity with the onset of semistable arc volcanism typically calc alkaline as the subduction zone matures and stabilizes and the intrusion of quartz diorite and eruption of silicic lavas and 4 death which is the sudden demise of active spreading and ophiolite related volcanism which in many cases is linked to collision with an active spreading center and the onset of shallow underthrusting of the buoyant spreading axis expressed as dikes and lavas with oceanic basalt compositions that crosscut or overlie rocks of the older suites 5 resurrection with emplacement by obduction onto a passive margin or accretionary uplift with continued subduction the early stages 1 3 may be diachronous and each stage may overlap in both time and space the existence of this consistent progression implies that ophiolite formation is not a stochastic event but is a natural consequence of the ssz tectonic setting

a recipe for microcontinent formation accreted slivers of continental margins are common in the geologic record but the processes that lead to their formation are poorly understood we observe an association of plume related microcontinent isolation and subsequent long term asymmetries in oceanic crustal accretion based on four recent examples the seychelles in the indian ocean jan mayen in the norwegian greenland sea and the east tasman plateau and the gilbert seamount complex in the tasman sea these microcontinents formed by rerifting of a young continental margin 25 m y old in the vicinity of a mantle plume stem followed by asymmetric seafloor spreading two dimensional numerical stochastic basin modeling suggests that a yield strength minimum along the landward edge of a rifted margin thermally enhanced by heating from a mantle plume may cause a spreading ridge to jump onto this zone of weakness this action isolates a passive margin segment the association of large igneous provinces and microcontinents should be useful for identifying similar events in the geological record

new paleomagnetic results of paleocene volcanic rocks from the lhasa block tectonic implications for the collision of india and asia abstract a combined geochronological and paleomagnetic investigation has been performed on paleocene volcanic sequences in the lhasa block near the localities of mendui 30 1 n 90 9 e a total of 15 sites have been sampled from rhyolitic tuffs stepwise thermal demagnetizations successfully isolated high unblocking temperature characteristic directions the tilt corrected mean direction is d i 359 0 26 1 with 95 9 2 and n 14 sites corresponding to a paleopole at 73 6 n 274 3 e with a 95 7 3 positive fold tests suggest a primary origin for the characteristic remanence in order to provide a more accurate pole we propose to combine site mean directions from this study and achache et al s 1984 the combined average palaeomagnetic direction from early paleocene volcanic rocks is d 356 6 i 25 9 21 7 95 6 8 after tilt correction n 22 sites corresponding to a pole at 73 2 n 282 4 e with a 95 5 4 the paleomagnetic results yield a paleolatitude of 13 6 5 4 n for the southern margin of eurasia at 55 ma compared with expected paleomagnetic directions from the stable india and eurasia blocks significant crustal shortening of 1400 600 km and 2000 550 km respectively may have occurred between the southern margin of eurasia and the stable india and within eurasia since the collision of india and eurasia

tibetan tectonic evolution inferred from spatial and temporal variations in post collisional magmatism cenozoic magmatism on the tibetan plateau shows systematic variations in space and time that must be considered in models concerning tibetan tectonic evolution after the india asia collision which started in the early tertiary and terminated the gangdese arc magmatism in the lhasa terrane made of the southern tibetan plateau widespread potassium rich lavas and subordinate sodium rich basalts were generated from f50 to 30 ma in the qiangtang terrane of northern tibet subsequent post collisional magmatism migrated southwards producing ultrapotassic and adakitic lavas coevally between f26 and 10 ma in the lhasa terrane then potassic volcanism was renewed to the north and has become extensive and semicontinuous since f13 ma in the western qiangtang and songpan ganze terranes such spatial temporal variations enable us to elaborate a geodynamic evolution model which depicts when and how the indian continental lithospheric mantle started thrusting under asia by involving rollback and breakoff of the subducted neo tethyan slab followed by removal of the thickened lhasa lithospheric root we propose that only after the lithospheric removal which occurred at f26 ma could the indian mantle lithosphere have commenced its northward underthrusting and henceforth served as a pivotal control to the himalayan tibetan orogenesis consequently the tibetan plateau is suggested to have risen diachronously from south to north whereas the southern part of the plateau may have been created and maintained since the late oligocene the northern plateau would have not attained its present day elevation and size until the mid miocene when the lower part of the western qiangtang and songpan ganze lithospheres began to founder owing to the push of the underthrust indian mantle lithosphere d 2004 elsevier b v all rights reserved

observation of anomalous elastic behavior in rock at low temperatures we have measured qualitative elastic properties ofbasaltandbereasandstonefromroomtemperaturedown to 4 k using resonant ultrasound spectroscopy rus a simple elastic solid should show a monotonic increase in the elastic constantsas temperature decreases the basalt sam ples show this gross behavior while the sandstone shows a veryunexpectedanomalous regimebetween40kand200k wheretheelasticconstantsdecreasewithdecreasingtemper ature bothrocks showtemperature dependentstructurein both the modulus and internal friction and also signicant hysteresis indicatinghistoryandrate dependentproperties this data provides insight into the time and energy scales of dynamical eects observed in sandstones the low tem perature range 4k 300k is of interest due to the renewed interest in mars and lunar rocks shown to display hysteretic dynamical properties tencate and shankland 1996 and a los alamos basalt rus data allows us to determine the qualitative behavior of the elas tic constants and the internal friction of a sample the data show a remarkable dierence between these two rock samples and previous work guyer et al 1999 also sug gests that the sandstone has exceptional elastic properties we observe evidence of characteristic temperatures in both rocks and a region of anomalous temperature coecient in thesandstone thecoolingandwarmingdatashowhystere sis in both samples with a signicant and repeatable loop apparent in the sandstone we believe the processes asso ciated with the characteristic temperatures evident in these data are related to the processes that produce the observed tencate and shankland 1996 strain driven slow dynamics at room temperatures

cenozoic volcanism and tectonic evolution of the tibetan plateau abstract cenozoic volcanism on the tibetan plateau which shows systematic variations in space and time is the volcanic response to the india asia continental collision the volcanism gradually changed from na rich k rich to potassic ultrapotassic adakitic compositions along with the india asia collision shifting from contact collision i e soft collision or syn collision to all sided collision i e hard collision the sodium rich and potasium rich lavas with ages of 65 40 ma distribute mainly in the lhasa terrane of southern tibet and subordinately in the qiangtang terrane of central tibet the widespread potassic ultrapotassic lavas and subordinate adakites were generated from 45 to 26 ma in the qiangtang terrane of central tibet subsequent post collisional volcanism migrated southwards producing ultrapotassic and adakitic lavas coevally between 26 and 8 ma in the lhasa terrane then potassic and minor adakitic volcanism was renewed to the north and has become extensive and semicontinuous since 20 ma in the western qiangtang and songpan ganze terranes such spatial temporal variations provide important constraints on the geodynamic processes that evolved at depth to form the tibetan plateau these processes involve roll back and break off of the subducted neo tethyan slab followed by removal of the thickened lhasa lithospheric root and consequently northward underthrusting of the indian lithosphere the tibetan plateau is suggested to have risen diachronously from south to north whereas the southern part of the plateau may have been created and maintained since the late oligocene the northern plateau would have not attained its present day elevation and size until the mid miocene when the lower part of the western qiangtang and songpan ganze lithospheres began to founder and detach owing to the persistently northward push of the underthrust indian lithosphere

tidal and residual circulations in coupled restricted and leaky lagoons abstract finite element numerical modelling based on field data is used to study the tidal and tidally induced residual circulation dynamics of a coupled restricted and leaky coastal lagoon system located in the magdalen islands gulf of saint lawrence havre aux maisons lagoon hml is of a restricted nature with a neutral inlet in terms of tidal asymmetry grande entree lagoon gel is of a leaky nature with a marked ebb dominance at the inlet due to direct interactions between the main astronomical tidal constituents the imbalance caused by the different tidal filtering characteristics of both inlets combines with the internal morphological asymmetries of the system to produce a residual throughflow from hml to gel the residual circulation is also characterized by strongest values at both inlets very weak residual currents in hml deep basin and a dipole of residual eddies over the deeper areas of gel further investigations including numerical tracer experiments will be necessary to achieve a full understanding of the long term circulation of this lagoonal system

resonant ultrasound spectroscopy studies of berea sandstone at high temperature resonant ultrasound spectroscopy was used to determine the elastic moduli of berea sandstone from room temperature to 478 k sandstone is a common component of oil reservoirs and the temperature range was chosen to be representative of typical downhole conditions down to about 8 km in agreement with previous works berea sandstone was found to be relatively soft with a bulk modulus of approximately 6 gpa as compared to 37 5 gpa for quartz at room temperature and pressure it was found that berea sandstone undergoes a 17 softening in bulk modulus between room temperature and 385 k followed by an abnormal behavior of similar stiffening between 385 k and 478 k

a new view of near shore dynamics based on observations from hf radar abstract since 1984 the oscr hf radar system has been used in over 50 deployments to measure near shore surface currents for both scientific and engineering applications the enhanced scope resolution and accuracy of these measurements have yielded new insights into the tidal wind and density driven dynamics of the near shore zone tidal current ellipses obtained from these radar measurements have been shown to be in good aggrement with values calculated by numerical models both for the predominant constituents and also for higher harmonics coherent patterns of wind forced currents ahve been determined with strong evidence of a slab like surface response in one deployment with offshore winds blowing over relatively deep water this slab rotated clockwise at near inertial frequency strong up to 20cm s 1 persistent surface residual currents are commonly observed these are almost certainly generated by small horizontal density gradients these observed surface residuals provide ideal data for rigorous testing of 3 d numerical models with a threatened rise in sea level hf radar is well suited for observing the expected changes in the dynamics of near shore regions continuing development of these radar systems offers exciting prospects of remote sensing of both surface waves and currents future applications may extend beyond the near shore region to measurements along the shelf edge in oceanic gyres and for beach processes

the northern suture pakistan margin of a cretaceous island arc the northern suture is a fault separating the cretaceous kohistan island arc terrain northwest himalayas from palaeozoic sediments of the asian plate to the north the kohistan arc includes volcanic and sedimentary rocks andesitic lavas tuffs volcaniclastics slates and limestones metamorphosed to greenschist facies and intruded by the two phase kohistan batholith asian continental margin sediments are mainly of shelf type are variably metamorphosed and intruded by the karakoram batholith the northern suture is a zone of melange from 150 m to 4 km wide and contains blocks of volcanic greenstone limestone red shale conglomerate quartzite and serpentinite in a slate matrix it has a strong planar fabric but in many places bedding is preserved in blocks and matrix and depositional rather than tectonic contacts are seen between the two the melange is inferred to be an olistostrome largely derived from the kohistan arc formed in a small back arc basin between kohistan and asia limestone blocks in the melange are dated as aptian albian post tectonic intrusions yield radiometric ages from 111 to 62 ma the northern suture therefore probably formed in the early late cretaceous during closure of the back arc basin the tethys ocean lay south of kohistan where the main mantle thrust represents the westward continuation of the indus tsangpo suture

estuarine salt conservation a lagrangian approach abstract understanding tidal and mean flow transport mechanisms that maintain the estuarine salt balance against seaward transport by river flow is one of the fundamental questions of estuarine oceanography previous salt transport studies have failed to resolve this problem for two principal methodological reasons in addition to the inherent variability of estuaries first salt transport expansions used to represent the salt balance have included a large number of terms that could not be related to any theory of estuarine circulation and whose physical meaning was thus unclear second it has proven difficult to measure small mean velocities in the presence of much larger tidal variations a new analysis method that expresses lagrangian scalar property transport in terms of time and space variations of eulerian variables is much simpler than expansions previously used and yields new insights into estuarine transport processes this approach like previous expansions is valid for narrow systems in which cross channel transport processes are weak and the ratio of tidal amplitude to mean depth is small a surface boundary condition is used to eliminate one class of troublesome terms the equivalence of the tidal cycle average transport of tidal variations in salinity and the stokes transport of mean salinity is then employed to show that the vertically integrated tidal flow plays no role in the overall salt balance that is seaward transport of salt caused by the river flow can only be balanced by landward transport resulting from correlations between shear and stratification in the mean tidal and overtide flows in a hypothetical inviscid estuary tide and overtide flows would be vertically uniform and outward fluvial salt transport could only be balanced by shear and stratification in the mean gravitational circulation in a more realistic example with strong friction the gravitational circulation would be severely damped and inward transport could only be accomplished by correlations of shear and stratification in the tide and overtide flows

bainang terrane yarlung tsangpo suture southern tibet xizang china a record of intra neotethyan subduction accretion processes preserved on the roof of the world the bainang terrane an intra oceanic island arc subduction complex into which tethyan oceanic rocks were accreted during the cretaceous is preserved within the yarlung tsangpo suture zone of tibet the lithostratigraphic succession established from field mapping records a long history of sedimentation in different portions of the central tethyan domain from late triassic to mid cretaceous time these rocks are preserved within a south verging imbricate thrust stack of thin 1 km thick northward younging tectonic slices five lithotectonic units were mapped in the terrane and these units are assigned to two distinct tracts the northern tract which accumulated on the north side of neotethys was probably separated from its southern counterpart by a mid ocean ridge detailed radiolarian biostratigraphy is used to constrain the timing of depositional events within each tract oceanic plate stratigraphy of the northern tract records its northward travel and mid cretaceous late aptian approach towards a south facing intra oceanic subduction zone rocks in the southern tract developed closer to the indian subcontinent and experienced thermotectonic subsidence and mid jurassic basic alkaline intraplate magmatism they were probably accreted late in the cretaceous variations in structural style across the terrane indicate deformation at different depths and vertical growth of the wedge rather than lateral accretion the overall tectonostratigraphy of the terrane reflects its development in a remote intra oceanic setting

subduction history of the tethyan region derived from seismic tomography and tectonic reconstructions 1 in the mantle underneath the tethyan suture zone large volumes of positive velocity anomalies have been imaged by seismic tomography and interpreted as the present day signature of subducted tethyan lithosphere we investigate the mesozoic cenozoic subduction history of the region by integrating independent information from mantle tomography and tectonic reconstructions three different subduction scenarios for the tethyan oceanic lithosphere representative for the available tectonic reconstructions are used to predict the present thermally anomalous volumes associated with the lithospheric surface subducted since the late mesozoic next these predicted thermal volumes and their expected positions are compared to the relevant anomalous volumes derived from seismic tomographic images in this analysis we include among others the possible effects of ridge subduction and slab detachment after the cenozoic continental collisions absolute plate motion and slab thickening in the mantle our preferred subduction model comprises the opening of large back arc oceanic basins within the eurasian margin the model points to slab thickening by a factor of 3 in the mantle in which case the estimated volumes allow for active oceanic spreading 1 2 5 cm yr in the tethyan lithosphere during convergence our results further indicate the occurrence of early oligocene slab detachment underneath the northern zagros suture zone followed by both westward and eastward propagation of the slab tear and diachronous eocene to miocene slab detachment below the eastern to western himalayas free sinking rates of the detached material of 2 cm yr in the lower mantle provide the best fit between the tomographic mantle structure and our tethyan subduction model

gondwana to asia plate tectonics paleogeography and the biological connectivity of the indian sub continent from the middle jurassic through latest eocene 166 35 ma abstract using the most up to the date information available we present a considerably revised plate tectonic and paleogeographic model for the indian ocean bordering continents from gondwana s middle jurassic break up through to india s collision with asia in the middle cenozoic the landmass framework is then used to explore the sometimes complex and occasionally counter intuitive patterns that have been observed in the fossil and extant biological records of india madagascar africa and eastern eurasia as well those of the more distal continents although the paleogeographic model confirms the traditional view that india became progressively more isolated from the major landmasses during the cretaceous and paleocene it is likely that at various times minor physiographic features principally ocean islands provided causeways and or stepping stone trails along which land animals could have migrated to from the sub continent aside from a likely link albeit broken by several marine gaps to africa for much of this time it is notable that the present day recent biota of madagascar indicates that the ancestors of five land mammal orders plus bats crossed the 400 km wide mozambique channel at different times in the cenozoic it is possible that the kerguelen plateau connected india and australia antarctica in the mid cretaceous approximately 115 90 ma later the seychelles mascarene plateau and nearby elevated sea floor areas could have allowed faunas to pass between southern india and madagascar in the late cretaceous from around 85 65 ma with an early cenozoic extension to this path forming as a result of the reunion hot spot trace islands growing on the ocean floor to the ssw of india the modelling also suggests that india s northward passage towards asia with eventual collision at 35 ma involved the ne corner of the sub continent making a glancing contact with sumatra followed by burma from 57 ma late paleocene onwards a scenario which is compatible with the fossil record indicating that india asia faunal exchanges began occurring at about this time finally we contend that a number of biologically based direct terrestrial migration routes that have been proposed for last 15 m y of the cretaceous asia to india antarctica to madagascar and or india can probably be dismissed because the marine barriers likely varying from 1000 up to 2500 km were simply too wide

plate tectonics and net lithosphere rotation over the past 150 my we have developed an improved model of global digital palaeo plate boundaries and plate motion to describe the distribution and history of plates since the late jurassic from this history we computed net lithospheric rotation nr through time confirming the so called westward drift but only for the past 30 myr the nr has significantly smaller magnitudes 0 13 my past 5 my than for some other plate models it averages to 0 11 0 03 my for the past 50 my with a small but systematic increase toward the present the westward drift seen only for the past 30 my is attributed to the increased dominance of a steadily growing and accelerating pacific plate nr shows peaks with time but only an early tertiary peak of 0 33 my when the indian plate was undergoing the largest known acceleration deceleration can be interpreted with some confidence we find a linear decreasing trend in net rotation over the past 150 my but attribute this trend to increasing reconstruction uncertainties back in time as subduction consumed more than half of the oceanic crust since the jurassic after removing a linear time trend we find a nr average of about 0 12 my for the past 150 my

on tide induced lagrangian residual current and residual transport 2 residual transport with application in south san francisco bay california the transports of solutes and other tracers are fundamental to estuarine processes the apparent transport mechanisms are convection by tidal current and current induced shear effect dispersion for processes which take place in a time period of the order of a tidal cycle however as emphasis is shifted toward the effects of intertidal processes the net transport is mainly determined by tide induced residual circulation and by residual circulation due to other processes the commonly used intertidal conservation equation takes the form of a convection dispersion equation in which the convective velocity is the eulerian residual current and the dispersion terms are often referred to as the phase effect dispersion or sometimes as the tidal dispersion the presence of these dispersion terms is merely the result of a fickian type hypothesis since the actual processes are not fickian thus a fickian hypothesis obscures the physical significance of this equation recent research results on residual circulation have suggested that long term transport phenomena are closely related to the lagrangian residual current or the lagrangian residual transport in this paper a new formulation of an intertidal conservation equation is presented and examined in detail in a weakly nonlinear tidal estuary the resultant intertidal transport equation also takes the form of a convection dispersion equation without the ad hoc introduction of phase effect dispersion in a form of dispersion tensor the convective velocity in the resultant equation is the first order lagrangian residual current the sum of the eulerian residual current and the stokes drift the remaining dispersion terms are important only in higher order solutions they are due to shear effect dispersion and turbulent mixing there exists a dispersion boundary layer adjacent to shoreline boundaries an order of magnitude estimate of the properties in the dispersion boundary layer is given the present treatment of intertidal transport processes is illustrated by an analytical solution for an amphidromic system and by a numerical application in south san francisco bay california the present formulation reveals that the mechanism for long term transport of solutes is mainly convection due to the lagrangian residual current in the interior of a tidal estuary this result also points out the weakness in the tidal dispersion formulation and explains the large variability of the observed values for tidal dispersion coefficients further research on properties of the dispersion boundary layer is needed

restoration of cenozoic deformation in asia and the size of greater india a long standing problem in the geological evolution of the india asia collision zone is how and where convergence between india and asia was accommodated since collision proposed collision ages vary from 65 to 35 ma although most data sets are consistent with collision being underway by 50 ma plate reconstructions show that since 50 ma 2400 3200 km west to east of india asia convergence occurred much more than the 450 900 km of documented himalayan shortening current models therefore suggest that most post 50 ma convergence was accommodated north of the indus yarlung suture zone we review kinematic data and construct an updated restoration of cenozoic asian deformation to test this assumption we show that geologic studies have documented 600 750 km of n s cenozoic shortening across and north of the tibetan plateau the pamir hindu kush region accommodated 1050 km of n s convergence geological evidence from tibet is inconsistent with models that propose 750 1250 km of eastward extrusion of indochina approximately 250 km of indochinese extrusion from 30 to 20 ma of indochina suggested by se asia reconstructions can be reconciled by dextral transpression in eastern tibet we use our reconstruction to calculate the required size of greater india as a function of collision age even with a 35 ma collision age the size of greater india is 2 3 times larger than himalayan shortening for a 50 ma collision the size of greater india from west to east is 1350 2600 km consistent with robust paleomagnetic data from upper cretaceous paleocene tethyan himalayan strata these estimates for the size of greater india far exceed documented shortening in the himalaya we conclude that most of greater india was consumed by subduction or underthrusting without leaving a geological record that has been recognized at the surface

resonance and elastic nonlinear phenomena in rock in a great variety of laboratory experiments over large intervals in stress strain and frequency rocks display pronounced nonlinear elastic behavior here we describe nonlinear response in rock from resonance experiments two important features of nonlinear resonant behavior are a shift in resonant frequency away from the linear resonant frequency as the amplitude of the disturbance is increased and the harmonics in the time signal that accompany this shift we have conducted young s mode resonance experiments using bars of a variety of rock types limestone sandstone marble chalk and of varying diameters and lengths typically samples with resonant frequencies of approximately 0 5 1 5 khz display resonant frequency shifts of 10 or more over strain intervals of 10 7 to 10 6 and under a variety of saturation conditions and ambient pressure conditions correspondingly rich harmonic spectra measured from the time signal progressively develop with increasing drive level in our experiments to date the resonant peak is observed to always shift downward if indeed the peak shifts indicating a net softening of the modulus with drive level this observation is in agreement with our pulse mode and static test observations and those of other researchers resonant peak shift is not always observed even at large drive levels however harmonics are always observed even in the absence of peak shift when detected strain levels exceed 10 7 or so this is an unexpected result important implications for the classical perturbation model approach to resonance results from this work observations imply that stress strain hysteresis and discrete memory may play an important role in dynamic measurements and should be included in modeling this work also illustrates that measurement of linear modulus and q must be undertaken with great caution when using resonance

a cretaceous dike swarm provides evidence of a spreading axis in the back arc basin of the kohistan paleo island arc northwestern himalaya pakistan abstract a 100 4000 m wide and 15 km long dike swarm consisting of basalt and dolerite occurs at the base of the thelichi formation in the kohistan paleo island arc terrane north pakistan the dikes contain hornblende altered from diopsidic augite diopsidic augite relics ophitic to subophitic texture chlorite epidote sphene apatite zircon ilmenite titanomagnetite and magnetite the geochemistry reveals two groups of dikes 1 higher tio 2 2 74 3 50 wt na 2 o fe 2 o 3 and lower al 2 o 3 12 65 14 16 wt and mgo 3 73 5 04 wt 2 lower tio 2 1 24 2 05 wt na 2 o fe 2 o 3 and higher al 2 o 3 14 02 16 52 wt and mgo 3 98 7 52 wt the mgo contents 3 73 7 52 wt show a variation in the dikes from relatively primitive to more evolved compositions the dikes contain high amounts of both lile and hfse the major trace and rare earth elements data confirm the morb affinity and the back arc basin origin of the dike swarm the nw se orientation of the dike swarm and its 134 3 ma k ar age suggest the spreading axis of the back arc basin in the early cretaceous

indian and african plate motions driven by the push force of the reunion plume head mantle plumes are thought to play an important part in the earth s tectonics yet it has been difficult to isolate the effect that plumes have on plate motions here we analyse the plate motions involved in two apparently disparate events the unusually rapid motion of india between 67 and 52 million years ago and a contemporaneous transitory slowing of africa s motion and show that the events are coupled with the common element being the position of the indian and african plates relative to the location of the reunion plume head the synchroneity of these events suggests that they were both driven by the force of the reunion plume head the recognition of this plume force has substantial tectonic implications the speed up and slowdown of india the possible cessation of convergence between africa and eurasia in the palaeocene epoch and the enigmatic bends of the fracture zones on the southwest indian ridge can all be attributed to the reunion plume an analysis of the plate motions involved in two apparently disparate tectonic events the unusually rapid motion of india between 67 million and 52 million years ago and a transitory slowing of africa s motion at around the same time has uncovered a possible common cause for the two the reunion volcanic hotspot beneath the indian ocean has been active for about 67 million years and the close coupling of the motion of india and africa strongly suggests that it is the waxing and waning of the reunion plume head that is driving these motions

survival times of subducted slab remnants in numerical models of mantle flow abstract we examine the rate of sinking of truncated subducted slabs in a series of idealized three dimensional and two dimensional simulations using a 3d cartesian numerical convection model our goal is to account for the inference from seismic tomographic studies that some slab remnants have persisted at mid mantle depths for 150 myr or more while others have descended to the base of the mantle where they have spread laterally for great distances along the core mantle boundary we consider the influence of initial slab dimensions mantle viscosity stratification slab dip angle and three dimensionality on the slab survival time the time required to clear the upper half of the mantle of all measurable traces of a previous subduction zone we find slab survival time is determined by two principal factors the initial slab length and the viscosity stratification of the mantle survival time decreases with increasing slab length and increases with increasing viscosity stratification for a viscosity jump by a factor of 100 at a depth of 660 km shallow slab remnants remain in the upper 1500 km of the mantle for times well in excess of 150 myr while deep slabs draw their upper portions out of the upper half of the mantle in about 150 myr the dependence of the age of slab remnants on both the total depth extent of the original slab and the degree of viscosity stratification has not been studied previously it can account for the different rates of descent of slabs lying in close geographic proximity

palaeomagnetism of late cretaceous sediments from southern tibet evidence for the consistent palaeolatitudes of the southern margin of eurasia prior to the collision with india abstract twenty sites were drilled in the late cretaceous shexing formation for palaeomagnetic studies in the lhasa terrane near the locality of maxiang 29 9 n 90 7 e the stepwise thermal demagnetizations successfully isolated high unblocking temperature characteristic directions the tilt corrected mean direction is d i 350 8 32 1 with 95 8 1 and n 20 sites corresponding to a paleopole at 75 0 n 306 7 e with a 95 6 8 positive fold tests indicate a primary origin for the characteristic remanence based on previous cretaceous data mainly from the takena formation and paleocene data from the linzizong volcanic rocks near the city of lhasa the latitude of the southern margin of asia is located at about 15 n and yields a stable position of the lhasa terrane during cretaceous and paleocene compared with expected paleomagnetic directions from the stable india and eurasia blocks the collision palaeolatitude further implies the total latitudinal convergence was accommodated by 1700 800 km 16 2 7 6 between southern tibet and stable eurasia and 1500 830 km 14 4 7 9 between southern tibet and stable india since the collision of india and eurasia a collision age between c 54 and 47 ma was determined using the results for the southern margin of eurasia according to our new data and the extent of greater india

using seismic waves to image earth s internal structure seismic waves generated in earth s interior provide images that help us to better understand the pattern of mantle convection that drives plate motions

towards absolute plate motions constrained by lower mantle slab remnants the movement of tectonic plates before the cretaceous period is poorly understood a global mantle tomography model suggests that the longitude of oceanic subduction zones up to 300 million years ago was offset by up to 18 compared with tectonic reconstructions for the same period

new paleomagnetic results from the lhasa block revised estimation of latitudinal shortening across tibet and implications for dating the india asia collision knowing how and when the tibetan plateau formed is of global concern as high topography modifies climate circulation patterns and alters species migration paths previous paleomagnetic work suggested that intracontinental shortening from deformation related to the india asia collision played a major role in building the tibetan plateau however these studies were derived from red bed facies which are often afflicted by inclination shallowing leading to overestimates of continental shortening we report new paleomagnetic data from the lhasa block the southern continental margin of asia before the india asia collision a high temperature characteristic remanent magnetization chrm was isolated from 43 sites of late cretaceous red beds 32 sites of late cretaceous lava flows and nine sites of eocene tuffs all the chrm directions are significantly different from the present day geomagnetic field and the mean directions of red beds lava flows and tuffs pass the fold test the distribution patterns of chrm directions derived from volcanic rocks are consistent with a secular variation model suggesting that secular variation was likely averaged out in contrast the red bed chrm directions exhibit an east west elongated distribution pattern indicative of inclination flattening the mean chrm inclination of the red beds 23 5 degrees 2 5 degrees is consistent with previous results but significantly shallower than that of 41 9 degrees 4 4 degrees carried by the coeval lava flows an elongation inclination analysis of the red bed chrm directions yielded an unflattened mean inclination of 39 7 degrees with 95 confidence limits between 37 7 degrees and 41 9 degrees the paleolatitudes of the late cretaceous volcanic rocks and inclination corrected sediments suggest about 10 degrees less shortening occurred between lhasa and siberia than previous estimates with total shortening lying between 3 3 degrees 2 6 degrees and 7 3 degrees 5 4 degrees since the late cretaceous the mean chrm inclination derived from the eocene tuffs 51 8 degrees 5 2 degrees is indistinguishable from the expected reference direction our new results imply that india asia convergence was accommodated mainly through consumption of the india plate and yield a much younger age ca 43 ma for the initial india asia collision than previously estimated by paleomagnetic data c 2010 elsevier b v all rights reserved

4 d evolution of se asia s mantle from geological reconstructions and seismic tomography abstract how the collision between india and asia is related to processes deeper in the mantle is unclear here we compare geological reconstructions of block motions within asia since 50 ma with the tomographically imaged three dimensional 3 d morphology of subducted lithosphere to obtain insight into the spatiotemporal evolution of mantle structure past positions of the convergent margin show remarkable similarities with slab geometry at specific depths the striking change in slab geometry from a linear structure beneath 1100 km to an increasingly distorted shape at depths of less than 700 km results from collision the slab contours match the progressive deformation of asia s margin including india s indentation and sundaland s extrusion ever since the onset of collision the indian plate appears to have overridden its own sinking mantle and it does not seem at present to underthrust tibet significantly north of the zangbo suture if correct this observation would provide further evidence against models of plateau build up involving indian lithosphere the tomographic images beneath india confirm that asian deformation has absorbed at least 1500 km of convergence since collision began from the match between the southeastward motion of sundaland between 40 and 20 ma and the principal change in slab structure between 700 and 1100 km depths we infer that lateral advection in the mantle is small and that the sinking rate beneath sunda was 2 cm yr in the lower mantle and 5 cm yr above the transition zone

inducing in situ nonlinear soil response applying an active source 1 it is well known that soil sites have a profound effect on ground motion during large earthquakes the complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source to our knowledge we are the first group to propose and test such an approach employing a large surface vibrator as a source we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra we apply a homodyne algorithm for measuring spectral amplitudes which provides robust signal to noise ratios at the frequencies of interest spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source providing the means to separate source and near source nonlinearity from pervasive nonlinearity in the soil column we find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios corresponding to material softening with amplitude observed across the array as the source amplitude is increased the observed peak shifts are consistent with laboratory measurements of soil nonlinearity our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes

reconstruction of the deformed collision zone between india and asia by backward motion of lithospheric blocks 1 on the basis of a synthesis of tectonic data available on the india asia collision we present a first attempt to reconstruct the evolution of the collision zone assuming that the deformation of the lithosphere is localized along narrow shear zones and that the interiors of mantle blocks in between remain relatively undeformed we define block contours from the fault pattern and move back the blocks along their boundary faults along convergent or extensional boundaries the crust is assumed to shorten or stretch coherently step by step we go backward in time to finally reach the collision onset for each time step we find a solution compatible with the data set available and the position of the adjacent blocks for each block the search for compatibility at the scale of the entire collision zone allows for solving the kinematics of regions with fewer data and suggests plausible scenarios for regions where data is lacking for each step we calculate large scale displacement maps and determine euler poles for each block for the most recent time step the map proposed is compared to gps motions the deformation budget implies that extrusion absorbed 30 of the convergence between india and siberia during the entire collision span but varied with time accounting for as little as 3 or as much as 60 of this convergence at different epochs

seismic evidence of nonlinear crustal deformation during a large slow slip event in mexico repeated cross correlations of ambient seismic noise indicate a long term seismic velocity change associated with the 2006 m7 5 slow slip event sse in the guerrero region mexico because the sse does not radiate seismic waves the measured velocity change cannot be associated with the response of superficial soil layers to strong shaking as observed for regular earthquakes the perturbation observed maximized at periods between 7 s and 17 s which correspond to surface waves with sensitivity to the upper and middle crust the amplitude of the relative velocity change 10 3 was much larger than the volumetric deformation 10 6 at the depths probed 5 20 km moreover the time dependence of the velocity perturbation indicated that it was related to the strain rate rather than the strain itself this suggests that during strong slow slip events the deformation of the overlying crust shows significant nonlinear elastic behavior

insights on the kinematics of the india eurasia collision from global geodynamic models the eocene india eurasia collision is a first order tectonic event whose nature and chronology remains controversial we test two end member collision scenarios using coupled global plate motion subduction models the first conventional model invokes a continental collision soon after 60 ma between a maximum extent greater india and an andean style eurasian margin the alternative scenario involves a collision between a minimum extent greater india and a neotethyan back arc at 60 ma that is subsequently subducted along southern lhasa at an andean style margin culminating with continent continent contact at 40 ma our numerical models suggest the conventional scenario does not adequately reproduce mantle structure related to tethyan convergence the alternative scenario better reproduces the discrete slab volumes and their lateral and vertical distribution in the mantle and is also supported by the distribution of ophiolites indicative of tethyan intraoceanic subduction magmatic gaps along southern lhasa and a two stage slowdown of india our models show a strong component of southward mantle return flow for the tethyan region suggesting that the common assumption of near vertical slab sinking is an oversimplification with significant consequences for interpretations of seismic tomography in the context of subduction reference frames

local atomic structure of fontainebleau sandstone evidence for an amorphous phase 1 rocks are traditionally viewed as crystallite aggregates in general the treatment of the composition and mechanical properties never includes solid non crystalline components such glass like materials are difficult to detect in small quantities by standard techniques e g thin section polarimetry x ray diffraction but if they exist at critical locations i e grain contacts they could affect the behavior of rocks considerably neutron scattering measurements on a solid sample of fontainebleau sandstone have shown clear evidence for the presence of an unexpected glass like component atomic pair distribution function pdf analysis shows significant local structural deviations from pure quartz these deviations appear as an excess of 5 10 of nearest neighbor nn si o and o o bonds which is consistent with a 5 10 volume fraction of vitreous silica these measurements may provide significant information about the still unexplained causes for the peculiar mechanics and dynamics of sedimentary rocks

episodicity in back arc tectonic regimes abstract the evolution of back arc basins is tied to the development of the dynamics of the subduction system they are a part of we present a study of back arc basins and model their development by implementing 3d time dependant computer models of subduction including an overriding plate we define three types of episodicity pseudo quasi and hyper episodicity and find evidence of these in nature observations of back arc basin ages histories of spreading quiescence and compression in the overriding plate give us an understanding of the time development of these subduction zones and back arc basins across the globe today a number of trenches are advancing the izu bonin trench the mariana trench the japan trench the java sunda trench and the central portion of the peru chile trench the andes subduction zone the izu bonin mariana and japan all have established back arc basins while the others have documented episodes of spreading quiescence compression or a combination of these the combination of advancing and retreating trench motion places these subduction zones in the category of hyper episodicity quasi episodicity in which the back arc shifts between phases of rifting spreading and quiescence is the dominant form of episodic back arc development in the present we find this type of episodicity in models for which the system is dynamically consistent that we have allowed the subducting plate s velocity to be determined by the sinking slabs buoyancy quasi and hyper episodicity are only found in subduction zones with relatively high subducting plate velocities between 6 and 9 cm year finally those subduction zones for which the subducting plate is moving slowly such as in the mediterranean or the scotia sea experience only pseudo episodicity where the spreading moves linearly towards the trench but often does so in discrete ridge jump events

multicyclic history of the northern india continental margin northwestern himalaya 1 the geologic evolution of northern india is best recorded in the stratigraphic succession of the zanskar range northwestern himalaya which represents the most complete transect through this ancient continental margin sedimentary history began in the late proterozoic and recorded a late pan african orogenic event around the cambrian ordovician boundary when the gondwana supercontinent was eventually assembled the following long period of epicontinental deposition in shallow seas linked to palaeo tethys lasted until the early permian when a neo tethyan rift began to open between paleo india and the cimmerian microcontinents neo tethyan history can be subdivided into two sedimentary megasequences both recording a major tectonic and magmatic event in the lower part the first one began with breakup in the late permian and lasted until the end of the jurassic the second one started in the early cretaceous with the final detachment of india from gondwana and the opening of the indian ocean and ended with the india eurasia collision in the early eocene the two megasequences can be in turn subdivided into six transgressive regressive supersequences bounded by tectonically enhanced unconformities basal sandstone units of early permian late permian norian callovian early cretaceous and paleocene age are invariably associated with oolitic ironstones or reworked glauco phosphorites and mark the transgre sive part of each supersequence next condensed nodular carbonates or shales with pelagic fauna are typically overlain by thick shallowing upward marly units capped by regressive platformal carbonates the six tectono eustatic supercycles reflect successive rifting episodes which punctuated the progressive separation of india from the rest of gondwana and document the combination of plate microplate reorganizations and eustatic climatic and oceanographic changes in the tethyan domain after the onset of collision between india and asia close to the paleocene eocene boundary obduction of the remnants of the neo tethys ocean floor onto the indian margin began and the latter underwent multiphase deformation with fold thrust shortening followed by heating and extension after the main metamorphic event ophiolitic nappes were re thrusted and finally emplaced with their sedimentary sole on top of the passive margin succession

seismic velocity changes strain rate and non volcanic tremors during the 2009 2010 slow slip event in guerrero mexico summary we use ambient noise cross correlations to monitor small but reliable changes in seismic velocities and to analyse non volcanic tremor nvt intensities during the slow slip event sse that occurred in 2009 and 2010 in guerrero we test the sensitivity of the seismic velocity to strain variations in absence of strong motions the 2009 2010 sse presents a complex slip sequence with two subevents occurring in two different portions of the fault from a seismic array of 59 seismometers installed in small antennas we detect a velocity drop with maximum amplitude at the time of thefirst subevent we analyse the velocity change at different period bands and observe that the velocity perturbation associated with the sse maximizes for periods longer than 12 s then a linearized inversion of the velocity change measured at different period bands is applied in order to determine the depth of the portion of the crust affected by this perturbation no velocity change in the first 10km is detected below the velocity perturbation increases with depth affecting the middle and lower crust finally we compute the transient deformation produced by the sse in an elastic model using the slip evolution recovered from the inversion of continuous gps the comparison between the velocity changes and the deformation suggests that the velocity change is correlated with the strain rate rather than with the strain this result is similar to what was observed during the 2006 sse in the same region and suggests a non linear behaviour of the crust the velocity changes can be interpreted together with other observables such as nvts during the 2009 2010 sse we measure nvt activity using continuous seismic records filtered between 2 and 8 hz we observe a correlation between velocity changes for period band greater than 14 s and tremor activity whereas no correlation exists between velocity changes and seismic noise energy measured at long periods these observations suggest that both seismic velocity change and nvt can be used as indication of transient deformation at depth

memory relaxation and microfracturing in dilatant rock when palisades diabase and westerly granite are made dilatant by cyclic differential stress not exceeding 85 of fracture strength they behave as nonelastic energy absorbing media however in some ways their properties show reversibility more typical of elastic media the quasi elastic behavior is characterized by the property memory memory in dilatant rock is manifested by the closure of material property hysteresis loops when the sample is returned to certain previous differential stress states simple examples of memory have been observed before in westerly granite the present work confirms the existence of memory for strain and compressional wave velocity in westerly granite and palisades diabase at confining pressures between 250 and 1500 bars memory of maxima or minima of the stress difference can be developed with multiple points remembered simultaneously a theory of memory based on reversible tensile griffith cracks was developed which predicts a variety of phenomena many of which have been observed a property that is not predicted is the healing which occurs when unloading is stopped and the load is held constant a time dependent increase in velocity and decrease in sample volume is observed consistent with the healing of cracks relaxation ceases after approximately 2 min after relaxation a further decrease in load produces an initial elastic response followed by a return to the previous nonelastic behavior perfect memory is lost when microfracturing occurs as a result of either exceeding the previous absolute maximum stress or extensive cycling acoustic emissions are observed for conditions where memory is not perfect and are not observed when memory is perfect using travel time data a new low value was found for the dependence of dilatancy onset on pressure

closing the gap between regional and global travel time tomography recent global travel time tomography studies by zhou 1996 and van der hilst et al 1997 have been performed with cell parameterizations of the order of those frequently used in regional tomography studies i e with cell sizes of 1 2 these new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations 5 cells the inferred structures are however of larger scale than is usually obtained in regional models and it is not clear where and if individual cells are actually resolved this study aims at resolving lateral heterogeneity on scales as small as 0 6 in the upper mantle and 1 2 3 in the lower mantle this allows for the adequate mapping of expected small scale structures induced by for example lithosphere subduction deep mantle upwellings and mid ocean ridges there are three major contributions that allow for this advancement first we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the earth s mantle by seismic waves spakman and bijwaard 1998 second we exploit the global data set of engdahl et al 1998 which is a reprocessed version of the global data set of the international seismological centre their reprocessing included hypocenter redetermination and phase reidentification finally we combine all data used p pp and pwp phases into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal to noise ratio is improved in the approximate solution of the huge inverse problem we obtain a variance reduction of 57 1 synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells 0 6 or 1 2 in the shallow parts of subduction zones decreasing to approximately 2 3 resolution in well sampled regions in the lower mantle vertical resolution can be worse up to several hundreds of kilometers in subduction zones with rays predominantly pointing along dip important features of the solution are as follows 100 200 km thick high velocity slabs beneath all major subduction zones sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle large high velocity anomalies in the lower mantle that have been attributed to subduction of the tethys ocean and the farallon plate and low velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under iceland east africa the canary islands yellowstone and the society islands our findings corroborate that the 660 km boundary may resist but not prevent present day large scale mass transfer from upper to lower mantle or vice versa this observation confirms the results of previous global mantle studies that employed coarser parameterizations

geochronology of the initiation and displacement of the altyn strike slip fault western china abstract 40ar 39ar dating studies have been carried out along the dangjin pass transect across the altyn strike slip fault assf the samples gave ages of 445 2 454 3 ma in the northern belt 164 3 178 4 ma in the mesozoic shear zone and 26 3 36 4 ma in the cenozoic shear zone using the piercing point of the bashikaogong fault and the cangma heihe fault an offset of 350 400 km along the assf has been estimated the 40ar 39ar dating of the syntectonic growth or syntectonic resetting minerals from the samples within the assf belt and offset estimations from different age piercing points suggest that the assf should be initiated in the middle jurassic 178 4 160 ma combined with previously reported ages our studies show that the assf is characterized by multi phase re activation during 85 100 25 40 and 8 10 ma following its initiation in the middle jurassic in the regional tectonic setting of convergence between the indian and eurasian continents

acceleration and deceleration of india asia convergence since the cretaceous roles of mantle plumes and continental collision 1 a strong 50 35 ma decrease in india asia convergence is generally ascribed to continent continent collision however a convergence rate increase of similar magnitude occurred between 65 50 ma an earlier increase occurred at 90 ma both episodes of accelerated convergence followed upon arrival of a mantle plume below and emplacement of a large igneous province lip on the indian plate we here first confirm these convergence rate trends reassessing the indo atlantic plate circuits then using two different numerical models we assess whether plume head arrival and its lateral asthenospheric flow may explain the plate velocity increases and whether decreased plume flux and increasing continent plume distance may explain deceleration even without continental collision the results show that plume head arrival can indeed lead to absolute indian plate motion accelerations on the order of several cm yr followed by decelerations on timescales similar to the reconstructed fluctuations the 90 ma increase could potentially be explained as response to the morondova mantle plume alone the 65 50 ma convergence rate increase however is larger than can be explained by plume head spreading alone we concur with previous hypotheses that plume induced weakening of the indian continental lithosphere asthenosphere coupling and an increased slab pull and ridge push efficiency are the most likely explanations for the large convergence rate increase the post 50 ma decrease is best explained by orogeny related increased trench resistivity decreased slab pull due to continental subduction and possibly restrengthening of lithosphere asthenosphere coupling upon plume demise

a geodynamic model of mantle density heterogeneity using cenozoic and mesozoic plate motion reconstructions we derive a model of present day mantle density heterogeneity under the assumption that subducted slabs sink vertically into the mantle the thermal buoyancy of these slabs is estimated from the observed thermal subsidence cooling of oceanic lithosphere slab velocities in the upper mantle are computed from the local convergence rate we assume that slabs cross the upper lower mantle interface and continue sinking into the lower mantle with a reduced velocity for a velocity reduction factor between 2 and 5 our slab heterogeneity model is as correlated with current tomographic models as these models are correlated with each other we have also computed a synthetic geoid from our density model for a viscosity increase of about a factor of 40 from the upper to lower mantle our model predicts the first 8 spherical harmonic degrees of the geoid with statistical confidence larger than 95 and explains 84 of the observed geoid assuming that the model c21 and s21 terms are absent due to a long relaxation time for earth s rotational bulge otherwise 73 of the geoid variance is explained the viscosity increase is consistent with our velocity reduction factor for slabs entering the lower mantle since downwelling velocities are expected to scale roughly as the logarithm of viscosity loge 40 3 7 these results show that the history of plate tectonics can explain the main features of the present day structure of the mantle the dynamic topography induced by this heterogeneity structure consists mainly of about 1 km amplitude lows concentrated along the active continental margins of the pacific basin our model can also be used to predict the time variation of mantle heterogeneity and the gravity field we find that the age of the geoid defined as the time in the past before which the geoid becomes uncorrelated with the present geoid is about 50 m y our model for the history of the degree 2 geoid which is equivalent to the history of the inertia tensor should give us a tool to study the variations in earth s rotation pole indicated in paleomagiietic studies

the closing of tethys and the tectonics of the himalaya recent geological and geophysical data from southern tibet allow refinement of models for the closing of southern neo tethys and formation of the himalaya shelf sediments of the indian passive continental margin which pass northward into deep sea tethyan sediments of the indus tsangpo suture zone were deposited in the late cretaceous an andean type margin with a 2 500 km long trans himalayan kohistan ladakh gangdese granitoid batholith formed parallel to the southern margin of the lhasa block together with extensive andesites rhyolites and ignimbrites lingzizong formation the southern part of the lhasa block was uplifted deformed and eroded between the cenomanian and the eocene in the western himalaya the kohistan island arc became accreted to the northern plate at this time the northern part of the lhasa block was affected by jurassic metamorphism and plutonism associated with the mid jurassic closure of the bangong nujiang suture zone to the north the timing of collision between the two continental plates ca 50 40 ma marking the closing of tethys is shown by 1 the change from marine flysch like to continental molasse like sedimentation in the indus tsangpo suture zone 2 the end of gangdese i type granitoid injection 3 eocene s type anatectic granites and migmatites in the lhasa block and 4 the start of compressional tectonics in the tibetan tethys and indus tsangpo suture zone south facing folds south directed thrusts after the eocene closure of tethys deformation spread southward across the tibetan tethys zone to the high himalaya deep crustal thrusting barrovian metamorphism migmatization and generation of oligocene miocene leucogranites were accompanied by south verging recumbent nappes inverting metamorphic isograds and by south directed intracontinental shear zones associated with the main central thrust continued convergence in the late tertiary resulted in large scale north directed backthrusting along the indus tsangpo suture zone more than 500 km shortening is recorded in the foreland thrust zones of the indian plate south of the suture and 150 km shortening is recorded across the indian shelf zanskar range and the indus suture in ladakh there was also large scale shortening of the karakoram and tibetan microplates north of the suture as much as 1 000 km shortening occurred in tibet the more recent deformation however involved the spreading of this thickened crust and the lateral motion of the tibetan block along major approximately east west trending strike slip fault zones

temporal changes of shallow seismic velocity around the karadere duzce branch of the north anatolian fault and strong ground motion we analyze temporal variations of seismic velocity along the karadere duzce branch of the north anatolian fault using seismograms generated by repeating earthquake clusters in the aftershock zones of the 1999 mw7 4 izmit and mw7 1 duzce earthquakes the analysis employs 36 sets of highly repeating earthquakes each containing 4 18 events the events in each cluster are relocated by detailed multi step analysis and are likely to rupture approximately the same fault patch at different times the decay rates of the repeating events in individual clusters are compatible with the omori s law for the decay rate of regional aftershocks a sliding window waveform cross correlation technique is used to measure travel time differences and evolving decorrelation in waveforms generated by each set of the repeating events we find clear step like delays in the direct s and early s coda waves sharp seismic velocity reduction immediately after the duzce main shock followed by gradual logarithmic type recoveries a gradual increase of seismic velocities is also observed before the duzce main shock probably reflecting post seismic recovery from the earlier izmit main shock the temporal behavior is similar at each station for clusters at various source locations indicating that the temporal changes of material properties occur in the top most portion of the crust the effects are most prominent at stations situated in the immediate vicinity of the recently ruptured fault zones and generally decrease with normal distance from the fault a strong correlation between the co seismic delays and intensities of the strong ground motion generated by the duzce main shock implies that the radiated seismic waves produced the velocity reductions in the shallow material

nonlinear dynamical triggering of slow slip on simulated earthquake faults with implications to earth 1 among the most fascinating recent discoveries in seismology are the phenomena of dynamically triggered fault slip including earthquakes tremor slow and silent slip during which little seismic energy is radiated and low frequency earthquakes dynamic triggering refers to the initiation of fault slip by a transient deformation perturbation most often in the form of passing seismic waves determining the frictional constitutive laws and the physical mechanism s governing triggered faulting is extremely challenging because slip nucleation depths for tectonic faults cannot be probed directly of the spectrum of slip behaviors triggered slow slip is particularly difficult to characterize due to the absence of significant seismic radiation implying mechanical conditions different from triggered earthquakes slow slip is often accompanied by nonvolcanic tremor in close spatial and temporal proximity the causal relationship between them has implications for the properties and physics governing the fault slip behavior we are characterizing the physical controls of triggered slow slip via laboratory experiments using sheared granular media to simulate fault gouge granular rock and glass beads are sheared under constant normal stress while subjected to transient stress perturbation by acoustic waves here we describe experiments with glass beads showing that slow and silent slip can be dynamically triggered on laboratory faults by ultrasonic waves the laboratory triggering may take place during stable sliding constant friction and slip velocity and or early in the slip cycle during unstable sliding stick slip experimental evidence indicates that the nonlinear dynamical response of the gouge material is responsible for the triggered slow slip

the tidally averaged momentum balance in a partially and periodically stratified estuary abstract observations of turbulent stresses and mean velocities over an entire spring neap cycle are used to evaluate the dynamics of tidally averaged flows in a partially stratified estuarine channel in a depth averaged sense the net flow in this channel is up estuary due to interaction of tidal forcing with the geometry of the larger basin the depth variable tidally averaged flow has the form of an estuarine exchange flow downstream at the surface upstream at depth and varies in response to the neap spring transition the weakening of the tidally averaged exchange during the spring tides appears to be a result of decreased stratification on the tidal time scale rather than changes in bed stress the dynamics of the estuarine exchange flow are defined by a balance between the vertical divergence of the tidally averaged turbulent stress and the tidally averaged pressure gradient in the lower water column in the upper water column tidal stresses are important contributors particularly during the n

latest marine horizon north of qomolangma mt everest implications for closure of tethys seaway and collision tectonics the newly discovered marine horizon at zhepure syncline north of qomolangma revises timing of the closure of the eastern neo tethys seaway in the central himalayas the marine pengqu formation conformably overlies a lutetian shallow water carbonate platform and comprises shales interbedded with sandstones deposited in a neritic shelf environment the strata are dated by nannofossils and foraminifera as late early lutetian to late priabonian age np15 np20 deposited 47 34 ma indicating that the final closure of the tethys seaway in this region occurred at 34 ma the newly discovered strata provide evidence about emplacement of a major thrust sheet south of southern tibet prior to 37 ma which affected the regional climate

structure and evolution of the himalaya tibet orogenic belt the 1981 french chinese expedition to tibet focused on the lhasa block extending earlier coverage 400 km north of the tsangpo suture the lhasa block stood between 10 and 15 n latitude over most of the upper cretaceous and eocene and if gondwanian in origin had detached from gondwana by early permian seismic profiles reveal a complex moho topography resulting both from multiple continental thrusting and large scale strike slip faulting subduction related granitoids representing mixtures of mantle and crustal components and anatectic granitoids have been analysed and dated this study emphasizes the role of smaller blocks in the accretion of the continental mosaic

on lagrangian residual currents with applications in south san francisco bay california the lagrangian residual circulation has often been introduced as the sum of the eulerian residual circulation and the stokes drift unfortunately this definition of the lagrangian residual circulation is conceptually incorrect because both the eulerian residual circulation and the stokes drift are eulerian variables in this paper a classification of various residual variables are reviewed and properly defined the lagrangian residual circulation is then studied by means of a two stage formulation of a computer model the tidal circulation is first computed in a conventional eulerian way and then the lagrangian residual circulation is determined by a method patterned after the method of markers and cells to demonstrate properties of the lagrangian residual circulation application of this approach in south san francisco bay california is considered with the aid of the model results properties of the eulerian and lagrangian residual circulation are examined it can be concluded that estimation of the lagrangian residual circulation from eulerian data may lead to unacceptable error particularly in a tidal estuary where the tidal excursion is of the same order of magnitude as the length scale of the basin a direction calculation of the lagrangian residual circulation must be made and has been shown to be feasible

constraints to the timing of india eurasia collision a re evaluation of evidence from the indus basin sedimentary rocks of the indus tsangpo suture zone ladakh india deposited within the indus tsangpo suture zone the cenozoic indus basin sedimentary rocks have been interpreted to hold evidence that may constrain the timing of india eurasia collision a conclusion challenged by data presented here the eurasian derived 50 8 51 ma chogdo formation was previously considered to overlie indian plate marine sedimentary rocks in sedimentary contact thus constraining the timing of collision as having occurred by this time using isotopic analysis u pb dating on detrital zircons ar ar dating on detrital white mica sm nd analyses on detrital apatite sandstone and conglomerate petrography mudstone geochemistry facies analysis and geological mapping to characterize and correlate the formations of the indus basin sedimentary rocks we review the nature of these contacts and the identification and correlation of the formations our results reveal that previously interpreted stratigraphic contacts identifying chogdo formation unconformably overlying indian plate sedimentary rocks are incorrect rather we suggest that the inaccuracy of previous interpretations is most likely a result of formation misidentification and thus cannot be used to constrain the timing of india asia collision

plate tectonic evolution of the southern margin of eurasia in the mesozoic and cenozoic abstract thirteen time interval maps were constructed which depict the triassic to neogene plate tectonic configuration paleogeography and general lithofacies of the southern margin of eurasia the aim of this paper is to provide an outline of the geodynamic evolution and position of the major tectonic elements of the area within a global framework the hercynian orogeny was completed by the collision of gondwana and laurussia whereas the tethys ocean formed the embayment between the eurasian and gondwanian branches of pangea during late triassic early jurassic times several microplates were sutured to the eurasian margin closing the paleotethys ocean a jurassic cretaceous north dipping subduction boundary was developed along this new continental margin south of the pontides transcaucasus and iranian plates the subduction zone trench pulling effect caused rifting creating the back arc basin of the greater caucasus proto south caspian sea which achieved its maximum width during the late cretaceous in the western tethys separation of eurasia from gondwana resulted in the formation of the ligurian penninic pieniny magura ocean alpine tethys as an extension of middle atlantic system and a part of the pangean breakup tectonic system during late jurassic early cretaceous times the outer carpathian rift developed the opening of the western black sea occurred by rifting and drifting of the western central pontides away from the moesian and scythian platforms of eurasia during the early cretaceous cenomanian the latest cretaceous paleogene was the time of the closure of the ligurian pieniny ocean adria alcapa terranes continued their northward movement during eocene early miocene times their oblique collision with the north european plate led to the development of the accretionary wedge of the outer carpathians and its foreland basin the formation of the west carpathian thrusts was completed by the miocene the thrust front was still propagating eastwards in the eastern carpathians during the late cretaceous the lesser caucasus sanandaj sirjan and makran plates were sutured to the iranian afghanistan plates in the caucasus caspian sea area a north dipping subduction zone jumped during paleogene to the scythian turan platform the shatski terrane moved northward closing the greater caucasus basin and opening the eastern black sea the south caspian underwent reorganization during oligocene neogene times the southwestern part of the south caspian basin was reopened while the northwestern part was gradually reduced in size the collision of india and the lut plate with eurasia caused the deformation of central asia and created a system of nw se wrench faults the remnants of jurassic cretaceous back arc systems oceanic and attenuated crust as well as tertiary oceanic and attenuated crust were locked between adjacent continental plates and orogenic systems

sinking slabs below fossil subduction zones abstract the fate of sinking lithospheric slabs following the cessation of subduction is examined using a two dimensional numerical model of mantle convection the model incorporates rigidly moving oceanic plates which recycle through the underlying viscous mantle and rigid continental plates which cannot be subducted due to their inherent buoyancy the model is used to study the effects of the cessation of subduction due to the accretion of a 4000 km wide block of continental crust to the margin of a large stationary continent e g the india asia collision we model a continental block representative of greater india that is carried along with an oceanic plate towards a subduction zone at the leading edge of a much larger stationary continent which represents asia prior to continental collision subducted oceanic plate material under rides the margin of the larger continent for several hundred kilometres before peeling away from the surface boundary layer and descending into the mantle as a cold narrow slab consequently once a moving continental block arrives at the continental margin and closes over the original subduction zone the surface location of the continental suture lies several hundred kilometres seaward of the location of the cold remnant of the sinking slab the ongoing convective flow in the subcontinental mantle adjusts to accommodate the change in mechanical boundary condition at the upper surface a step increase in viscosity at 660 km depth slows the rate of slab descent substantially but cannot account for mid mantle thermal anomalies attributed to subduction zones that became extinct more than 140 my ago the horizontal displacement of the slab remnants relative to the continental crustal suture depends on the relative motion of the fossil plate boundary at the surface and the underlying mantle over riding of an intra oceanic subduction zone or the indenting of a smaller subcontinent into a larger continental mass can cause the tomographically observed cold anomalies to be displaced laterally by thousands of kilometers from the location of the former converging plate boundary

geoframework coupling multiple models of mantle convection within a computational framework solver coupling can extend the capability of existing modeling software and provide a new venue to address previously intractable problems a software package has been developed to couple geophysical solvers demonstrating a method to accurately and efficiently solve multiscale geophysical problems with reengineered software using a computational framework pyre pyre is a modeling framework capable of handling all aspects of the specification and launching of numerical investigations we restructured and ported citcoms a finite element code for mantle convection into the pyre framework two citcoms solvers are coupled to investigate the interaction of a plume at high resolution with global mantle flow at low resolution a comparison of the coupled models with parameterized models demonstrates the accuracy and efficiency of the coupled models and illustrates the limitations and utility of parameterized models

joint seismic geodynamic and mineral physical constraints on three dimensional mantle heterogeneity implications for the relative importance of thermal versus compositional heterogeneity summary the joint interpretation of seismic and geodynamic data requires mineral physical parameters linking seismic velocity to density perturbations in the earth s mantle the most common approach is to link velocity and density through relative scaling or conversion factors r s dln dlnv s however the range of possible r s values remains large even when only considering thermal effects we directly test the validity of several proposed depthdependent conversion profiles developed from mineral physics studies for thermally varying properties of mantle materials the tests are conducted by simultaneously inverting shear wave traveltime data and a diverse suite of convection related constraints interpreted with viscous flow response functions of the mantle these geodynamic constraints are represented by surface spherical harmonic basis functions up to harmonic degree 16 and they consist of the global free air gravity field tectonic plate divergences dynamic surface topography and the excess ellipticity of the core mantle boundary cmb the tests yield an optimum 1 d thermal r s profile that is generally compatible with all considered data with the exception of the dynamic surface topography that is most sensitive to the shallow upper mantle this result is not surprising given that cratonic roots are known to be compositionally distinct from the surrounding ambient mantle moreover it is expected that the temperature dependence of the thermal r s in the upper mantle is significant due to the temperature dependence of seismic attenuation or q therefore a simple 1 d density velocity relationship is insufficient to address this problem we obtained independent estimates of the first order correction factors to the selected r s profile within the cratonic roots and in the ambient thermal upper mantle these correction factors defined as r s lnv s are highly negative within the cratons signifying considerable compositional buoyancy this result confirms that the negative buoyancy associated with the low temperatures in the cratons is significantly counteracted by the composition induced positive buoyancy resulting in near neutral buoyancy of the cratonic roots within the ambient upper mantle the correction factors are found to be positive which is consistent with the expected behaviour of the r s relationship in thermally varying upper mantle material we obtain a significantly greater reconciliation of the global gravity anomalies and dynamic surface topography when applying these laterally varying corrections compared to a simple 1 d r s relationship inversion for a fully 3 d r s relationship reveals secondary effects including additional compositional variation within the cratonic roots and the deep mantle superplume structures we estimate the relative magnitude of the thermal and compositional non thermal contributions to mantle density anomalies and conclude that thermal effects dominate shear wave and density heterogeneity throughout the non cratonic mantle we also demonstrate the potential pitfalls of scaling a purely seismically derived model to obtain density rather than performing a true joint inversion to obtain velocity and density simultaneously

reconstructing early himalayan tectonic evolution and paleogeography from tertiary foreland basin sedimentary rocks northern india the latest paleocene middle eocene subathu formation and the oligocene miocene dagshai and kasauli formations of the indian foreland basin record the early evolution of the himalayan fold thrust belt sandstone petrography of the subathu formation shows a predominantly recycled sedimentary source with a distinct ophiolitic and volcanic influence that was drastically reduced by the time of deposition of the dagshai formation sandstones in the dagshai and kasauli formations consist predominantly of metapelitic detritus the metamorphic grade of metapelitic lithic grains increases with time from dominantly very low grade at the base of the dagshai formation to dominantly low grade in the kasauli sandstones mudstone geochemistry documents the presence of a mafic ultramafic source during the time of deposition of the subathu formation that becomes significantly less important by the time of deposition of the dagshai formation compositions of subathu formation detrital spinels show they were either derived from both mid ocean ridge basalt type and arc type ophiolites or from an ophiolite of composite origin and that southerly derivation from the deccan trap continental flood basalts is unlikely kasauli formation garnet compositions suggest derivation from medium grade metamorphic rocks to amphibolite facies subathu formation composition reflects provenance influence from the indus suture zone during latest paleocene middle eocene time indicating initiation of continent continent collision and development of the foreland basin by this time suture zone influence was drastically reduced by the time of deposition of the dagshai formation when the embryonic thrust belt provided a barrier sufficient to partially separate the suture zone from the himalayan foreland basin the first appearance of himalayan metamorphic detritus occurs in the dagshai formation at the close of oligocene time whereas the kasauli formation records erosion to deeper metamorphic levels during earliest miocene time the occurrence of garnet and higher grade metamorphic lithic grains during early neogene time is coincident with the timing of displacement along the main central thrust and south tibetan detachment zone composition of early foreland basin sediments from pakistan balakot murree and kamlial formations to nepal bhainskati and dumri formations and bangladesh kopili barail formations surma group indicates diachronous arrival of ophiolitic to low grade metamorphic detritus derived respectively from the indus suture and early himalayan thrust sheets in the north this is consistent with progressively later closure of neotethys along the suture from latest paleocene time in the west to eocene time or even later in the east

resonant acoustic spectroscopy of microfracture in a westerly granite sample 1 the experimental setup for precision acoustic measurements of the elastic tensor of structurally inhomogeneous materials has been designed and fabricated the achieved accuracy of measurements is high enough for studying a fine structure of rocks it is shown that resonant acoustic spectroscopy is a useful tool for measuring the small content of oriented microcracks that can be used for diagnostics of early stages of fracture when oriented microcracking occurs elastic characteristics of a westerly granite sample have been investigated the elastic tensor has been found to be hexagonal oriented with respect to the rift plane the obtained data have shown the presence of oriented microcracks and have been explained sensibly by a model of microcracks

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