climatic control on rapid exhumation along the southern himalayan front abstract along the southern himalayan front shf areas with concentrated precipitation coincide with rapid exhumation as indicated by young mineral cooling ages twenty new young 80 of the annual precipitation aft cooling ages reveal a coincidence between rapid erosion and exhumation that is focused in a 50 70 km wide sector of the himalaya rather than encompassing the entire orogen assuming simplified constant exhumation rates the rocks of two age vs elevation transects were exhumed at 1 4 0 2 and 1 1 0 4 mm a with an average cooling rate of 40 50 c ma during pliocene quarternary time following other recently published hypotheses regarding the relation between tectonics and climate in the himalaya we suggest that this concentrated loss of material was accommodated by motion along a back stepping thrust to the south and a normal fault zone to the north as part of an extruding wedge climatically controlled erosional processes focus on this wedge and suggest that climatically controlled surface processes determine tectonic deformation in the internal part of the himalaya

ten million year history of a thrust sheet the final deformed state of a fold and thrust belt may be reached by an infinite number of kinematic paths two end member kinematic paths are due to continuous or discontinuous rates of deformation we have used a new magnetostratigraphic section from the spanish pyrenees to calibrate the emplacement history over 10 m y of a major thrust sheet the sierras marginales thrust sheet and the deformation of both its hanging wall and footwall six time windows from before 36 5 ma until after 24 7 ma were recognized on the basis of structural and stratigraphic relationships between syntectonic strata and major structures footwall deformation of the sierras marginales thrust sheet occurred continuously on a detachment within the foreland basin sequence 4 5 km net shortening shortening rates steadily increased from 0 14 to 1 5 mm yr with time although the detachments at the base of the foreland and the detachment at the base of the sierras marginales thrust sheet were active coevally the rate of displacement on the sierras marginales detachment decayed with time and shows considerable variability emplacement of the thrust sheet may be divided into three distinct periods a rapid 13 8 km translation from 37 0 to 36 5 ma 27 6 mm yr a gradual climb of the toe of the thrust sheet up a 4 3 km long ramp across the foreland basin succession from 36 5 to 32 0 ma 0 95 mm yr and a final 8 9 km translation from 32 0 to 29 5 ma 3 56 mm yr internal deformation of the thrust sheet occurred only after it reached its present position at some time before 29 5 ma shortening rates steadily decreased from 0 6 mm yr between 29 5 and 27 8 ma to 0 26 mm yr from 27 8 until after 24 7 ma only 1 km of shortening accumulated during each of the two periods the last shortening localized on the most northerly thrust in the study area folding and subsequent faulting above the detachment beneath the foreland suggest that slip was transmitted to its tip point continuously throughout the deformation in contrast translation followed by internal deformation on the hinterland side of the toe of the sierras marginales thrust sheet indicates a successive deactivation of the southern parts of the detachment with time rate of deformation on both detachments was discontinuous and shows substantial variability about the mean in general the spatial and temporal pattern of deformation was distributed and continuous in detail however shortening was spatially and temporally discontinuous above each detachment and structures related to each one display distinctly different deformational patterns rates and styles

numerical simulations of glacial valley longitudinal profile evolution glaciers shape alpine landscapes they broaden valley bottoms enhance local valley relief generate multiple steps overdeepen valley floors and cause tributary valleys to hang these distinctive glacial signatures result from 10 4 10 5 yr of erosion during which swings in climate drive advances and retreats of alpine glaciers we use a numerical model of glacial erosion to explore the development of the longitudinal profiles of glaciated valleys the model is driven by the past 400 k y of variable climate because both sliding speed which dictates abrasion rate and water pressure fluctuations which strongly modulate quarrying rate should peak at the equilibrium line altitude ela we expect the locus of most rapid erosion to follow the transient ela simulations of a single glacial valley show rapid flattening of the longitudinal profile inclusion of a tributary glacier creates a step immediately downvalley of the tributary junction that persists over multiple glaciations and commonly leaves the tributary valley hanging steps and overdeepenings result from an increase in ice discharge immediately below the tributary junction which is accommodated primarily by increased ice thickness and hence sliding rate the size of the step increases with the ratio of tributary to trunk ice discharge while the height of a hanging valley reflects the difference in the time integrated ice discharge in tributary and trunk valleys and therefore increases as the discharge ratio decreases

