retrograde net transfer reaction insurance for pressure temperature estimates retrograde net transfer reactions significantly affect compositions of metamorphic minerals yet are rarely considered when determining pressure temperature p t conditions two natural amphibolite facies metapelites from the central himalaya of nepal exhibit extremely common compositional patterns including increases in mn and fe fe mg at the rims of garnets which are the result of retrograde garnet dissolution and fe mg exchange with biotite however typical thermobarometric approaches for these rocks result in errors of hundreds of degrees and 3 6 kbar compared with thermobarometry of nearby rocks and petrogenetic grids these large errors result because dissolution of high fe garnet has strongly affected the fe mg ratio of matrix biotite x ray maps help evaluate the extent and chemical effects of retrograde reactions in these samples by identifying mineral regions that retain highest t compositions or through a new data processing approach by permitting correction of mineral compositions to original high t values these approaches ensure against retrograde net transfer reactions and should be applied routinely in thermobarometric studies they ultimately yield p t estimates that are more petrologically reasonable and permit rapid screening of samples for those least affected by retrograde reactions reconsideration of thermobarometry in the central and eastern himalaya indicates that retrograde net transfer reactions are extremely common therefore previous thermobarometric studies based on garnet major element compositions from that region should be reevaluated

a model for garnet and plagioclase growth in pelitic schists implications for thermobarometry and p t path determinations numerical models of the progressive evolution of pelitic schists in the ncmnkfmash system with the assemblage garnet biotite chlorite staurolite plagioclase muscovite quartz h2o are presented with the goal of predicting compositional changes in garnet and plagioclase along different p t paths the numerical models support several conclusions that should prove useful for interpreting the p t paths of natural parageneses i garnet may grow along p t vectors ranging from heating with decompression to cooling with compression p t paths deduced from garnet zoning that are inconsistent with these growth vectors are self contradictory ii there is a systematic relation between garnet and plagioclase composition and growth such that for most p t paths garnet growth requires plagioclase consumption furthermore mass balance in a closed system requires that as plagioclase is consumed the remaining plagioclase becomes increasingly albitic inclusions of plagioclase in the core of garnet should be more anorthitic than those near the rim and zoned matrix plagioclase should have rims that are more albitic than the cores complex plagioclase textures may arise from the local variability of growth and precipitation kinetics iii a decrease of fe fe mg in a garnet zoning profile is a reliable indicator of increasing temperature for the assemblage modelled however there is no single reliable p monitor and inferences about p can only be made by considering plagioclase and garnet together iv consumption of garnet during the production of staurolite removes material from the outer shell of a garnet and may make recovery of peak metamorphic compositions and p t conditions impossible low peak temperatures typically recorded by staurolite bearing assemblages may reflect this phenomenon v diffusional homogenization of garnet affects the computed p t path and results in a clockwise rotation of the computed p t vector relative to the true p t path

on the interpretation of peak metamorphic temperatures in light of garnet diffusion during cooling abstract finite difference models of fe mg diffusion in garnet undergoing cooling from metamorphic peak conditions are used to infer the significance of temperatures calculated using garnet biotite fe mg exchange thermometry for rocks cooled from high grades where the garnet was initially homogeneous the calculated temperature tcalc using garnet core and matrix biotite depends on the size of the garnet the ratio of garnet to biotite in the rock vgarnet vbiotite and the cooling rate for garnets with radii of 1 mm and vgarnet vbiotite 1 tcalc is 633 700 and 777 c for cooling rates of 1 10 and 100 c ma for vgarnet vbiotite 1 and 4 and a cooling rate of 10 c ma tcalc is approximately 660 and 610 c respectively smaller and larger garnets have lower and higher tcalc respectively these results suggest that peak metamorphic temperatures may be reliably attained from rocks crystallized at conditions below tcalc of the garnet core provided that vgarnet vbiotite is sufficiently small 0 1 and that the composition of the biotite at the metamorphic peak has not been altered during cooling numerical experiments on amphibolite facies garnets with nominal peak temperatures of 550 600 c generate a well in fe fe mg near the rim during cooling maximum calculated temperatures for the assemblage garnet chlorite biotite muscovite plagioclase quartz using the fe fe mg at the bottom of the well with matrix biotite range from 23 43 c to 5 12 c below the peak metamorphic temperature for cooling rates of 1 and 100 c ma respectively maximum calculated temperatures for the assemblage garnet staurolite biotite muscovite plagioclase quartz are approximately 70 c below the peak metamorphic temperature and are not strongly dependent on cooling rate the results of this study indicate that it may be very difficult to calculate peak metamorphic temperatures using garnet biotite fe mg exchange thermometry on amphibolite facies rocks tmax 550 c because the rim composition of the garnet which is required to calculate the peak temperature is that most easily destroyed by diffusion

