aftsolve a program for multi kinetic modeling of apatite fission track data aftsolve is a computer program for deriving thermal history information from apatite fission track data it implements a new fission track annealing model that takes into account the known kinetic variability among different apatite species to fully utilize this model a fission track worker must obtain data that can be used to infer the kinetic characteristics of each apatite grain from which a measurement was taken such data can consist of etch figure lengths or chemical composition the benefit of this overall approach is that it allows useful information to be derived from previously unusable analyses extends the practical range of geological temperatures constrained by fission track analyses and increases overall confidence in model predictions aftsolve also incorporates the effects of fission track orientation relative to the apatite crystallographic c axis variation in initial track length and the biasing effect of 2 5 2 cf irradiation for enhancing confined horizontal track length detection aftsolve is written for windows operating systems and has a graphical interface that allows interactive input of thermal histories and real time generation of estimates for fission track length distributions and ages for up to six simultaneously modeled kinetic populations it also includes procedures for estimating the range of time temperature histories that are statistically consistent with a data set and constraints entered by the user

indian punch rifts tibet the strikes of late cenozoic normal faults in tibet fan systematically from a northwesterly to northeasterly direction from west to east they exhibit an axis of symmetry that bisects the himalayan arc and trends subparallel to india asia motion the two dimensional stress field inferred from the spatial pattern of normal fault orientations can be explained solely by collisional stresses localized along a southern part of the himalayan arc continued insertion of indian crust into tibet can explain the persistence of extension and may be increasing plateau crustal thickness and elevation coeval with rifting

normal fault thermal regimes and the interpretation of low temperature thermochronometers abstract exhumation rates inferred from thermochronometers are dependent on spatial and temporal variations in temperature in active extensional mountain belts the temperature field is complicated by tectonic and surface processes including 1 lateral heat flow across large range bounding normal faults due to the juxtaposition of a cool hanging wall and a relatively warmer footwall 2 uplift and erosion of the footwall 3 sedimentation and burial of the hanging wall 4 lateral heat refraction around low thermal conductivity sediments deposited in the hanging wall basin and 5 3d temperature variations due to high relief topography developed on the footwall we explore these mechanisms through a series of 2d conductive thermal models designed to investigate the effect of tectonics and topography on apatite fission track aft and u th he thermochronometer data models were tuned to the geometry and kinematics of the wasatch mountains utah usa the principal parameters in our model are exhumation and burial rates ranging from 0 2 to 5 7 mm per year at the range front and decreasing with distance from the fault surface morphology taken from usgs digital elevation models and basin geometries inferred from seismic and gravity surveys predicted aft and u th he ages were generated using cooling rate dependent annealing and diffusion kinetic models results indicate after 10 million years of exhumation footwall u th he and aft closure temperature isotherms within 10 km of the fault are advected upward 500 and 1000 m respectively from their initial position the upward advection of isotherms and the 2d nature of the thermal regime can result in erroneous exhumation rates calculated from plots of sample elevation versus age using 1d thermal models for simulations with a uniform vertical uplift rate and canyon and ridge topography 1d and 2d exhumation rate differences were 20 70 for u th he and 10 for aft data samples collected perpendicular to fault strike up the range front are sensitive to the exhumation rate and footwall tilt differences between 1d and 2d range front exhumation rates were 10 95 for u th he and 10 40 for aft data furthermore 1d thermal models are incapable of deciphering footwall tilt at simulated exhumation rates of 3 4 and 5 mm per year predicted u th he and aft ages form closely spaced near vertical lines on a plot of elevation versus sample age therefore suggesting a decreased sensitivity of data to topography and surface processes at high exhumation rates despite large errors in previously collected aft data from the wasatch mountains our model suggests a revised fault adjacent exhumation rate of 0 5 0 6 mm per year compared to previous estimates of 0 7 mm per year based on 1d thermal models

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

lunar volatile depletion due to incomplete accretion within an impact generated disk the moon may have accreted from a disk of debris after a giant impact simulations suggest that part of the moon derives from volatile poor melt in the hot inner disk with most of the volatile elements condensing later and accreting to earth

