tectonic control on 10 be derived erosion rates in the garhwal himalaya india erosion in the himalaya is responsible for one of the greatest mass redistributions on earth and has fueled models of feedback loops between climate and tectonics although the general trends of erosion across the himalaya are reasonably well known the relative importance of factors controlling erosion is less well constrained here we present 25 10 be derived catchment averaged erosion rates from the yamuna catchment in the garhwal himalaya northern india tributary erosion rates range between 0 1 and 0 5 mm yr 1 in the lesser himalaya and 1 and 2 mm yr 1 in the high himalaya despite uniform hillslope angles the erosion rate data correlate with catchment averaged values of 5 km radius relief channel steepness indices and specific stream power but to varying degrees of nonlinearity similar nonlinear relationships and coefficients of determination suggest that topographic steepness is the major control on the spatial variability of erosion and that twofold to threefold differences in annual runoff are of minor importance in this area instead the spatial distribution of erosion in the study area is consistent with a tectonic model in which the rock uplift pattern is largely controlled by the shortening rate and the geometry of the main himalayan thrust fault mht our data support a shallow dip of the mht underneath the lesser himalaya followed by a midcrustal ramp underneath the high himalaya as indicated by geophysical data finally analysis of sample results from larger main stem rivers indicates significant variability of 10 be derived erosion rates possibly related to nonproportional sediment supply from different tributaries and incomplete mixing in main stem channels

increased late pleistocene erosion rates during fluvial aggradation in the garhwal himalaya northern india the response of surface processes to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution in the himalaya most large rivers feature prominent fill terraces that record an imbalance between sediment supply and transport capacity presumably due to past fluctuations in monsoon precipitation and or effects of glaciation at high elevation here we present volume estimates chronological constraints and 10 be derived paleo erosion rates from a prominent valley fill in the yamuna catchment garhwal himalaya to elucidate the coupled response of rivers and hillslopes to pleistocene climate change although precise age control is complicated due to methodological problems the new data support formation of the valley fill during the late pleistocene and its incision during the holocene we interpret this timing to indicate that changes in discharge and river transport capacity were major controls compared to the present day late pleistocene hillslope erosion rates were higher by a factor of 2 4 but appear to have decreased during valley aggradation the higher late pleistocene erosion rates are largely unrelated to glacial erosion and could be explained by enhanced sediment production on steep hillslopes due to increased periglacial activity that declined as temperatures increased alternatively erosion rates that decrease during valley aggradation are also consistent with reduced landsliding from threshold hillslopes as a result of rising base levels in that case the similarity of paleo erosion rates near the end of the aggradation period with modern erosion rates might imply that channels and hillslopes are not yet fully coupled everywhere and that present day hillslope erosion rates may underrepresent long term incision rates

modeling mountain building and the seismic cycle in the himalaya of nepal a host of information is now available regarding the geological and thermal structure as well as deformation rate across the himalaya of central nepal these data are reconciled in a two dimensional mechanical model that incorporates the rheological layering of the crust which depends on the local temperature and surface processes over geological timescale 5 ma the 20 mm yr estimated shortening rate across the range is accommodated by localized thrust faulting along the main himalayan thrust fault mht the mht reaches the surface along the foothills where it is called the main frontal thrust fault mft the mht flattens beneath the lesser himalaya and forms a midcrustal ramp at the front of the higher himalaya consistent with the river incision and the anticlinal structure of the lesser himalaya farther northward the mht roots into a subhorizontal shear zone that coincides with a midcrustal seismic reflector aseismic slip along this shear zone is accommodated in the interseismic period by elastic straining of the upper crust increasing the coulomb stress beneath the front of the higher himalaya where most of the microseismic activity clusters negligible deformation of the hanging wall requires a low apparent friction coefficient less than 0 3 on the flat portion of the mht on the ramp might be as high as 0 6 sensitivity tests show that a rather compliant quartz rich rheology and a high radioactive heat production in the upper crust of 2 5 w m 3 is required erosion affects the thermal structure and interplays with crustal deformation a dynamic equilibrium is obtained in which erosion balances tectonic uplift maintaining steady state thermal structure topography and deformation field using a linear diffusion model of erosion we constrain the value of the mass diffusivity coefficient to 0 5 1 6 l0 4 m 2 yr this study demonstrates that the data are internally consistent and compatible with current understanding of the mechanics of crustal deformation and highlight the role of viscous flow in the lower crust and of surface erosion in orogeny processes on the long term as well as during interseismic period

isotopic constraints on the structural relationships between the lesser himalayan series and the high himalayan crystalline series garhwal himalaya nd and sr isotope systematics may provide important constraints on the location of major thrust systems that separate lithologically similar sedimentary sequences the potential of the technique is illustrated by this isotopic study of the main central thrust system of the himalaya nd isotope data from the garhwal himalaya indicate that metasedimentary rocks from the vaikrita group nd 14 to 19 correlate closely with those from the high himalayan crystalline series which constitutes the hanging wall lithologies of the main central thrust in contrast metasedimentary rocks from the munsiari group nd 23 to 28 show marked similarities to the lesser himalayan series in the footwall of the main central thrust sr isotopes support the correlations in that the vaikrita group shows partial reequilibration at 500 ma whereas the munsiari group has not undergone sr isotope homogenization since 1800 ma thus the vaikrita thrust that juxtaposes these two formations is recognized as the main central thrust in garhwal himalaya the thrust coincides approximately with the location of the kyanite isograd confirming that inverted metamorphism is characteristic of both hanging wall and footwall of the main central thrust along the tons thrust known locally as the srinagar thrust 50 km south of the main central thrust low grade quartzarenites with nd sr isotope and trace element characteristics typical of lesser himalayan formations have been emplaced on phyllites and siltstones with geochemical characteristics of the high himalayan crystalline series the field relationships most probably result from out of sequence thrusting in which lesser himalayan series rocks to the north were emplaced over low grade equivalents of the high himalayan crystalline series preserved in the external part of the orogen this study establishes the value of isotope data for lithostratigraphic correlations within orogenic belts

topography relief and trmm derived rainfall variations along the himalaya 1 along the southern himalayan topographic front the indian summer monsoon modulates erosive processes and rates to investigate the influence of topography and relief on rainfall generation and resultant erosion we processed satellite rainfall amounts for the last 8 years 1998 2005 from the tropical rainfall measurement mission trmm based upon a spatial resolution of 5 5 km for the himalaya we identify 1 the spatial distribution of rainfall and 2 the large scale relationships between topography relief and rainfall locations our results show two distinct rainfall maxima along strike in the himalaya the first outer rainfall peak occurs along the southern margin of the lesser himalaya within a narrow band of mean elevation 0 9 0 4 km and mean relief 1 2 0 2 km the second discontinuous inner band typically occurs along the southern flank of the greater himalaya elevation and relief both 2 1 0 3 km

on the inference of denudation rates from cooling ages of minerals abstract the inference of denudation rates from the temperature time histories of rocks is not straightforward because the geothermal gradient evolves with time as denudation advects heat towards the land surface we show that the curvature in the temperature time path recorded by an exhumed rock is strongly dependent upon denudation rate and upon the temperature of the rock before denudation began but only weakly dependent upon the initial geothermal gradient the original depth of burial of the rock or the duration of the denudation for the case in which denudation rate is constant a systematic search of parameter space yields robust estimates of denudation rate and of the initial temperature if there is prior knowledge of the initial geothermal gradient the original depth of burial of the rock or the duration of the denudation this may be used further to limit the acceptable ranges of the other two parameters application of this procedure to two geological examples suggests that denudation rates in the nanga parbat massif are lower and less variable than previously suggested and that the inference drawn from accelerated cooling rates of accelerating denudation in the andes is incorrect

cosmic ray labeling of erosion surfaces in situ nuclide production rates and erosion models a number of in situ cosmogenic radionuclides and stable nuclides have been measured in natural exposed rock surfaces with a view to study their in situ production and rock erosion rates 1 the in situ radionuclides can be used for a high resolution tomography of the erosional history of an exposed surface two stable nuclides 3he 21ne and five radionuclides 10be 26al 36cl 14c 39ar having half lives in the range of 300 1 5 106 yr half life are measurable in many rock types a prerequisite for the application of the in situ nuclides for the study of erosional histories of surfaces is a knowledge of their production rates under different irradiation conditions altitude latitude irradiation geometry and shielding relative nuclide production rates can be determined fairly accurately using the extensive available data on cosmic ray neutrons 2 absolute nuclide production rates cannot generally be predicted with any accuracy because of lack of data on excitation functions of nuclides unless some normalization is possible as was done in the case of several cosmic ray produced isotopes in the atmosphere 3 based on a recent natural calibration experiment in which erosion free surfaces exposed to cosmic radiation for 11 000 yrs were sampled the absolute production rates of 10be and 26al in quartz have been accurately estimated for mountain altitudes in sierra nevada 4 the absolute production rates of 10be and 26al in quartz can therefore be estimated fairly accurately for any given latitude and altitude some measurements of 14c in rocks of low erosion rate 5 similarly allow an estimate of its production rate attempts made to measure the in situ production rates of 3he in rocks have not yet led to a convergent production rate in view of the importance of knowing the production rates of isotopes of he ne and ar i present here theoretical estimates of their production rates based on available cross section data i discuss the information that can be extracted from the study of the in situ nuclides in rocks useful parameters characterizing the exposure history of a rock surface are 1 the effective surface exposure age and 2 the time averaged erosion rate the implications of these parameters for single and multiple nuclide studies are discussed in terms of the erosion models considered

thermotectonic history of the chiplakot crystalline belt in the lesser himalaya kumaon india constraints from apatite fission track thermochronology abstract fission track ages and confined track length distribution of apatite samples separated from the chiplakot crystalline belt ccb of the lesser himalayan crystalline lhc zone located to the south of the main central thrust mct munsiari thrust mt in kumaon india have been determined ages from the ccb along the kali and darma valleys fall in two distinct groups in the northern part of the ccb the ages range from 9 8 0 6 to 7 6 0 6 ma with a weighted mean of 9 6 0 1 ma while in the southern part the ages vary from 17 9 0 9 to 12 9 1 1 ma with a weighted mean of 14 1 0 1 ma the bimodal distribution of track lengths indicates that the ages are mixed ages rather than simple cooling ages the apatite fission track aft ages and already published structural data of the ccb suggest a complex erosional denudation history within the upper 3 4 km of the crust of the ccb the ages further indicate that the ccb was thrust into place earlier than the middle miocene i e at the time of development of the mct since then these rocks have remained within the upper 3 km of crust and were affected by only moderate to slow erosion and exhumation these results have important implications for the tectonic evolution of the lhc zone to the south of the mct mt the exhumation of the lhc zone in different parts of the himalaya was not uniform in the kumaon himalaya it was not controlled as in the himachal himalaya by any major tectonic event since it was thrust over the lesser himalayan meta sedimentary lhms zone and underwent moderate to slow erosion and exhumation

