crustal reworking at nanga parbat pakistan metamorphic consequences of thermal mechanical coupling facilitated by erosion within the syntaxial bends of the india asia collision the himalaya terminate abruptly in a pair of metamorphic massifs nanga parbat in the west and namche barwa in the east are actively deforming antiformal domes which expose quaternary metamorphic rocks and granites the massifs are transected by major himalayan rivers indus and tsangpo and are loci of deep and rapid exhumation on the basis of velocity and attenuation tomography and microseismic magnetotelluric geochronological petrological structural and geomorphic data we have collected at nanga parbat we propose a model in which this intense metamorphic and structural reworking of crustal lithosphere is a consequence of strain focusing caused by significant erosion within deep gorges cut by the indus and tsangpo as these rivers turn sharply toward the foreland and exit their host syntaxes the localization of this phenomenon at the terminations of the himalayan arc owes its origin to both regional and local feedbacks between erosion and tectonics

kinematics of active deformation in the sulaiman lobe and range pakistan the western margin of the indian plate is highly oblique to the direction of convergence between india and asia and represents an excellent example of large scale oblique continent continent collision determining the strain field in western pakistan and how it relates to the plate motion and plate margin geometry affords an exceptional opportunity for understanding oblique margin processes in general through the inversion of regional and teleseismic body waves we have determined the source parameters of 10 moderate sized earthquakes that occurred between 1964 and 1985 in and around the sulaiman range pakistan the earthquakes are dominantly thrust events with slip vectors that are approximately perpendicular to the lobate sulaiman mountain front slip vector orientations rotate 60 70 from a n s to a wnw ese direction of compression consistent with the geometries of the complex festoon shaped mountain belts of this margin we have estimated the spatial variation of the horizontal strain rate and velocity fields within sulaiman using vertically averaged models that accommodate plate motion constraints within a deforming layer the most important factors determining the style of strain rotation in the sulaiman lobe and range are the presence of pure strike slip motion along the chaman fault and the relatively rigid and undeformed katawaz basin that is therefore allowed to translate obliquely relative to india this same conclusion is obtained using either a three dimensional frictional analogue model with significant basal tractions or a thin sheet viscous numerical model without basal tractions thrusting in a predominantly nw se direction in the sulaiman range accommodates 5 14 mm yr of n s motion between india eurasia and 3 6 mm yr of e w shortening seismic moment release this century within the india eurasia plate boundary zone west of the western himalayan syntaxis constitutes roughly 40 of the expected total seismic moment release for this time period particularly significant moment rate deficits exist within the sulaiman range and along the chaman fault

the late neogene 87sr 86sr record of lowland himalayan rivers fossil shells and paleosol carbonate from ancestral himalayan river deposits provide a 87 sr 86 sr record of lowland himalayan river water during the late neogene reconstructed 87 sr 86 sr river values increased sharply in the late miocene probably marking the beginning of exhumation of high 87 sr 86 sr metalimestones more in the central than in the western himalayas these results imply that the marine 87 sr 86 sr record may not be a proxy for silicate weathering or consumption of atmospheric co 2 resulting from that weathering

possible thermal buffering by crustal anatexis in collisional orogens thermobarometric evidence from the nepalese himalaya one dimensional thermal models suggest that large thrust faults may exert first order control on the thermal structure of collisional orogenic belts and an increasing body of petrologic data may be interpreted as consistent with these models we present new thermobarometric data from the central nepalese himalaya which imply that large sections of the himalayan orogen may have been at roughly uniform temperatures during active thrusting a phenomenon that is not predicted by existing models we attribute this phenomenon to thermal buffering as a consequence of widespread in situ anatexis if this process is common in the intermediate and lower crust of active orogens and if it is localized spatially one dimensional thermal models may produce an extremely oversimplified picture of heat transfer processes during collisional orogenesis