the steady state thermal structure of eroding orogenic belts and accretionary prisms an analytical expression derived in this paper for the steady state thermal structure of thrust belts and subduction systems shows that surface erosion rate basal accretion rate and upper plate heat production are extremely important in controlling the temperature and metamorphic history of orogenic belts and accretionary prisms especially where upper plate thickness h is large for example even at rapid erosion and accretion rates and with large upper plate heat productions temperatures near the toe of a thrust wedge e g at h 5 km are nearly the same as those that would be present in the absence of erosion accretion and upper plate heat production in thicker parts of the same wedge e g at h 30 km the combined effects of erosion accretion and upper plate heat production can be dramatic for sufficiently large magnitudes of upper plate heat production or sufficiently rapid erosion rates inverted negative geothermal gradients are produced in the lower part of the upper plate and the upper part of the lower plate although it is difficult to produce inverted geotherms where h 20 km the potential importance of erosion and accretion can be illustrated through two examples 1 beneath the himalayas increasing the upper plate radiogenic heat production from 0 to 2 5 mw m3 and increasing the rate of surface erosion from 0 to 1 mm yr raises the computed steady state temperatures on the main central thrust at h 30 km from 265 c to 600 c and produces a local temperature maximum of 635 c located 6 km above the fault and 2 within the upper plate of the peru trench increasing the rate of surface erosion and basal accretion from 0 to 2 4 mm yr or increasing the upper plate radiogenic heat production from 0 to 1 3 mw m3 increases the computed steady state surface heat flow at h 30 km from 10 to 33 mw m3 but only increases temperatures at the fault surface from 125 c to 170 c shear heating along thrust faults or subduction boundaries is not necessary to produce the high metamorphic temperatures or surface heat flows observed in these two examples and the generation of large amounts of heat by shear stresses on the faults more than a few tens of bars appears to be inconsistent with some of the observations made in these areas this study indicates that many of the major discrepancies between existing thermal models of overthrust belts and observed temperatures as inferred from field data can be attributed to neglect of horizontal and vertical advection of heat caused by surface erosion and by accretion of material from lower plate to upper plate across the fault surface

effects of diffusional modification of garnet growth zoning on p t path calculations the effect of intragranular diffusion on chemical zoning in garnet and on p t paths calculated from that zoning was evaluated using a numerical model of multicomponent diffusion in combination with simulations of garnet growth syn and post growth diffusion of mg fe mn ca species in garnet was calculated for a model pelitic assemblage over a range of temperatures from 485 to 635 c compositions from zoned garnet as modified by diffusion hypothetical inclusions of plagioclase within garnet and matrix phases were used to reconstruct pressure temperature p t paths from isobaric and polybaric model histories p t path calculations based on numerical simulations conducted over an input isobaric heating path that reached peak temperatures between 585 and 635 c show that relaxation of garnet compositional gradients by diffusion can induce modest to appreciable curvature in the inferred paths retrieved paths also indicated somewhat smaller overall temperature changes relative to the actual temperature difference of the input path the magnitude of these distortions is shown to depend upon the heating and cooling rate and garnet crystal size as well as the actual peak temperature condition the effect of diffusion on path trajectories in simulations with thermal histories that also included cooling were comparable to heating only models that reached peak temperatures approximately 15 30 c higher compositions of garnets with radii less than 1 mm that reached actual peak temperatures of 605 c along temperature time histories characteristic of regional metamorphism experienced sufficient diffusional relaxation to introduce errors of hundreds of bars to in excess of one kilobar in path trajectories path distortions were significant at heating cooling rates less than 10 c ma but rapidly diminished for rates faster than this in polybaric simulations diffusion effects were least noticeable when the actual pressure temperature conditions changed in a clockwise sense i e convex to higher p and higher t but apprecciable modification was seen in path models that underwent counterclockwise changes in p and t reequilibration of garnet rim compositions occurred during cooling on all paths and temperature maxima obtained from garnet biotite geothermometry underestimated actual peak conditions by 40 to 70 c calculations suggest that p t path trajectories calculated from garnets of at least 1 mm size and that experienced actual thermal maxima below 585 c are not likely to be distorted by diffusional effects during regional metamorphism however p t path reconstructions based on garnet zonation with smaller grains or higher temperatures may lead to misinterpretation of crystallization history the partitioning record of peak metamorphic temperatures may be destroyed by diffusional reequilibration of garnet rim compositions during cooling seriously complicating the task of quantitatively estimating diffusion effects on path calculations