the north pyrenean aquitaine basin france evolution and hydrocarbons the mesozoic tertiary aquitaine basin overlaps the aquitaine block and the northern edge of the iberian margin both units are situated in the southwesternmost part of the european continent the aquitaine shelf is a stable platform overlying a relatively thin crustal segment it underwent extensional block faulting and many of its folded structures are related to salt tectonism the iberian block is a thicker lithospheric unit that acted as a buttress at the junction of these two structural domains the south aquitaine or north pyrenean area developed where crustal stretching rifting wrenching and orogenic compression were maximal from the mesozoic to the tertiary the history of the basin has been a suite of rifting attempts in the context of the progressive opening of the bay of biscay which never completely succeeded the relative plate motions and the later convergence led nevertheless to the formation of the aquitaine basin and to the emergence of the pyrenean fold belt the mesozoic tertiary sedimentary infill of the basin is mostly marine with thick evaporites stable platform carbonates subsiding platform shaly calcareous deposits and a characteristic diachronous turbiditic flysch and molasse complex events during the basin forming extensional and translational periods contributed more to the generation of the hydrocarbon accumulations than did the effects of the subsequent compressive structural regime the latter however may have enhanced thermal flow and thus maturation of hydrocarbon source rocks

river incision into bedrock mechanics and relative efficacy of plucking abrasion and cavitation improved formulation of bedrock erosion laws requires knowledge of the actual processes operative at the bed we present qualitative field evidence from a wide range of settings that the relative efficacy of the various processes of fluvial erosion e g plucking abrasion cavitation solution is a strong function of substrate lithology and that joint spacing fractures and bedding planes exert the most direct control the relative importance of the various processes and the nature of the interplay between them are inferred from detailed observations of the morphology of erosional forms on channel bed and banks and their spatial distributions we find that plucking dominates wherever rocks are well jointed on a submeter scale hydraulic wedging of small clasts into cracks bashing and abrasion by bedload and chemical and physical weathering all contribute to the loosening and removal of joint blocks in more massive rocks abrasion by suspended sand appears to be rate limiting in the systems studied here concentration of erosion on downstream sides of obstacles and tight coupling between fluid flow patterns and fine scale morphology of erosion forms testify to the importance of abrasion by suspended load rather than bedload particles mechanical analyses indicate that erosion by suspended load abrasion is considerably more nonlinear in shear stress than erosion by plucking in addition a new analysis indicates that cavitation is more likely to occur in natural systems than previously argued cavitation must be considered a viable process in many actively incising bedrock channels and may contribute to the fluting and potholing of massive unjointed rocks that is otherwise attributed to suspended load abrasion direct field evidence of cavitation erosion is however lacking in terms of the well known shear stress or stream power erosion law erosion by plucking is consistent with a slope exponent n of 2 3 to 1 whereas erosion by suspended load abrasion is more consistent with a slope exponent of 5 3 given that substrate lithology appears to dictate the dominant erosion process this finding has important implications for long term landscape evolution and the models used to study it

orogeny and orography the effects of erosion on the structure of mountain belts a numerical model of the coupled processes of tectonic deformation and surface erosion in convergent orogens is developed to investigate the nature of the interaction between these processes crustal deformation is calculated by a two dimensional finite element model of deformation in response to subduction and accretion of continental crust erosion operates on the uplifted surface of this model through fluvial incision which is taken to be proportional to stream power the relative importance of the tectonic and erosion processes is given by a dimensionless erosion number relating convergence velocity rock erodibility and precipitation rate this number determines the time required for a system to reach steady state and the final topographic shape and size of a mountain belt fundamental characteristics of the model orogens include asymmetric topography with shallower slopes facing the subducting plate and an asymmetric pattern of exhumation with the deepest levels of exhumation opposite to subduction these characteristics are modified when the regional climate exhibits a dominant wind direction and orographically enhanced precipitation on one side of the mountain belt the two possible cases are dominant wind in the direction of motion of the subducting plate and dominant wind direction in the opposite direction of the subducting plate velocity models of the former case predict a broad zone of exhumation with maximum exhumation in the orogen interior models of the latter case predict a focused zone of exhumation at the margin of the orogen and at high erosion number a reversal in the topographic asymmetry natural examples of these two cases are presented the southern alps of new zealand exhibits the climate and exhumation asymmetry characteristic of wind in the direction opposite to motion of the subducting plate the asymmetry of topography suggests that erosion is not efficient enough to have reversed the topographic asymmetry the contrasting example of dominant wind in the direction of subduction motion is provided by the olympic mountains of washington state in this case exhumation of deep levels of the cascadia accretionary wedge shows a broad domal pattern consistent with the observed orographic precipitation

