post oligocene river incision southern sierra madre occidental mexico the sierra madre occidental of western mexico consists of a granitic basement covered by oligocene ignimbrites that define a reference surface from which to estimate late cenozoic river incision a 90 m grid digital elevation model was used to characterize contemporary topography and interpolate the late oligocene surface of the ignimbrite plateau from a surface fit to the highest points in the relatively undissected uplands between major river valleys long term river incision rates calculated from the difference between this reference surface and longitudinal profiles of 11 rivers that flow toward the tepic zacoalco rift zone range from about 0 01 to 0 2 mm year 1 river profiles of this region also show evidence of river capture driven by flexural uplift along the flank of the rift zone river profile concavity values in the sierra madre occidental range from 0 22 to 0 63 a range similar to that reported previously for a wide range of environments in contrast the empirically constrained ratio of exponents in the stream power model of river incision m n ranges from 0 44 to 0 52 close to the expected theoretical value of 0 5 the wider range of observed values may illustrate how can differ from the driving values of m n in non steady state bedrock river systems

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

dynamic equilibrium between an active uplift and the sevier river utah in a zone of active tectonic uplift one expects a river to change morphology simply because valley slope deforms but by cutting through a transverse uplift the sevier river in southern utah also adapts to the tectonics by eroding and depositing these localized changes of transport are caused by varying channel geometry pattern sedimentology and slope in the uplift the river is narrow efficiently shaped coarser grained and expends a greater proportion of energy on bed particles it deposits sediments downstream in exaggerated meandering reaches with wider less efficient channels localized erosion and deposition create a dynamic equilibrium with the structural growth but this equilibrium includes nonlinear responses like knickpoints and depositional lobes in response to the required erosion and deposition the river so dramatically alters channel grain size and channel geometry that the increase of sinuosity expected in the tectonically steepened valley does not occur

mathematical modeling of graded river profiles numerical modeling of the longitudinal profiles of rivers at grade is accomplished using the basic equations of open channel flow sediment transport equations and empirical relations for downstream variation in flow discharge sediment discharge sediment caliber and channel width only in some cases are the computed stream profiles fit exactly by any one of the commonly supported mathematical function analogs to graded profile form exponential logarithmic or power function but in most cases any of these functions can provide a fit with a degree of error smaller than would be noted in treating field data profiles dominated by spatial change in fluid and sediment discharge are distinctly power functions while profiles dominated by sediment size reduction are not necessarily exponential in form other important controls on profile shape are the degree of downstream width change in response to increasing discharge and the general range of sediment size a dynamic model of a river s approach to grade

experimental channel response to tectonic uplift 1 we have performed a series of experiments to investigate the relationship between channel geometry and tectonic forcing in steady state landscapes at various uplift rates the experimental setup consists of uniformly uplifted silica paste eroded by artificial rainfall in this setup erosional landscapes evolve by growth and amalgamation of incisions which organize into a drainage network whose dynamics are a function of the interaction between vertical channel cutting hillslope erosion and sediment transport high precision 0 5 mm pixel size digital elevation models were constructed using a stereogrammetric camera system from this data channel bed slope was found to be independent of discharge and position on the experimental surface and to increase linearly with uplift rate geometric parameters of the flow such as channel width or depth could not be observed directly during the experiments using laminar flow equations for the mean flow velocity these parameters were back calculated and relationships to the imposed substrate uplift rate were derived channel width cross sectional area and wetted perimeter decrease with increasing uplift rate to a limit value hydraulic radius and flow depth increase slightly and flow velocity increases approximately linearly with increasing uplift rate these results are qualitatively consistent with recent field surveys and highlight the importance of channel width and slope variations in accommodating channel response to variable rock uplift rates

channel changes in badlands stripping of the vegetation and soil from a 13 hectare site in virginia underlain by coastal plain sediments created a rapidly evolving badland topography two types of channels developed 1 sand bed alluvial channels were graded to transport the bed material load supplied from slope erosion with available runoff but they also generally eroded their beds slowly and 2 steeper bedrock floored channels incised rapidly in bedrock channels the erosion rate was proportional to the 4 9ths power of drainage area and the 2 3rds power of gradient these exponents are consistent with a model in which the erosion rate is proportional to the bed shear during high flows due to rapid mass wasting and reduced runoff the alluvial channels became as much as 50 steeper during the winter than the summer with an attendant yearly cycle of winter aggradation and summer entrenchment the gradients their seasonal variability and their downstream hydraulic geometry were consistent with the predictions of total load transport formulas for sand beds and high loads the hydraulic geometry of alluvial channels in the virginia badlands were similar to that on the morrison formation in the western united states

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

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

