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

laboratory experiments simulating the geomorphic response to tectonic uplift 1 we present the results of an experimental study of topography dynamics under conditions of constant precipitation and uplift rate the experiment is designed to develop a complete drainage network by the growth and propagation of erosion instabilities in response to tectonic perturbations the quantitative analysis of topographic evolution is made possible by using telemetric lasers that perform elevation measurements at an excellent level of precision we focus our study on the effect of initial surface organization and of uplift rate on both the transient dynamics and the steady state forms of topography we show that the transient phase is strongly dependent on the initial internally drained area which is found to decrease exponentially with time the topography always reaches a steady state whose mean elevation depends linearly on uplift rate with a strictly positive value when uplift is zero steady state surfaces are characterized by a well defined slope area power law with a constant exponent of 0 12 and an amplitude that depends linearly on uplift rate with a strictly positive value when uplift is zero these results are consistent with a stream power law erosion model that includes a nonnegligible threshold for particle detachment uncertainty regarding the sediment transport length is resolved by calibrating the transient dynamics with a surface process model reappraising published results on the linear dependency between mean elevation or relief and denudation rate we suggest that an erosion threshold is worth considering for large scale systems

landsliding and the evolution of normal fault bounded mountains much of the tectonic and climatic history in high relief regions such as the mountains of the western u s basin and range province is contained in the morphology of hillslopes drainage networks and other landforms that range in scale from 10 1 to 101km to understand how these landforms evolve we have developed a numerical landscape evolution model that combines a detailed tectonic displacement field with a set of physically based geomorphic rules bedrock landsliding long recognized as a significant geomorphic process in mountainous topography is for the first time explicitly included in the rule set in a series of numerical experiments we generate synthetic landscapes that closely resemble mountainous topography observed in the basin and range the production of realistic landscapes depends critically on the presence of bedrock landslides and landsliding yields rates of long term erosion that are comparable in magnitude to those of fluvial erosion the erosive efficiency of bedrock landsliding implies that hillslopes may respond very quickly to changes in local base level and that fluvial erosion is the rate limiting process in steady state experimental landscapes our experiments generate power law distributions of landslide sizes somewhat similar to both field and laboratory observations thus even a simple model of bedrock landsliding is capable of quantitatively reproducing mountainous topography and landslide distributions and represents a significant step forward in our understanding of the evolution of normal fault bounded ranges

coupling and decoupling of crust and mantle in convergent orogens implications for strain partitioning in the crust topographic relief strain partitioning and crust mantle coupling within convergent orogens may be systematically related to the distribution of strength as a function of depth within continental crust this result is obtained by solving analytically for basally driven three dimensional flow within an idealized crust with a laterally invariant depth dependent newtonian viscosity two end member flow states are shown to occur in the first crustal flow is directly coupled to the underlying mantle and is controlled by stresses transmitted vertically through the crust in the second upper crustal to midcrustal flow is decoupled from the underlying mantle and is controlled by stresses transmitted horizontally through the crust this model is applied to deformation of continental crust above plate like mantle that undergoes oblique convergence across a narrow suture if no significant low viscosity zone develops within the crust crustal deformation is coupled to motion of the underlying mantle and a triangular mountain range develops surface shortening occurs at the foot of the mountains normal to the trend of the range while strike slip and extension are partitioned along the axis of the orogen if a low viscosity zone is initially absent but develops during crustal thickening a steep sided flat topped plateau forms crustal deformation is coupled to mantle motions on the plateau flanks but may be strongly decoupled beneath the plateau surface shortening normal to the plateau occurs at the foot of the plateau while surface extension normal to the plateau occurs at the plateau edge if a low viscosity zone is present in the lower crust prior to convergence a wide orogen with low topographic relief develops crustal deformation is decoupled from the mantle except at the edges of the deforming region and surface shortening is partitioned uniformly across the belt this study shows that horizontal and vertical partitioning of strain occur even in the absence of lateral strength variations and that shortening and extension may exist simultaneously in different parts of an orogen similar morphology and strain partitioning can be observed in many real orogenic belts including the eastern alps the tibetan plateau and the lachlan tasman fold belt

supply and removal of sediment in a landslide dominated mountain belt central range taiwan abstract a strong coupling between hillslope and valley systems is often inferred for mountain landscapes dominated by bedrock landsliding we reveal the nature of this link using data sets on landsliding and sediment transport from two montane catchments draining the eastern central range of taiwan here the magnitude frequency distribution of landslides can be modeled by a robust power law but this scale invariance is not mirrored in the sediment discharge at the mountain front instead downstream sediment loads reflect a complex response to both sediment supply and ambient hydraulic conditions the rivers do not transport significant amounts of sediment unless it is provided by hillslope mass wasting in the catchment removal of landslide debris is a function of the transport capacity of the stream at the site of entry thus there is a dual supply and transport control on sediment loads in bedrock floored streams over a monitoring period of 25 yr the bulk of the sediment leaving the mountain bel

controls on the strength of coupling among climate erosion and deformation in two sided frictional orogenic wedges at steady state 1 many important insights regarding the coupling among climate erosion and tectonics have come from numerical simulations using coupled tectonic and surface process models however analyses to date have left the strength of the coupling between climate and tectonics uncertain and many questions unanswered we present an approximate analytical solution for two sided orogenic wedges obeying a frictional rheology and in a condition of flux steady state that makes explicit the nature and sensitivity of the coupling between climate and deformation we make the simplifying assumption that the wedge grows in a self similar fashion consistent with airy isostasy such that topographic taper is invariant with orogen width tectonic influx rate and climate we illustrate first how and why the form of the erosion rule matters to orogen evolution and then derive a physically based orogen scale erosion rule we show that steady state orogen width crest elevation and crustal thickness are controlled by the ratio of accretionary flux to erosional efficiency to a power dictated by the erosion process remarkably we show that for most combinations of parameters in the erosion law rock uplift rate is more strongly controlled by erosional efficiency than it is by the accretionary flux further assuming frontal accretion with no underplating the spatial distribution of erosional efficiency dictates the relative rock uplift rates on the pro wedge and retro wedge and the time averaged trajectories of rocks through the orogen the restriction to invariant frictional properties is conservative in these respects systems subject to positive feedback between erosion and rheology will exhibit even stronger coupling among climate erosion and deformation than shown here