thermochronology constraints for the propagation sequence of the south pyrenean basement thrust system france spain 1 in this work we combined apatite fission track and biotite k feldspar 40ar 39ar ages with tectonic data in the west central part of the axial zone of the pyrenees we discuss the exhumation ages and rates of the neouvielle bordere louron and bielsa variscan granites and their relationships with the timing and sequence of south vergent basement thrusting within the pyrenean orogenic prism the 40ar 39ar ages on k feldspars from the neouvielle massif sample nv7 seem to indicate tectonic movements on the eaux chaudes thrust during the early middle eocene fission track results suggest that the exhumation of the neouvielle massif occurred around 35 ma and exhumation of the bordere louron massif around 32 ma in relation to thrusting on the gavarnie thrust the bielsa massif was exhumed from around 19 ma by out of sequence movements on the bielsa thrust we thus show that whereas most of the pyrenean basement thrust faults here the eaux chaudes gavarnie and guarga thrusts were active in sequence toward the southern foreland from the early eocene to the earliest miocene some of them here the bielsa thrust were activated out of sequence in the hinterland later than the generally accepted aquitanian age for the end of the pyrenean compression finally the apatite fission track modeling indicate a last cooling episode starting around 5 ma which is most certainly related to the pliocene reexcavation of the southern and northern flanks of the pyrenees

implications of the shear stress river incision model for the timescale of postorogenic decay of topography 1 the reason for the survival of mountainous topography in ancient orogenic belts is a long standing problem in geomorphology and geodynamics we explore the geomorphologic controls on the timescale for the postorogenic decay of topography to address the question of whether there is a viable geomorphologic explanation for the persistence of topography in ancient orogenic belts or whether alternative geodynamic explanations must be sought using both approximate analytical solutions and numerical simulations we show that the standard detachment limited stream power river incision model predicts for reasonable initial topographies relatively short postorogenic decay times of 1 10 myr additional layers of complexity are introduced to this simplest model including isostasy and flexural strength a transition to transport limited conditions during decline and the incorporation of the combined effects of a critical threshold for erosion and the stochastic variability of flood magnitudes each of these additional factors acts to lengthen the decay timescale pure airy isostatic rebound of a thick crustal root results in decay times that are at most a factor of 6 longer than the detachment limited result the transition to transport limited conditions involves partial protection of the bed by a thin layer of alluvium which inhibits erosion and therefore increases decay times by a factor of 2 3 in our analysis finally including critical shear stress results increases decay time by approximately a factor of 20 for the parameter values in our example calculation more importantly however a significant portion of topography remains after predicted lowering rates have dropped to values less than measured rates of denudation in weathering limited low relief environments thus a model combining isostasy a transition to transport limited conditions and a critical shear stress for erosion could account for the presence of residual topography for hundreds of million years all three of these factors can be expected to play a role in natural settings

statistical modelling of annealing kinetics of fission tracks in zircon reassessment of laboratory experiments abstract we remodel the annealing kinetics of fission tracks fts in zircon to fit data from 43 laboratory experiments in three earlier publications using spontaneous ft lengths in zircon from the nisatai dacite heating temperature and time are in the ranges of 350 912 c and 10 3 104 h respectively the data sets illustrate a transition in annealing process at a certain temperature zone in order to express the annealing behavior involving the possible transitional process with improved goodness of fit two types of models are constructed the hybrid linear model that consists of the fanning linear model at high temperature and the parallel linear model at low temperature connected with a transitional temperature zone and the parallel curvilinear model that approximates smoothly the transition in the annealing process the validity of the extrapolation of the model to geological time scale is confirmed by comparison with some deep borehole core data the length estimates by the hybrid linear model agree well with the borehole data whilst the estimates by the parallel curvilinear model show a significant disagreement with highly annealed borehole samples although this model provides good agreement with the experimental data for the hybrid linear model estimates of the partial annealing zone corresponding to isothermal heating for 106 and 107 years are 281 352 c and 262 330 c where the top and bottom of the partial annealing zone are defined by the mean length reduction ratio of 0 4 and 0 8 respectively