knickpoint initiation and distribution within fluvial networks 236 waterfalls in the waipaoa river north island new zealand if knickpoints transmit signals of base level fall in river networks then improvements in our understanding of their retreat rate and basin wide distribution helps constrain the transient response following perturbation many studies of knickpoint retreat focus on the response of trunk streams to base level fall here we examine the response of an entire fluvial network as recorded by 236 active knickpoints distributed within the waipaoa river on the north island of new zealand base level fall within the waipaoa catchment initiated 18 000 years ago in response to a climatically triggered and tectonically exacerbated pulse of incision using observations from field work aerial photo analysis and a digital elevation model dem we study the knickpoint positions within the network we find that 70 of the knickpoints are located at drainage areas between 1 10 5 m 2 and 1 10 6 m 2 and more than half are b1 km upstream of a large change in drainage area for the knickpoints b1 km upstream of large tributary junctions we find that the retreat distances were well correlated with the tributaries drainage areas in order to determine how a pulse of incision distributes itself throughout a fluvial network we develop two simple end member models and compare their behavior to the observed knickpoint distribution in the waipaoa in the first model we propose that a knickpoint initiated at the basin outlet retreats upstream and distributes the signal throughout the network at a rate that is a power law function of drainage area in the second model we propose that knickpoints form near a threshold drainage area below which channels cannot incise with the same efficiency as possible in downstream reaches though neither model addresses the along stream variability in substrate or knickpoint form the misfit between the modeled and the observed knickpoints along stream positions are surprisingly low b1 km 1 25 of the total stream length for knickpoints with drainage areas b1 10 6 m 2 large misfits up to 3 5 km are observed for knickpoints with present day drainage areas greater than 1 10 6 m 2 the large single step in channel elevation that characterizes knickpoints presently observed in tributaries of the waipaoa river may not characterize the base level fall signal that propagated through the trunk streams evidence for progressive rather than instantaneous incision in the trunk streams the knickpoints vicinities to tributary junctions and the equivalent success of the two end member models lead us to conclude that the present positions of the 236 observed knickpoints are largely a consequence of thresholds in channel incision at low drainage areas

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

palaeomagnetic correlation and dating of mineralization and intrusive events abstract in the mining district of santa eulalia chihuahua northern mexico the mineralization consisting of simple massive concentrations of lead zinc iron sulphides and oxides confined in largely unaltered cretaceous limestones has traditionally been related to several intrusive emplacement events that occurred during the tertiary however the relative age relationships of the intrusives volcanics and mineralization have proved difficult to establish in this work we report palaeomagnetic results which permit definition of the mean relative age relationships for the igneous units and their correlation with the mineralizing event the mean directions and pole positions for the intrusives mineralization and volcano sedimentary cover capping series indicate a tertiary age eocene miocene the intrusives show normal intermediate and reverse polarities suggesting several emplacement events or a long period of intrusion the direction and polarity of the remanent magnetization for the mineralization correlates with those for the felsites which occur as a complex set of sills and dykes and they are therefore penecontemporaneous with the mineralization

the role of sediment in controlling steady state bedrock channel slope implications of the saltation abrasion incision model abstract the saltation abrasion model is a mechanistic model for river incision into bedrock by saltating bedload which we have previously derived and used experimental data to constrain all parameter values here we develop a method for applying the saltation abrasion model at a landscape scale and use the model as a reference for evaluating the behavior of a wide range of alternative incision models in order to consider the implications of the saltation abrasion model as well as other models for predicting topographic steady state channel slope to determine the single valued discharge that best represents the effects of the full discharge distribution in transporting sediment and wearing bedrock we assume all runoff can be partitioned between a low flow and a high flow discharge in which all bedload sediment transport occurs during high flow we then use the gauged discharge record and measurements of channel characteristics at a reference field site and find that the optimum discharge has a moderate magnitude and frequency due to the constraints of the threshold of grain motion and bed alluviation by high relative sediment supply incision models can be classified according to which of the effects of sediment on bedrock incision are accounted for using the predictions of the saltation abrasion model as a reference we find that the threshold of motion is the most important effect that should be represented explicitly followed in order of decreasing importance by the cover effect the tools effect and the threshold of suspension effect models that lack the threshold of motion over predict incision rate for low shear stresses and under predict the steady state channel slope for low to moderate rock uplift rates and rock strengths models that lack the cover effect over predict incision rate for high sediment supply rates and fail to represent the degree of freedom in slope adjustment provided by partial bed coverage models that lack the tools effect over predict incision rate for low sediment supply rates and do not allow for the possibility that incision rate can decline for increases in shear stress above a peak value overall the saltation abrasion model predicts that steady state channel slope is most sensitive to changes in grain size such that the effect of variations in rock uplift rate and rock strength may affect slope indirectly through their possible but as yet poorly understood influence on the size distribution of sediments delivered to channel networks by hillslopes