transtensional deformation in the central himalaya and its role in accommodating growth of the himalayan orogen 1 field mapping structural analysis and geochronologic data from northwestern nepal reveal major normal right slip motion along a previously unrecognized west northwest striking system of shear zones that we term the gurla mandhata humla fault system gmh the gmh obliquely cuts across the greater himalayan crystalline sequence and into the lesser himalayan imbricate thrust belt via two right step over structures the average slip direction on the gmh parallels the strike of the himalayan orogen motion along this fault system has resulted in an apparent left separation of the south tibet detachment main central thrust zone and lesser himalayan imbricate thrust belt along a north striking segment of the fault system we estimate a minimum of 21 km of net slip on the southern branch of the gmh by restoring the trace of the main central thrust zone parallel to the average slip direction on the fault taking into account slip estimates from the northern branch of the gmh yields a minimum net slip estimate of 24 4 to 32 4 km for the gmh the 232th 208pb ion microprobe monazite ages from leucogranite bodies indicate that motion on the gmh occurred after 15 ma its initiation immediately followed crustal thickening between the main central thrust zone and indus yalu suture zone motion on the gmh is contemporaneous with arc normal contraction in the southernmost himalayan orogen these observations can be explained by a model that involves foreland propagating structural systems facilitating arc normal contraction in the foreland and arc parallel extension in the hinterland that work together to maintain the arcuate shape of the himalayan orogen

exhumation of the inyo mountains california implications for the timing of extension along the western boundary of the basin and range province and distribution of dextral fault slip rates across the eastern california shear zone 1 new geologic mapping tectonic geomorphologic 10be terrestrial cosmogenic nuclide and u th he zircon and apatite thermochronometric data provide the first numerical constraints on late cretaceous to late quaternary exhumation of the inyo mountains and vertical slip and horizontal extension rates across the eastern inyo fault zone california the east dipping eastern inyo fault zone bounds the eastern flank of the inyo mountains a prominent geomorphic feature within the western basin and range province and eastern california shear zone u th he zircon and apatite thermochronometry yield age patterns across the range that are interpreted as indicating 1 two episodes of moderate to rapid exhumation associated with laramide deformation during the late cretaceous early tertiary 2 development of a slowly eroding surface during a prolonged period from early eocene to middle miocene 3 rapid cooling exhumation and initiation of normal slip along the eastern inyo fault zone accommodated by westward tilting of the inyo mountains block at 15 6 ma and 4 rapid cooling exhumation and renewed normal slip along the eastern inyo fault zone at 2 8 ma fault slip continues today as indicated by fault scarps that cut late pleistocene alluvial fan surfaces the second episode of normal slip at 2 8 ma also signals onset of dextral slip along the hunter mountain fault yielding a pliocene dextral slip rate of 3 3 1 0 mm a where a is years summing this dextral slip rate with estimated dextral slip rates along the owens valley death valley and stateline faults yields a net geologic dextral slip rate across the eastern california shear zone of 9 3 2 2 1 4 to 9 8 1 4 1 0 mm a

high times on the tibetan plateau paleoelevation of the thakkhola graben nepal east west extension in the tibetan plateau is generally assumed to have resulted from gravitational collapse following thickening and uplift on the basis of this assumption several studies have dated east west extensional structures to determine when the plateau attained its current high elevation however independent estimates of elevation are needed to determine whether extension occurred before during or after the plateau achieved its current elevation because the isotopic composition of meteoric water decreases with increasing elevation significant change in local elevation throughout the thakkhola graben depositional history should be recorded by change in 18 o values of fluvial and lacustrine carbonates the 18 o values of 16 to 23 of thakkhola graben carbonates reflect meteoric water values similar to modern values and suggest that the southern tibetan plateau attained its current elevation prior to eastwest extension initiation of thakkhola graben extension is constrained between 10 and 11 ma based on magnetostratigraphy of the older tetang formation the 13 c values of soil carbonates suggest an age younger than 8 ma for the base of the thakkhola formation