the age and rate of displacement along the main central thrust in the western bhutan himalaya abstract in western bhutan the main central thrust mct is broadly folded creating multiple exposures of the fault surface over a 70 km across strike distance this unusual map pattern presents a unique opportunity to map the mct and document both the magnitude and age of displacement in situ th pb sims and la icp ms geochronology of metamorphic monazite from the immediate hanging wall of the mct indicates that prograde monazite growth in greater himalayan gh rocks continued until 20 8 1 1 ma whereas crystallization of in situ melts characterized by high y monazite overgrowths occurred during cooling from ca 15 10 ma prograde monazite growth at 15 ma in lesser himalaya lh rocks in the immediate footwall requires that lh footwall strata began to be buried at this time and the mct had reached its southernmost exposed extent by combining prograde monazite ages in the immediate hanging wall and footwall the duration of mct displacement is bracketed between 20 8 1 1 and 15 0 2 4 ma immediately north of our study area a published estimate of shearing along the outer south tibetan detachment std argues for displacement between 20 and 15 ma coeval with the age range for mct displacement that we document in this study however retrograde monazite grains as young as 10 ma suggest that gh rocks were cooling until 10 ma 5 myr later than motion on the outer std immediately to the north this cooling was either the result of continued displacement on the mct or growth of a duplex that passively folded the mct using a sequential reconstruction we estimate a total displacement of 230 km which is the sum of displacements on the mct and the structurally lower paro thrust over a duration of 5 8 2 6 myr this indicates a horizontal shortening rate of 4 0 3 2 1 3 cm yr which exceeds present rates estimated from geodetic measurements across the himalaya and mct displacement rates c 2 cm yr inferred from petrologic and thermal models in central nepal but is indistinguishable from plate convergence rates calculated for eastern bhutan between 23 and 20 ma 3 3 0 7 cm yr our study highlights that displacement on the mct alone achieved plate velocity rates in western bhutan and that the age and rate of mct displacement varied significantly across the himalaya

centrifuge modelling of the influence of crustal fabrics on the development of transfer zones insights into the mechanics of continental rifting architecture abstract centrifuge analogue experiments are used to model the reactivation of pre existing crustal fabrics during extension the models reproduced a weakness zone in the lower crust whose geometry was varied in order to investigate its role in controlling the architecture of rift segments and related transfer zones the typical rift system geometry was characterised by two offset rift segments connected by a major transfer zone in which boundary faults were oblique to the extension vector and displayed a significant transcurrent component of movement the transfer zone was also characterised by cross basin faults with both trend and strike slip component of movement opposite to that displayed by the master faults typically different structural patterns were obtained by changing the offset angle between the rift segments supporting that the structural pattern at transfer zones is strongly influenced by the orientation of pre existing discontinuities with respect to the stretching vector in the models the aspect ratio ratio of length vs width of the transfer zone shows a positive correlation with the offset angle i e the more the inherited fabric is parallel to the extension direction the longer and narrower the transfer zones in case of staircase offset of the rift segments 90 the structural pattern was characterised by two isolated rift depressions linked by a narrow transfer zone in which border faults with alternating polarity overlapped prominent rise of the ductile lower crust was also observed at the transfer zone many of these geometrical features display striking similarities with natural rift systems the results of the current experiments provide useful insights into the mechanics of continental rift architecture supporting that rift propagation width and along axis segmentation may be strongly controlled by the reactivation of pre existing pervasive crustal fabrics