sediment transport mechanisms on soil mantled hillslopes landscape evolution is modeled widely using a simple creep law for complex processes of sediment transport here field data show how a new transport model combined with an exponential soil production law better captures spatial variations of soil thick ness on hillslopes we combine parameterizations of simple and depth dependent creep with overland flow to predict soil thickness and suggest how soil distribution evolves in response to climatic and tectonic forcing we present an empirical expression for the response time of the system to external forcing that shows strong dependence on relief and is independent of soil production rate we suggest that this parameterization may be used to quantify upland carbon storage and removal and predict impacts of defor estation or rapid climatic changes

late cenozoic uplift of mountain ranges and global climate change chicken or egg the high altitudes of most mountain ranges have commonly been ascribed to late cenozoic uplift without reference to when crustal thickening and other tectonic processes occurred deep incision and recent denudation of these mountain ranges abundant late cenozoic coarse sediment near them and palaeobotanical evidence for warmer climates where high mountain climates today are relatively cold have traditionally been interpreted as evidence for recent uplift an alternative cause of these phenomena is late cenozoic global climate change towards lower temperatures increased alpine glaciation a stormier climate and perturbations to humidity vegetative cover and precipitation

mantle structure beneath the western united states and its implications for convection processes 1 we present tomographic images of the mantle structure beneath the western united states our dynamic north america models of p and s velocity structure dna07 p and dna07 s use teleseismic body waves recorded at 600 seismic stations provided by the earthscope transportable array and regional networks dna07 p and s benefit from the unprecedented aperture of the network while maintaining a dense station distribution providing high resolution body wave imaging of features through the transition zone and into the lower mantle the main features imaged include 1 the juan de fuca subduction system that bottoms at 400 km beneath oregon implying interaction with the yellowstone anomaly 2 a low velocity conduit beneath yellowstone national park that bottoms at 500 km and dips toward the northwest 3 shallow low velocity anomalies upper 200 km beneath the eastern snake river plain esrp and the high lava plains and a deep low velocity anomaly 600 km beneath the esrp but not newberry 4 a low velocity slab gap to 400 km depth immediately south of the mendocino triple junction and south of the gorda slab and 5 high velocity drips beneath the transverse ranges the southern central valley sierra nevada and central nevada these observations reveal extremely heterogeneous mantle structure for the western united states and suggest that we are only just beginning to image the complex interactions between geologic objects the transportable array allows for analysis of the relationships between these anomalies in an internally consistent single tomographic model the dna velocity models are available for download and slicing at http dna berkeley edu

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

fluvial landsculpting and the fractal dimension of topography abstract quantitative models of landform development can help us to understand the evolution of mountains and regional topography and the effects of tectonic motions and climate on landscape including its fractal geometry this paper presents a general and powerful three dimensional model of fluvial erosion and deposition at hill to mountain range scale the model works by accumulating the effects of randomly seeded storms or floods precipitons that cause diffusional smoothing then move downslope on digital topography grids that erode portions of elevation differences that transport a slope limited amount of eroded material and that deposit alluvium when their sediment carrying capacity is exceeded the iteration of these simple and almost linear rules produce very complicated simulated landscapes demonstrating that complex landscapes do not require complex laws each process implemented in the model is affected differently by changes in horizontal scale erosion a scale free process roughens topography at all wavelengths this roughening is balanced by diffusive processes scaling as 1 l 2 at short wavelengths and deposition scaling scaling as 1 l 2 at long wavelenghts such a mixture of scale free and scale dependent processes can produce multifractal behavior in the models the fractal dimension of the model topography is much more sensitive to climatic variables than to tectonic uplift landscape evolution may be fractal but it does not seem to be chaotic analysis of topography of areas in southern arizona using variograms shows approximately fractal behavior with mean fractal dimension around 2 2 2 3 departures from an exact fractal relatioship imply that the topography is in detail multifractal the fractal dimension at short wavelengths is less than that at long wavelengths this variation could either be caused by the relative strengths of diffusive and erosional processes shaping the topography or a result of changes in climatic or tectonic conditions still preserved in the landscape

dating topography of the sierra nevada california using apatite u th he ages the upward motion of rock masses relative to the earth s surface has been documented for most of the main mountain belts using thermochronological and petrological techniques more fundamental to the physical processes of mountain building however is the motion of the earth s surface itself which remains elusive here we describe a technique for estimating the age of topographic relief by mapping the low temperature thermal structure imparted by river incision using the ages of apatites determined from their uranium thorium and helium contents the technique exploits horizontal variations in temperature in the shallow crust that result from range normal river drainages because cooling beneath ancient river valleys occurs earlier than beneath intervening ridges our results from the sierra nevada california indicate that two of the modern transverse drainages the kings and the san joaquin had developed deep canyons by the late cretaceous period suggesting that the high topography of the range is 50 60 million years older than generally thought

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