uplift and convergence along the himalayan frontal thrust of india along the himalayan thrust front in northwestern india terrace deposits exposed 20 to 30 m above modern stream level are interpreted to have been uplifted by displacement on the underlying himalayan frontal thrust a radiocarbon age limits the age of the terrace to 1665 215 calendar bc 3663 215 radiocarbon years before present yielding a vertical uplift rate of 6 9 1 8 mm yr in combination with published studies constraining the dip of the himalayan frontal thrust fault to about 30 in the study area the observed uplift rate equates to horizontal shortening across the himalayan frontal thrust of 11 9 3 1 mm yr and the slip rate of the himalayan frontal thrust of 13 8 3 6 mm yr this is similar to previously reported rate estimates along the himalayan arc based on displacement of older plio miocene age rocks or the much shorter records of geodesy and historical seismicity the similarity is consistent with the idea that convergence across the himalayan front has occurred at a relatively steady rate through time the seismic expression of this deformation includes several great m 8 historical earthquakes which due to lack of surface rupture during those events have been attributed to their occurrence on blind thrusts yet the occurrence of a possible fault scarp in the field area indicates that past earthquakes have been sufficiently large to rupture to the surface and produce coseismic scarps these observations suggest a potential for earthquakes along the himalayan frontal thrust larger than those observed historically

effects of bedrock landslides on cosmogenically determined erosion rates abstract the successful quantification of long term erosion rates underpins our understanding of landscape formation the topographic evolution of mountain ranges and the mass balance within active orogens the measurement of in situ produced cosmogenic radionuclides crns in fluvial and alluvial sediments is perhaps the method with the greatest ability to provide such long term erosion rates although deep seated bedrock landsliding is an important erosional process in active orogens its effect on crn derived erosion rates is largely unquantified we present a numerical simulation of cosmogenic nuclide production and distribution in landslide dominated catchments to address the effect of bedrock landsliding on erosion rates results of the simulation indicate that the temporal stability of erosion rates determined from crn concentrations in sediment decreases with increased ratios of landsliding to bedrock weathering rates within a given catchment area and that as the frequency of landsliding increases larger catchment areas must be sampled in order to accurately evaluate long term erosion rates in addition results of this simulation suggest that sediment sampling for crns is the appropriate method for determining long term erosion rates in regions dominated by mass wasting processes whereas bedrock surface sampling for crns may generally underestimate long term erosion rates response times of crn concentrations to changes in erosion rate indicate that climatically driven cycles of erosion may be detected but that complete equilibration of crn concentrations to new erosional conditions may take tens of thousands of years comparison of simulated crn derived erosion rates with a new data set of such rates from the nepalese himalaya underscore the conclusions drawn from the model results

beyond threshold hillslopes channel adjustment to base level fall in tectonically active mountain ranges numerous empirical and model based studies argue that in general hillslopes and river channels increase their gradients to accommodate high rates of base level fall to date however few data sets show the dynamic range of both these relationships needed to test theoretical models of hillslope evolution and river incision here we utilize concentrations of 10 be in quartz extracted from river sand on the eastern margin of the tibetan plateau to explore relationships among short term 10 2 10 5 a erosion rate hillslope gradient and channel steepness our data illustrate nonlinear behavior and a threshold in the relationship between erosion rate and mean hillslope gradient confi rming the generalization that hillslopes around the world are limited by slope stability and cease to provide a metric for erosion at high rates 0 2 mm a the relationship between channel steepness index and erosion rate is also nonlinear but channels continue to steepen beyond the point where threshold hillslopes emerge up to at least 0 6 mm a demonstrating that channel steepness is a more reliable topographic metric than mean hillslope gradient for erosion rate and that channels ultimately drive landscape adjustment to increasing rates of base level fall in tectonically active settings

large scale erosion rates from in situ produced cosmogenic nuclides in european river sediments abstract we have calculated long term erosion rates of 20 100 mm kyr from quartz contained 10be in the bedload of middle european rivers for catchments ranging from 102 to 105 km2 these rates average over 10 40 kyr and agree broadly with rock uplift incision and exhumation rates historic soil erosion rates and erosion rates calculated from the measured sediment loads of the same rivers moreover our new erosion rate estimates correlate well with lithology and relief however in the regen neckar loire and meuse catchments cosmogenic nuclide derived erosion rates are consistently 1 5 4 times greater than the equivalent rates derived from measured river loads this may be due to the systematic under representation of high magnitude low frequency transport events in the gauging records which cover less than a century alternatively the discrepancy may derive from spatially non uniform erosion and preferential tapping of deeper sections of the irradiation profile a third explanation relates the high cosmogenic nuclide derived erosion rates to inheritance of an elevated pleistocene erosion signal uncertainties associated with the cosmogenic nuclide derived erosion rate estimates are not greater than the potential errors in conventional estimates therefore the cosmogenic nuclide approach is an effective tool for rapid catchment wide assessment of time integrated rates of bedrock weathering and erosion and we anticipate its fruitful application to the quaternary sedimentary record

landscape form and millennial erosion rates in the san gabriel mountains ca abstract it has been long hypothesized that topography as well as climate and rock strength exert first order controls on erosion rates here we use detrital cosmogenic 10 be from 50 basins ranging in size from 1 to 150 km 2 to measure millennial erosion rates across the san gabriel mountains in southern california where a strong e w gradient in relief compared to weak variation in precipitation and lithology allow us to isolate the relationship between topographic form and erosion rate our erosion rates range from 35 to 1100 m ma and generally agree with both decadal sediment fluxes and long term exhumation rates inferred from low temperature thermochronometry catchment mean hillslope angle increases with erosion rate until 300 m ma at which point slopes become invariant with erosion rate although this sort of relation has been offered as support for non linear models of soil transport we use 1 d analytical hillslope profiles derived from existing soil transport laws to show that a model with soil flux linear in slope but including a slope stability threshold is indistinguishable from a non linear law within the scatter of our data catchment mean normalized channel steepness index increases monotonically though non linearly with erosion rate throughout the san gabriel mountains even where catchment mean hillslope angles have reached a threshold this non linearity can be mostly accounted for by a stochastic threshold incision model though additional factors likely contribute to the observed relationship between channel steepness and erosion rate these findings substantiate the claim that the normalized channel steepness index is an important topographic metric in active ranges

fission track and 40ar39ar evidence for episodic denudation of the gangotri granites in the garhwal higher himalaya india abstract in order to document quantitatively the cooling and denudation history of the higher himalayan granites bordering the tethyan sedimentary zone fission track ft apatite and 40 ar 39 ar mica ages have been determined on the gangotri leucogranites and biotite granites in the garhwal region of india gangotri is the source area of the ganges river and lies in the midst of the highest himalayan peaks in india a total of 15 apatite ages from a vertical profile 2580 4370 m on the gangotri granites yields ft ages in the range of 1 5 0 6 to 2 4 0 5 ma indicating that the rock column with a relief of 1800 m cooled through 130 10 c within only 1 million years during the late pliocene an average denudation rate of 2 mm yr is estimated for the past 2 4 million years from the gangotri granites we also report a muscovite 40 ar 39 ar age of 17 9 0 1 ma and a biotite age of 18 0 0 1 ma these reflect cooling of the rocks through 300 350 c probably related to an early miocene pulse of denudation caused by a basement cover detachment the martoli normal fault above the leucogranites time temperature pathways indicate that the cooling of the rocks in the late pliocene quaternary was five to six times the magnitude of cooling between 18 and 2 ma indicating a distinct pulse of rapid denudation in the late pliocene quaternary we interpret these young apatite ages and fast denudation as a geomorphic response increased erosion and cooling of the rocks to a major pulse of tectonic uplift in the higher himalaya shortly before 2 4 ma the effect of climatic cooling on this denudation is considered secondary to the role of tectonic forcing and indeed produced a positive feedback to the primary cause although our study is confined to the garhwal region it is probable that other granitic bodies of the higher himalaya bordering the tethyan sedimentary rocks and forming the loftiest summits in the himalaya have also experienced episodic denudation one major pulse in the early miocene which was mainly tectonic denudation and another in the late pliocene quaternary which was mainly erosional the latter is well recorded by apatite ft data and is consistent with the hypothesis that rapid uplift and denudation of the himalayan rocks may have influenced the initiation of the ice ages in the northern hemisphere

transient fluvial incision in the headwaters of the yellow river northeastern tibet china 1 we utilize topographic analysis of channel profiles combined with field measurements of erosion rates to explore the distribution of channel incision in the anyemaqen shan a broad mountainous region in the northeastern tibetan plateau tributary channels to the yellow river display systematic downstream increases in channel gradient associated with convex upward longitudinal profiles steep lower reaches of channels are associated with rapid 1 m ka incision rates along the yellow river while upstream reaches are associated with relatively slow 0 05 0 1 m ka erosion of soil mantled uplands covariance between erosion rates and channel steepness indices suggest that channels are adjusted to match long wavelength differential rock uplift across the range geologic constraints indicate that rapid incision downstream of the range is associated with excavation of basin fill driven by changes in relative base level farther downstream the upstream limit of this wave of transient incision is marked by a series of knickpoints that are found at nearly the same elevation throughout the watershed consistent with knickpoint migration as a kinematic rather than diffusional wave tributary channel gradients downstream of knickpoints however display a progressive adjustment to increased incision rates that may reflect the influence of increased sediment flux comparison of observed channel profiles to a stream power model of fluvial behavior reveals that the rate of knickpoint propagation can only be explained if the erosional efficiency coefficient k increases during incision our results thus highlight the utility of channel profile analysis to reconstruct the fluvial response to both active tectonism and external changes in base level

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

structure and shortening of the kangra and dehra dun reentrants sub himalaya india surface geology oil well seismic reflection and magnetostratigraphic data are integrated to evaluate the structural style and the shortening rate at the himalayan front sub himalaya of northwest india the sub himalaya between the main boundary thrust and the himalayan frontal fault is the primary surface expression of shortening between the himalaya and the indian plate at certain locations the himalayan frontal fault is a blind thrust beneath anticlines of siwalik tertiary molasse parallel to the himalayan orogen the main boundary thrust is sinuous so the width of the sub himalaya ranges from 30 to 80 km where the sub himalaya is narrow nahan salient tertiary rocks are exposed in imbricate thrust sheets where the sub himalaya is broad kangra and dehra dun reentrants alluvium fills wide synclinal valleys duns seismic reflection data reveal that surface anticlines form in association with south vergent thrusts that root in a decollement at the base of the tertiary section reflection profiles and well data also indicate that the basement lithology changes northward from precambrian crystalline rocks beneath the indo gangetic plains to precambrian and cambrian metasedimentary rocks beneath the sub himalaya the sub himalayan decollement dips 2 5 northward beneath the kangra reentrant but it is steeper 6 beneath the dehra dun reentrant the kangra and dehra dun reentrants display fault propagation folds having steep limbs in the north and fault propagation and fault bend folds that have gently north dipping limbs in the south a balanced cross section of the kangra reentrant shows that a minimum of 23 km shortening has occurred since 1 9 1 5 ma yielding a shortening rate of 14 2 mm yr shortening has occurred at a rate of 6 16 mm yr across the dehra dun reentrant these data are similar to other published shortening rates and indicate that approximately 25 of the total india eurasia convergence at this longitude is accommodated within the sub himalaya given continued convergence and the presence of overpressured wells in the kangra reentrant the region is likely at risk from moderate and or great earthquakes in the future