how flat is tibet high resolution digital topography three arc second grid for most of tibet provides new information to characterize the relief of the highest and largest plateau on earth the arid to semiarid central and northern part of the plateau interior has low relief average slopes of over 250 m windows and a mean elevation of 5023 m above sea level at moderate wavelengths of m relief is or less for most of tibet as opposed to the much higher relief of up to 6 km on the plateau edges where glacial and fluvial dissection is greater because of higher levels of precipitation the only faults manifesting significant topographic relief are the relatively small scale generally north trending graben systems primarily in southern tibet and several large scale fault systems near the edges of tibet the flatness of tibet implies that 1 there has been little deformation especially shortening of the uppermost crust north of the graben systems during the late cenozoic and 2 shallow crustal isostatic compensation has been acting to level the surface of the plateau

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

causes and consequences of protracted melting of the mid crust exposed in the north himalayan antiform gneiss domes exposed by the north himalayan antiform in southern tibet provide evidence for protracted melting of the mid crust during convergent tectonics anatexis was initiated during underthrusting and recurred during the southwards extrusion of a low viscosity crustal sheet sr nd isotope systematics establish that the cores of the domes dominated by cambrian gneiss form a window of the high himalayan crystalline series hhcs indian plate combined trace element and isotopic analyses from 4 two mica granites that intrude both core and the metasedimentary mantle of the domes indicate an origin from melting of the high himalayan crystalline series by the incongruent breakdown of muscovite and to a lesser extent biotite at temperatures from 680 to 760 c u pb dating of accessory phases extends the known period of granite magmatism from the late oligocene to mid miocene a ca rich garnet bearing leucogranite of oligocene age 27 5 ma provides a rare insight into the pre miocene history of the underthrust indian plate reflecting prograde melting at greater depths during a period of crustal thickening the origin of the miocene granites 15 10 ma can be modeled from the exhumation of an extruding channel at a rate of ca 1 mm year constrained by ar isotope data melting over a protracted time scale along the leading edge of the indian plate less than 50 km from the suture accelerated horizontal shortening within the collision zone and provided a positive feedback mechanism between magmagenesis and deformation

middle to late miocene extremely rapid exhumation and thermal reequilibration in the kung co rift southern tibet 1 the kung co rift is an approximately nnw striking wsw dipping normal fault exposed in southern tibet and is part of an extensive network of active approximately ns striking normal faults exposed across the tibetan plateau detailed new and published u th he zircon and apatite thermochronometric data from the footwall of the early miocene kung co granite provide constraints on the middle miocene to present day exhumation history of the footwall to the kung co fault inverse modeling of thermochronometric data yield age patterns that are interpreted as indicating 1 initiation of normal fault slip at 12 13 ma and rapid exhumation of the footwall between 13 and 10 ma 2 acceleration of normal fault slip at rates of 21 9 6 9 mm yr at 10 ma 3 rapid thermal reequilibration between 10 and 9 ma and 4 slow exhumation and or quiescence from 9 ma to the present day hanging glacial valleys in the footwall and fault scarps that cut late quaternary till and moraine deposits indicate that fault slip continues today middle to late miocene initiation of extension across the kung co rift is broadly the same as the documented initiation of ew extension across the south central tibetan plateau eastward flow of middle or lower crust from beneath tibet accommodated by northward underthrusting of indian crust beneath tibet provides a plausible explanation for the onset of ew extension across the tibetan plateau

role of oblique convergence in the active deformation of the himalayas and southern tibet plateau noting similarities with subduction along curved oceanic trenches and using a simple block model we show that radial vergence evident in earthquake slip vectors along the himalayan deformation front east west extension on north trending normal faults in the himalayas and southern tibet and right lateral strike slip on the karakorum jiali fault zone can all result from basal shear caused by the indian plate sliding obliquely beneath tibet along a gently dipping arcuate plate boundary within the framework of this mechanism the normal faults in the himalayas and southern tibet are not proxies for the uplift history of tibet the distribution and style of the faults in the himalayas and southern tibet suggest that the basal drag from the underthrusting indian lithosphere extends northward beneath most of southern tibet