effect of alpha damage annealing on apatite u th he thermochronology recent data on he diffusion challenge the temperature sensitivity of apatite u th he thermochronology the damage induced by recoil of u and th decay series during emission of particles recoil damage has been proposed to modify he diffusion properties through time however we propose that annealing of these irradiation defects may be an important phenomenon and may be significant in case of slowly cooled or reheated basement rocks to test this hypothesis we developed a quantitative model including an explicit treatment of recoil damage annealing and their effect on he diffusion kinetics and calibrate it against literature data our model is based on two hypotheses 1 helium is in equilibrium between an apatite crystal and its defects and 2 alpha recoil damage annealing can be described analogously to fission track annealing this model has been embedded into a monte carlo simulation of helium production ejection diffusion and applied to data from the french massif central a complex slowly cooled terrain with burial reheating where the thermal history has been constrained by previous fission track ft data including ft length distributions u th he ages are close to the ft ages from the same samples and are generally reproducible among replicates but some samples present he age dispersion that is not correlated with crystal size our model reproduces the massif central data very well except for three samples where he ages are older than corresponding ft ages we show that annealing of irradiation damage has an important impact on retentivity of helium and that the he content he is only a rough approximation of the damage level in particular our results show that independence of he ages on crystal sizes in case of reheated samples is a clear indication of the higher he retentivity induced by recoil defects and that an explicit treatment of defect annealing is required for a correct interpretation of u th he ages in such a case more generally a correlation with the crystal size can bring information on the thermal path only if the age of defects well represented by the fission track age is available due to the dependence of the partial retention zone on damages conversely in case of rapid cooling or for samples having low u and th contents damage effects can be ignored without significant effects on he ages

lithological control on the post orogenic topography and erosion history of the pyrenees abstract numerous studies on active mountain ranges have demonstrated the interaction between tectonics and climate in shaping topography here we explore how variations in rock types have affected the topographic development of the pyrenees since cessation of orogenesis ca 20 ma our study is based on a multidisciplinary approach and integrates topographic analyses rock strength measurements and thermal modelling of low temperature thermochronological data published across the central pyrenees results indicate a strong influence of rock strength in determining the post orogenic morphology of the pyrenees we observe a correlation between rock strength and the normalized channel steepness index k s n of the different lithologies moreover the highest topography is dominated by the variscan plutonic massifs which have highest rock strength consequently the drainage divide appears to track the position of these massifs abrupt deceleration of exhumation recorded in inverse modelling of low temperature thermochronologic data suggests that the exhumation of the variscan massifs also played role in lowering in erosion rates over the massifs during orogenesis

the sis palaeovalley a record of proximal fluvial sedimentation and drainage basin development in response to pyrenean mountain building the sis conglomerate body represents the middle eocene to oligocene transfer zone trunk palaeovalley fill of the sis fluvial system a drainage system established within the pyrenees during late palaeocene times the spatial stability of the fluvial transfer zone active for at least 38 my and hence the longevity of its aggradational palaeovalley component 19 5 my was controlled by its location between long lived pre existing structures coarse grained fluvial facies dominate the palaeovalley fill with alluvial fan facies shed from its defining marginal structures the detailed sedimentology of very proximal fluvial facies deposited within the dominantly erosional realm of an active mountain belt has rarely been documented before because of their poor preservation potential the sis conglomerate body contains a robust internal stratigraphy with stratigraphic units defined by distinct bounding surfaces across which there are pronounced changes in facies and provenance these mark the reorganization of the headward portions of the sis fluvial system during the evolution of the pyrenean axial zone antiformal stack major changes in discharge resulted demonstrating the highly variable nature of even mountain belt scale fluvial systems when viewed on timescales of several to tens of millions of years provenance details indicate that initial unroofing of hercynian granitoids situated within the pyrenean axial zone occurred around 54 5 ma early ypresian immediately before the first significant exhumation event within the drainage basin of the sis fluvial system this is earlier than previously constrained by apatite fission track studies rock uplift accelerated in the lutetian and bartonian with the initial aggradation of the palaeovalley fill the cajigar and cornudella formations and sis one and two members this became marked in the priabonian sis three and four members with significant activity on local structures including the morreres backthrust an increase in basement derived clasts and a headwater decapitation event also indicate pronounced axial zone antiformal stack development at this time axial zone development intensified further in the oligocene with the deposition of the collegats formation and the switch in the main depositional loci of the system from the tremp graus thrust sheet top basin to the ebro basin to the south