line of sight displacement from alos 2 interferometry mw 7 8 gorkha earthquake and mw 7 3 aftershock interferometric synthetic aperture radar insar is a key tool for the analysis of displacement and stress changes caused by large crustal earthquakes particularly in remote areas a challenge for traditional insar has been its limited spatial and temporal coverage especially for very large events whose dimensions exceed the typical swath width of 70 100 km this problem is addressed by the alos 2 satellite whose palsar 2 instrument operates in scansar mode enabling a repeat time of 2 weeks and a swath width of 350 km here we present insar line of sight displacement data from alos 2 palsar 2 observations covering the mw 7 8 gorkha nepal earthquake and its mw 7 3 aftershock that were acquired within 1 week of each event the data are made freely available and we encourage their use in models of the fault slip and associated stress changes the mw 7 3 aftershock not only extended the rupture area of the main shock toward the east but also left a 20 km gap where the fault has little or no coseismic slip we estimate this unslipped fault patch has the potential to generate a mw 6 9 event

on steady states in mountain belts the dynamic system of tectonics and erosion contains important feedback mechanisms such that orogenic systems tend toward a steady state this concept is often invoked but the nature of the steady state is commonly not specified we identify four types of steady state that characterize the orogenic system and illustrate these cases by using numerical model results and natural examples these types are 1 flux steady state 2 topographic steady state 3 thermal steady state and 4 exhumational steady state they refer to the erosional flux the topography the subsurface temperature field and the spatial pattern of cooling ages respectively models suggest that the topography will reach a steady mean form at the scale of an orogenic belt but perfect topographic steady state is unlikely to be achieved at shorter length scales thermal steady state is a precondition for exhumational steady state and in the case of temperature dependent deformation topographic steady state exhumational steady state is characterized by reset age zones spatially nested according to closure temperature as illustrated in natural systems from new zealand the cascadia accretionary margin and taiwan

erosion rates driven by channel network incision in the bolivian andes the bolivian andes flank one of earth s major topographic features and dominate sediment input into the amazon basin millennial scale erosion rates and dominant controls on erosion patterns in this range are poorly known to define these patterns we present 48 erosion rate estimates derived from analysis of in situ 10 be in quartz bearing alluvium collected from the upper beni river basin erosion rates corrected for the non uniform distribution of quartz in the sample basins range from 0 04 mm a 1 to 1 35 mm a 1 and thus integrate over 10 2 10 4 years mean and modal values are 0 42 standard deviation 0 29 and 0 2 0 4 mm a 1 respectively within the range of long term average erosion rates in this area derived from apatite fission track thermochronology 0 1 0 6 mm a 1 hence our data do not record any significant variation in erosion rate over the last several million years mean and modal short term erosion rates for the andes are an order of magnitude lower than rates in the ganges river headwaters in the high himalaya and an order of magnitude greater than rates typical of the european alps in the upper beni river region of the bolivian andes short term basin averaged erosion rates correlate with normalized channel steepness index a metric of relative channel gradient corrected for drainage area neither normalized channel steepness index nor basin averaged erosion rate shows strong correlation with mean basin hillslope gradient or mean basin local relief because many hillslopes in the upper beni river region are at threshold values of slope and local relief patterns of normalized channel steepness index appear primarily to reflect tectonic patterns and transient adjustment to those patterns by channel networks climate and lithology do not appear to exert first order controls on patterns of basin averaged erosion rates in the bolivian andes copyright 2005 john wiley sons ltd

seismotectonics of the himalayan collision zone geometry of the underthrusting indian plate beneath the himalaya fault plane solutions and well determined focal depths of medium sized earthquakes topography and landsat imagery in conjunction with seismicity maps cross sections and available geological information are used to investigate the present tectonics of the himalayan continental collision zone most of the accurately located epicenters of events along the himalayan arc 78 e 95 e that occurred between 1961 and 1981 are concentrated in a narrow zone about 50 km wide lying between the northerly dipping main boundary mbt and main central mct thrusts most of these events are located just south of the mct though the epicenters of the events are in general well located their depths as determined by teleseismic travel time data are very unreliable events with accurately determined depths obtained from identification of surface reflected phases define a simple planar zone from about 10 km and 20 km depth with an apparent dip of about 15 this result is all the more remarkable considering that the events used were located along about an 1800 km length of the himalyan arc except for one all available focal mechanisms of events within this zone indicate shallow 30 north dipping thrusts this shallow north dipping zone apparently defines a part of the detachment that separates the underthrusting indian plate from the lesser himalayan crustal block the spatial extent and the geometry of this interplate thrust zone strongly indicate that the mbt and nearby subsidiary surface and blind thrusts rather than the mct are currently the most active structures of the himalayan arc we suggest that the great himalayan earthquakes m 8 occur along the same detachment surface as defined by the thrust type medium sized events events located to the south of the mbt and beneath the ganges foredeep show normal faulting with t axes perpendicular to the himalayan trend the above results suggest that the indian continental plate is underthrusting the himalayan crustal blocks in a relatively coherent and simple geometry and that this geometry is not much different from that observed along oceanic subduction zones the november 19 1980 earthquake that occurred near the mct near 88 5 e shows a predominantly strike slip focal mechanism one of the nodal planes of this mechanism is transverse to the himalayan structural grain and moreover this plane has a trend similar to that of the recently mapped yadong gulu rift in the tethyan himalaya and in southern tibet just northeast of the earthquake we interpret this predominantly left lateral strike slip mechanism to indicate a possible genetic relationship between transverse structural features in the underthrusting indian plate the kishangang basement fault and the upper himalayan blocks and tibet

seismotectonics and crustal stress field in the kumaon garhwal himalaya abstract we present fault plane solutions of 94 well located small to moderate sized 1 5 ml 5 4 earthquakes which occurred in the kumaon garhwal himalaya during 2005 2008 using p wave polarity and body wave amplitudes these earthquakes show a mixture of thrust normal and strike slip type mechanism with a majority of thrust type most of the thrust earthquakes occur at a depth of 8 22 km in the main central thrust mct zone and the lower himalaya the spatial distribution of these earthquakes suggest that the strain resulting from the ongoing collision of the indian plate with the eurasian plate is being consumed by thrust fault movement mainly on the north dipping munsiari thrust and south dipping tons thrust the strike slip earthquakes are mainly observed in the lower himalaya as well as around the munsiari region in the mct zone the normal earthquakes are also observed in different parts of the kumaon garhwal himalaya and the gangetic plain their occurrence is attributed to the local structure s as well as the flexure of the indian plate stress tensor inversion of the calculated fault plane solutions indicates that the maximum compressive stress in the gangetic plain is n s directed and near vertical whereas in the kumaon garhwal himalaya it is near horizontal and nne ssw directed and correlating with the prevailing stress condition due to northward movement of indian plate

one dimensional reference velocity model and precise locations of earthquake hypocenters in the kumaon garhwal himalaya we report an optimum 1d velocity model of the kumaon garhwal himalaya based on analysis of 5000 seismograms generated at 50 digital broadband seismic stations in the region during april 2005 and june 2008 these include data from 385 local earthquakes with azimuth gaps of 180 each of which had a minimum of seven p and five s phase readings this velocity model was then used to map the locations of 1150 earthquakes 1 m 5 that occurred in the region a majority 75 of these earthquakes occurred in the upper 20 km of the crust largely above the main himalayan thrust and forms an 10 km wide band along the surface trace of the munsiari thrust however we also find another parallel band of earthquakes about 70 km to its southwest through the tehri dam in the lesser himalaya and a significant number 15 both in the tethys himalaya and within the flexed crust of the underthrusting indian plate beneath the ganga basin furthermore unlike the reported absence of seismicity in the lower crust of the nepal himalaya and its reappearance in the shallow mantle the kumaon garhwal himalayan crust is found to be seismically active throughout with no event reliably recorded in the shallow mantle of this region during the epoch of our seismic recording online material figures show constrains and distribution of data used in joint inversion and table shows p and s wave corrections obtained from travel time inversion

fashion and phases of late pleistocene aggradation and incision in the alaknanda river valley western himalaya india abstract we study the aggradation and incision of the alaknanda river valley during the late pleistocene and holocene the morphostratigraphy in the river valley at deoprayag shows the active riverbed a cut terrace and a fill terrace the sedimentary fabric of the fill terrace comprises four lithofacies representing 1 riverbed accretion 2 locally derived debris fan 3 the deposits of waning floods and 4 palaeoflood records the sedimentation style coupled with geochemical analysis and optically stimulated luminescence osl dating indicate that this terrace formed in a drier climate and the river valley aggraded in two phases during 21 18 ka and 13 9 ka during these periods sediment supply was relatively higher incision began after 10 ka in response to a strengthened monsoon and aided by increase of the tectonic gradient the cut terrace formed at 5 ka during a phase of stable climate and tectonic quiescence the palaeoflood records suggest wetter climate 200 300 yr ago when the floods originated in the upper catchment of the higher himalaya and in the relatively drier climate 1 2 ka when locally derived sediments from the lesser himalaya dominated flood deposits maximum and minimum limits of bedrock incision rate at deoprayag are 2 3 mm a and 1 4 mm a

convergence across the northwest himalaya from gps measurements 1 horizontal velocities of 26 global positioning system gps stations in the northwest himalayan region provide new constraints on the partitioning of india eurasia convergence and elastic strain accumulation about the locked main frontal thrust mft the northwest striking karakorum fault slips at 11 4 mm yr and contributes to east west extension of southern tibet and westward motion of the northwest himalaya towards nanga parbat rather than playing a role in eastward extrusion of tibet crustal shortening across the himalaya occurs within a zone centered about 100 km north of the siwalik foothills and the mft model inversions of the gps data indicate that the mft is locked over a width of 100 km comparison with geologic mft slip rate estimates suggests that this zone is building up a slip deficit at a rate of 14 1 mm yr and will eventually fail in future great earthquakes

fission track dating of phosphate minerals and the thermochronology of apatite several phosphate minerals have been investigated for their usefulness in geochronology and thermochronology by the fission track method of these apatite ca5 po4 3 f cl oh has proved pre eminently suitable for this purpose for reasons discussed below apatite was one of the first minerals amongst many others to be investigated for fission track dating by fleischer and price 1964 and has subsequently become by far the most important of all the minerals used for dating by this method the usefulness of apatites for fission track dating arises from their near universal tendency to concentrate uranium within their structure at the time of crystallization and their widespread occurrence in all of the major rock groups price and walker 1963 first recognized the possibility that the accumulation of radiation damage tracks in natural minerals from the spontaneous nuclear fission of 238u within their lattices could be used for geological dating they also demonstrated that a simple chemical etching procedure served to enlarge these fission tracks to optical dimensions so that they could be observed and measured under an ordinary optical microscope price and walker 1962b 1963 this simple procedure quickly opened the way for a wide range of nuclear particle track studies in natural minerals and glasses e g fleischer et al 1975 the most important of which was fission track dating fleischer et al 1964 also showed that spontaneous or fossil fission tracks provided an explanation for various anomalous etch pits in apatite which had puzzled crystallographers over many years the essence of a fission track age determination involves measuring the number of tracks that have accumulated over the lifetime of the mineral along with an estimate of the amount of uranium that is present knowing the rate of spontaneous fission decay a geological age can be calculated detailed