a new approach to stable isotope based paleoaltimetry implications for paleoaltimetry and paleohypsometry of the high himalaya since the late miocene the change in oxygen isotopic composition of precipitation is modeled using equilibrium fractionation during rayleigh distillation linked to the thermodynamics of atmospheric ascent and water vapor condensation the primary controllers of the vertical variation in isotopic composition with elevation are the low elevation temperature and relative humidity as these control the vertical distribution of condensation an empirical fit of precipitation versus model condensation based on alpine stations is derived this fit is represented in the model as the weighted mean composition of condensation within a 1000 m thick air parcel 1500 500 m above the ground surface and is used for all other regions comparison of model versus observed modern precipitation reveals a close fit particularly of more highly elevated sites comparison of modern waters in the himalayas and southern tibet with model predictions particularly as revealed by comparison of observed and predicted hypsometry provides additional support to the validity of the model finally application of this model to estimates of paleo waters in the himalayas and southern tibet reveals that this region had already achieved its present hypsometry by the late miocene about 10 ma ago fl 2001 published by elsevier science b v

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

modeling of the temperature sensitivity of the apatite u th he thermochronometer apatite u th he apparent ages will generally reflect residence for extended periods at temperatures where helium is neither quantitatively retained nor lost by diffusion to characterize the response of apatite he ages to thermal histories involving partial he retention we explored solutions to the he production diffusion equation under thermally static conditions the analytical solution to this equation coupled with published diffusivity data demonstrates that the zone of partial he retention extends from about 40 c to 85 c this zone lies at temperatures 35 c cooler than the analogous fission track partial annealing zone he ages within the partial retention zone ultimately achieve a balance between he production and loss yielding a steady state age both the ultimate age and the time it takes to achieve this age are temperature dependent for example an apatite held at 75 c equilibrates to an age of 2 ma after 17 myr regardless of whether equilibrium is approached from a higher or a lower initial he age for representative dynamic thermal histories we evaluated apatite he ages using a numerical solution to the ingrowth diffusion equation the results illustrate the sensitivity of he ages to various geologic histories and are useful for understanding he age elevation relationships and for testing time temperature paths derived from apatite fission track length distributions in addition although he diffuses rapidly from apatite at shallow crustal temperatures modeling of ambient temperature fluctuations indicates that he ages are nearly unaffected by surficial processes

when did the roof collapse late miocene north south extension in the high himalaya revealed by th pb monazite dating of the khula kangri granite th pb ion microprobe measurements made on 12 monazite grains from the khula kangri granite tibet bhutan frontier are interpreted to indicate that crystallization occurred at 12 5 0 4 ma the leucogranite is cut by the gonto la detachment part of the southern tibet detachment system that has allowed upper level north directed extension of the himalayan orogen significant orogen normal extension in southern tibet appears to have continued 8 10 m y later than previously recognized this is the first reported crystallization age for a leucogranite east of the yadong cross structure an apparent 70 km offset of the high himalaya and southern tibet detachment west of the yadong cross structure reliable ages for high himalaya events major main central thrust slip granite generation and emplacement attainment of critical topography and major detachment extension group between ca 24 and 19 ma we interpret the west to east change across the yadong cross structure to be due to either 1 an abrupt 10 m y younging of principal high himalayan events or 2 a deeper thus younger exposed part of the footwall of the southern tibet detachment near khula kangri the southern tibet detachment is cut by the highly oblique yadong gulu rift a manifestation of tibet plateau east west extension integrated estimates of magnitude and rate of detachment displacement suggest that the observed postcrystallization north directed extension lasted for 1 3 m y after which time the yadong gulu rift formed this interpretation is consistent with initiation of east west extension of tibet at ca 8 ma