quantitative testing of bedrock incision models for the clearwater river nw washington state 1 the incision of rivers in bedrock is thought to be an important factor that influences the evolution of relief in tectonically active orogens at present there are at least six competing models for incision of bedrock rivers but these models have received little quantitative testing we statistically evaluate these models using observations from the clearwater river in northwestern washington state which crosses the actively rising forearc high of the cascadia margin a previous study has used fluvial terraces along the clearwater to estimate bedrock incision rates over the last 150 kyr they show that incision rates have been steady over the long term 50 kyr consistent with other evidence based on isotopic cooling ages for steady long term 1 myr erosion rates the steady state character of the river allows us to use the relatively simple time invariant solutions for the various incision models and also to estimate long term sediment discharge along the river which is a critical variable for some incision models an interesting feature of the clearwater river is that it has a downstream decrease in the rate of incision from 0 9 mm yr in the headwater to 0 1 mm yr at the coast none of the incision models including the shear stress model successfully accounts for this relationship this result may be due to the simple way in which these models are used commonly without consideration for the distribution of discharge with time and the variable capacity of the river channel to contain peak flows along its course we suggest some general improvements for the incision models and also guidelines for selecting those rivers that will allow good discrimination between competing models

a new surface processes model combining glacial and fluvial erosion we have developed a new surface processes model incorporating large scale fluvial processes local hill slope processes and glacial erosion ice thickness and velocity are calculated under a shallow ice approximation simulation experiments in fast growing orogens comparing the efficiencies of fluvial and glacial erosion where the two are operating simultaneously over several glacial cycles show that glacial landscapes can support greater ice masses than fluvial landscapes glacial valley and lake shapes create a disequilibrium between landform and land forming process that leads to pulses of high erosion at the end of glacial periods glacial erosion rates can reach a constant value in a uniformly growing orogen and glacial erosion is capable of eroding drainage divides when the ice is moderately thick

footwall topographic development during continental extension 1 we examine the progressive development of footwall topography associated with a set of active normal faults in the northeastern basin and range province of the western united states fault length and displacement increase monotonically from northeast to southwest in the study area allowing us to track both variations in footwall morphology with increasing displacement and along strike changes in morphology on a single fault we show that patterns of catchment area footwall relief and catchment outlet spacing vary predictably and are related to the growth of the range bounding normal fault array in this semiarid region full parsing of footwall drainage area and removal of antecedent topography do not occur until fault arrays grow beyond two crustal scale segments multiple segment faults with lengths of up to 150 km have footwall relief that is limited to 1000 m in the center of the footwall and that decays to zero at the fault tips over a length scale of 15 km we hypothesize that this relatively uniform footwall relief is erosionally limited and reflects the efficacy of surface processes in removing footwall material in the center of the footwall if the fault array grows by relatively steady propagation of the tips we suggest that the 15 km length scale required to reach uniform relief is related to a timescale of relief generation by the fault tip propagation rate while such propagation rates are poorly known an average rate of 10 mm yr 1 would imply footwall relief generation over a timescale of 1 myr

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

topographic evolution and morphology of surfaces evolving in response to coupled fluvial and hillslope sediment transport 1 this paper quantifies how the ratio of sediment transport on hillslopes to sediment transport in channels influences surface and channel network morphologies and the dynamics of topographic evolution this problem is investigated by development and investigation of a simple deterministic model incorporating mass balance of sediment and runoff coupled with a law combining dispersive and concentrative sediment transport processes our analysis includes the identification of a new nondimensional parameter de that is a function of rainfall system size rock type and hydraulic regime and that is a measure of the relative importance of fluvial and hillslope sediment transport we show that de has an important influence on the surface morphology e g total exposed surface area and interface width which reflects surface roughness and relief channel network form e g channel sinuosity channel spacing and timescale of surface evolution surface and channel network morphologies are also strongly influenced by the overall surface slope relative to the magnitude of initial topographic roughness topographic evolution occurs in distinct phases of relief growth and decay the transition between which is controlled by a saturation phenomenon related to the growth of spatial correlations the scaling behavior of simulated topography with respect to both time and space is obtained and is shown to be independent of de roughness exponents are found to be independent of de but dependent on the magnitude of initial roughness interface width is shown to grow and decay as a logarithm of time

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

geomorphic limits to climate induced increases in topographic relief recognition of the potential for strong dynamic coupling between atmospheric and tectonic processes has sparked intense cross disciplinary investigation and debate on the question of whether tectonics have driven long term climate change or vice versa it has been proposed that climate change might have driven the uplift of mountain summits through an isostatic response to valley incision because isostasy acts to compensate mean elevations the debate hinges on the question of whether climate change can significantly increase topographic relief or more precisely increase the volume of missing mass between summits and ridges here we show that in tectonically active mountain ranges geomorphic constraints allow only a relatively small increase in topographic relief in response to climate change thus although climate change may cause significant increases in denudation rates potentially establishing an important feedback between surficial and crustal processes neither fluvial nor glacial erosion is likely to induce significant isostatic peak uplift

river transport effects on compressional belts first results from an integrated analogue numerical model 1 erosion deposition and tectonic processes interact to define the style of deformation in compressional belts this paper introduces a method of integrating a numerical model for surface transport with analogue sandbox models of compressional belts to further improve our understanding of this interplay the analogue model provides the numerical model with the rising topography whereas the numerical model calculates the river network and the resulting amount of erosion and sediment transport compared to previous analogue models with uniform two dimensional erosion this method provides a tool to incorporate process based rates of erosion and deposition with spatial distributions that depend on the river network developed on the evolving topography rivers running parallel to the evolving mountain belt induce along strike changes in erosion and deposition that in turn result in asymmetric tectonic deformation the additional load of sediments controls the propagation of thrust as well as the direction of vergence of the main thrust high erosion rates do not always lead to fewer thrusts within the compressional belts but generally steepen and lengthen the lifetime of the affected faults