landscape response to tectonic forcing digital elevation model analysis of stream profiles in the mendocino triple junction region northern california the topographic evolution of orogens is fundamentally dictated by rates and patterns of bedrock channel incision quantitative field assessments of process based laws are needed to accurately describe landscape uplift and denudation in response to tectonics and climate we evaluate and calibrate the shear stress or similar unit stream power bedrock incision model by studying stream profiles in a tectonically active mountain range previous work on emergent marine terraces in the mendocino triple junction region of northern california provides spatial and temporal control on rock uplift rates digital elevation models and field data are used to quantify differences in landscape morphology associated with along strike northwest to southeast changes in tectonic and climatic conditions analysis of longitudinal profiles supports the hypothesis that the study area channels are in equilibrium with current uplift and climatic conditions consistent with theoretical calculations of system response time based on the shear stress model within uncertainty the profile concavity of the trunk streams is constant throughout the study area 0 43 as predicted by the model channel steepness correlates with uplift rate these data help constrain the two key unknown model parameters the coefficient of erosion k and the exponent associated with channel gradient n this analysis shows that k cannot be treated as a constant throughout the study area despite generally homogeneous substrate properties for a reasonable range of slope exponent values n best fit values of k are positively correlated with uplift rate this correlation has important implications for landscape evolution models and likely reflects dynamic adjustment of k to tectonic changes due to variations in orographic precipitation and perhaps channel width sediment load and frequency of debris flows the apparent variation in k makes a unique value of n impossible to constrain with present data

widespread aggradation in the mountainous catchment of the alaknanda ganga river system timescales and implications to hinterland foreland relationships abstract cut and fill type fluvial terraces are ubiquitous in the lesser himalayan zone of the alaknanda ganga ganges rivers system which flows perpendicular to himalayan lithotectonic units and traverses a steep climatic gradient lithofacies analysis of the sedimentary sequences of cut and fill terraces indicated that the valley aggradation took place via 1 channel bar development and excess sediment supply 2 debris flows composed of mixed rounded to sub rounded lithoclasts resulting from episodic high intensity rainfalls in the upper catchment or 3 debris flows or rockfalls generated by local landslides the luminescence chronology indicates that valley aggradation took place in two phases 49 25 ka and 18 11 ka the incision of the fill started soon after 11 ka paleoclimatic records from marine sediments indicate that the aggradation and incision in the alaknanda ganga river has oscillated in phase with global climatic variations glaciation deglaciation processes in the upper catchment produced huge amounts of sediment between 63 and 11 ka which was fluvially transferred to the lower valley via several cycles of erosion and deposition leading to extensive aggradation the climatic amelioration at 11 ka and the completion of deglaciation processes led to increased fluvial discharge and decreased sediment supply conditions conducive for incision of the alluvial fills records from the indo gangetic plain and the ganga delta demonstrate that the phase of aggradation was regional but that incision in the foreland started at least 2 3 ka later after 7 ka bedrock incision rates as calculated from dated alluvial covers of terraces that are separated by bedrock steps are spatially variable and fall within the range of rates reported from across the himalaya these estimated rates however are higher than the basin average erosion rates calculated using isotopic mass balance in riverbed sediments this study suggests that during the last 50 ka river dynamics in the himalayas were dominated by monsoon variability and the role of tectonic activity was limited to bedrock incision in few reaches only

has focused denudation sustained active thrusting at the himalayan topographic front the geomorphic character of major river drainages in the himalayan foothills of central nepal suggests the existence of a discrete west northwest trending break in rock uplift rates that does not correspond to previously mapped faults the 40 ar 39 ar thermochronologic data from detrital muscovites with provenance from both sides of the discontinuity indicate that this geomorphic break also corresponds to a major discontinuity in cooling ages samples to the south are proterozoic to paleozoic whereas those to the north are miocene and younger combined these observations virtually require recent pliocene holocene motion on a thrust sense shear zone in the central nepal himalaya 20 30 km south of the main central thrust field observations are consistent with motion on a broad shear zone subparallel to the fabric of the lesser himalayan lithotectonic sequence the results suggest that motion on thrusts in the toe of the himalayan wedge may be synchronous with deeper exhumation on more hinterland structures in central nepal we speculate that this continued exhumation in the hinterland may be related to intense sustained erosion driven by focused orographic precipitation at the foot of the high himalaya

dominance of tectonics over climate in himalayan denudation landscape denudation in actively deforming mountain ranges is controlled by a combination of rock uplift and surface runoff induced by precipitation whereas the relative contribution of these factors is important to our understanding of the evolution of orogenic topography no consensus currently exists concerning their respective infl uences to address this question denudation rates at centennial to millennial time scales were deduced from 10be concentrations in detrital sediments derived from 30 small basins 10 600 km2 in an 200 km wide region in central nepal along a northward strike perpendicular transect average denudation rates sharply increase from 0 5 mm yr in the lesser himalayas to 1 mm yr when crossing the physiographic transition and then accelerate to 2 3 mm yr on the southern flank of the high peaks in the greater himalayas despite a more than five fold increase in denudation rate between the southern and northern parts of this transect the corresponding areas display similar precipitation rates the primary parameter that presents a signifi cant co variation with denudation is the long term rockuplift rate that is interpreted to result from the ramp fl at transition along the main himalayan thrust we propose that in this rapidly uplifting mountain range landscapes adjust quickly to changing climatic conditions such that denudation is mainly limited by the rate at which material is pushed upward by tectonic processes and made available for removal by surface processes in this particular context variations in precipitation appear to have only a second order role in modulating the denudation signal that is primarily set by the background rock uplift rate 2014 geological society of america

dynamics of the 2015 m7 8 nepal earthquake the 2015 m7 8 nepal earthquake ruptured part of the main himalayan thrust beneath kathmandu to study the dynamics of this event we compute p wave spectra of the main shock and of two large aftershocks to estimate stress drop and radiated energy we find that surface reflections depth phases of these shallow earthquakes produce interference that severely biases spectral measurements unless corrections are applied measures of earthquake dynamics for the main shock are within the range of estimates from global and regional earthquakes we explore the azimuthal and temporal variations of radiated energy and highlight unique aspects of the m7 8 rupture the beginning of the earthquake likely experienced a dynamic weakening mechanism immediately followed by an abrupt change in fault geometry correlation of backprojection results with frequency dependent variations in the radiated energy rate and with the suggested geometry of the main himalayan thrust yields new constraints on dynamic ruptures through geometrical barriers

aftershocks of the march 1999 chamoli earthquake and seismotectonic structure of the garhwal himalaya the aftershocks of the chamoli earthquake m b 6 3 of 28 march 1999 are analyzed to understand the seismotectonic structure of the garhwal himalaya the permanent as well as temporary microearthquake network data are used the mainshock occurred on the basement thrust by thrust faulting the aftershocks on the other hand are triggered at the seismogenic faults to the south of the main central thrust mct by thrust faulting as well as by strike slip faulting and all the events occurred above the plane of detachment basement thrust the mct is not the seismogenic fault

southward extrusion of tibetan crust and its effect on himalayan tectonics the tibetan plateau is a storehouse of excess gravitational potential energy accumulated through crustal thickening during india asia collision and the contrast in potential energy between the plateau and its surroundings strongly influences the modern tectonics of south asia the distribution of potential energy anomalies across the region derived from geopotential models indicates that the himalayan front is the optimal location for focused dissipation of excess energy stored in the plateau the modern pattern of deformation and erosion in the himalaya provides an efficient mechanism for such dissipation and a review of the neogene geological evolution of southern tibet and the himalaya shows that this mechanism has been operational for at least the past 20 million years this persistence of deformational and erosional style suggests to us that orogens like other complex systems can evolve toward steady state configurations maintained by the continuous flow of energy the capacity of erogenic systems to self organize into temporally persistent structural and erosional patterns suggests that the tectonic history of a mountain range may depend on local energetics as much as it does on far field plate interactions

seismotectonics of the nepal himalaya from a local seismic network the national seismological network of nepal consists of 17 short period seismic stations operated since 1994 it provides an exceptional view of the microseismic activity over nearly one third of the himalayan arc including the only segment between longitudes 78 e and 85 e that has not produced any m 8 earthquakes over the last century it shows a belt of seismicity that follows approximately the front of the higher himalaya with most of the seismic moment being released at depths between 10 and 20 km this belt of seismicity is interpreted to reflect interseismic stress accumulation in the upper crust associated with creep in the lower crust beneath the higher himalaya the seismic activity is more intense around 82 e in far western nepal and around 87 e in eastern nepal western nepal between 82 5 and 85 e is characterized by a particularly low level of seismic activity we propose that these lateral variations are related to segmentation of the main himalayan thrust fault the major junctions between the different segments would thus lie at about 87 e and 82 e with possibly an intermediate one at about 85 e these junctions seem to coincide with some of the active normal faults in southern tibet lateral variation of seismic activity is also found to correlate with lateral variations of geological structures suggesting that segmentation is a long lived feature we infer four 250 400 km long segments that could produce earthquakes comparable to the m 8 4 bihar nepal earthquake that struck eastern nepal in 1934 assuming the model of the characteristic earthquake the recurrence interval between two such earthquakes on a given segment is between 130 and 260 years

assessing quaternary reactivation of the main central thrust zone central nepal himalaya new thermochronologic data and numerical modeling we study the recent dynamics of the central nepal himalaya focusing on possible reactivation of the footwall of the main central thrust which is marked by an abrupt topographic transition different tectonic mechanisms such as overthrusting of a major crustal ramp underplating or out of sequence thrusting have been suggested to explain the morphology and exhumation patterns in this area we present 25 new apatite fission track ages collected along a north south transect in central nepal as well as two age elevation profiles ages are consistently younger than 3 ma old in the main central thrust zone and increase continuously to 4 6 ma old in the south no jump in apatite fission track ages is observed across the topographic transition apparent exhumation rates from age elevation relationships vary from 0 46 0 13 0 09 km ma in the palung granite south of kathmandu to 4 4 4 8 1 5 km ma in the main central thrust zone the latter rate is probably overestimated by a factor of two due to topographic effects as shown by a new numerical model these strongly varying exhumation rates can be explained by overthrusting of a crustal ramp which exerts a primary control on age patterns and do not require out of sequence reactivation of thrusts in the main central thrust zone