river profiles along the himalayan arc as indicators of active tectonics abstract longitudinal profiles along sixteen major transverse himalayan rivers add important constraints to models of active continental subduction and its evolution these profiles are characterized by a zone of relatively high gradient that cannot be associated with differential resistence to erosion in all cases the base of the zone of increased gradients correlates with 1 the topographic front between the lesser and high himalayas 2 the narrow belt of intermediate magnitude thrust earthquakes 3 the main central thrust zone mct these features define a small circle in the central portion of the himalayan arc these correlations suggest that the discontinuity in the river profiles and the other features are controlled by a major tectonic boundary between the rising high himalayas and the lesser himalayas no sharp increases in gradient are observed near the main boundary thrust mbt except on a few rivers such as the jhelum or kundar where the mbt lies close to both the mct and the seismic belt thus it is unlikely that the mbt is a major tectonic boundary the diversion of river courses along the mbt and around anticlines in the sub himalayas has probably been caused by aggradation near the rosion deposition boundary upstream of uplifts in the mahabharat range and sub himalayas a parallel is drawn between the himalayas and new guinea based on the hypothesis that continent arc collision of the type occurring in northern australia preceded continent continent collision in the himalayas the present sedimentary tectonic phase in new guinea resembles the subathu paleocene eocene phase in the himalayas incipient counterparts of the major himalayan structures including the mct and the mbt are recognized in new guinea the drainage patterns in the himalayas and in new guinea bear a similar relation to major structures this suggests that 1 the tectonic evolution of the himalayas has been rather uniform since early stages of collision and 2 the himalayan drainage was also formed at these early stages and is therefore antecedent to the rise of the high himalayas

isostatic rebound due to tectonic denudation a viscous flow model of a layered lithosphere a four layer model of the upper 150 km of the earth is used to calculate the viscous response of continental crust and the underlying mantle to tectonic denudation the model comprises a strong upper crustal layer a weak lower crustal layer a very strong mantle lithosphere layer and a weak mantle asthenosphere layer which is in accord with experimental constraints on strength depth profiles for continental lithosphere the strength of each layer is represented by its effective viscosity flow in the crust and mantle is driven by buoyancy forces which arise from the unloading of an allochthon along a detachment fault by a series of instantaneous displacements earthquakes or rapid creep events numerical solutions obtained by using a finite element technique predict footwall uplift moho deflection and surface topography that are consistent with observations from the basin and range province of the western united states the calculated curvature of the footwall uplift is also similar to that observed and is sensitive to the geometry of the detachment fault such bending need not be elastically controlled hence the curvatures of footwall domes do not clearly place limits on the effective elastic thickness of the extending crust the upward deflection of the moho and the surface topography are sensitive to the viscosity structure and enable us to bound the range of the various viscosities by matching observations from the basin and range province which indicate no moho deflection and low magnitude of surface topography 3 5 km we estimate the upper crustal lower crustal and mantle lithospheric viscosities in the ranges 1021 1023 pa s 1019 1021 pa s and 1021 1023 pa s respectively

geochronologic constraints on syntaxial development in the nanga parbat region pakistan 40ar 39ar data hornblende biotite muscovite and k feldspar and u pb data zircons were obtained from the nanga parbat haramosh massif nphm nw pakistan along three transects in the southern regions of the nphm we have based our interpretations on our new data as well as geochronologic dates from previous studies in the northern regions of the massif geochronologic data show that the nphm has experienced exceptionally high denudation and cooling rates over the past 10 m y u pb ages determined through sensitive high resolution ion microprobe shrimp depth profiling experiments on metamorphic zircons and conventional u pb monazite dates suggest that the timing of metamorphism varied across the massif in addition we have documented that the massif has experienced postmetamorphic differential cooling both along and across strike thermochronologic data on currently exposed surface rocks suggest that cooling occurred more recently and at greater rates in the south central regions of the massif representing deeper crustal levels than along the margins and northern regions of the massif within the tato region cooling following peak metamorphic temperatures of 600 700 c was as high as 140 c m y following partial melting of pelitic units biotites from this area record plateau ages of 0 9 0 1 ma along the astor and indus gorges cooling was less rapid approximately 70 80 c m y following peak metamorphism as indicated by u pb monazite ages of 6 8 ma and 40ar 39ar muscovite cooling ages of 2 2 3 4 ma cooling over the last 3 m y occurred at rates of 100 140 c m y the overall cooling age pattern within the massif is interpreted syntaxial growth through the development of north plunging antiforms prior to 3 ma followed by reverse faulting along east dipping fault zones along the raikot river transect the biotite cooling age pattern is consistent with the folding of isotherms during folding of the foliation surfaces the age pattern was disrupted at 1 ma due to faulting along the raikot and tato faults an electronic supplement of tables a1 a2 and a3 may be obtained on a diskette or anonymous ftp from kosmos agu org login to agu s ftp account using anonymous as the username and guest as the password go to the right directory by typing cd apend type ls to see what files are available type guest and the name of the file to get it finally type exit to leave the system paper 95tc00032 geochronologic constraints on syntaxial development in the nanga parbat region pakistan david m winslow peter k zeitler c page chamberlain and ian s williams diskette may be ordered from american geophysical union 2000 florida avenue n w washington dc 20009 15 00 payment must accompany order

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