apatite u th he age constraints on the development of the great escarpment on the southeastern australian passive margin the southeast australian margin like other high elevation passive margins is characterised by a steep escarpment that separates a dissected coastal plain from a low relief inland plateau quantitative constraints on the generation of escarpments can be provided by apatite u th he ages here we use a coast perpendicular traverse across the coastal lowlands escarpment and plateau to test the three prevailing models of se australian escarpment formation namely retreat into a downwarped rift shoulder escarpment retreat and down wearing on high elevation rift shoulder with flexural rebound apatites from the coastal plain have he ages of between 87 and 112 ma suggesting that the coastal lowlands developed very rapidly after rifting and continental break up at 85 100 ma the he age data are inconsistent with the erosion of a downwarped rift margin and cannot be explained by a constant post break up rate of lateral escarpment retreat across the coastal plain or by constant down wearing the data require either rapid escarpment retreat or rapid in place excavation of the escarpment soon after break up in response to rifting and the lowering of base levels on the margin of the new continent at break up followed by a period of landscape stability and low erosion combined with the existing apatite fission track record the he data are consistent with erosion of 3 4 km within a maximum of 28 myr of break up at a minimum vertical erosion rate of 130 m myr 1 along the coast the rapid denudation period across the coastal plain in this region took less than 48 myr from the coast to the escarpment base which corresponds to an average vertical erosion rate of 45 m myr 1 this is equivalent to a mean escarpment retreat rate of 5 10 km myr 1 apatite he ages from the plateau 183 247 ma indicate that the highlands remained stable throughout continental break up experiencing average erosion rates of less than 10 m myr 1 since the late palaeozoic early mesozoic

steady state erosion of critical coulomb wedges with applications to taiwan and the himalaya 1 orogenic structure appears to be partially controlled by the addition to and removal of material from the mountain belt by tectonic accretion and geomorphic erosion respectively we developed a coupled erosion deformation model for orogenic wedges that are in erosional steady state and deform at their coulomb failure limit erosional steady state is reached when all material introduced into the wedge is removed by erosion that is limited by the rate at which rivers erode through bedrock we found that the ultimate form of a wedge is controlled by the wedge mechanical properties sole out depth of the basal decollement erosional exponents basin geometry and the ratio of the added material flux to the erosional constant as this latter ratio is increased wedge width and surface slopes increase we applied these models to the taiwan and himalayan orogenic wedges and found that despite a higher flux of material entering the former the inferred ratio was larger for the latter calculated values for the erodibility of each wedge showed at least an order of magnitude lower value for the himalaya relative to taiwan these values are consistent with the lower precipitation regime in the himalaya relative to taiwan and the exposure of crystalline rocks within the himalayan orogenic wedge independently determined rock erodibility estimates are consistent with the accretionary wedge sediments and metasediments and the crystalline and high grade metamorphic rocks exposed within taiwan and the himalaya respectively therefore differences in rock type and climate apparently lead to key differences in the erosion and hence orogenic structure of these two mountain belts

the tectonic evolution of the southern alps new zealand insights from fully thermally coupled dynamical modelling summary we present here the results of numerical models of the southern alps in which the thermal and dynamical developments of the orogen are fully coupled to one another this interlinkage allows the use of a wide range of thermal and physical features of the orogen as constraints on the applicability of our models in particular the thermal aspect of this method enables the prediction of the distribution of thermochronological ages at the surface of the model perturbation of the geothermal structure due to rapid uplift and exhumation causes an intrinsic weakening and concentration of strain along the equivalent of the alpine fault zone in these models this thermal weakening of the crust also produces a zone of high strain antithetic to the alpine fault in the upper crust that is comparable to the main divide fault zone within the southern alps the introduction of orographic rainfall and erosional processes into the model leads to the development of high asymmetric topography comparable to that of the southern alps this topographic profile results from the capture of available precipitation by the windward side of the orogen and the resultant rain shadow on the leeward side due to the time required to accumulate suycient topography for this rain capture evect to become significant the establishment of this high asymmetric topography lags the initiation of surface uplift by several million years comparison of the observed thermal history of the southern alps with that seen in models based on divering hypotheses of the tectonic evolution of the south island shows that the present tectonic regime of the orogen most probably developed in a single rapid reorganization of plate motions at approximately 5 ma with relative stability of the tectonic regime since that time the variation in isotopic ages along the southern alps is consistent with that expected from variation in accumulated uplift and exhumation along the orogen arising from the obliquity of convergence of the australian and pacific plates

erosional dynamics flexural isostasy and long lived escarpments a numerical modeling study erosional escarpments are common features of high elevation rifted continents fission track data suggest that these escarpments form by base level lowering and or marginal uplift during rifting followed by lateral retreat of an erosion front across tens to hundreds of kilometers previous modeling studies have shown that this characteristic pattern of denudation can have a profound impact upon marginal isostatic uplift and the evolution of offshore sedimentary basins yet at present there is only a rudimentary understanding of the geomorphic mechanisms capable of driving such prolonged escarpment retreat in this study we present a nonlinear two dimensional landscape evolution model that is used to assess the necessary and sufficient conditions for long term retreat of a rift generated escarpment the model represents topography as a grid of cells with drainage networks evolving as water flows across the grid in the direction of steepest descent the model accounts for sediment production by weathering fluvial sediment transport bedrock channel erosion and hillslope sediment transport by diffusive mechanisms and by mass failure numerical experiments presented explore the effects of different combinations of erosion processes and of dynamic coupling between denudation and flexural isostatic uplift model results suggest that the necessary and sufficient conditions for long term escarpment retreat are 1 incising bedrock channels in which the erosion rate increases with increasing drainage area so that the channels steepen and propagate headward 2 a low rate of sediment production relative to sediment transport efficiency which promotes relief generating processes over diffusive ones 3 high continental elevation which allows greater freedom for fluvial dissection and 4 any process including flexural isostatic uplift that helps to maintain a drainage divide near an escarpment crest flexural isostatic uplift also facilitates escarpment retreat by elevating topography in the vicinity of an eroding escarpment thereby increasing channel gradients and accelerating erosion which in turn generates additional isostatic uplift of all the above conditions high continental elevation is common to most rift margin escarpments and may ultimately be the most important factor