constraints on the tectonic and landscape evolution of the bhutan himalaya from thermochronometry the observed geomorphology and calculated thermal histories of the bhutan himalaya provide an excellent platform to test ideas regarding the influence of tectonics and climate on the evolution of a convergentmountain range however little consensus has been reached regarding the late cenozoic history of the bhutan himalaya some researchers have argued that observed geologic relationships show slowing deformation rates such that the range is decaying from a geomorphic perspective while others see the range as growing and steepening we suggest that a better understanding is possible through the integrated interpretation of geomorphic and thermochronometric data from the comparison of predictions from models of landscape evolution and thermal kinematic models of orogenic systems new thermochronometric data throughout bhutan aremost consistent with a significant decrease in erosion rates from2 to 3 km ma down to 0 1 0 3 km ma around 6 4ma we interpret this pattern as a decrease in rock uplift rates due to the activation of contractional structures of the shillong plateau an uplifted region approximately 100 km south of bhutan however low relief fluvial landscapes throughout the bhutanese hinterland record a late pulse of surface uplift likely due to a recent increase in rock uplift rates constraints from our youngest thermochronometers suggest that this increase in rock uplift and surface uplift occurred within the last 1 75ma these results imply that the dynamics of the bhutan himalaya and shillong plateau have been linked during the late cenozoic with structural elements of both regions active in variable ways and times over that interval

the status of central seismic gap a perspective based on the spatial and temporal aspects of the large himalayan earthquakes abstract the central himalaya considered as a prominent seismic gap is generally believed to be the most vulnerable segment due for a great plate boundary earthquake m 8 two significant historical earthquakes are known to have occurred in this region in ad 1505 and 1803 interpreted by some workers as great earthquakes occurrences of these events question the validity of central himalaya being a seismic gap here we analyze the macroseismal data from the central himalaya and the gangetic plains and suggest that neither of these earthquakes can be qualified as a plate boundary type event and the size of at least the 1803 does not exceed mw 7 7 our studies further indicate that the central himalayan frontal thrusts may be undergoing a quiescence of 1000 years in terms of generation of great earthquakes assuming a steady elastic strain accumulation due to plate convergence and its simple translation as periodic great earthquakes such long quiescence may appear intriguing we suggest that the long term quiescence and the spatial and temporal clustering that characterize the distribution of large earthquakes in the himalaya is a manifestation of the tectonic deformation associated with a highly evolved fold and thrust belt our analyses also indicate that the large moderate earthquakes are mostly concentrated on higher level thrusts these structures probably accommodate most of the seismogenic slip in the himalaya and the rest of the strain could possibly be expended on passive folding consequently great earthquakes m 8 are relatively infrequent as in the case of central himalaya apparently leading to long periods of quiescence

thrust geometries and deep structure of the outer and lesser himalaya kumaon and garhwal india implications for evolution of the himalayan fold and thrust belt the kumaon garhwal region of the himalaya lies near the center of the himalayan fold and thrust belt we have drawn two balanced cross sections 100 km apart through the outer and lesser himalaya the cross sections incorporate all the surface well log and earthquake seismic data currently available from the region two branch line maps showing trailing and leading branch lines and cutoff lines of the major thrusts in the region are also drawn the three dimensional deep structure of the outer and lesser himalaya is interpreted based on the balanced cross sections and the branch line maps deep structure of the higher and tethyan himalaya is extrapolated based on surface geology and is subject to revision as more surface and seismic data become available from these areas a sequential evolutionary model for the kumaon himalaya along the eastern pindari section is proposed according to this model the kumaon himalaya evolved by an overall forelandward progression of thrusting with some reactivation along the munsiari thrust mt the main boundary thrust mbt and the main central thrust mct we use structural stratigraphic and radiometric criteria to place time constraints on the motion of these thrusts earliest motion along the mbt may have occurred in early middle paleocene but the main episode probably started in late eocene and may still be continuing emplacement of the mt had occurred by middle late eocene whereas the mct shows activity around 20 ma thus exhibiting break back thrusting with respect to the mt shortening estimates are obtained from the pindari section minimum shortening in the sedimentary thrust sheets of the outer and lesser himalaya is 161 km or 65 as a first approximation we have also restored the crystalline sheets in order to obtain shortening estimates for the entire himalaya minimum shortening for the himalaya after restoring the mct sheet varies from 354 76 to 421 km 79 these estimates were further combined with published data from the area between the mct sheet and the indus tsangpo suture zone itsz minimum shortening between the indo ganga foreland and the itsz thus obtained lies in the range 687 754 km or 69 72 we compare our shortening estimates with those available from the pakistan and nepal himalaya

erosion and exhumation in the himalaya from cosmogenic isotope inventories of river sediments the outward erosional flux is a key factor in the tectonic evolution of mountain belts and there is much debate about the feedbacks between tectonics erosion and climate here we use cosmogenic nuclides 10be and 26al analysed in quartz from river sediments from the upper ganges catchment to make the first direct measurements of large scale erosion rates in a rapidly uplifting mountain belt the erosion rates are highest in the high himalaya at 2 7 0 3 mm yr 1 errors fall to 1 2 0 1 mm yr on the southern edge of the tibetan plateau and are 0 8 0 3 to 0 6 mm yr in the foothills to the south of the high mountains these relative estimates are corroborated by the nd isotopic mass balance of the river sediment analysis of sediment from an abandoned terrace suggests that similar erosion rates have been maintained for at least the last few thousand years the data presented here along with data recently published for european river catchments demonstrate that a log linear relationship between relief and erosion rate holds over three orders of magnitude variation in erosion rate and between very different climatic and tectonic regimes the erosion rate estimates from cosmogenic nuclides correlate well with exhumation rates calculated from previously published apatite fission track ages in the indian himalaya this confirms that much of the exhumation in the himalayan mountain chain is now balanced by erosion however exhumation rates calculated from high blocking temperature systems such as 40ar 39ar in muscovite imply lower exhumation rates rocks presently at the surface must have undergone a three to six fold increase in exhumation rate within the last few million years we show how this could be explained either by climatic forcing of erosion rate changes or by tectonics published evidence for equally rapid changes of exhumation rate in the past and the probable diachroneity in the time at which the present exhumation rates accelerated imply that tectonics has moderated at least some of the change in exhumation rates

assessment of the martian gravity field at short wavelength with mars express 1 the gravity part of the mars express radio science experiment consists in measuring gravity data near pericenter above selected target areas of geophysical interest the low altitude of the mars express at pericenter 263 329 km makes it a very sensitive gravity sensor at small wavelengths which can give new constraints on the local structure of the crust and lithosphere mars express gravity data can also be used to check the quality of the existing global gravity solutions in this paper we show that mars express gravity data confirm in the tharsis area the validity of existing global gravity solutions in particular the resolution of the most recent global gravity solution is excellent up to harmonic degree 73 and the amplitude at short wavelength is not biased toward zero by the power law regularization up to the same degree

thrusting and back thrusting as post emplacement kinematics of the almora klippe insights from low temperature thermochronology abstract crystalline klippen over the lesser himalayan sequence lhs in the kumaon and garhwal regions of nw himalaya are the representative of southern portion of the main central thrust mct hanging wall these were tectonically transported over the juxtaposed thrust sheets berinag tons and ramgarh of the lhs zone along the mct these klippen comprise of nw se trending synformal folded thrust sheet bounded by thrusts in the south and north in the present study the exhumation histories of two well known klippen namely almora and baijnath and the ramgarh thrust sheet in the kumaon and garhwal regions vis a vis himalayan orogeny have been investigated using apatite fission track aft ages along a 60 km long orogen perpendicular transect across the almora klippe and the ramgarh thrust sheet 16 aft cooling ages from the almora klippe and 2 from the ramgarh thrust sheet have been found to range from 3 7 0 8 to 13 2 2 7 ma and 6 3 0 8 to 7 2 1 0 ma respectively from lhs meta sedimentary rocks only a single aft age of 3 6 0 8 ma could be obtained three aft ages from the baijnath klippe range between 4 7 0 5 and 6 6 0 8 ma aft ages and exhumation rates of different klippen show a dynamic coupling between tectonic and erosion processes in the kumaon and garhwal regions of nw himalaya however the tectonic processes play a dominant role in controlling the exhumation thrusting and back thrusting within the almora klippe and ramgarh thrust sheet are the post emplacement kinematics that controlled the exhumation of the almora klippe combining these results with the already published aft ages from the crystalline klippen and the higher himalayan crystalline hhc the kinematics of emplacement of the klippen over the lhs and exhumation pattern across the mct in the kumaon and garhwal regions of nw himalaya have been investigated

decoupling of erosion and precipitation in the himalayas the hypothesis that abrupt spatial gradients in erosion can cause high strain rates in active orogens has been supported by numerical models that couple erosional processes with lithospheric deformation via gravitational feedbacks1 2 3 most such models invoke a stream power rule in which either increased discharge or steeper channel slopes cause higher erosion rates spatial variations in precipitation and slopes are therefore predicted to correlate with gradients in both erosion rates and crustal strain here we combine observations from a meteorological network across the greater himalaya nepal along with estimates of erosion rates at geologic timescales greater than 100 000 yr from low temperature thermochronometry across a zone of about 20 km length spanning the himalayan crest and encompassing a more than fivefold difference in monsoon precipitation significant spatial variations in geologic erosion rates are not detectable decreased rainfall is not balanced by steeper channels instead additional factors that influence river incision rates such as channel width and sediment concentrations must compensate for decreasing precipitation overall spatially constant erosion is a response to uniform upward tectonic transport of greater himalayan rock above a crustal ramp

frictional slip of granite at hydrothermal conditions sliding on faults in much of the continental crust likely occurs at hydrothermal conditions i e at elevated temperature and elevated pressure of aqueous pore fluids yet there have been few relevant laboratory studies to measure the strength sliding behavior and friction constitutive properties of faults at hydrothermal conditions we slid laboratory granite faults containing a layer of granite powder simulated gouge velocity stepping experiments were performed at temperatures of 23 to 600 c pore fluid pressures ph2o of 0 dry and 100 mpa wet effective normal stress of 400 mpa and sliding velocities v of 0 01 to 1 m s 0 32 to 32 m yr conditions were similar to those in earlier tests on dry granite to 845 c by lockner et al 1986 the mechanical results define two regimes the first regime includes dry granite up to at least 845 and wet granite below 250 c in this regime the coefficient of friction is high 0 7 to 0 8 and depends only modestly on temperature slip rate and ph2o the second regime includes wet granite above 350 c in this regime friction decreases considerably with increasing temperature temperature weakening and with decreasing slip rate velocity strengthening these regimes correspond well to those identified in sliding tests on ultrafine quartz we infer that one or more fluid assisted deformation mechanisms are activated in the second hydrothermal regime and operate concurrently with cataclastic flow slip in the first cool and or dry regime is characterized by pervasive shearing and particle size reduction slip in the second hot and wet regime is localized primarily onto narrow shear bands adjacent to the gouge rock interfaces weakness of these boundary shears may result either from an abundance of phyllosilicates preferentially aligned for easy dislocation glide or from a dependence of strength on gouge particle size major features of the granite data set can be fit reasonably well by a rate and temperature dependent three regime friction constitutive model chester this issue we extrapolate the experimental data and model fit in order to estimate steady state shear strength versus depth along natural slipping faults for sliding rates as low as 31 mm yr we do this for two end member cases in the first case pore pressure is assumed hydrostatic at all depths shallow crustal strength in this case is similar to that calculated in previous work from room temperature friction data while at depths below about 9 13 km depending on slip rate strength becomes less sensitive to depth but sensitive to slip rate in the second case pore pressure is assumed to be near lithostatic at depths below 5 km strength is low at all depths in this case 20 mpa in agreement with observations of weak faults such as the san andreas the predicted depth of transition from velocity weakening to velocity strengthening lies at about 13 km depth for a slip rate of 31 mm yr in rough agreement with the seismic aseismic transition depth observed on mature continental faults these results highlight the importance of fluid assisted deformation processes active in faults at depth and the need for laboratory studies on the roles of additional factors such as fluid chemistry large displacements higher concentrations of phyllosilicates and time dependent fault healing