strong rocks sustain ancient postorogenic topography in southern africa the cape mountains of southern africa exhibit an alpine like topography in conjunction with some of the lowest denudation rates in the world this presents an exception to the often cited coupling of topography and denudation rates and suggests that steep slopes alone are not sufficient to incite the high denudation rates with which they are commonly associated within the cape mountains slope angles are often in excess of 30 degrees and relief frequently exceeds 1 km yet be 10 based catchment averaged denudation rates vary between 2 32 0 29 m m y and 7 95 0 90 m m y we attribute the maintenance of rugged topography and suppression of denudation rates primarily to the presence of physically robust and chemically inert quartzites that constitute the backbone of the mountains be 10 based bedrock denudation rates on the interfluves of the mountains vary between 1 98 0 23 m m y and 4 61 0 53 m m y the close agreement between the rates of catchment averaged and interfluve denudation indicates topography in steady state these low denudation rates in conjunction with the suggestion of geomorphic stability are in agreement with the low denudation rates 20 m m y estimated for southern africa during the late cenozoic by means of cosmogenic nuclide thermochronology and offshore sedimentation analyses accumulatively these data suggest that the coastal hinterland of the subcontinent may have experienced relative tectonic stability throughout the cenozoic

punctuated thrust deformation in the context of doubly vergent thrust wedges implications for the localization of uplift and exhumation understanding uplift and exhumation in orogenesis requires integrating the effects of both climatic and tectonic forcing punctuated thrust activity is widely documented from field studies in a range of mountain belts the potential for climatic forcing of thrust activity has been a subject of recent debate here we use a modeling approach to analyze the behavior of individual thrust units in the context of asymmetric doubly vergent thrust wedges the model predicts that rates of surface uplift frontal accretion and exhumation should be punctuated on a time scale linked to thrust sheet geometry and convergence rates this time scale ranges from 0 1 to 5 m y for various settings and should be calculated before external forcings such as climate are invoked

an ibexian early ordovician trilobite faunule from the type section of the rabbitkettle formation southern mackenzie mountains northwest territories the rabbitkettle formation is a widespread lower paleozoic unit in the northern cordillera a trilobite faunule collected from the middle part of the type section and described here documents an early ibexian early ordovician age for this interval and provides the first biostratigraphic reference point for the type section this faunule is dominated by the agnostoid micragnostus the leiostegiid leiostegium and the remopleuridid kainella these genera and rarer taxa have north american south american and austral asian affinities a new species kainella eleutherolfi is established

expression of active tectonics in erosional landscapes understanding the manner and degree to which topography in active mountain ranges reflects deformation of the earth s surface remains a first order goal of tectonic geomorphology a substantial body of research in the past decade demonstrates that incising channel systems play a central role in setting relationships among topographic relief differential rock uplift rate and climatically modulated erosional efficiency this review provides an introduction to the analysis and interpretation of channel profiles in erosional mountain ranges we show that existing data support theoretical expectations of positive monotonic relationships between channel steepness index a measure of channel gradient normalized for downstream increases in drainage area and erosion rate at equilibrium and that the transient response to perturbations away from equilibrium engenders specific spatial patterns in channel profiles that can be used to infer aspects of the forcing these aspects of channel behavior lay the foundation for a series of case studies that we use to illustrate how focused quantitative analysis of channel morphology can provide insight into the spatial and temporal dynamics of active deformation although the complexities of river response to climate lithology and uplift patterns mean that multiple interpretations of topographic data alone will always possible we show that application of stream profile analysis can be a powerful reconnaissance tool with which to interrogate the rates and patterns of deformation in active mountain belts