late quaternary erosion in southeastern australia a field example using cosmogenic nuclides late quaternary rates of apparent soil production bedrock incision and average erosion are determined for the southeastern highlands of australia using in situ produced cosmogenic nuclide concentrations of 10 be and 26 al apparent soil production rates define a steep inverse exponential function of soil depth with a maximum of 143 m ma 1 under zero soil depth there were no observed soil depths between about 25 cm and zero however such that the maximum observed rate is about 50 m ma 1 the bredbo river catchment average erosion rate is 157 1mm a 1 and is similar to the average hillslope erosion rate of 167 1mm a 1 bedrock incision rates average 9 mma 1 and suggest that the higher rates of hillslope erosion may be in response to a pulse of incision perhaps generated by knickpoint propagation bedrock erosion rates inferred froma tor profile average 3 8 mma 1 with higher rates on other more weathered tor tops an aboveground tor profile of nuclide concentrations is consistent with a simple model of rapid stripping of the surrounding saprolite supporting the view that at least one episodic period of increased denudation has affected the landscape evolution of the highlands we test this hypothesis by using a simple landscape evolution model to reasonably predict the spatial variation of soil depth as well as the emergence of tors r 2001 elsevier science ltd and inqua all rights reserved

soil creep problems raised by a 23 year study in australia in october 1965 and february 1966 55 young pits were installed in tropical northern territory nt and temperate new south wales nsw pits were monitored in 1968 1971 and 1974 also for the nt only in 1988 in each region half of the pits are on weathered granite and half on sandstone local relief is 30m or less and slopes are up to 20 annual rainfall is evenly distributed in the nsw sites 800 mm a 1 but is confined to the five to six month wet season in nt 1200 mm a 1 six pits suffered external disturbance and so were not analysed analysis of 160 rods in 49 undisturbed pits shows 1 vectorial movement generally not downslope parallel to the ground surface but dominated by a vertically downward component 2 significant uphill and vertically upward components of movement for many rods 3 a weak correlation between total movement and sine of slope 4 rapid movement during 1965 68 and slow movement thereafter 5 significantly higher creep rates on the nt granites than on all other sites perhaps because mound building termites are especially active there we conclude that our data do not support soil creep models which assume that all movement is downslope and slope parallel repeated long term measurements are essential to distinguish long term creep rates from the short term effects of disturbance copyright 1999 john wiley sons ltd

cosmogenic nuclides topography and the spatial variation of soil depth abstract if the rate of bedrock conversion to a mobile layer of soil depends on the local thickness of soil then hillslopes on uniform bedrock in a landscape approaching dynamic equilibrium should be mantled by a uniform thickness of soil conversely if the depth of soil varies across an actively eroding landscape then rates of soil production will also vary and consequently the landscape will not be in morphologic equilibrium the slow evolution of hillslopes relative to the tempo of climatic variations and tectonic adjustments would suggest that local morphologic disequilibrium may be expected in many landscapes here we explore this issue of equilibrium landscapes through a previously developed model that predicts the spatial variation in thickness of soil as a consequence of the local balance between soil production and erosion first we confirm the assumption in the model that soil production varies inversely with the thickness of soil using two independent methods one method uses the theoretical prediction that at local steady state soil production equals removal the depth of soil should vary inversely with hillslope curvature the second method relies on direct measurements of in situ produced concentrations of cosmogenic 10 be and 26 al in bedrock at the base of the soil column for our study site in northern california the two methods agree and yield the expression that the rate of soil production declines exponentially with the thickness of soil from 0 077 mm year with no soil mantle to 0 0077 mm year under 1 m of soil we then use this function of soil production in a coupled production and diffusive model of sediment transport to explore the controls on the spatial variation of the depth of soil on four separate spur ridges noses where we measured the data for the function of soil production model predictions are sensitive to boundary conditions grid scale and run time nonetheless we found good agreement between predicted and observed depths of soil as long as we used the observed function of soil production the four noses each have spatially varying curvature and consequently have varying depths of soil implying morphologic disequilibrium we suggest that our study site has been subjected to a wave of incision and varying intensities of erosion because of tectonic and climatic oscillations that have a frequency shorter than the morphologic response time of the landscape

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

the continental collision zone south island new zealand comparison of geodynamical models and observations the south island zone of oblique continent continent convergence occurs along a 400 km long section of the modern australia pacific plate boundary zone across which about 50 km of shortening has been accommodated since about 10 ma the orogen comprises a central mountain range southern alps flanked on both sides by what are interpreted to be foreland basins two essential features that characterize the orogen are 1 the degree of denudation that accompanied deformation and 2 a fundamental structural asymmetry the architectural asymmetry of the orogen can be explained by plane strain finite element models of continental convergence incorporating mantle subduction comparison of model and orogen polarity implies that pacific plate mantle subducts the models predict two crustal scale dipping shear zones that form above the point where the pacific mantle subducts the localized one more distant from the incoming plate retro step up shear zone corresponds to the alpine fault whereas its conjugate pro step up shear zone corresponds to the distributed strain and thrusting along the eastern margin of the mountain belt parameters that modify the model boundary conditions top surface degree of denudation basal zone subduction load crust mantle velocity discontinuity subduction of lower crust mantle retreat and distributed decrease in mantle velocity and the internal strength of the crust two layer crust with moderate coupling temperature distribution strain weakening are varied in a series of numerical model calculations that establish the combination of material properties and boundary conditions that lead to different cross sectional architectures of the modeled collision zone in turn these are compared with observations about the south island orogen the calculations show how the style and extent of deformation across the whole orogen depend on the rheological properties of the crustal layer and on the balance between its internal strength and the combined effects of the boundary and gravitational stresses north to south along strike differences in the width and two dimensional architecture of the orogen simulated in the experiments by varying the model parameters can be explained by a combination of southward increases in preconvergent crustal thickness geothermal gradient convergence and potentially subduction retreat with the added possibility of a southward decrease in the component of lower crustal subduction