north south extension within the convergent himalayan region recent work by burg et al 1984 indicates the presence of east west striking gently north dipping normal faults in the higher himalayas and southern tibet that formed during postcollisional convergence of india and tibet these faults extend for at least 800 km along strike we interpret these normal faults as probable late miocene extensional features with perhaps several tens of kilometres of down to the north displacement a simple elastic model suggests that these normal faults may have formed during gravitational collapse of the miocene topographic front between india and tibet in this interpretation gravitational collapse occurred by southward motion relative to india and tibet of a wedge of crustal rocks bounded above by gently north dipping normal faults and below by thrust faults that probably dip north north south extension produced in this way is probably confined to upper crustal levels only and does not reflect regional extension of the entire lithosphere

new insight into the south tibetan detachment system not a single progressive deformation 1 low angle normal faults lanf typically regarded as accommodating crustal or lithospheric extension may also form during lithospheric shortening the best studied system of syn contractional lanfs is the south tibetan detachment system a network of low angle normal sense faults and shear zones that formed coevally with and parallel to south vergent thrusts during lithospheric shortening accompanying development of the himalayan orogen in the eastern himalaya there are several across strike exposures of the south tibetan detachment system we present new structural and thermometry data from the eastern himalaya that demonstrate that the south tibetan detachment system cannot have formed as a single progressive structure we characterize and distinguish two distinct structural and tectonic components within the currently recognized system 1 an extensive diffuse sheared layer that formed the boundary between strong upper crust and weak southward flowing middle crust and 2 a network of brittle ductile lanfs that locally exhume partly excise and overprint the earlier mylonite zone at the topographic break between the himalayan orogen and the tibetan plateau the sheared layer not a lanf formed the boundary between upper and middle crust during ductile flow of the middle crust and is extensively exposed in the himalaya at the base of klippen of upper crustal rocks preserved in bhutan along the crest of the himalaya where it has been excised and exhumed by the brittle ductile extrusion lanfs and bounding the cores of the north himalayan gneiss domes citation kellett d a and d grujic 2012 new insight into the south tibetan detachment system not a single progressive deformation tectonics 31 tc2007 doi 10 1029 2011tc002957

crustal channel flows 1 numerical models with applications to the tectonics of the himalayan tibetan orogen 1 plane strain thermal mechanical numerical models are used to examine the development of midcrustal channel flows in large hot orogens in the models radioactive self heating reduces the viscosity of tectonically thickened crust and increases its susceptibility to large scale horizontal flow channels can be exhumed and exposed by denudation focused on the high relief transition between plateau and foreland we interpret the himalaya to have evolved in this manner channel flows are poorly developed if the channel has a ductile rheology based on wet quartz flow laws and well developed if there is an additional reduction in viscosity to 1019 pa s this reduction occurs from 700 c to 750 c in the models and is attributed to a small percentage of in situ partial melt melt weakening model ht1 provides an internally consistent explanation for the tectonic evolution of many features of the himalayan tibetan orogenic system erosional exhumation exposes the migmatitic channel equivalent to the greater himalayan sequence ghs between coeval normal and thrust sense ductile shear zones corresponding to the south tibetan detachment and the main central thrust systems outward flow of unstable upper crust rotates these shears to low dip angles in the model both the ghs and the lesser himalayan sequence are derived from indian crust with the latter from much farther south similar models exhibit a range of tectonic styles including the formation of domes resembling north himalayan gneiss domes model results are relatively insensitive to channel heterogeneities and to variations in the behavior of the mantle lithosphere beneath the model plateau

ductile extrusion of the higher himalayan crystalline in bhutan evidence from quartz microfabrics quartz textures measured from deformed quartz tectonites within the lesser himalaya and higher himalaya crystalline of bhutan show similar patterns orientation and distribution of the quartz crystallographic axes were used to confirm the regional shear sense the asymmetry of c axis and a axis patterns consistently indicates top to the south shearing the obliquity of the texture and the inferred finite strain plane strain to moderately constrictional suggest the strain regime had a combination of rotational and irrotational strain path in most of the samples from the bhutan himalaya the inferred deformation mechanisms suggest moderate to high temperature conditions of deformation that produced the observed crystallographic preferred orientation much higher temperature of deformation is indicated in the quartz veins from a leucogranite the observed ductile deformation is pervasively developed in the rocks throughout the investigated area the intensity of deformation increases only slightly in the vicinity of the main central thrust simultaneous southward shearing within a large part of the higher himalaya crystalline near and above the main central thrust and normal faulting across the south tibetan detachment is explained by the tectonically induced extrusion of a ductily deforming wedge the process of extrusive flow suggested here can be approximated quantitatively by channel flow models that have been used to describe subduction zone processes channel flow accounts for some observed phenomena in the himalayan orogen such as inverted metamorphic sequences near the main central thrust not related to an inversion of isotherms and the syntectonic emplacement of leucogranites into the extruding wedge locally leading to an inversion of isotherms due to heat advection