late oligocene early miocene initiation of shortening in the southwestern chinese tian shan implications for neogene shortening rate variations abstract the cenozoic tian shan is one of the preeminent examples of an intracontinental orogen however there remains a significant controversy over when deformation related to the india asia collision commenced and therefore how shortening within the mountains has been partitioned over time one approach has been to look at the modern shortening rate as measured by geodetic studies combined with estimates of the total shortening across the range and extrapolate backwards this approach suggests that the onset of range construction was ca 10 ma k y abdrakhmatov s a aldazhanov b h hager m w hamburger t a herring k b kalabaev k b kalabayev v i makarov p molnar s v panasyuk m t prilepin r e reilinger i s sadybakasov b j souter y a trapeznikov v y tsurkov a v zubovich relatively recent construction of the tien shan inferred from gps measurements of present day crustal deformation rates nature 384 6608 1996 450 453 an alternate method is to determine the age of the onset of exhumation using thermochronology we present 19 new apatite fission track aft results from the southwestern chinese portion of the belt this region represents the first area exhumed during the late tertiary along a transect at ca 76 e exhumation commenced at the oligocene miocene boundary 24 ma along the maidan and muziduke thrusts which bound the southern side of the kokshaal range subsequently deformation propagated ca 20 km south to the kashi basin bounding thrust kbt which was exhumed by no earlier than 18 9 3 3 ma three detrital aft samples from plio pleistocene strata deposited ca 20 km farther south contain fission track grain age peaks that young monotonically upwards from 20 9 7 0 5 3 ma to 15 9 5 4 4 0 ma with a fairly constant lag time of 16 to 18 ma these ages combined with structural data suggest that both the hanging wall and the footwall of the kbt experienced a renewed episode of exhumation during the latest cenozoic the discrepancy between the late oligocene miocene initiation of significant exhumation shown herein and the 10 ma initiation estimate from geodesy suggests that the tian shan has undergone a complex late cenozoic shortening history assuming that the present shortening rate could account for the total amount of cenozoic shortening in 10 ma and realizing that shortening initiated at least 15 myr earlier we conclude that the shortening rate must have varied over time possibly in pulsed southward migrating events and that the present rate may not reflect the average rate since initiation of range uplift

combined use of tracer approach and numerical simulation to estimate groundwater recharge in an alluvial aquifer system a case study of nasunogahara area central japan summary in this study we simulate the spatial and temporal distribution of groundwater recharge in an alluvial aquifer system in the nasunogahara area of japan natural stable isotopes 18 o d were considered as additional calibration targets in a numerical model the reliability of the model outputs was further validated by comparing the results from the numerical simulation and an independent tracer approach the results indicated that the calibrated model can effectively simulate the spatial and temporal characteristics of the contribution ratios of recharge sources to groundwater in the nasunogahara area however the tracer approach i e end member mixing analysis provided more reliable results at point scale particularly for the estimated contribution ratios of paddy field water the precipitation in the nasunogahara area is the major recharge source its mean contribution ratio is 58 for a one year period over the entire alluvial fan river seepage is significant in the upstream area of the alluvial fan and the contribution ratio of river waters along the river channels in the upstream area increases during the wet season paddy field water is a highly important recharge source in the midstream and downstream areas of the alluvial fan and the contribution ratio of paddy field water obviously increases from dry season to wet season because of irrigation this study demonstrates that combined use of the tracer approach and numerical simulation with stable isotopes as additional calibration targets can eliminate their respective limitations and can assist in better understanding the groundwater recharge mechanism in alluvial aquifer systems

syntectonic burial and post tectonic exhumation of the southern pyrenees foreland fold thrust belt in the oligocene to earliest miocene the entire southern pyrenean foreland fold and thrust belt was backfilled and buried in up to 3 km of syntectonic continental conglomerates from the range9s own erosional debris then starting in the mid to late miocene the pyrenees were exhumed by erosional excavation to their present relief a model is proposed for this unusual development whereby late eocene oligocene tectonic uplift of the margins of the ebro basin raised local base level and blocked normal dispersal of erosional debris to adjacent oceans thus causing the basin to fill then backfill northwards across the entire southern flank of the pyrenees subsequent miocene rifting of the catalan mediterranean margin and the messinian salinity crisis lowered base level allowing the ebro river to cut headward capture the ebro basin and re excavate the pyrenees