subsurface geology ancient hydrothermal systems and crater excavation processes beneath lake rotomahana evidence from lithic clasts of the 1886 ad rotomahana pyroclastics abstract the craters associated with the 1886 ad phreatomagmatic rotomahana eruption okataina volcanic centre new zealand and the near vent geology are now hidden beneath lake rotomahana and its post eruptive sediment fill lithic clasts from the near vent lithic lapilli ash deposits of the rotomahana pyroclastics are used in this study to trace geological and geothermal conditions before the eruption as well as vent excavation dynamics near vent deposit characteristics were described in the field representative lithic clasts were documented petrographically and unaltered clasts were analysed for major and trace element compositions the majority of the lithic clasts were rhyolites with subordinate ignimbrites and hydrothermally altered clasts and trace siltstone and silicified clasts the rhyolites were classified into four petrographic groups according to phenocryst content and assemblage and were more diverse with respect to geochemical compositions most of the rhyolite lithics in the rotomahana pyroclastics did not match the rhyolite domes exposed subaerially around the lake but did have affinities with the pre matahina caldera wairua rhyolite and potentially other older non exposed domes ignimbrites most likely correlated either to the matahina ignimbrite or older non exposed units hydrothermally altered rhyolite and ignimbrite lithic clasts are common and suggest that there has been a long lived hydrothermal system in this sector possibly dating back to early activity of the okataina volcanic centre the diversity in lithic types indicate a spatial variation in country rock lithology and strength which probably contributed to the vent position and morphology along the rotomahana fissure

east west extension and miocene environmental change in the southern tibetan plateau thakkhola graben central nepal east west extensional basins are distributed across the southern half of the tibetan plateau at an elevation of 4 km these basins have generated much interest because of their potential implications for the regional tectonics and force distribution in the plateau this study documents the sedimentology of the miocene pliocene thakkhola graben fill in order to reconstruct basin evolution and paleoenvironment analysis of depositional systems paleodrainage patterns and conglomerate clast provenance of the 1 km thick graben fill sets limits on the timing of activity of the basin bounding faults and the development of southward axial drainage in the basin during the deposition of the oldest basin fill tetang formation ca 11 9 6 ma probably in a restricted basin minor motion occurred on the basin bounding fault systems an angular unconformity separates the tetang and overlying thakkhola formations where this contact can be observed in the southern part of the basin southward axial drainage was established by ca 7 ma with the onset of deposition of the thakkhola formation several episodes of damming of this drainage system are recorded by widespread lacustrine deposits in the southern part of the basin facies distribution and the progressive rotation of strata in the thakkhola formation indicate that the dangardzong fault on the western edge of the basin was active at this time and drainage ponding may have been related to displacement on normal faults associated with the south tibetan detachment system to the south of thakkhola graben contrasts between deposits of the tetang and thakkhola formations provide evidence for environmental change in the basin in the tetang formation the abundance of lacustrine facies the pollen record and the absence of paleosol carbonate suggest that conditions were more humid than during subsequent deposition of the thakkhola formation environmental change in the thakkhola graben coincided with environmental change observed in the siwalik foreland basin sequence arabian sea and bay of bengal at ca 8 7 ma although this climate change event has been previously attributed to intensification of the asian monsoon in response to uplift of the tibetan plateau paleoaltimetry data indicate that this region had already attained a high elevation by ca 11 ma thus the thakkhola graben stratigraphic record suggests that uplift of the southern tibetan plateau and the onset of the asian monsoon as inferred from paleoclimatic indicators were not directly related in a simple way