sediment flux from a mountain belt derived by landslide mapping in humid uplands landsliding is the dominant mass wasting process in the western southern alps of new zealand landslides are scale invariant and have a power law magnitude frequency distribution independent studies from other regions suggest that this is a general property of landsliding this observation is of critical importance to the evaluation of the impact of events of different length scales over different time intervals on landscape evolution it is particularly useful when estimating regional geomorphic rates because it constrains the frequency and overall significance of extreme events which cannot otherwise be evaluated by integrating the complete response of the system we estimate the regional denudation rate due to landsliding to be 9 4 mm yr 1 sediment discharge from the western southern alps is dominated by landslide derived material

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

fluvial response to horizontal shortening and glaciations a study in the southern alps of new zealand 1 it has been postulated that a steady state between erosional and tectonic processes may develop in continental collision however it is not clear whether steady state conditions can be reached for all components of the landscape here we show using landscape evolution models and field evidence that a true geomorphic steady state may never be reached in the southern alps of new zealand the strong asymmetries in tectonic uplift and tectonic advection and the onset of glaciations constantly interact to prevent the landscape from reaching a topographic steady state evidence suggests that the first order geomorphology on the western side of the southern alps is controlled by orographic precipitation combined with extreme rates of tectonic uplift whereas the development of deep glacial valleys on the eastern side is initiated by differential uplift along large faults we also develop a first order equation governing the dynamics of the main divide to show that both tectonic advection and fluvial erosion efficiency control the position and the height of the main drainage divide using a two dimensional landscape evolution model we demonstrate that the transition from glacial to fluvial conditions at the end of the last glaciation led to substantial modifications of the landscape while the main trunk channels get slowly uplifted ridges are leveled down causing the relief to decrease hillslopes appear to be affected by fluvial processes which seem to be driven by incision of river tributaries this reduction of relief will probably never reach a steady state since warmer interglacial periods are substantially shorter than glacial periods

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

a dynamic model to investigate coupling between erosion deposition and three dimensional thin plate deformation 1 this paper presents a mathematical model developed to investigate fully coupled interactions between erosion deposition and dynamic three dimensional deformation within fold belt settings the mechanical part of the model is an elastic plastic thin plate overlying an inviscid substrate this plate can be compressed from the side and responds by the development of buckle folds the surface topography is advected along with the underlying deforming plate and can be modified through time by a combination of fluvial and hillslope sediment transport the resulting mass redistribution creates loads on the plate which then feed back and influence subsequent deformation erosion and deposition preliminary results calculated with the model indicate that erosion and deposition have two main effects on deformation the first is to decrease the amount of shortening required to initiate folding while the second is to increase the fold wavelength both effects result from the ability of erosion and deposition to reduce the influence of gravity on deformation in the extreme case of efficient erosion and deposition deformation becomes localized on a single large scale mega anticline as opposed to being distributed on a train of anticlines and synclines this study indicates that erosion and deposition play a major role in governing the deformation and topography of fold belts and the sediment routing system acting within such a setting

three dimensional landscape process modelling the effect of dem resolution many landscape models have been developed over the past decades however relatively little is known about handling the effects of changing spatial and temporal resolutions therefore resolution effects remain a factor of uncertainty in many hydrological and geomorphological modelling approaches in this paper we present an experimental multi scale study of landscape process modelling emphasis was laid on quantifying the effect of changing the spatial resolution upon modelling the processes of erosion and sedimentation a simple single process model was constructed and equal boundary conditions were created using artificial digital elevation models dems eliminated effects of landscape representation the only variable factors were dem resolution and the method of flow routing both steepest descent and multiple flow directions our experiments revealed an important dependency of modelled erosion and sedimentation rates on these main variables the general trend is an increase of erosion predictions with coarser resolutions an artificial mathematical overestimation of erosion and a realistic natural modelling effect of underestimating resedimentation cause this increasing the spatial extent eliminates the artificial effect while at the same time the realistic effect is enhanced both effects can be quantified and are expected to increase within natural landscapes the modelling of landscape processes will benefit from integrating these types of results at different resolutions

sediment transport capacity of sheet and rill flow application of unit stream power theory theoretical equations for calculating the unit stream power of both sheet and rill flow were developed and used to predict the sediment transport capacity of such flows independent data sets from three sources representing both finely aggregated clay soils and coarse textured nonaggregated soils sheet rill and composite sheet rill flow systems and a range of slopes were used to test the utility of the method the results were very good and demonstrated the simplicity and robustness of the method for shallow overland flow the best results were obtained when the critical unit stream power at incipient sediment motion was treated as a constant value that was independent of slope the results also suggest that a unique value of critical unit stream power for rill initiation exists that is independent of soil type for noncohesive loams or fine sands and finely aggregated clay soils the sediment transport capacity can be accurately predicted from a knowledge of the physical characteristics of the soil or bed material alone for aggregated clay soils this requires information on the aggregate size distribution and the effects of soil particle size differentiation as flow rates and unit stream powers increase with the transition from sheet to rill flow

growth and lateral propagation of fault related folds in the siwaliks of western nepal rates mechanisms and geomorphic signature 1 we study the controls on drainage development in tectonically active regions using a numerical tectonic surface processes model combining tectonic uplift caused by fault related folding with erosion by fluvial incision hillslope diffusion and landsliding our model shows the fundamental control exerted by the dip of the detachment underlying the folds on drainage evolution when the detachment is horizontal the relative rates of tectonic uplift and fluvial incision control the evolution for a nonzero dip in contrast the lateral displacement gradient associated with fold propagation sets up a lateral slope behind the active structure which deflects the stream network we also demonstrate the importance of landsliding for the attainment of realistic and steady state topography in contrast discontinuous tectonic movements do not appear to influence the system we apply our model to the dundwa fault related fold ridge siwalik nepal we estimate a remarkably low mean propagation rate for different segments of the structure this finding together with the structure and morphology of the ridge leads us to propose that the ridge developed by linkage of several component segments the drainage evolution predicted when modeling this scenario compares favorably with field observations our models provide insights into the dynamics of fault related fold propagation in particular as a result of fault segment linkage estimates of propagation rate of these structures may be strongly scale dependent and the observed morphology and drainage patterns on fault related fold ridges may be controlled by fault geometry rather than by the relative rates of tectonics and surface processes