geochemistry of a dry steam geothermal zone formed during rapid uplift of nanga parbat northern pakistan abstract natural dry steam zones vapour only are relatively rare most geothermal systems contain both liquid and vapour and typically follow a boiling point depth bpd relationship the nanga parbat uplift driven conductive thermal anomaly results in a geothermal system which follows a bpd relationship at shallow levels but below about 3 km fluid inclusions show that the hydrothermal fluid is dry steam with fluid densities from 0 36 to as low as 0 07 g cm 3 this dry steam zone may persist down to the brittle ductile transition the dry steam has salinities up to 5 wt dissolved salts and up to 22 mole dissolved co 2 the dry steam originated as meteoric water high on the slopes of nanga parbat with 18 o as low as 16 oxygen isotopic exchange with the host rock was facilitated by high temperatures 340 to 450 c and low fluid densities so that the fluid meteoric isotopic signature was completely obliterated hence quartz veins formed by the migrating dry steam have 18 o between 9 and 15 a range which is indistinguishable from quartz in the host rocks quartz vein precipitation from dry steam requires 3 to 5 orders of magnitude greater volume of fluid than typical hydrothermal fluids the dry steam zone at nanga parbat has formed due to near isothermal depressurization of very hot fluid during rapid tectonic uplift at rates 3 6 mm year

earthquake strain rates and instantaneous relative motions within central and eastern asia summary we use the spatial distribution of moment tensors of earthquakes in this century to estimate the velocity field in asia within a eurasian reference frame in a least squares inversion strain rates on the surface of the earth are matched with continuous spline functions in order to recover the velocity gradient tensor associated with the seismic moment release in asia earthquakes account for 40 60 per cent of the expected motion of india relative to eurasia with the missing component of strain rate equivalent to about 20 mm yr 1 of n s shortening between siberia and india in this solution south china rotates counterclockwise and moves eastwards relative to siberia using rigid plate constraints we next investigate the characteristics of the complete horizontal strain field in asia that accommodates plate motions our strain rate solutions are analogous to the response of a newtonian thin viscous sheet in which the rate of work done by the straining medium in accommodating the velocity boundary conditions is a minimum in these solutions the euler pole for india relative to eurasia is constrained nuvel 1a demets et al 1994 but in the process of fitting the vlbi velocity at shanghai china ward 1994 heki et al 1995 the euler pole for south china is determined in the inversion a solution that both fits the velocity at shanghai china and yields a strain rate field consistent with the earthquake mechanisms is one where the south china block has a motion relative to siberia described by the pole at 51 n 131 e 0 3 deg myr 1 comparison of the complete strain field that accommodates plate motion with the seismic strains indicates that earthquake moment release rates in this century within mongolia are about a factor of 4 larger than the long term rate within gansu ningxia the earthquake moment rates have been about a factor of 2 higher than the long term rate the strike slip faulting within mongolia gansu ningxia western sichuan and yunnan is possibly a direct result of velocity boundary conditions imposed on the south china block by forces unrelated to continental collision of india and eurasia such as forces associated with subduction along the margins of south eastern asia verification of this requires a better understanding of the role of pre existing zones of weakness within the asian continental lithosphere