diffusion of helium in zircon and apatite abstract diffusion of helium has been characterized in natural zircon and apatite polished slabs of zircon and apatite oriented either normal or parallel to c were implanted with 100 kev 3he at a dose of 5 1015 3 he cm2 diffusion experiments on implanted zircon and apatite were run in pt capsules in 1 atm furnaces 3he distributions following experiments were measured with nuclear reaction analysis using the reaction 3he d p 4he for diffusion in zircon we obtain the following arrhenius relations d c 2 3 10 7 exp 146 11 kj mol 1 r t m 2 s 1 d c 1 7 10 5 exp 148 17 kj mol 1 r t m 2 s 1 although activation energies for diffusion normal and parallel to c are comparable there is marked diffusional anisotropy with diffusion parallel to c nearly 2 orders of magnitude faster than transport normal to c these diffusivities bracket the range of values determined for he diffusion in zircon in bulk release experiments although the role of anisotropy could not be directly evaluated in those measurements in apatite the following arrhenius relation was obtained over the temperature range of 148 449 c for diffusion normal to c d 2 10 10 6 exp 117 6 kj mol 1 r t m 2 s 1 in contrast to zircon apatite shows little evidence of anisotropy he diffusivities obtained in this study fall about an order of magnitude lower than diffusivities measured through bulk release of he through step heating and within an order of magnitude of determinations where ion implantation was used to introduce helium and he distributions measured with elastic recoil detection since the diffusion of he in zircon exhibits such pronounced anisotropy helium diffusional loss and closure cannot be modeled with simple spherical geometries and the assumption of isotropic diffusion a finite element code cylmod has recently been created to simulate diffusion in cylindrical geometry with differing radial and axial diffusion coefficients we present some applications of the code in evaluating helium lost from zircon grains as a function of grain size and length to diameter ratios and consider the effects of shape anisotropy where diffusion is isotropic as in the case of apatite but shapes of crystal grains or fragments may depart significantly from spherical geometry

the role of crustal strength variations in shaping orogenic plateaus with application to tibet 1 the tibetan plateau is the type example of a large orogenic plateau formed as a result of continent continent collision the morphology of the plateau and its margins suggests that preexisting variations in crustal strength have influenced the growth of the plateau we have developed a three dimensional numerical model of deformation in a viscous crust in order to investigate the effects of lateral heterogeneities on plateau growth the model includes a two layer crust and allows for lateral variation of viscosity in both the upper and lower layers model results indicate that crustal strength variations have a dramatic effect on the morphology and dynamics of a developing plateau a region of strong crust is characterized by a very steep plateau margin that propagates extremely slowly does not accommodate significant shortening strain and is subparallel to local upper and lower crustal velocities a weak crustal region develops a gently sloping margin uplift propagates rapidly across the weak zone and crustal material within the plateau is diverted toward the low strength region with a relatively simple distribution of strength variations corresponding to strong tarim and sichuan basin crust and a weak southeastern corner our model produces a plateau with many similarities to the tibetan plateau including the overall morphology rotation around the eastern syntaxis and e w extension analysis of model results suggests that e w extension of the central plateau may be related to the rapid eastward flow of crustal material into a weak zone in the se corner without significant change in plateau elevation

hydrogen and oxygen isotope variations in the high himalaya peraluminous manaslu leucogranite evidence for heterogeneous sedimentary source abstract the manaslu granite belongs to the high himalaya leucogranitic belt which was produced by melting of the crust during postcollisional thrusting d and 18o values have been determined for whole rock and coexisting minerals from the 8 km thick manaslu massif and its 50 km long dyke sheet its country rocks and the formation 1 f1 paragneisses which are the source of the granite for the granite dmusc range from 70 to 85 and 18o w r from 10 9 to 12 8 h and o isotope fractionations among minerals are consistent with high temperature equilibrium and for oxygen closed system evolution a few samples coming mainly from the dyke swarm have very low d values down to 188 and biotite muscovite h isotope fractionations indicative of disequilibrium the d h ratios of associated magmatic tourmaline are essentially unmodified from the distribution of d values in the granite and its country rocks circulation of very low deuterium meteoric hydrothermal waters was extremely localized because these waters are depleted in deuterium by up to 50 relative to modern meteoric waters the manaslu area was either at an altitude substantially higher than that of today 2500 6000 m for the analyzed samples or a mountain chain once existed to the south the f1 gneisses have 18oquartz between 12 and 14 3 which confirms that the granite was generated from f1 but d values are 20 higher than in the granite such a difference can be a result of degassing of the magma and or introduction of fluid in the melting zone infiltration of low d fluid 90 into the hot but dry f1 probably triggered partial melting these fluids could have come from the dehydration of the midlands sediments which are separated from the overlying f1 by the main central thrust the correlations among 18o 87 sr 86 sr 20 ma and nd values in both f1 and the granite indicate that the variations of these isotopic ratios in the manaslu are inherited from those in f1 at the time of melting in turn these ratios in f1 are related to the proportion of quartz and phyllosilicates for the isotopic ratios of nd and o and to the quantity of radiogenic sr generated within the sediment which is a function of age and rb content amount of phyllosilicate some other himalayan leucogranites require either other crustal source rocks or the 18o and 87 sr 86 sr ratios of f1 vary along the himalaya