implications of the saltation abrasion bedrock incision model for steady state river longitudinal profile relief and concavity the saltation abrasion model predicts rates of river incision into bedrock as an explicit function of sediment supply grain size boundary shear stress and rock strength here we use this experimentally calibrated model to explore the controls on river longitudinal profile concavity and relief for the simple but illustrative case of steady state topography over a wide range of rock uplift rates we find a characteristic downstream trend in which upstream reaches are close to the threshold of sediment motion with large extents of bedrock exposure in the channel bed while downstream reaches have higher excess shear stresses and lesser extents of bedrock exposure profile concavity is most sensitive to spatial gradients in runoff and the rate of downstream sediment fining concavity is also sensitive to the supply rate of coarse sediment which varies with rock uplift rate and with the fraction of the total sediment load in the bedload size class variations in rock strength have little influence on profile concavity profile relief is most sensitive to grain size and amount of runoff rock uplift rate and rock strength influence relief most strongly for high rates of rock uplift analysis of potential covariation of grain size with rock uplift rate and rock strength suggests that the influence of these variables on profile form could occur in large part through their influence on grain size similarly covariation between grain size and the fraction of sediment load in the bedload size class provides another indirect avenue for rock uplift and strength to influence profile form copyright 2008 john wiley sons ltd

hanging valleys in fluvial systems controls on occurrence and implications for landscape evolution 1 we document and characterize hanging valleys in a fluvially eroded landscape in eastern taiwan our conceptual model for the initiation of hanging valleys builds on a recently proposed model of bedrock incision in which erosion actually becomes less efficient on very steep channel gradients if a pulse of incision in the main stem outpaces the tributary response the gradients at tributary mouths may therefore pass a threshold value beyond which erosional efficiency declines giving rise to a mismatch between trunk and tributary erosion rates this mismatch is expected at junctions with small tributaries where a step function decrease in drainage area also leads to sharp contrasts in water and sediment flux between trunk and tributary channels the occurrence of hanging valleys in actively uplifting fluvial landscapes such as the central range of taiwan suggests that the most common parameterizations of bedrock erosion which typically assume a monotonic positive correlation between channel gradient and incision rate may be violated in very steep channels in addition hanging valleys could greatly increase the response time of landscapes to tectonic perturbations since catchments above these tributary mouths will be insulated from these perturbations until a new suite of processes e g weathering and rock mass failure wear through the hanging valley lip the results of this study underscore the need for a more complete understanding of bedrock erosion processes and the incorporation of process transitions and threshold conditions into landscape evolution models

equilibrium and time scales in geomorphology application to sand bed alluvial streams equilibrium is defined as a single valued temporally invariant functional relationship between the values of an output variable and the values of the input variable s in a geomorphic system disequilibrium occurs if the output deviates from the functional relationship by more than a consensual degree natural geomorphic variables are characterized by a relaxation time output variables are insensitive to cyclical inputs with frequencies much greater than the relaxation time but can respond completely for sufficiently low frequencies rapid trends recent step changes or pulse inputs and intermediate frequency inputs can cause disequilibrium the gradient of sand bed alluvial channels the output variable is determined by sediment and water delivery from slopes the input variables and changes in this hydraulic regime require regrading by erosion and deposition initial stages of adjustment to changed regime in a long unbranched channel with sediment and water delivery only at the upstream end propagate downstream but later stages of adjustment occur simultaneously throughout the reach in a dendritic channel network the gradient responds rather uniformly throughout the network to changes in regime during all stages of adjustment the time scale of adjustment to changed regime depends upon the size of the channel network or stream length the sediment and water discharges and to a lesser degree upon the magnitude of the change grade as defined by mackin 1948 is synonymous with equilibrium as used in this paper if a period of years is replaced by a time period commensurate with the relaxation time of the gradient the use of the term grade is best restricted to a single valued relationship between channel gradient and the hydraulic regime