a late miocene pliocene origin for the central himalayan inverted metamorphism abstract perhaps the best known occurrence of an inverted metamorphic sequence is that found immediately beneath the himalayan main central thrust mct generally thought to have been active during the early miocene however in situ 208 pb 232 th dating of monazite inclusions in garnet indicates that peak metamorphic recrystallization of the mct footwall occurred in this portion of the central himalaya at only ca 6 ma the apparent inverted metamorphism appears to have resulted from activation of a broad shear zone beneath the mct which tectonically telescoped the young metamorphic sequence this explanation may resolve some outstanding problems in himalayan tectonics such why the mct and not the more recently initiated thrusts marks the break in slope of the present day mountain range it also renders unnecessary the need for exceptional physical conditions e g high shear stress to explain the generation of the himalayan leucogranites

timing of initial seafloor spreading in the newfoundland iberia rift broad areas of subcontinental lithospheric mantle are commonly exposed along ocean continent transition zones in magma poor rifts and are thought to be exhumed along lithospheric scale detachment faults during the final stages of rifting however the nature of the transition from final rifting to seafloor spreading is controversial we present the first high precision u pb zircon geochronologic and hf isotopic data from gabbros that intrude exhumed mantle at ocean drilling program odp sites 1070 and 1277 in the newfoundland iberia rift north atlantic the sites are conjugate to one another within crust that is commonly considered to have been emplaced during early seafloor spreading magnetic data suggest that crustal accretion occurred at both sites during magnetic polarity chrons m3 m0 130 126 ma however our data indicate that asthenospheric melts were emplaced over brief intervals 1 m y prior to or coeval with mantle exhumation at 124 ma at odp site 1070 and 115 ma at odp site 1277 we suggest that this discrepancy is the result of continued mantle exhumation along large west dipping detachment faults until lithospheric breakup the breakup location is likely coincident with the large amplitude magnetic j anomaly and our 115 ma date for magmatism within this anomaly provides the best available age constraint for breakup along the studied transect

cooling history of the nw himalaya pakistan fission track and 40ar 39ar cooling ages indicate that the late tertiary cooling history of the himalayan ranges of northern pakistan is largely a function of uplift and erosion interpretation of cooling ages which range from under 0 5 ma to over 80 ma suggests that during the late tertiary long term uplift rates at least doubled from under 0 2 mm yr to in some cases well over 0 5 mm yr uplift rates show strong and systematic regional variations as well which reflect the greater uplift of eastern and northern regions the association of very rapid uplift and erosion with the nanga parbat haramosh massif can be explained by a locally vigorous collision of india with eurasia near a promontory of indian crust the resultant rapid uplift of the nanga parbat haramosh massif reactivated the main mantle thrust melange zone with a reversed sense of motion discontinuities in the coolingage distribution along the main mantle thrust in the southern swat hazara region may be the result of the thermal effects of overthrusting

quantitative determination of the effect of temperature on mudstone decay during wet dry cycles a case study of purple mudstone from south western china although temperature is considered a key factor influencing rock decay little is known about the quantitative relationship between a rock s decay rate and temperature which makes it difficult to quantitatively predict the rate of rock decay under varied temperature during wetting drying cycle conditions thus the aim of the present paper was to observe the decay rate of various purple mudstones and to develop a model with which to calculate their decay rates under varied temperature conditions three types of purple mudstone were sampled from the tuodian group j 3 t matoushan group k 2 m and lufeng group j 1 l located in the chuxiong district of yunnan province south western china all samples were manually cut into cubes of 50 x 50 x 50 mm before being subjected to 39 wetting drying cycles in a 2 mm sieve under artificial constant temperature natural daily temperature variation and artificial temperature alternation minimum and maximum treatments in the laboratory decay rates were calculated by weighing the mass remaining in the sieve after each treatment cycle the results showed that the decay rates of the tested rocks rose with an increase in temperature and temperature difference maximum minimum with the rank order of decay rate being j 3 t k 2 m j 1 l quantitative analysis revealed rock decay rate to be significantly related to temperature by a power function and exponentially related to temperature difference our results suggested that temperature difference has a stronger effect on purple mudstone decay than that of temperature magnitude during the wetting drying processes c 2015 elsevier b v all rights reserved