slip model of the 2015 mw 7 8 gorkha nepal earthquake from inversions of alos 2 and gps data we use surface deformation measurements including interferometric synthetic aperture radar data acquired by the alos 2 mission of the japanese aerospace exploration agency and global positioning system gps data to invert for the fault geometry and coseismic slip distribution of the 2015 mw 7 8 gorkha earthquake in nepal assuming that the ruptured fault connects to the surface trace of the main frontal thrust mft fault between 84 34 e and 86 19 e the best fitting model suggests a dip angle of 7 the moment calculated from the slip model is 6 08 1020 nm corresponding to the moment magnitude of 7 79 the rupture of the 2015 gorkha earthquake was dominated by thrust motion that was primarily concentrated in a 150 km long zone 50 to 100 km northward from the surface trace of the main frontal thrust mft with maximum slip of 5 8 m at a depth of 8 km data thus indicate that the 2015 gorkha earthquake ruptured a deep part of the seismogenic zone in contrast to the 1934 bihar nepal earthquake which had ruptured a shallow part of the adjacent fault segment to the east

doppler ducting of short period gravity waves by midlatitude tidal wind structure 1 multiwavelength airglow image data depicting a short period 4 9 min atmospheric gravity wave characterized by a sharp leading front have been analyzed together with synoptic meteor radar wind data recorded simultaneously from bear lake observatory utah 41 6 n 111 6 w the wind data suggest the presence of a semidiurnal tide with horizontal winds peaking at around 60 m s along the sse direction of motion 170 from north of this short period wave it was found that the gravity wave was most probably ducted because of the doppler shift imposed by this wind structure a marked 180 phase shift was observed between the near infrared oh and the oi 557 7 nm emissions numerical simulation results for similar ducted waves excited by idealized model sources suggest that the phase shift between the wave modulated airglow intensities may be explained simply by chemical processes rather than by wave dynamics phase velocities of simulated waves however appear higher than those of observed waves suggesting the importance of tidal thermal structure in determining the doppler ducted wave characteristics

analysis and characterization of the vertical accuracy of digital elevation models from the shuttle radar topography mission 1 the first near global high resolution digital elevation model dem of the earth has recently been released following the successful shuttle radar topography mission srtm of 2000 this data set will have applications in a wide range of fields and will be especially valuable in the earth sciences prior to widespread dissemination and use it is important to acquire knowledge regarding the accuracy characteristics in this work a comprehensive analysis of the vertical errors present in the data set and the assessment of their effects on different hydrogeomorphic products is performed in particular the work consisted of 1 measuring the vertical accuracy of the data set in two areas with different topographic characteristics 2 characterizing the error structure by comparing elevation residuals with terrain attributes 3 assessing a wavelet based filter for removing speckle and 4 assessing the effects of vertical errors on hydrogeomorphic products and on slope stability modeling the results indicate that in the two sites relief has a strong effect on the vertical accuracy of the srtm dem in the high relief terrain large errors and data voids are frequent and their location is strongly influenced by topography while in the low to medium relief site errors are smaller although the hilly terrain still produces an effect on the sign of the errors speckling generates deviations in the drainage network in one of the investigated areas but the application of a wavelet filter proved to be an effective tool for removing vertical noise although further fine tuning is necessary vertical errors cause differences in automatically extracted hydrogeomorphic products that range between 4 and 1090

a signature of transience in bedrock river incision rates over timescales of 10 4 10 7 years measured rates of river incision into bedrock are commonly interpreted as proxies for rates of rock uplift see refs 1 and 2 for example and indices of the strength of climatic forcing of erosion over time see refs 3 and 4 for example this approach implicitly assumes that river incision rates are in equilibrium with external forcings over a wide range of timescales here we directly test this assumption by examining the temporal scaling of bedrock river incision from 155 independent measurements of river incision compiled from 14 sites of these sites 11 exhibit a negative power law dependence of bedrock river incision rate on measurement interval a relationship that is apparent over timescales of 10 4 10 7 years and is independent of tectonic and geomorphic setting thus like rates of sediment accumulation rates of river incision into bedrock exhibit non steady state behaviour even over very long measurement intervals non steady state behaviour can be explained by episodic hiatuses in river incision triggered by alluvial deposition if such hiatuses have a heavy tailed length distribution regardless of its cause the dependence of incision rate on measurement interval complicates efforts to infer tectonic or climatic forcing from changes in rates of river incision over time or from comparison of rates computed over different timescales

current strain accumulation in the hinterland of the northwest himalaya constrained by landscape analyses basin wide denudation rates and low temperature thermochronology abstract rupture associated with the 25 april 2015 m w 7 8 gorkha nepal earthquake highlighted our incomplete understanding of the structural architecture and seismic cycle processes that lead to himalayan mountain building in central nepal in this paper we investigate the style and kinematics of active mountain building in the himalayan hinterland of northwest india approximately 400 km to the west of the hypocenter of the nepal earthquake via a combination of landscape metrics and long ma and short term ka erosion rate estimates from low temperature thermochronometry and basin wide denudation rate estimates from 10 be concentrations we focus our analysis on the area straddling the pt 2 the physiographic transition between the lesser and high himalaya that has yielded important insights into the nature of hinterland deformation across much of the himalaya our results from northwest india reveal a distinctive pt 2 that separates a lesser himalaya region with moderate relief 1000 m and relatively slow erosion 2 has persisted in the same relative position since at least the past 1 5 ma we interpret these observations to suggest that strain accumulation in this hinterland region throughout at least the past 1 5 ma has been accomplished both by crustal thickening via duplexing and overthrusting along transient emergent faults despite the 400 km distance between them similar spatiotemporal patterns of erosion and deformation observed in northwest india and central nepal suggest both regions experience similar styles of active strain accumulation and both are susceptible to large seismic events

impact of glacial erosion on 10 be concentrations in fluvial sediments of the marsyandi catchment central nepal several processes contribute to denudation in high mountain environments of these glacial erosion is particularly difficult to constrain despite its critical importance in the evolution of many mountain ranges in this study we present a new data set of 10be concentrations in fluvial sediments sampled along the marsyandi river and its main tributaries in central nepal we interpret the 10be concentrations as being significantly impacted by glacially derived sediments along the marsyandi river such additions complicate conventional interpretations of 10be derived catchment scale denudation rates using a simple linear mass conservation formulation we invert our data set in order to separate the different denudational contributions to the observed signal as well as to constrain their magnitude and spatial distribution our results suggest significant variations in glacial erosion both in space and magnitude within the marsyandi catchment

characterizing the main himalayan thrust in the garhwal himalaya india with receiver function ccp stacking abstract we use common conversion point ccp stacking of ps receiver functions to image the crustal structure and moho of the garhwal himalaya of india our seismic array of 21 broadband seismometers spanned the himalayan thrust wedge at 79 80 e between the main frontal thrust and the south tibet detachment in 2005 2006 our ccp image shows the main himalayan thrust mht the detachment at the base of the himalayan thrust wedge with a flat ramp flat geometry seismic impedance contrasts inferred from geologic cross sections in garhwal imply a negative impedance contrast velocity decreasing downward for the upper flat located beneath the lower himalaya and a positive impedance contrast velocity increasing downward for the ramp located beneath the surface trace of the munsiari thrust or mct i at the lower flat located beneath the higher himalaya spatially coincident measurements of very high electrical conductivities require the presence of free fluids and we infer a negative impedance contrast on the mht caused by ponding of these fluids beneath the detachment our seismic image indicates that the upper flat of the mht is 10 km below sea level and dips north at 2 connecting to a mid crustal ramp which is 10 km high and dips at 16 the lower flat is 20 25 km below sea level and dips at 4 the main central thrust mct appears as a negative impedance contrast dipping at 16 the moho is nearly horizontal at 35 45 km depth beneath the sub himalaya and lower himalaya deepening to 50 km or more beneath the higher himalaya this depth is 10 25 km shallower than in the nw indian himalaya and 5 10 km shallower than in central nepal requiring significant along strike variations in crustal thickness the observed thickness of subducted indian crust in garwhal is 20 28 km

discharge discharge variability and the bedrock channel profile long term bedrock incision is driven by daily discharge events of variable magnitude and frequency with ineffective events below an incision threshold we explore theoretically how this short term stochastic behavior controls long term steady state incision rates and bedrock channel profiles combining a realistic frequency magnitude distribution of discharge with a deterministic detachment limited incision model in which incision rate is a power function of basal shear stress above a critical shear stress our model predicts a power law relationship between steady state slope and drainage area consistent with observations the exponent of this power law is independent of discharge mean and variability while the amplitude factor which controls mountain belt relief is a power law function of mean runoff with an exponent of 0 5 and a complex function of runoff variability in accordance with evidence that incision occurs between 6 and 20 of time in rapidly incising rivers 1 mm yr our model predicts that channel steepness is virtually insensitive to runoff variability runoff variability can only decrease channel steepness for very slow incision rates and or weak lithologies the relationship between channel steepness and incision rate is always a power law whose exponent depends on the channel cross sectional geometry and runoff variability this contradicts models neglecting discharge stochasticity in which the steepness incision scaling is set by the incision law exponent our results suggest that changes in climate variability cannot explain an increase in bedrock incision rates during the late cenozoic within the context of a detachment limited model

spatially averaged long term erosion rates measured from in situ produced cosmogenic nuclides in alluvial sediment spatially averaged erosion rates of small catchments can be accurately inferred from the concentrations of cosmogenic nuclides in stream sediment here we test this technique at two catchments by comparing erosion rates inferred from cosmogenic nuclides with rates of alluvial fan deposition over the past 16 000 years these two independent estimates agree within one standard error demonstrating that cosmogenic nuclide signatures of stream sediment can be used to measure spatially averaged long term erosion rates using this technique we show that long term erosion rates are an exponential function of average hillslope gradient at these sites

earthquake recurrence and rupture dynamics of himalayan frontal thrust india the black mango fault is a structural discontinuity that transforms motion between two segments of the active himalayan frontal thrust hft in northwestern india the black mango fault displays evidence of two large surface rupture earthquakes during the past 650 years subsequent to 1294 a d and 1423 a d and possibly another rupture at about 260 a d displacement during the last two earthquakes was at minimum 4 6 meters and 2 4 to 4 0 meters respectively and possibly larger for the 260 a d event abandoned terraces of the adjacent markanda river record uplift due to slip on the underlying hft of 4 8 0 9 millimeters per year or greater since the mid holocene the uplift rate is equivalent to rates of fault slip and crustal shortening of 9 6 3 5 7 0 millimeters per year and 8 4 3 6 7 3 millimeters per year respectively when it is assumed that the hft dips 30 10