the use of a landscape simulator in the validation of the siberia landscape evolution model transient landforms siberia is a physically based model for the geomorphic evolution of landforms it is essential that the siberia model be tested or validated against controlled landform development previous studies have demonstrated that siberia is able to simulate declining equilibrium landforms and in this paper we examine siberia s ability to simulate landforms as they evolve to their declining equilibrium form these landscapes are termed transient landforms landscapes generated by siberia were compared to those produced by a physical model experimental model landforms at stages of evolution comparison of the experimental landscapes with the simulated landscapes using total mass hypsometric curve width function cumulative area distribution and area slope demonstrate that siberia can simulate the experimental model landscape during development i e transient landscapes copyright 2002 john wiley sons ltd

longitudinal profiles in simple alluvial systems theoretical considerations suggest that exponential or quadratic functions are apt for describing the longitudinal profiles of aggrading alluvial systems that are unaffected by significant lateral inputs of water or sediment a new set of field data confirms this for individual sedimentary links along a wandering gravel bed river in british columbia canada each link is viewed as a simple alluvial system without major water or sediment inputs within which grain size typically fines downstream in a relatively systematic manner such homogeneous reaches are suitable for the investigation of simple profile form it is found that quadratic approximations are the most flexible descriptor for link longitudinal profiles overall the river forms an irregularly cuspate long profile structured by these fundamental length scale units both link gradients and link fining rate reflect size selective aggradation associated with the bounding lateral inputs

calculating bed load transport in steep boulder bed channels 1 steep rough channels occupy a large fraction of the total channel length in mountainous regions most sediment mobilized on hillslopes must pass through these streams before reaching lower gradient channels steep channels have wide grain size distributions that are composed of finer more mobile sediment and large rarely mobile grains the large grains can bear a significant portion of the total shear stress and thereby reduce the stress available to move the finer sediment conventional bed load transport equations often overpredict the sediment flux in steep channels by several orders of magnitude we hypothesize that sediment transport equations overpredict the sediment flux because they do not 1 account for the stress borne by rarely mobile grains 2 differentiate between highly and rarely mobile sediment and 3 account for the limited availability of mobile sediment here we modify a conventional bed load transport equation to include these three effects we use measurements of the flow bed properties and sediment flux in a small steep flume to test this equation we supply gravel at a constant rate through fields of regularly spaced immobile spheres and measure the bed coverage by gravel and sphere protrusion the percent of the sphere that protrudes above the gravel deposit for a given sphere spacing the proportion of the bed covered by gravel increases and the sphere protrusion decreases with greater sediment supply thus bed coverage and immobile grain protrusion may serve as proxies for sediment availability in steep rough streams unlike most transport equations that we tested our modified bed load equation predicts sediment fluxes to within an order of magnitude of the measured values our results demonstrate that accurately predicting bed load transport in steep rough streams may require accounting for the effects of local sediment availability coverage by mobile sediment and drag due to rarely mobile particles

grain size along two gravel bed rivers statistical variation spatial pattern and sedimentary links a new set of field data facilitates a detailed analysis of variations in bed material grain size within two confluent gravel bed rivers in northeastern british columbia canada a preliminary assessment of grain size variability establishes a basis for examination of the spatial pattern of grain size change standard anova techniques are inappropriate because individual samples have unequal variances and are not normally distributed alternative tests for homoscedasticity and comparison of means are therefore utilized within site between sample variability is not significant the grain size distributions that were obtained at individual sites are therefore representative of the depositional environments that were sampled in both rivers mean grain size does vary significantly between sites and there is therefore a basis for examining the data for spatial patterns such as downstream fining textural variations along the two rivers studied here are complex and show negligible overall fining in over 100 km this is the consequence of a large number of tributary inputs and non alluvial sediment sources which are the legacy of late pleistocene glaciation the identification of lateral sources like these is fundamental for understanding textural changes within rivers the sedimentary link a channel reach between significant lateral sediment inputs provides a means of isolating fluvial maturation processes abrasion and sorting from contingent lateral inputs strong fining trends are apparent in most links and classification of grain size measurements according to their location within particular links greatly improves the statistical explanation of textural variation identification of sedimentary links provides a means of applying models of fluvial fining processes so isolation of link networks will aid the development of basin scale models of textural variation 1998 john wiley sons ltd