mechanical model for the tectonics of doubly vergent compressional orogens a mechanical model of crustal shortening and deformation driven by the relative convergence of rigid underlying mantle plates explains many features of convergent orogens results based on numerical models and supported by sandbox models show that a coulomb crustal layer subject to basal velocity boundary conditions corresponding to asymmetric detachment and subduction of the underlying mantle passes through three stages of orogenic growth 1 block uplift bounded by step up shear zones 2 development of a low taper wedge over the underthrusting mantle plate and 3 development of a low taperwedge overlying the overthrusting mantle plate and verging in the opposite direction when modified by isostasy basal viscous flow surface erosion and denudation and sedimentation the resultant model orogens exhibit a variety of styles with characteristics in common with small rapidly denuded orogens large orogens with plateaus and extensional characteristics and active subduction margins with doubly vergent accretionary wedges and deformed fore arc basins

three dimensional critical wedges tectonics and topography in oblique collisional orogens in orogens deforming in oblique compression both the long wavelength mean topography of averaged slopes and the higher frequency topography first harmonic of valleys and ridges are related to the mechanics of deformation in three dimensions and cannot be adequately modeled in two dimensions a method of analysis of the principal stresses failure planes and topography is presented for three dimensional critical wedges in a critical coulomb orogen which is at failure throughout the components of the stress tensor are interrelated through the failure condition defined in terms of the principal stresses variation of any of the cartesian components must be accommodated by adjustment of one or more of the other components so that the orogen remains at failure this coupling of cartesian components to maintain the wedge at failure permits exploration of the behavior of a three dimensional critical orogen in which shear stresses arise outside of the plane of convergence internal adjustment of shear stresses to externally imposed stress conditions results in the formation of characteristic orogen shapes for an oblique collision zone such as the southern alps of new zealand the resultant mountain belt is one of an obliquely deforming outboard fold thrust belt with subdued relief rising abruptly into a steep backslope to the main divide which falls steeply to the indentor forming the inboard slope the outboard topographic slope is reduced in relation to that of the two dimensional orogen the rigid indentor has a fixed strike and therefore imposes a failure orientation adjacent to the indentor that differs from that of the ideal three dimensional critical wedge solution this produces several different slip orientations and minor partitioning of strain near the indentor first harmonic topography of valleys and ridges arises due to nonlinear heterogeneous behavior of both erosional and mechanical processes and is characterized by uplift controlled ridges in the outboard and incision controlled valleys in the inboard

petrogenetic and tectonic significance of young leucogranites from the northwestern himalaya pakistan in the himalaya of northwestern pakistan small discordant bodies of leucogranite intrude metamorphosed basement of the indian crust as the first step in assessing the tectonic significance of these granites we have used the shrimp ion microprobe to obtain u pb dates on zircon zircons of igneous appearance are scarce in the leucogranites and with one exception are associated with xenocrystic components which occur as discrete grains and as cores rimmed by young zircon zircons of igneous appearance are high in u 5000 to 40 000 ppm low in th u 0 01 to 0 05 and frequently colored a distinct blue xenocrystic components give a range of proterozoic ages consistent with ages previously reported for indian crust in sharp contrast zircons of igneous appearance yield clusters of concordant tertiary ages two leucogranite dikes from southern localities yield rather old emplacement ages of about 35 ma swat and about 50 ma naran together with previously determined 40ar 39ar and fission track cooling ages these ages document immediate and rapid metamorphism and denudation following eocene collision three leucogranite dikes from northern localities within the nanga parbat haramosh massif are very young having emplacement ages of about 2 3 ma 5 ma and 7 ma the formation and emplacement of these leucogranites during the rapid late tertiary denudation of the nanga parbat region suggests to us that decompression melting may be a viable mechanism for leucogranite genesis