natural controls of fluvial denudation rates in major world drainage basins we present a new compilation of estimates of modern rates of mechanical and chemical denudation for externally drained basins exceeding 5 105 km2 in area these estimates are based on sediment and solute load data selected in order to represent natural rates as far as possible chemical denudation rates have been calculated by deducting the nondenudational component of solute load mechanical denudation rates range from 1 mm kyr 1 for the st lawrence and dnepr basins to 670 mm kyr 1 for the brahmaputra basin chemical denudation rates vary from 1 mm kyr 1 kolyma niger nile and rio grande basins to 27 mm kyr 1 chiang jiang basin the kolyma basin has the lowest 4 mm kyr 1 and the brahmaputra basin the highest overall rate of denudation 688 mm kyr 1 relationships between denudation rates and a range of morphometric hydrologic and climatic variables are investigated through correlation and regression analysis morphometric variables such as mean local relief are accurately calculated for large basins for the first time by using the national geophysical data center 10 minute topographic database variables expressing basin relief characteristics and runoff are found to be most strongly associated with both mechanical and chemical denudation rates with more than 60 of the variance in total denudation being accounted for by basin relief ratio and runoff basin area runoff variability and mean temperature however are only weakly associated with rates of denudation although direct comparisons cannot be made it appears that rates of basin denudation derived from present day mass flux estimates are not overall significantly different from estimates of long term rates based on sediment volume and thermochronologic data it therefore appears that the key factors identified as controlling denudation rates here are also applicable to the geological time spans relevant to the interaction between tectonic and denudational processes

a relict landscape in the centre of fennoscandian glaciation cosmogenic radionuclide evidence of tors preserved through multiple glacial cycles the presence of well developed tors boulder fields and weathering mantles in the parkajoki area of northeastern sweden near the centre of fennoscandian glaciation has been used to suggest that these landscapes were preserved during all glacial cycles since ice sheet initiation in the late cainozoic this implies that all successive large scale glaciations must have had frozen bed conditions across this area to allow for subglacial landscape preservation cosmogenic 10be and 26al data from three tors and a meltwater channel in the parkajoki area were used to test this hypothesis of landscape preservation through multiple glacial cycles apparent exposure ages of tor summit bedrock surfaces ranging between 79 and 37 ka in an area deglaciated at 11 ka are consistent with the interpretation of these features as relict landforms that have survived glaciation with little or no erosion single nuclide minimum exposure age data require that the tors have survived at least two complete glacial cycles this estimate is based on i the approximate duration of periods of ice sheet cover versus ice free conditions as deduced from the dsdp 607 marine benthic foraminifer oxygen isotope record in conjunction with ii a record of fennoscandian ice sheet flow traces and hence ice sheet extent and iii noting that cosmogenic nuclides are accumulated only during ice free periods in addition mean cosmogenic 26al 10be concentration ratios from two of the sites indicate a minimum model total history of 605 ka and a maximum erosion rate of 1 6 m ma 1 thus the numerical ages confirm the overall qualitative interpretation of landscape preservation through multiple glacial cycles

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

geologic constraints on bedrock river incision using the stream power law denudation rate in unextended terranes is limited by the rate of bedrock channel incision often modeled as work rate on the channel bed by water and sediment or stream power the latter can be generalized as ka m s n where k represents the channel bed s resistance to lowenng whose variation with lithology is unknown a is drainage area a surrogate for discharge s is local slope and m and n are exponents whose values are debated we address these uncertainties by simulating the lowering of ancient river profiles using the finite difference method we vary m n and k to match the evolved profile as closely as possible to the corresponding modern river profile over a time period constrained by the age of the mapped paleoprofiles we find at least two end member incision laws ka 0 3 0 5 s for australian rivers with stable base levels and k f a 0 1 0 2 s n for rivers in kauai subject to abrupt base level change the long term lowering rate on the latter expression is a function of the frequency and magnitude of knickpoint erosion characterized by k f incision patterns from japan and california could follow either expression if they follow the first expression with m 0 4 k varies from 10 7 10 6 m 0 2 yr for granite and metamorphic rocks to 10 5 10 4 m 0 2 yr for volcaniclastic rocks and 10 4 10 2 m 0 2 yr for mudstones this potentially large variation in k with lithology could drive strong variability in the rate of long term landscape change including denudation rate and sediment yield

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

escarpment evolution on high elevation rifted margins insights derived from a surface processes model that combines diffusion advection and reaction experiments with a surface processes model of large scale 1 1000 km long term 1 100 m y erosional denudation are used to establish the controls on the evolution of a model escarpment that is related to the rifting of a continent the model describes changes in topographic form as a result of simultaneous short and long range mass transport representing hillslope diffusive processes and fluvial transport advection respectively fluvial entrainment is modeled as a first order kinetic reaction which reflects the credibility of the substrate and therefore the fluvial system is not necessarily carrying at capacity one dimensional and planform models demonstrate that the principal controls on the evolution of an initially steep model escarpment are 1 antecedent topography drainage 2 the timescale or equivalently a length scale in the fluvial entrainment reaction 3 the flexural response of the lithosphere to denudation and 4 the relative efficiencies of the short and long range transport processes when rainfall and substrate lithology are uniform a significant amount of discharge draining over the escarpment top causes it to degrade only when the top of the model escarpment coincides with a drainage divide can escarpment retreat occur for these conditions an additional requirement for retreat of a model escarpment without decline is a long reaction time scale for fluvial entrainment this corresponds to a substrate that is hard to detach by fluvial erosion and therefore to fluvial erosion that is not transport limited continuous backtilting of an escarpment due to flexural isostatic uplift in response to denudational unloading helps maintain the scarp top as a divide it is essential if the escarpment gradient is to be preserved during retreat in a uniform lithology low flexural rigidities promote steep and slowly retreating escarpments for given rainfall and substrate conditions the morphology of a retreating model escarpment is determined by the ratio of the short range diffusive and long range advective transport efficiencies a low ratio which is interpreted to correspond to a relatively arid climate and weathering limited conditions promotes steep sharp topped escarpments with straight main slopes and escarpment retreat occurs over a wide range of height scales a high ratio interpreted to correspond to a more humid temperate climate produces a convex upper slope and concave lower slope morphology and only major escarpments are predicted to preserve a high scarp gradient lithological contrasts in the model produce more complex morphologies and predict the formation of scarps crowned by an erosionally resistant caprock however resistant caprocks are not an essential requirement for model scarps to retreat we conclude that the inferred controls and model behavior are both consistent with the present day morphology of rifted continental margins and with modern conceptual models of landscape evolution