late cenozoic uplift of southeastern tibet the age of surface uplift in southeastern tibet is currently unknown but the initiation of major river incision can be used as a proxy for the timing of initial uplift the topographically high eastern plateau and gently dipping southeastern plateau margin are mantled by an elevated low relief relict landscape that formed at a time of slow erosion at low elevation and low tectonic uplift rates prior to uplift of the eastern tibetan plateau thermochronology from deep river gorges that are cut into the relict landscape shows slow cooling between ca 100 and ca 10 20 ma and a change to rapid cooling after ca 13 ma with initiation of rapid river incision at 0 25 0 5 mm yr between 9 and 13 ma a rapid increase in mean elevation of eastern tibet beginning at this time supports tectonic climate models that correlate the lateral eastern expansion of high topography in tibet with the late miocene intensification of the indian and east asian monsoons

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

implications of shortening in the himalayan fold thrust belt for uplift of the tibetan plateau 1 recent research in the himalayan fold thrust belt provides two new sets of observations that are crucial to understanding the evolution of the himalayan tibetan orogenic system first u pb zircon ages and sm nd isotopic studies demonstrate that 1 greater himalayan medium to high grade metasedimentary rocks are much younger than true indian cratonic basement and 2 these rocks were tectonically mobilized and consolidated with the northern margin of gondwana during early paleozoic orogenic activity these observations require that greater himalayan rocks be treated as supracrustal material in restorations of the himalayan fold thrust belt rather than as indian cratonic basement in turn this implies the existence of greater himalayan lower crust that is not exposed anywhere in the fold thrust belt second a regional compilation of shortening estimates along the himalayan arc from pakistan to sikkim reveals that 1 total minimum shortening in the fold thrust belt is up to 670 km 2 total shortening is greatest in western nepal and northern india near the apex of the himalayan salient and 3 the amount of himalayan shortening is equal to the present width of the tibetan plateau measured in an arc normal direction north of the indus yalu suture zone this information suggests that a slab of greater indian lower crust composed of both indian cratonic lower crust and greater himalayan lower crust with a north south length of 700 km may have been inserted beneath the tibetan crust during the cenozoic orogeny we present a modified version of the crustal underthrusting model for himalayan tibetan orogenesis that integrates surface geological data recent results of mantle tomographic studies and broadband seismic studies of the crust and upper mantle beneath the tibetan plateau previous studies have shown that incremental mesozoic and early cenozoic shortening had probably thickened tibetan crust to 45 55 km before the onset of the main cenozoic orogenic event thus the insertion of a slab of greater indian lower crust with maximum thickness of 20 km tapering northward could explain the cenozoic uplift of the tibetan plateau the need for tibetan crust to stretch laterally as the greater indian lower crust was inserted may explain the widespread east west extension in the southern half of the plateau our reconstruction of the himalayan fold thrust belt suggests that indian cratonic lower crust of presumed mafic composition and high strength should extend approximately halfway across the tibetan plateau to the banggong suture from there northward we predict that the tibetan plateau is underlain by more felsic and therefore weaker lower crust of greater himalayan affinity two slab break off events are predicted by the model the first involved neotethyan oceanic lithosphere that foundered 45 35 ma and the second consisted of greater indian lithosphere most likely composed of greater himalayan material that delaminated and foundered 20 10 ma asthenospheric upwelling associated with the break off events may explain patterns of cenozoic volcanism on the tibetan plateau although the model predicts a northward migrating topographic front due solely to insertion of greater indian lower crust the actual uplift history of the plateau was complicated by early middle tertiary shortening of tibetan crust