temperatures heat flux and frictional stress near major thrust faults we examine the two dimensional advective and conductive transport of heat in a region of thrust faulting both simple theoretical considerations and numerical experiments show that the steady state temperatures near the fault are reduced by a divisor s 1 b z vsin below what they would be with the same heat sources but in the absence of advection in this expression z is the depth to the fault v is the slip rate is the component of the dip of the fault in the direction of underthrusting is thermal diffusivity and b is a dimensionless factor that is essentially equal to one for most forms of heating initial changes in temperatures near the fault are given by t z t t z 0 12vsin t z 0 z these simple formulae are successful because of the neglible influence of lateral conduction of heat at least for slip rates of a few mm yr or more two time constants govern the transition from the initial change in temperature to steady state t1 u v where u is the distance along the fault in the direction of underthrusting from the surface to the depth in question and t2 z 2 2 where z is the depth to the fault when elapsed times exceed the sum of these two time constants temperatures differ from their steady state values by only about 10 percent the numerical experiments indicate that the simple formulae are sufficiently accurate that sophisticiated numerical modeling of temperatures in specific regions is unwarranted putain an application of these simple formulae to measurements of conductive heat flow at island arcs implies that shear stresses at island arcs approach 100 mpa and are greater than 30 mpa calculations of temperatures appropriate for the himalaya suggest that shear stresses of 100 mpa on the main central thrust probably are required to account for the tertiary granites of the region if melting took place after slip began on the thrust similarly the cut off in seismicity at a depth of about 15 km in the himalaya if due to temperatures exceeding 350 to 450 c implies a deviatoric stress close to 100 mpa

constraints on the tectonic evolution of the northwestern himalaya from geochronologic and petrologic studies of babusar pass pakistan petrologic geochronologic and structural studies of indian plate rocks in the footwall of the main mantle thrust indicate that peak metamorphism occurred in thrust imbricates well before the final assembly of the kohistan island arc with the indian plate as the northern margin of the indian plate was thickened by this early thrusting the metamorphic rocks evolved along increasing pressure temperature paths and in some regions reached temperatures sufficient for partial melting it is difficult to constrain the timing of this early deformation and metamorphism but 40 ar 39 ar hornblende ages and u pb ages of crosscutting leucogranite dikes suggest that in most areas peak metamorphism was earlier than 40 ma and in some areas earlier than 50 ma after peak metamorphism in the interval from 40 to 25 ma rocks of the indian plate continued to be thrust southward creating sharp metamorphic discontinuities at thrust boundaries from 25 to 20 ma rocks of the indian plate underwent rapid cooling

metamorphism and melting of the lithosphere due to rapid denudation nanga parbat massif himalaya the nanga parbat haramosh massif of the western himalaya has undergone a complex metamorphic and denudational history over the past 1 8 ga metamorphism in the tato region is characterized by a barrovian type metamorphism followed by partial melting and recrystallization along shear zones during nearly isothermal decompression while the massif was still at temperatures 550 c mineralogical and textural evidence for this path includes garnet zoning patterns late stage migmatization abundant cordierite both in schists and leucogranite dikes replacement of high pressure assemblages by low pressure assemblages and abundant low density inclusions thermobarometry on the nonmigmatized gneisses reveals conditions of 540 740 c at 7 1 13 1 kbar thermobarometry in the migmat ized rocks reveals final equilibration at 608 675 c at 3 9 6 8 kbar early fluid inclusions occur in quartz inclusions within garnet porphyroblasts that grew during decompression the last fluid inclusions to be trapped occur in microfrac

uplift of tibetan plateau a history of the elevation and crustal thickness of the tibetan plateau since the continental collision at 45 ma has been computed using a kinematic model based on plate tectonic reconstructions and conservation of crustal volume in our reconstructions the indochina block was attached to the southern margin of eurasia before being extruded to the southeast as evidenced by the opening of the south china sea this opening affected the mass input into tibet so that tibet did not uplift to an elevation higher than 1 km until early miocene when spreading of the south china sea ceased the model then gives a steady rise of tibet from middle miocene to its present elevation