volumetric budget and grain size fractionation of a geological sediment routing system eocene escanilla formation south central pyrenees the supply of sediment and its characteristic grain size mix are key controls on depositional facies and stratigraphic architectures in sedimentary basins consequently constraints on sediment caliber budgets and fl uxes are a prerequisite for effective stratigraphic prediction here we investigate a mid to late eocene 41 6 33 9 ma sediment routing system in the spanish pyrenees we derive a full volumetric sediment budget weighted for grain size fractions partitioned between terrestrial and marine depositional sectors and we quantify sediment fl uxes between depocenters the paleo sediment routing system was controlled by syndepositional thrust tectonics and consisted of two major feeder systems eroding the high pyre nees that supplied a river system draining parallel to the regional tectonic strike and that ultimately exported sediment to coastal shallow marine and deep marine depo centers we show signifi cant changes in both the volume and grain size distribution of sediment eroded from the pyrenean mountain belt during three different time intervals in the mid to late eocene which controlled the characteristics of stratigraphy preserved in a series of wedge top basins the time averaged sediment discharge from source areas increased from 250 km 3 m y to 700 km 3 m y over the 7 7 m y interval investigated this temporal increase in sediment supply caused major westward progradation of facies belts and led to substantial sediment bypass through the terrestrial routing system to the initially marine jaca basin the grain size mix measured as size fractions of gravel sand and fi ner than sand also changed over the three time intervals integration of volumetric and grain size information from source to sink provides an estimate of the long term grain size distribution of the sediment supply comprising 9 gravel 24 sand and 67 fi ner than sand the techniques and concepts used in the escanilla study can profi tably be applied to paleo sediment routing systems in other tectonic and climatic settings and to catchments with a range of bedrock lithology and vegetation this will promote a better generic understanding of the dynamics of source tosink systems and provide a powerful tool for forward stratigraphic modeling the sediment routing system approach has the potential to contribute strongly to new models of sequence stratigraphy

himalayan tectonics explained by extrusion of a low viscosity crustal channel coupled to focused surface denudation recent interpretations of himalayan tibetan tectonics have proposed that channel flow in the middle to lower crust can explain outward growth of the tibetan plateau1 2 3 and that ductile extrusion of high grade metamorphic rocks between coeval normal and thrust sense shear zones can explain exhumation of the greater himalayan sequence4 5 6 7 here we use coupled thermal mechanical numerical models to show that these two processes channel flow and ductile extrusion may be dynamically linked through the effects of surface denudation focused at the edge of a plateau that is underlain by low viscosity material our models provide an internally self consistent explanation for many observed features of the himalayan tibetan system8 9 10

two sided orogen collision and erosion from the sandbox to the southern alps new zealand in continental collision zones the indentor is of approximately the same vertical dimension as the material being deformed thereby allowing material to move over the indentor and out of the deforming orogen this simple constraint on the height of the indentor differs from that imposed by the bulldozer model of wedge deformation and results in the formation of an orogen with two opposite dipping topographic slopes these two wedges are mechanically coupled and the mechanics are strongly influenced by the erosion conditions imposed on the wedge where the growing mountains perturb prevailing moisture laden winds the erosion pattern is highly asymmetric erosional work being concentrated at the toe of the steep wedge facing the indentor the inboard wedge concentration of erosional energy results in concentration of mechanical energy in the same region and produces the distinctive geological and geophysical patterns associated with continental collision these include rapid uplift and exposure of deep crustal assemblages adjacent to the indentor accompanied by high crustal heat flow and anomalously low seismicity the outcrop pattern along the inboard wedge evolves from that of a crustal cross section currently exposed in the southern alps to one dominated by repetitive nappe structures consisting predominantly of lower crustal units as in the western alps of europe the other wedge the outboard wedge undergoes relatively little erosion is broad with a more gentle surface slope and maintains an outcrop pattern consisting of upper crustal material

sediment and rock strength controls on river incision into bedrock recent theoretical investigations suggest that the rate of river incision into bedrock depends nonlinearly on sediment supply challenging the common assumption that incision rate is simply proportional to stream power our measurements from laboratory abrasion mills support the hypothesis that sediment promotes erosion at low supply rates by providing tools for abrasion but inhibits erosion at high supply rates by burying underlying bedrock beneath transient deposits maximum erosion rates occur at a critical level of coarse grained sediment supply where the bedrock is only partially exposed fine grained sediments provide poor abrasive tools for lowering bedrock river beds because they tend to travel in suspension experiments also reveal that rock resistance to fluvial erosion scales with the square of rock tensile strength our results suggest that spatial and temporal variations in the extent of bedrock exposure provide incising rivers with a previously unrecognized degree of freedom in adjusting to changes in rock uplift rate and climate furthermore we conclude that the grain size distribution of sediment supplied by hillslopes to the channel network is a fundamental control on bedrock channel gradients and topographic relief