longitudinal profile development into bedrock an analysis of hawaiian channels analysis of topographic maps of rivers incised into dated hawaiian lava flows shows that the long term average bedrock erosion rate along certain reaches is linearly related to stream power field observations suggest that two processes may control hawaiian channel downcutting 1 stream power dependent erosion including abrasion of the channel bed by transported particles and 2 step wise lowering caused by knickpoint propagation modeling results indicate that a simple stream power dependent erosion law predicts the straight to weakly convex longitudinal profiles characteristic of kauai channels but is insufficient to predict two other characteristic features the upslope propagation of knickpoints and the straight 5 8 channel slopes below the knickpoints thus more than this single transport law is apparently required to model bedrock channel incision field surveys also indicate that significant portions of the channel lengths below the knickpoints are mantled with large boulders we propose that

differentiating the roles of lithology and uplift in the steepening of bedrock river long profiles an example from southeastern australia hack determined that the equilibrium river long profile on a single lithology is a straight line when plotted semi logarithmically deviation from this straight line may be an equilibrium response to lithological variation or disequilibrium following base level change the amount of deviation d ast at any point downstream of the point of departure from the semi log straight line can be measured by projection of the upstream equilibrium profile and subtraction of the actual stream elevation from the projected elevation comparisons between the values of d ast for a trunk and a tributary at their junction can be used to differentiate steepening as a result of base level change from lithologically induced steepening in the former case d ast represents the incision induced by the base level change and will be equal for the trunk and the tributary it is very unlikely that values of d ast will be equal in the case of lithologically induced steepening because d ast in an

implications of sediment flux dependent river incision models for landscape evolution 1 developing a quantitative understanding of the factors that control the rate of river incision into bedrock is critical to studies of landscape evolution and the linkages between climate erosion and tectonics current models of long term river network incision differ significantly in their treatment of the role of sediment flux we analyze the implications of various sediment flux dependent incision models for large scale topography in an attempt 1 to identify quantifiable and diagnostic differences between models that could be detected from topographic data or from the transient responses of perturbed systems and 2 to explain the apparent ubiquity of mixed bedrock alluvial channels in active orogens although certain forms of the various models can be discarded as inconsistent with morphological data we find that the relative intrinsic concavity indices of detachment and transport limited systems defined herein largely dictate whether the various models can be tied to distinctive steady state morphologies preliminary data suggest that no such diagnostic differences may exist and other methods must be developed to test models accordingly we develop and explore differences in the scaling behavior of topographic relief and the extent of detachment versus transport limited channels as a function of rock uplift rate that may allow discrimination among various models further we explore potentially diagnostic differences in the rates and patterns of transient channel response to changes in rock uplift rate in addition to general differences between detachment and transport limited systems our analysis identifies an interesting hysteresis in landscape evolution hybrid channels at the threshold between detachment and transport limited conditions are expected to act as detachment limited systems in response to an increase in rock uplift rate or base level fall and as transport limited systems in response to a decrease in rock uplift rate especially during postorogenic topographic decline the analyses presented set the stage for field studies designed to test quantitatively the various river incision models that have been proposed

cenozoic river profile development in the upper lachlan catchment se australia as a test of quantitative fluvial incision models 1 we have used early miocene valley filling basalts to reconstruct fluvial long profiles in the upper lachlan catchment se australia in order to use these as well constrained initial conditions in a forward model of fluvial incision many different fluvial incision algorithms have been proposed and it is not clear at present which one of these best captures the behavior of bedrock rivers we test five different formulations the ability of these models to reproduce the observed present day stream profiles and amounts of incision is assessed using a weighted mean misfit criterion as well as the structure of the misfit function the results show that for all models parameter combinations can be found that reproduce the amounts of incision reasonably well however for some models these best fit parameter combinations do not seem to have a physical significance whereas for some others best fit parameter combinations are such that the models tend to mimic the behavior of other models overall best fit model predictions are obtained for a detachment limited stream power model or an undercapacity model that includes a river width term that varies as a function of drainage area the uncertainty in initial conditions does not have a strong impact on model outcomes the model results suggest however that lithological variation may be responsible for variations in parameter values of a factor of 3 5