exhumation of the main central thrust from lower crustal depths eastern bhutan himalaya geothermometry and mineral assemblages show an increase of temperature structurally upwards across the main central thrust mct however peak metamorphic pressures are similar across the boundary and correspond to depths of 35 45 km garnet bearing samples from the uppermost lesser himalayan sequence lhs yield metamorphic conditions of 650 675 c and 9 13 kbar staurolite kyanite schists about 30 m above the mct yield p t conditions near 650 c 8 10 kbar kyanite bearing migmatites from the greater himalayan sequence ghs yield pressures of 10 14 kbar at 750 800 c top to the south shearing is synchronous with and postdates peak metamorphic mineral growth metamorphic monazite from a deformed and metamorphosed proterozoic gneiss within the upper lhs yield u pb ages of 20 18 ma staurolite kyanite schists within the ghs a few metres above the mct yield monazite ages of c 22 1 ma we interpret these ages to reflect that prograde meta morphism and deformation within the main central thrust zone mctz was underway by c 23 ma u pb crystallization ages of monazite and xenotime in a deformed kyanite bearing leucogranite and kyanite garnet migmatites about 2 km above the mct suggest crystallization of partial melts at 18 16 ma higher in the hanging wall south verging shear bands filled with leucogranite and pegmatite yield u pb crystallization ages for monazite and xenotime of 14 15 ma and a 1 2 km thick leuco granite sill is 13 4 0 2 ma thus metamorphism plutonism and deformation within the ghs con tinued until at least 13 ma p t conditions at this time are estimated to be 500 600 c and near 5 kbar from these data we infer that the exhumation of the mct zone from 35 to 45 km to around 18 km occurred from 18 to 16 to c 13 ma yielding an average exhumation rate of 3 9 mm year 1 this process of exhumation may reflect the ductile extrusion by channel flow of the mctz from between the overlying tibetan plateau and the underthrusting indian plate coupled with rapid erosion

a short duration pulse of ductile normal shear on the outer south tibetan detachment in bhutan alternating channel flow and critical taper mechanics of the eastern himalaya in easternmost bhutan the south tibetan detachment std is a ductile shear zone that juxtaposes the radi or sakteng klippe of the tethyan sedimentary series from underlying high grade greater himalayan rocks in situ la icpms u th pb analysis of metamorphic monazite from the immediate footwall and hanging wall of the std within the shear zone at the base of the klippe constrains north vergent normal shear to between 25 and 20 ma coeval thrusting on the main central thrust during this time supports a phase of channel flow viscous wedge model activity lasting only ca 3 ma geochronologic data from the eastern himalaya indicate alternating mechanisms for extrusion of the metamorphic core of the orogen from the late oligocene through to the late miocene switching from channel flow viscous wedge behavior to critical taper frictional wedge behavior each phase lasting approximately only 2 to 5 ma the tectonic evolution of the eastern himalaya is comparable to central and western himalayan tectonics during the early miocene but during the middle miocene metamorphism and magmatism in the eastern himalaya migrated towards the orogenic hinterland a process not widely documented elsewhere in the himalaya thus highlighting the need for an orogenic model in three spatial dimensions