geochemical evaluation of fenghuoshan group lacustrine carbonates north central tibet implications for the paleoaltimetry of the eocene tibetan plateau abstract a sedimentologic petrologic and geochemical evaluation of lacustrine carbonates from the eocene oligocene fenghuoshan group of north central tibet provides insight into the paleoenvironmental and paleolimnological setting of the hoh xil basin during the mid tertiary fenghuoshan lacustrine rocks consist primarily of carbonate mudstones and wackestones these carbonates are generally less than 3 m in thickness and are intercalated with siliciclastic lacustrine and fluvial alluvial plain deposits individual limestone beds coarsen and contain increasing amounts of siliciclastic material upward fenghuoshan carbonates also contain abundant ostracode bivalve and charaphyte fossils sedimentologic evidence suggests that carbonate deposition occurred in shallow relatively short lived lacustrine systems that were subsequently filled by alluvial deposits of coeval fluvial systems stable c and o isotopic analysis of fenghuoshan carbonates show 18o values ranging from 11 7 to 10 3 vienna pee de

underplating in the himalaya tibet collision zone revealed by the hi climb experiment we studied the formation of the himalayan mountain range and the tibetan plateau by investigating their lithospheric structure using an 800 kilometer long densely spaced seismic array we have constructed an image of the crust and upper mantle beneath the himalayas and the southern tibetan plateau the image reveals in a continuous fashion the main himalayan thrust fault as it extends from a shallow depth under nepal to the mid crust under southern tibet indian crust can be traced to 31 n the crust mantle interface beneath tibet is anisotropic indicating shearing during its formation the dipping mantle fabric suggests that the indian mantle is subducting in a diffuse fashion along several evolving subparallel structures

normal faulting in central tibet since at least 13 5 myr ago tectonic models for the evolution of the tibetan plateau interpret observed east west thinning of the upper crust to be the result of either increased potential energy of elevated crust1 or geodynamic processes that may be unrelated to plateau formation2 3 4 5 6 a key piece of information needed to evaluate these models is the timing of deformation within the plateau the onset of normal faulting has been estimated to have commenced in southern tibet between about 14 myr ago7 and about 8 myr ago8 and in central tibet about 4 myr ago9 here however we report a minimum age of approximately 13 5 myr for the onset of graben formation in central tibet based on mineralization ages determined with rb sr and 40ar 39ar data that post date a major graben bounding normal fault these data along with evidence for prolonged activity of normal faulting in this and other tibetan grabens support models that relate normal faulting to processes occurring beneath the plateau thinning of the upper crust is most plausibly the result of potential energy increases resulting from spatially and temporally heterogeneous changes in thermal structure and density distribution within the crust and upper mantle beneath tibet this is supported by recent geophysical and geological data10 11 12 13 14 15 16 17 which indicate that spatial heterogeneity exists in both the tibetan crust and lithospheric mantle

dislocation creep regimes in quartz aggregates abstract using optical and tem microscopy we have determined that three regimes of dislocation creep occur in experimentally deformed quartz aggregates depending on the relative rates of grain boundary migration dislocation climb and dislocation production within each regime a distinctive microstructure is produced due primarily to the operation of different mechanisms of dynamic recrystallization at lower temperatures and faster strain rates the rate of dislocation production is too great for diffusion controlled dislocation climb to be an effective recovery mechanism in this regime recovery is accommodated by strain induced grain boundary migration recrystallization with an increase in temperature or decrease in strain rate the rate of dislocation climb becomes sufficiently rapid to accommodate recovery in this regime dynamic recrystallization occurs by progressive subgrain rotation with a further increase in temperature or decrease in strain rate dislocation climb remains sufficiently rapid to accommodate recovery however in this regime grain boundary migration is rapid thus recrystallization occurs by both grain boundary migration and progressive subgrain rotation the identification of the three regimes of dislocation creep may have important implications for the determination of flow law parameters and the calibration of recrystallized grain size piezometers in addition the identification of a particular dislocation creep regime could be useful in helping to constrain the conditions at which a given natural deformation has occurred