flexural models of continental lithosphere based on the long term erosional decay of topography summary the most easily observed geophysical characteristic of continents is their surface topography topographic power in the spectral range 1000 100 km has been computed for north american geological provinces and is found to decrease systematically as the tectonic age of the underlying crust increases these data have been quantitatively modelled in terms of topographic erosion and stress relaxation within the continental lithosphere through geological time spans the lithosphere is treated as a thin maxwell viscoelastic plate which deforms flexurally in response to the vertical surface loads of topography and its erosion effective rheological model parameters are flexural rigidity d a measure of the lithosphere s elastic strength and viscous relaxation time constant a measure of its viscosity erosion is assumed to occur at a rate linearly proportional to the height of topography at any given time the proportionality being determined by the erosion time constant the models are constructed in space and time frequency domains using simple linear systems theory two loading schemes are considered model 1 assumes that a flexurally competent lithosphere is suddenly loaded by topography which subsequently erodes whereas model 2 ignores the emplacement and erosion of topography during orogenesis when the lithosphere is thermally activated and considers only the load effected by erosion beginning some time later when the lithosphere has cooled and strengthened both models predict that the remnant topography of very old regions such as the canadian shield can be maintained in fact is favoured by a viscoelastic lithosphere with as small as 1 10ma the continued existence of such topography does not necessarily imply significant long term lithospheric strength model 1 results indicate low values of the flexural rigidity d 1021 1022nm which are interpreted to characterize not so much presentday d but that which was effective when topography was formed model 2 more successfully reproduces the data and constrains the product d to values 1024 1025nm ma which compare very favourably to those found from studies of individual loads on continental lithosphere the model results show that the erosion time constant is probably wavenumber dependent and falls in the range 200 400ma for the topographic wavelengths under consideration

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

bed load transport i mechanical characteristics this paper deals with the mechanical characteristics of saltation of bed load transport in flowing water experiments have been carried out by means of high speed photography and advanced data processing techniques taking saltation as the mode of bed load transport systematic analysis has been made with regard to the forces acting on grains and the variation of their magnitudes along the saltation trajectories to obtain the kinematic and dynamic properties of these grains moveover the relations between magnitudes of saltation parameters and flow intensity are attained from the measured data some parameters in bagnold s theory of bed load transport are also discussed

climate exposed source rock lithologies crustal uplift and surface erosion a theoretical analysis calibrated with data from the alps north alpine foreland basin system paleofloristic data imply that paleoclimate changed in the swiss alps at the oligocene miocene boundary from humid and hot conditions toward a climate with high temperature and low humidity the aridization is associated with a change in depositional pattern from alluvial fans to lakes and floodplains suggesting decreasing sediment discharge a further 25 40 decrease of sediment discharge occurred at ca 20 ma when the orogenic core of the alps became exposed to the surface we applied a surface processes model to explore potential controls on the pattern of sediment discharge and on the evolution of the alpine drainage basin the model is based on the presumption that the rates of fluvial incision into bedrock are proportional to shear stress exerted by the flowing water the model results imply that the paleoclimate change resulted in an instantaneous decrease of sediment discharge and a vertical topographic growth until steady state conditions between erosional and crustal mass flux are established however exposure of the crystalline core of the alps at ca 20 ma is likely to have resulted in the 25 40 decrease of sediment discharge and the reorganization of the drainage pattern from an orogen normal to an orogen parallel orientation of dispersion

causes of concavity in longitudinal profiles of rivers a smooth concave configuration represents the most common quasi equilibrium profile of alluvial rivers in the present paper a generalized unified set of models for the longitudinal profile of rivers is presented various independent effects that drive concavity in longitudinal rivers are considered and classified the driving mechanisms considered here are 1 horizontal wavelike progradation of the river profile 2 abrasion of bed particles 3 aggradation of the river profile balancing subsidence at a constant speed and 4 the effect of tributaries adding sediment and water to the main stem of the river for simplicity it is assumed that the long term supply of water and sediment is constant hence a single discharge that is bank full discharge is used as an approximate characterization of the hydrologic regime for each point along the river sediment sorting has been studied elsewhere and is excluded from this analysis the river is assumed to have uniform bed material although the mean size may vary downstream it is concluded that the quasi equilibrium longitudinal profile created by each of the first three driving mechanisms mentioned above is concave the degree of concavity varies however depending on a set of physical conditions that can be identified by a set of dimensionless numbers a comparison of these numbers reveals the relative importance of the driving mechanisms for concavity a condition for concavity driven by the fourth mechanism that is tributary input is delineated