preliminary stratigraphic and structural architecture of bhutan implications for the along strike architecture of the himalayan system article i nfo preliminary mapping and stratigraphic correlation of lesser himalayan rock in eastern bhutan using field characteristics u pb detrital zircon dating and nd geochemistry define the first order stratigraphic architecture of the indian passive margin sequence in the eastern himalaya we use this new image of the lateral and vertical relationships of the original stratigraphy to determine the structural framework of the eastern himalayan fold thrust belt in bhutan we propose that lesser himalayan rock in bhutan can be divided into lower lesser himalayan rocks with a paleoproterozoic detrital zircon signal and upper lesser himalayan rocks with detrital zircon signals of 1000 500 ma the 500 ma detrital grains are from rocks in the frontal portions of the fold thrust belt north of the main boundary thrust as well as directly in the footwall of the main central thrust in the hinterland our preliminary stratigraphic study coupled with mapping allows us to construct a composite balanced cross section through the kuru chu valley in eastern bhutan which provides the first image of the geometry and amount of shortening through bhutan the main frontal thrust tilts a 6 km of thick section of neogene foreland basin deposits these deposits are separated from folded and faulted upper lesser himalayan rocks by the main boundary thrust north of the main boundary thrust we propose two duplex systems the southernmost duplex contains 9 repeated sections of upper lesser himalayan units including the permian gondwana sequence and the cambrian or younger baxa group the northernmost duplex is located in the footwall of the main central thrust and is comprised of two repeated sections of the proterozoic shumar and daling formations the southern boundary of this northern duplex is the shumar thrust which acts as a roof thrust for the southern duplex system and may correlate to the ramgarh thrust in nepal and india we propose that the development of the duplex systems passively folded overlying tethyan and greater himalayan rocks minimum shortening for this part of the himalayan fold thrust belt is 359 km all of which has occurred from 22 ma to present defining a long term shortening rate of 16 mm yr

flattening the bhutan himalaya a detailed thermobarometric transect of 35 samples across the greater and tethyan himalayan sequences in central bhutan demonstrates a tectonostratigraphically intact section with uniform apparent thermal and baric field gradients of 2072 1c km and 0 5770 08 kbar km pressure temperature paths determined from chemically zoned garnets in 6 samples demonstrate that these p t conditions correspond with maximum pressures the super lithostatic baric gradient cannot be explained by pre to syn metamorphic tectonic processes or by extension within an inclined slab instead the data imply 50 post peakmetamorphic flattening of the himalayan metamorphic core accommodated by distributed top to thenorth shear consistent with microstructural analysis orogenic flattening best explains the development of the south tibetan detachment system as a strain incompatibility feature rather than a structure bounding the top of a tectonically inserted wedge and helps reconcile debate attributing first order himalayan structural features to either wedge failure critical taper or pipe flow channel flow

preserved zr temperatures and u pb ages in high grade metamorphic titanite evidence for a static hot channel in the himalayan orogen abstract titanite grains from 800 c gneisses of the c 6 km thick greater himalayan sequence of central nepal were analyzed for zr in titanite temperatures and u pb ages to investigate the formation and evolution of a former weak mid crustal channel in the himalaya laser ablation icp ms depth profiles and spot analyses of titanite zr and u pb isotopic compositions are inconsistent with diffusional reequilibration on scales greater than 1 2 m and instead data appear to faithfully record both temperature and age titanite records protracted heating from 700 750 c at 37 ma to 775 c at 24 ma and possibly slight cooling to 765 c at 20 ma such temperatures exceed partial melting reactions in associated metapelites and imply profoundly weak rheologies i e a hot channel yet predate initiation of the bounding main central thrust and south tibetan detachment systems by as much as 15 myr evidently thick but weak crustal channels may remain stationary with respect to the convergence direction for 10 myr even during one of earth s biggest continent continent collisions

thermal models of dyke intrusion during development of continent ocean transition a consensus has emerged in recent years from a variety of geoscientific disciplines that extension during continental rifting is achieved only partly by plate stretching dyke intrusion also plays an important role magma intrusion can accommodate extension at lower yield stresses than are required to extend thick strong unmodified continental lithosphere mechanically thereby aiding the breakup process dyke intrusion is also expected to heat and thereby weaken the plate but the spatial extent of heating and the effect of different rates of magmatic extension on the timescales over which heating occurs are poorly understood to address this issue a numerical solution to the heat flow equation is developed here to quantify the thermal effects of dyke intrusion on the continental crust during rifting the thermal models are benchmarked against a priori constraints on crustal structure and dyke intrusion episodes in ethiopia finite difference models demonstrate that magmatic extension rate exerts a first order control on the crustal thermal structure once dyke intrusion supersedes faulting and stretching as the principal extensional mechanism the crust will heat and weaken rapidly less than 1 ma in the main ethiopian rift mer the majority of present day extension is focused on 20 km wide quaternary recent axial magmatic segments that are mostly seismogenic to mid crustal depths and show p wave seismic velocities characteristic of heavily intruded continental crust when reviewed in light of our models these observations require that no more than half of the mer s extension since 2 ma has been achieved by dyke intrusion magmatic heating and weakening of the crust would have rendered it aseismic if dyke intrusion accounted for the entire 6 mm yr extension rate in the older faster extending 16 mm yr red sea rift rsr in afar dyke intrusion is expected to have had a more dramatic impact on crustal rheology accordingly effective elastic plate thickness and moho depth in the danakil region of northernmost afar are markedly reduced and seismicity is shallower than in the mer thermally driven variations in crustal rheology over time in response to dyke intrusion thus play an important role in the development of continent ocean transition