dating fluvial terraces with and profiles application to the wind river wyoming abstract fluvial strath terraces provide a record of river incision and the timing of climatic perturbations to the fluvial system dating depositional surfaces like terraces that are older than the range of 14 c however is difficult we employ a cosmogenic radionuclide crn profile technique that addresses a major problem of crn dating on such surfaces nuclide inheritance by measuring 10 be and 26 al profiles we constrain the exposure age and the mean crn inheritance for the deposit the crn profile also yields a self check on the assumptions underlying the method we report our attempts to date terraces along the wind river wy like many sequences of western north american fluvial terraces these are inferred to reflect oscillation between glacial and interglacial conditions in the headwaters previous dating of some of these terraces and the associated terraces and glacial deposits makes this a unique location to compare dating methods dates from five sites along the bull lake glacial correlative terrace wr 3 are 118 125 ka which agrees with dates on bull lake age moraines and independent age estimates on the terrace and is consistent with the model of terrace glacial relationship crn inheritance is significant and highly variable requiring it be considered despite the additional sampling complexity assuming all inheritance in wr 3 deposits arises during exhumation in the headwaters we obtain minimum mean rates of exhumation of 13 130 m my for the source rocks alternatively assuming the crns are inherited during clast transport the time of fluvial transport from source to terrace is 10 ka it increases downstream and is lower for sand than cobbles the crn ages for older terraces wr 7 300 ka and wr 15 510 ka are lower by 50 than previous estimates based on tephrochronology the most plausible explanation is eolian deflation of a once thicker loess cover on the terrace surfaces mean thicknesses of loess of 0 5 1 5 m are required to reconcile these concentrations of crn with the previous estimates of age difficulty in dating the older terraces emphasizes that geologic caution independent estimates of age and multiple sample sites should still be part of dating depositional surfaces with crns even when employing the inheritance correction technique

p and s waveform modeling of continental sub lithospheric detachment at the eastern edge of the rio grande rift 1 east of the rio grande rift tomographic images of teleseismic data have revealed a se dipping slab like structure underneath the western edge of the great plains in the southwestern united states however finite difference synthetics require an amplified tomographic model to reproduce the waveform distortions as observed in broadband waveform data recorded along the la ristra transect in addition to travel time anomalies song and helmberger 2007 demonstrated how to use s waveforms and their amplitude patterns to further constrain the magnitude of the anomalous structure their preferred s velocity model suggests that the slab like structure is about 4 fast 120 km thick and dipping 70 75 to the se to about a depth of 600 km we adapt the preferred s wave model from song and helmberger 2007 and scale the p wave model using a suite of scaling factors sf dlnvs dlnvp we find that synthetics from the p model generated by sf 1 25 can distort p waveforms and fit the amplitude pattern best such a low sf indicates that the slab like anomaly is not only cold but also compositionally distinct we make use of sf and the s wave anomaly simultaneously to obtain corresponding vp vs anomaly and infer the origin of the slab like anomaly our result suggests that the observed sub lithospheric detachment is 310 20 c colder and more depleted than the adjoining mantle asthenosphere by 3 units of mg mg mg fe 100 it is negatively buoyant and is geodynamically consistent with the observed foundering of the continental lithosphere at the eastern edge of the rio grande rift in short we demonstrate the importance of seismic waveform anomalies in geochemical and geodynamic inferences

a process based model for colluvial soil depth and shallow landsliding using digital elevation data a model is proposed for predicting the spatial variation in colluvial soil depth the results of which are used in a separate model to examine the effects of root strength and vertically varying saturated conductivity on slope stability the soil depth model solves for the mass balance between soil production from underlying bedrock and the divergence of diffusive soil transport this model is applied using high resolution digital elevation data of a well studied site in northern california and the evolving soil depth is solved using a finite difference model under varying initial conditions the field data support an exponential decline of soil production with increasing soil depth and a diffusivity of about 50 cm 2 yr the predicted pattern of thick and thin colluvium corresponds well with field observations soil thickness on ridges rapidly obtain an equilibrium depth which suggests that detailed field observations relating soil depth to local topographic curvature could further test this model bedrock emerges where the curvature causes divergent transport to exceed the soil production rate hence the spatial pattern of bedrock outcrops places constraints on the production law the infinite slope stability model uses the predicted soil depth to estimate the effects of root cohesion and vertically varying saturated conductivity low cohesion soils overlying low conductivity bedrock are shown to be least stable the model may be most useful in analyses of slope instability associated with vegetation changes from either land use or climate change although practical applications may be limited by the need to assign values to several spatially varying parameters although both the soil depth and slope stability models offer local mechanistic predictions that can be applied to large areas representation of the finest scale valleys in the digital terrain model significantly influences local model predictions this argues for preserving fine scale topographic detail and using relatively fine grid sizes even in analyses of large catchments