wasp waisted hysteresis loops mineral magnetic characteristics and discrimination of components in mixed magnetic systems rock magnetic studies of complex systems that contain mixtures of magnetic minerals or mixed grain size distributions have demonstrated the need for a better method of distinguishing between different magnetic components in geological materials hysteresis loops that are constricted in the middle section but are wider above and below the middle section are commonly observed in mixed magnetic assemblages such wasp waisted hysteresis loops have been widely documented particularly with respect to rare earth permanent magnets basaltic lava flows remagnetized paleozoic carbonate rocks and an increasingly wide range of other rocks our modelling combined with a review of previous work indicates that there are several conditions that give rise to as well as magnetic properties that are characteristic of wasp waisted hysteresis loops first at least two magnetic components with strongly contrasting coercivities must coexist this condition can arise from either mixtures of grain sizes of a single magnetic mineral or a combination of magnetic minerals with contrasting cocrcivities or a combination of these two situations second materials that give rise to wasp waisted hysteresis loops will have relatively high ratios of the coercivity of remanence to coercive force b cr b c because b0 is controlled by the soft low coercivity component whereas bcris controlled by the hard high coercivity component third values of b cr b c 10 usually only occur for strongly wasp waisted loops when the low coercivity component comprises an overwhelmingly large fraction of the total volume of magnetic grains fourth a given mixture of superparamagnetic and single domain sd grains is more likely to give rise to wasp waisted hysteresis loops than an equivalent mixture of sd and multidomain grains fifth our results provide empirical confirmation that the total magnetization of a material is the sum of the weighted contributions of each component in the absence of significant magnetic interaction between particles thus to contribute significantly to wasp waisted behavior a mineral magnetic component must give rise to a significant portion of the total magnetization of the rock as a result minerals with weak magnetic moments such as hematite need to occur in large concentrations to cause wasp waistedness in materials that also contain ferrimagnetic minerals we outline a method for determining the magnetic components that can give rise to wasp waisted hysteresis loops this method is based on high and low temperature magnetic measurements that are used to identify the dominant remanence bearing mineral s and on mineral magnetic techniques that are used to discriminate between different magnetic domain states the method is illustrated with several examples from archaeological geological and synthetic materials

estimating the return times of great himalayan earthquakes in eastern nepal evidence from the patu and bardibas strands of the main frontal thrust the return times of large himalayan earthquakes are poorly constrained despite historical devastation of cities along the mountain range definitive links between events and specific segments of the main frontal thrust mft are not established and paleoseismological records have not documented the occurrence of several similar events at the same location in east central nepal however recently discovered primary surface ruptures of that megathrust in the a d 1255 and 1934 earthquakes are associated with flights of tectonically uplifted terraces we present here a refined longer slip history of the mft s two overlapping strands patu and bardibas thrusts in that region based on updated geomorphic neotectonic mapping of active faulting two 1 3 km long shallow seismic profiles and logging of two river cut cliffs three paleoseismological trenches and several pits with constraints from 74 detrital charcoals and 14 cosmogenic nuclide ages the amount of hanging wall uplift on the patu thrust since 3650 450 years requires three more events than the two aforementioned the uplift rate 8 5 1 5mm yr thrust dip 25 5 n and apparent characteristic behavior imply 12 17 5m of slip per event on the bardibas thrust discrete pulses of colluvial deposition resulting from the coseismic growth of a flexural fold scarp suggest the occurrence of six or seven paleo earthquakes in the last 4500 50 years the coeval rupture of both strands during great himalayan earthquakes implies that in eastern nepal the late holocene return times of such earthquakes probably ranged between 750 140 and 870 350 years

mechanical analysis of controls on strain partitioning in the himalayas of central nepal we present a mechanical analysis of the problem of slip partitioning between the major thrust systems in a collisional range we focus on two structures in the himalayas of central nepal the main himalayan thrust mht and the main central thrust mct we use finite element modeling to test the influence of various parameters such as friction coefficients and surface processes and we investigate how they affect the distribution of deformation between these two faults we observe that reproduction of the late quaternary kinematic pattern across the range with our model requires strict conditions on the friction coefficients such that the mht is very weak whereas the mct is significantly stronger the most important parameter that controls slip partitioning appears to be the dip angle of the mct with a gentler or steeper mct promoting or inhibiting slip respectively we also show that transient loading and unloading through focused glacial erosion in the higher part of the range can unclamp the mct and allow a significant increase in slip rates the results of this mechanical sensitivity investigation have important implications for the dynamics of the himalayan wedge and point toward along strike structural variations as a first order control on slip partitioning

fluvial incision and tectonic uplift across the himalayas of central nepal the pattern of fluvial incision across the himalayas of central nepal is estimated from the distribution of holocene and pleistocene terraces and from the geometry of modern channels along major rivers draining across the range the terraces provide good constraints on incision rates across the himalayan frontal folds sub himalaya or siwaliks hills where rivers are forced to cut down into rising anticlines and have abandoned numerous strath terraces farther north and upstream in the lesser himalaya prominent fill terraces were deposited probably during the late pleistocene and were subsequently incised the amount of bedrock incision beneath the fill deposits is generally small suggesting a slow rate of fluvial incision in the lesser himalaya the terrace record is lost in the high range where the rivers are cutting steep gorges to complement the terrace study fluvial incision was also estimated from the modern channel geometries using an estimate of the shear stress exerted by the flowing water at the bottom of the channel as a proxy for river incision rate this approach allows quantification of the effect of variations in channel slope width and discharge on the incision rate of a river the determination of incision rates requires an additional lithological calibration the two approaches are shown to yield consistent results when applied to the same reach or if incision profiles along nearby parallel reaches are compared in the sub himalaya river incision is rapid with values up to 10 15 mm yr it does not exceed a few millimeters per year in the lesser himalaya and rises abruptly at the front of the high range to reach values of 4 8 mm yr within a 50 km wide zone that coincides with the position of the highest himalayan peaks sediment yield derived from the measurement of suspended load in himalayan rivers suggests that fluvial incision drives hillslope denudation of the landscape at the scale of the whole range the observed pattern of erosion is found to closely mimic uplift as predicted by a mechanical model taking into account erosion and slip along the flat ramp flat geometry of the main himalayan thrust fault the morphology of the range reflects a dynamic equilibrium between present day tectonics and surface processes the sharp relief together with the high uplift rates in the higher himalaya reflects thrusting over the midcrustal ramp rather than the isostatic response to reincision of the tibetan plateau driven by late cenozoic climate change or late miocene reactivation of the main central thrust

a model for oxygen isotope transport in hydrothermal systems as groundwater passes through porous rocks exchange of oxygen between the fluid and the solid matrix causes a change in the oxygen isotope concentrations in both water and rock if the rate at which the exchange takes place can be estimated as a function of the isotope concentrations and temperature then the time taken for a rock water system to come to equilibrium with respect to isotope concentration might be calculated in this paper conservation laws are used to derive equations for the isotope transport in a porous medium these are combined with fluid flow equations and a simple equation describing the rate of oxygen isotope exchange between the water and the rock the resulting set of equations is solved analytically for the zero dimensional lumped parameter model and numerically using typical parameters for the one dimensional case

formation of gypsum and bassanite by cation exchange reactions in the absence of free liquid h2o implications for mars 1 smectites and hydrated mg sulfate minerals have been identified in close association at various locations on the martian surface the hydration states of sulfates and smectites are dependent on temperature and relative humidity rh and therefore these minerals have the potential to affect cycling and bioavailability of h2o on mars we have conducted x ray powder diffraction experiments to investigate cycling of h2o within mixtures of ca bearing smectites and hydrated mg sulfate minerals under conditions of varying rh similar to those that exist at or just beneath the martian surface our experiments show that under conditions of varying rh cation exchange reactions occur between these two potential components of the martian regolith producing gypsum caso4 2h2o and bassanite caso4 0 5h2o in the absence of free liquid h2o cation exchange reactions were accompanied by significant loss of porosity warping of the sample surface and in some cases volume expansion the formation of ca sulfate minerals in these experiments provides evidence for the development of thin films of h2o at mineral surfaces and suggests that similar processes may operate at the arid surface of mars humidity driven cation exchange reactions between smectites and hydrated mg sulfate minerals may therefore play a role in shaping the present day martian surface and could have provided a transient source of h2o and nutrients e g major and trace elements and possibly organic micro macronutrients for putative microorganisms

gps measurements of present day convergence across the nepal himalaya the high elevations of the himalaya and tibet result from the continuing collision between india and asia which started more than 60 million years ago1 4 from geological and seismic studies of the slip rate of faults in asia5 it is believed that approximately one third of the present day convergence rate between india and asia 58 4mmyr 1 is responsible for the shortening uplift and moderate seismicity of the himalaya great earthquakes also occur infrequently in this region releasing in minutes the elastic strain accumulated near the boundary zone over several centuries and accounting for most of the advance of the himalaya over the plains of india the recurrence time for these great earthquakes is determined by the rate of slip of india beneath tibet which has hitherto been estimated indirectly from global plate motions6 from the slip rates of faults in asia7 8 from seismic productivity9 and from the advance of sediments on the northern ganges plain10 here we report geodetic measurements using the global positioning system gps of the rate of contraction across the himalaya which we find to be 17 52 2 mm yr 1 from the form of the deformation field we estimate the rate of slip of india beneath tibet to be 20 5 2 mmyr 1 strain sufficient to drive one or more great himalayan earthquakes with slip similar to that accompanying the magnitude 8 1 bihar nepal 1934 earthquake may currently be available in western nepal

lower edge of locked main himalayan thrust unzipped by the 2015 gorkha earthquake faults are unlocked by earthquakes analysis of seismic data from the 2015 nepal earthquake shows that only part of the main himalayan thrust fault was unzipped by the quake leaving much of the fault locked and ready to slip in a future event

estimating rates of denudation using cosmogenic isotope abundances in sediment we propose as a testable hypothesis a basin scale approach for interpreting the abundance of in situ produced cosmogenic isotopes an approach which considers explicitly both the isotope and sediment flux through a drainage basin unlike most existing models which are appropriate for evaluating in situ produced cosmogenic isotope abundance at discrete points on earth s surface our model is designed for interpreting isotope abundance in sediment because sediment is a mixture of materials in favourable cases derived from throughout a drainage basin we suggest that measured isotope abundances may reflect spatially averaged rates of erosion we investigate the assumptions and behaviour of our model and conclude that it could provide geomorphologists with a relatively simple means by which to constrain the rate of landscape evolution if a basin is in isotopic steady state and if sampled sediments are well mixed