crustal channel flows 2 numerical models with implications for metamorphism in the himalayan tibetan orogen 1 results from a thermal mechanical model ht1 that includes midcrustal channel flow are compatible with many features of the himalayan tibetan system radioactive self heating and rheological weakening of thickened model orogenic crust lead to the formation of a hot low viscosity midcrustal channel and a broad plateau channel material corresponding to the greater himalayan sequence ghs flows outward from beneath the plateau in response to topographically induced differential pressure at the plateau flank it is exhumed by focused surface denudation and juxtaposed with cooler newly accreted material corresponding to the lesser himalayan sequence lhs the model channel is bounded by coeval thrust and normal sense ductile shear zones interpreted to represent the main central thrust mct zone and south tibetan detachment system respectively inverted metamorphism associated with the model mct zone results from distributed ductile shear along the mct and extrusion of the hot channel a variety of model p t t path styles resembling those observed in the ghs and lhs are produced for points traveling through contrasting tectonic regimes that coexist in different parts of the model predicted times of peak metamorphism cooling and erosion of metamorphic facies are generally compatible with observations although model ghs cooling ages are too young the times of m1 and m2 metamorphic events observed in the ghs correspond to model times of maximum burial and maximum heating respectively the results highlight the need to integrate tectonics and metamorphism in continental collision models and demonstrate the importance of lateral transport of both heat and material in large hot orogens

quantifying internal strain and deformation temperature in the eastern himalaya bhutan implications for the evolution of strain in thrust sheets abstract quartz microstructural analysis together with quantification of 3 d crystal plastic strain of quartz clasts normalized fry and rf methods constrain deformation temperature ranges and internal strain magnitude orientation within himalayan thrust sheets in bhutan lesser himalayan lh thrust sheets display an inverted deformation temperature gradient from 250 to 310 c at the main boundary thrust to 500 670 c at the main central thrust mct attributed primarily to stacking of lh thrust sheets deformed at progressively higher temperatures toward the hinterland the hot iron effect of the mct hanging wall only affects the upper few hundred meters of hinterland lh thrust sheets frontal thrust sheets exhibit 7 layer parallel shortening lps strain all other thrust sheets exhibit layer normal flattening lnf strain with 1 9 1 8 1 0 lh rocks and 2 1 1 8 1 0 greater and tethyan himalayan rocks median ellipsoids x parallel to lineation we propose that lps strain developed foreland ward of the thrust deformation front and that lnf strain resulted from later tectonic loading the lps to lnf transition occurs at minimum temperatures for quartz plasticity ca 250 270 c lnf ellipsoids are folded along with strata within thrust sheets indicating that internal strain preceded thrust imbrication and translation at the scale of 10 s of meters structural distance strain magnitude of hanging wall rocks does not increase near thrusts which supports the existence of discrete faults with large translations as generally depicted in balanced cross sections strain magnitudes in low grade greater and tethyan himalayan rocks in central bhutan indicate that the top to the north component of channel flow extrusion is most likely less than 19 km