a statistical framework to quantify spatial variation in channel gradients using the integral method of channel profile analysis 1 we present a statistical technique for analyzing longitudinal channel profiles our technique is based on the integral approach to channel analysis drainage area is integrated over flow distance to produce a transformed coordinate which has dimensions of length assuming that profile geometry is conditioned by the stream power law defined as e kamsn where e is erosion rate k is erodibility a is drainage area s is channel gradient and m and n are constants the slope of a transformed profile in elevation space should reflect the ratio of erosion rate to channel erodibility raised to a power 1 n this quantity is often referred to as the channel steepness and represents channel slope normalized for drainage area our technique tests all possible contiguous segments in the channel network to identify the most likely locations where channel steepness changes and also identifies the most likely m n ratio the technique identifies locations where either erodibility or erosion rates are most likely to be changing tests on a simulated landscape demonstrate that the technique can accurately retrieve both the m n ratio and the correct number and location of segments eroding at different rates where model assumptions apply tests on natural landscapes illustrate how the method can distinguish between spurious channel convexities due to incorrect selection of the m n ratio from those which are candidates for changing erodibility or erosion rates we also show how given erosion or uplift rate constraints the method can be used to constrain the slope exponent n

the influence of the evolution of the pyrenees on adjacent basins abstract important thicknesses of syn and post orogenic sediments are observed in the southern part of the aquitaine basin and in the north of the ebro basin their deposition cannot be explained only as a result of the extensional history of the early cretaceous thus the influence of the compression associated with the pyrenees formation upon the subsidence of these basins is studied europe and iberia are considered to behave as two semi infinite elastic plates the evolution through time of the sedimentary thicknesses is explained by the deflections of these two plates induced not only by the topographic load of the pyrenees but also by additional forces 10 12 n m that may have their origin in the formation of a cold lithospheric root during the shortening a series of calculations is presented for an approximately north south profile across the aquitaine basin the pyrenees and the ebro basin a good fit for the shape of the basins and the gravity anomalies was obtained with flexural rigidities of 1 5 10 23 nm for the northern plate and the outer part of the southern plate and 1 4 10 21 nm for the inner part of the southern plate

asymmetric growth of the pyrenees revealed through measurement and modeling of orogenic fluxes orogenic asymmetry is ubiquitous in terms of topography structure and metamorphism but the relative influences of climate and tectonics in determining such asymmetry remain contentious this study demonstrates that for the case of the european pyrenees orogenic asymmetry is not defined by inherited structural asymmetry nor by orographic asymmetry of precipitation but is the simple consequence of the inherent tectonic growth of a doubly vergent orogen the pyrenees record approximately 60 myr of upper crustal accretion erosion and sedimentation that developed from late cretaceous to early miocene times approximately 165 km of plate convergence was accommodated by four phases of orogenic growth that can be summarized in terms of their generic orogenic development as follows 1 early inversion of passive margin extensional faults initiating the pro wedge 2 growth by a combination of frontal accretion and underplating to the pro wedge symmetrical and steady record of erosion and foreland basin sedimentation 3 cessation of pro wedge frontal accretion with propagation of the retro wedge deformation front 4 accelerated growth of pro wedge antiformal stack induced by underplating accelerated erosion over antiformal stack and increase in sediment discharge to pro foreland basin stabilization of outer pro wedge by sediment blanketing discrete element model experiments are used to explore the natural tendency of the system towards asymmetry model output demonstrates that during active frontal accretion to the pro wedge of an orogen the maxima in erosional denudation is offset towards the pro wedge relative to the maximum depths of exhumation this asymmetry is fundamentally driven by the progressive reduction in advection velocities of material within the orogen from the pro wedge into the retro wedge hence for an efficient erosional regime a dynamic coupling between erosion and structural uplift of the pro wedge will reduce activity of the retro wedge and hence drive asymmetry in terms of low temperature thermochronometry erosion and sedimentation the models also demonstrate the tendency for pro wedge growth through frontal accretion to slow however the cessation of frontal accretion as recorded by the south pyrenean thrust belt phases 3 and 4 requires additional forcing it is proposed that the thick wedge top sedimentation recorded at this time increased taper angles and hence stabilized the wedge accretionary influx into the orogen at this time was accomplished through underplating these processes governing asymmetry do not require a spatially variable climate or a laterally heterogeneous underthrust plate