medieval pulse of great earthquakes in the central himalaya viewing past activities on the frontal thrust the himalaya has experienced three great earthquakes during the last century1934 nepal bihar 1950 upper assam and arguably the 1905 kangra focus here is on the central himalayan segment between the 1905 and the 1934 ruptures where previous studies have identified a great earthquake between thirteenth and sixteenth centuries historical data suggest damaging earthquakes in a d 1255 1344 1505 1803 and 1833 although their sources and magnitudes remain debated we present new evidence for a great earthquake from a trench across the base of a 13m high scarp near ramnagar at the himalayan frontal thrust the section exposed four south verging fault strands and a backthrust offsetting a broad spectrum of lithounits including colluvial deposits age data suggest that the last great earthquake in the central himalaya most likely occurred between a d 1259 and 1433 while evidence for this rupture is unmistakable the stratigraphic clues imply an earlier event which can most tentatively be placed between a d 1050 and 1250 the postulated existence of this earlier event however requires further validation if the two earthquake scenario is realistic then the successive ruptures may have occurred in close intervals and were sourced on adjacent segments that overlapped at the trench site rupture s identified in the trench closely correlate with two damaging earthquakes of 1255 and 1344 reported from nepal the present study suggests that the frontal thrust in central himalaya may have remained seismically inactive during the last similar to 700years considering this long elapsed time a great earthquake may be due in the region

how perfect is the himalayan arc the himalayan plate boundary because it is entirely subaerial is both the most dramatic and the most accessible to direct observation of all active convergent boundaries on earth the shape of this boundary can be described as a small circle of radius 1696 55 km centered at long 91 6 1 6 e and lat 42 4 2 1 n for the extent of the arc between long 77 2 and 92 1 e the pole of this small circle is consistent whether seismicity topography or stress state is used to define the position of the tectonic boundary the defined small circle also coincides with a peak in microseismicity the maximum horizontal strain rate and a peak in the vertical velocity field this quantitative definition of a stable curved tectonic boundary is a prerequisite to modeling the dynamics of curvature in convergent arcs and applying appropriate boundary conditions to other regional models

plio plistocene in sequence thrust propagation along the main central thrust zone kumaon garhwal himalaya india new thermochronological data abstract the kumaon and garhwal himalayas nw india have similar geologic tectonic topographic and precipitation pattern however low temperature thermochronological age pattern varies significantly along the strike as well as across the strike of major faults here we interpret the new set of thermochronological age data across the strike of major faults namely vaikrita thrust vt munsiari thrust mt berinag thrust bt and the baijnath nappe along pindari valley in the kumaon garhwal himalaya in the present work 18 new apatite fission track aft ages of samples collected along a north south transect have been reported ages in the hanging wall of vt which range from 2 1 0 2 to 4 2 0 7 ma have been found to be becoming younger linearly with distance from north to the vt this trend crosses the vt with significant jump in ages in the south of the vt ages are lying between 0 3 0 1 and 3 9 0 8 ma and show linear variation with distance from the vt to bt no significant jump in ages across the mt is observed and the linear trend of younging southward is continuing till the bt the break and displacement in age pattern with youngest ages close to thrusts possibly indicate progressive late thrust movement it reflects sequential uplift and cooling towards the south from the vt to bt which is consistent with an in sequence style of thrust propagation we compare our new data to the published thermochronological data from its adjacent traverses along the goriganga and dhauliganga valleys 25 km away from the present traverse towards east and west respectively in the kumaon garhwal himalaya in order to understand potential along strike variations in the kinematics our observations indicate a dynamic coupling between tectonic and surface processes in the garhwal kumaon himalaya and support that the geometry of crustal scale faults and their associated kinematics control exhumation pathways of rocks

effects of inherited cores and magmatic overgrowths on zircon u th he ages and age eu trends from greater himalayan sequence rocks mount everest region tibet previous constraints on the timing and rate of exhumation of the footwall of the south tibetan detachment system stds north of mount everest suggest rapid miocene cooling from 700 c to 120 c between 14 17 ma however 25 new single grain zircon he ages from leucogranites intruding greater himalayan sequence rocks in the footwall of the stds are between 9 9 and 15 ma with weighted mean ages between 10 and 12 ma zircon grains exhibit a positive correlation between age and effective uranium eu laser ablation zircon u pb geochronology detailed sem observations and laser ablation depth profiling of these zircons reveal low eu 0 5 2 5 ga inherited cores overgrown by high eu 17 22 ma rims this intragranular zonation produces ages as much as 32 too young when a standard alpha ejection correction assuming uniform eu distribution is applied modeling of the effects of varying rim thickness and rim eu concentration on the bulk grain eu and alpha ejection correction suggests that zonation also exerts the primary control on the form of the age eu correlation observed application of grain specific zonation dependent age corrections to our data yields zircon he ages between 14 and 17 ma in agreement with aft and 40ar 39ar ages growth of magmatic rims followed by cooling to 120 c within 1 6 million years supports rapid tectonic exhumation associated with slip along the stds in the miocene this study highlights the importance of characterizing parent nuclide zonation in zircon he studies which seek to understand the timing of exhumation along exhumed crustal blocks

scaling factors for production rates of in situ produced cosmogenic nuclides a critical reevaluation abstract new scaling factors are presented describing the altitude and latitudinal dependence of production rates for in situ produced cosmogenic nuclides the new factors incorporate the influence of the non dipole contributions to the geomagnetic field on the cosmic ray flux the currently used scaling factors of lal earth planet sci lett 104 1991 424 439 are based on the assumption of a dipole field the overall strategy used in deriving the new scaling factors is however very similar to that of lal and relies on the same and or similar nuclear data in this reevaluation data are only used if the effective geomagnetic parameters inclination horizontal field strength can be reconstructed for the time of measurement the absorption free pathlengths for cosmic rays selected for this study are based on observational data at altitudes relevant for exposure age dating sea level to 7000 m at sea level and latitudes between 20 and 40 the new factors are up to 18 lower and at high altitudes more than 30 higher than those of lal in addition to accounting for the influence of the effective geomagnetic field and providing more applicable estimates of the new factors also allow for correction for the significant deviations from the standard pressure altitude relationship that exist in the atmosphere 10

some analytical methods for converting thermochronometric age to erosion rate 1 an analytical method is presented for converting thermochronometric ages to surface erosion or equivalently exhumation rate the method incorporates the two most important thermal processes during cooling by erosion the dependence of closure temperature on cooling rate and the advection of heat by rock motion toward the earth s surface two thermal models are considered 1 a steady state model valid for low erosion rates and 2 an eroding half space model which has no steady state but captures the transient increase of geothermal gradient with erosion in each case it is assumed that data consist of one or more thermochronometric ages present day surface geothermal gradient and topographic information including the elevation at which the age was obtained analytical solutions are provided to derive the erosion rate from these data either as an explicit expression for the steady case or as a root finding problem for the transient case a graphical method for plotting age against erosion rate and geothermal gradient is presented as a method for solving the root finding problem and for tracking analytical errors in observations of age and surface geothermal gradient the graphical method is also appropriate for comparing data from different elevations or from different thermochronometric systems examples are provided using synthetic data or published data from the literature

calculation of the cosmogenic nuclide production topographic shielding scaling factor for large areas using dems the recent surge of applications using terrestrial cosmogenic nuclides tcns to calculate catchment averaged erosion rates from isotopic concentrations in fluvial sediment and the prospect of coupling tcn production functions with numerical surface process models spms necessitate a fast and accurate algorithm for the calculation of topographic shielding topographic shielding refers to the proportion of the incoming cosmic radiation that is shielded by the surrounding topography the scaling factor being defined as the ratio of the unshielded total minus shielded to the total or maximum cosmic ray flux i e the flux received by a horizontal unobstructed surface topography contributes to the reduction of tcn production by obstructing a certain proportion of the incoming flux and by modifying the angle of incidence available algorithms calculate the proportion of obstructed radiation by dividing the horizon as seen by the sample a grid cell in the case of a dem into arc segments usually of equal length for which the average obstruction heights expressed as zenith angles are calculated the use of these methods is feasible only when dealing with a small number of isolated samples since the identification of obstructions when dealing with an entire area is computationally very intensive this paper describes a method that uses a relief shadow modelling technique to identify those areas of a dem that are under shadow i e shielded and thus to account for the obstructed radiation this method produces results that are very similar to those obtained using a direct implementation of available methods maximum difference between results of c 0 1 the method based on relief shadow modelling is also faster than a direct implementation of any available method and can be readily implemented in any gis system with raster capabilities copyright 2006 john wiley sons ltd

interpretation of u th he single grain ages from slowly cooled crustal terranes a case study from the transantarctic mountains of southern victoria land low temperature thermochronologic techniques e g apatite fission track aft thermochronology and u th he dating constrain near surface t t paths and are often applied to uplift denudation and landscape evolution studies samples collected in vertical profiles from granitic walls on either side of the ferrar glacier southern victoria land antarctica were analyzed using aft thermochronology and apatite u th he dating to further constrain the lowest temperature thermal history of this portion of the transantarctic mountains aft central ages vary systematically with elevation and together with track length information define a multi stage cooling denudation history in the cretaceous and early tertiary apatite u th he single grain age variation with elevation is not as systematic with considerable intra sample age variation although many complicating factors e g u and th rich micro inclusions fluid inclusions variation in crystal size particle ejection correction zonation and particle ejection correction implantation of he into a crystal or impediment of he diffusion out of a crystal and 147sm derived particles may contribute to age dispersion we found that variation in single grain ages correlated with cooling rate samples that cooled relatively quickly have less variation in single grain ages whereas samples that cooled relatively slowly 3 c m y or resided within an u th he partial retention zone heprz prior to more rapid cooling have a comparatively greater variation in ages decay of u and th via particle emission creates a 4he concentration profile dependent upon the initial parent u th within a crystal variation of single grain ages for samples with non homogeneous u th distributions will be enhanced with long residence time in the partial retention zone i e slow cooling because of the relative importance of loss via volume diffusion and loss via particle ejection with respect to the u th zonation and the grain boundary correction of ages for particle ejection ft correction factor typically assumes uniform u and th distribution within the crystal and when applied to a population of crystals with different u and th distributions will enhance the variation in ages most complicating factors listed above for apatite u th he ages result in ages that are too old we propose that if considerable variation in u th he single grain ages exists that a weighted mean age is determined once outlier single crystal ages are excluded using the criterion of chauvenet or a similar approach we suggest that the true age or most representative age for that age population lies between the minimum u th he age and the weighted mean age we apply this approach coupled with composite age profiles to better constrain the t t history of the profiles along the ferrar glacier significant intra sample variation in single crystal apatite u th he ages and other minerals dated by the u th he method should be expected especially when the cooling rate is slow the variation of u th he single crystal ages is therefore another parameter that can be used to constrain low temperature thermal histories

interseismic coupling on the main himalayan thrust we determine the slip rate and pattern of interseismic coupling on the main himalayan thrust along the entire himalayan arc based on a compilation of geodetic interferometric synthetic aperture radar and microseismicity data we show that convergence is perpendicular to the arc and increases eastwards from 13 3 1 7 mm yr to 21 2 2 0 mm yr these rates are comparable to geological and geomorphic estimates indicating an essentially elastic geodetic surface strain the interseismic uplift rate predicted from the coupling model closely mimics the topography suggesting that a small percentage of the interseismic strain is permanent we find that the fault is fully locked along its complete length over about 100 km width we don t find any resolvable aseismic barrier that could affect the seismic segmentation of the arc and limit the along strike propagation of seismic ruptures the moment deficit builds up at a rate of 15 1 1 10 19 n m yr for the entire length of the himalaya