himalayan metamorphic sequence as an orogenic channel insight from bhutan the bhutan himalayas differ from the rest of the himalayas in two major ways i low grade sedimentary rocks lie above the greater himalayan sequence as klippen i e erosional remains of the south tibetan detachment and ii an out of sequence thrust the kakhtang thrust lies structurally above the klippen and it doubles the exposed thickness of the greater himalayan sequence our field observations and geochronological data constrain the main kinematic events in the bhutan himalayas crystallisation ages of leucogranite dykes deformed by the main central thrust and the south tibetan detachment indicate that these two structures operated together between 16 and 22 ma the out of sequence kakhtang thrust was active at 10 14 ma and was concurrent with reactivation of the south tibetan detachment restoration of the bhutan himalayas prior to the out of sequence thrusting shows that the greater himalayan sequence was the core of a long low viscosity crustal channel extending under the tibetan plateau we propose that the gravity driven southward extrusion of the channel material from underneath the tibetan plateau caused the inverted metamorphic sequence in the lesser himalayan sequence and in the greater himalayan sequence this process also led to occurrences of present day surface rocks that were derived from variable distances down dip but from similar crustal depths such an exhumation pattern can explain the similar peak pressures for the greater himalayan sequence along the length of the himalayas

variable shortening rates in the eastern himalayan thrust belt bhutan insights from multiple thermochronologic and geochronologic data sets tied to kinematic reconstructions 1 we present data on the burial displacement and exhumation history of the himalayan fold thrust belt in eastern bhutan these data document the magnitude and timing of displacement of large discrete structures and highlight temporal variability in shortening rates eight new40ar 39ar ages from white mica 32 new zircon u th he ages 7 new apatite fission track ages and 1 new u pb zircon la mc icp ms metamorphic rim growth age are combined with published cooling ages and deformation temperatures and incremental shortening magnitudes from restorations of two published balanced cross sections to illustrate the kinematic and temporal development of the bhutan thrust belt integrating these data from 23 ma to the present illustrates rapid horizontal shortening rates 28 35 mm yr between 23 20 ma and 15 10 ma separated by more moderate rates 10 23 mm yr shortening rates decrease significantly to 7 10 mm yr and possibly as low as 3 4 mm yr from 10 to 0 ma this decrease is interpreted to represent the onset of strain partitioning in the eastern part of the himalayan tibetan orogenic system between shortening in the bhutan thrust belt uplift of the shillong plateau and deformation and outward growth of the northern and eastern tibetan plateau within estimated error horizontal shortening rates during emplacement of the main central thrust sheet and during construction of the upper lesser himalayan duplex approached india asia tectonic velocities thus for periods of time between 23 20 ma and 15 10 ma the bhutan thrust belt may have absorbed nearly all india asian convergence at this longitude

strain deformation temperatures and vorticity of flow at the top of the greater himalayan slab everest massif tibet this paper presents quantitative data on strain deformation temperatures and vorticity of flow at the top of the greater himalayan slab the data were collected from the tibetan side of the everest massif where two low angle normal faults bound the upper surface of the greater himalayan slab the earlier and structurally lower lhotse detachment and the later and structurally higher qomolangma detachment greenschist to sillimanite grade quartz rich metasedimentary rocks exposed in the rongbuk to north col region of the everest massif are characterized by cross girdle quartz c axis fabrics indicating approximate plane strain conditions fabric opening angles progressively increase with depth beneath the overlying lhotse detachment and indicate progressively rising deformation temperatures of 525 625 50 c at depths of 300 600 m beneath the detachment deformation temperatures of c 450 c are indicated by fabric opening angles in epidote amphibolite facies mylonites located closer to the overlying detachment a top down to the north normal shear sense is indicated by the asymmetry of microstructures and c axis fabrics but the degree of asymmetry is low at distances greater than 400 m beneath the detachment and sillimanite grains are drawn into adjacent conjugate shear bands but still appear pristine indicating that deformation occurred at close to peak metamorphic temperatures these quenched fabrics and microstructures indicate rapid exhumation in agreement with previous isotopic dating studies mean kinematic vorticity numbers w m were independently calculated by three analytical methods calculated w m values range between 0 67 and 0 98 and indicate that although a simple shear component is generally dominant particularly in greenschist facies mylonites located between the lhotse and overlying qomolangma detachments there is also a major component of pure shear in samples located at 400 600 m beneath the lhotse detachment pure and simple shear make equal contributions at w k 0 71 our integrated strain and vorticity data indicate a shortening of 10 30 perpendicular to the upper surface of the greater himalayan slab and confirm that the upper surface of the slab is a stretching fault with estimated down dip stretches of 10 40 assuming plane strain deformation measured parallel to the flow plane transport direction

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