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

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

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

evidence for large earthquakes in metropolitan los angeles the sierra madre fault along the southern flank of the san gabriel mountains in the los angeles region has failed in magnitude 7 2 to 7 6 events at least twice in the past 15 000 years restoration of slip on the fault indicated a minimum of about 4 0 meters of slip from the most recent earthquake and suggests a total cumulative slip of about 10 5 meters for the past two prehistoric earthquakes large surface displacements and strong ground motions resulting from greater than magnitude 7 earthquakes within the los angeles region are not yet considered in most seismic hazard and risk assessments

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

downstream changes in the lithology and grain size of fluvial gravels the watarase river japan evidence of the role of abrasion in downstream fining two processes have been considered to explain downstream grain size reduction in gravel bed rivers hydraulic sorting by which finer grains are transported faster and farther than coarse grains and abrasion by which individual particles are diminished in size most studies on the cause of downstream fining in rivers emphasize hydraulic sorting as the dominant fining mechanism however broken boulders are often observed in the bed material of japanese rivers also grain size changes may be correlated with the difference in resistance to breakdown among lithologies the purpose of this study is to use the lithologic composition of bed material on the watarase river to demonstrate that some downstream fining cannot be explained by hydraulic sorting alone the study reach on the watarase river flows on a dissected alluvial fan which means that the river is no longer aggrading the study reach is influenced only slightly by input of gravel from tributaries bed gravel in the study reach consists of several lithologies andesite quartz porphyry sandstone hornfels and chert the grain size distribution and lithology of the river bed material were examined at nine sites there are two principal results 1 size distribution of the gravel is strongly related to lithology andesite boulders or large cobbles make up the framework sizes in the upstream part of the study reach while chert pebbles make up the framework sizes in the downstream part there are few andesite pebbles or chert boulders in the river bed 2 longitudinal changes in the composition of each size class show selective transportation by lithology in every gravel size because the mobility of gravel depends mainly on its size hydraulic sorting by lithology does not occur within the same size class these results indicate clearly that particle abrasion occurs in the watarase river and is responsible at least in part for the downstream decrease in particle size of bed material

fe mg interdiffusion rates in clinopyroxene experimental data and implications for fe mg exchange geothermometers chemical interdiffusion of fe mg along the c axis 001 in natural diopside crystals x di 0 93 was experimentally studied at ambient pressure at temperatures ranging from 800 to 1 200 c and oxygen fugacities from 10 11 to 10 17 bar diffusion couples were prepared by ablating an olivine x fo 0 3 target to deposit a thin film 20 100 nm onto a polished surface of a natural oriented diopside crystal using the pulsed laser deposition technique after diffusion anneals compositional depth profiles at the near surface region 400 nm were measured using rutherford backscattering spectroscopy in the experimental temperature and compositional range no strong dependence of d fe mg on composition of clinopyroxene fe mg ratio between di93 di65 or oxygen fugacity could be detected within the resolution of the study the lack of fo2 dependence may be related to the relatively high al content of the crystals used in this study diffusion coefficients d fe mg can be described by a single arrhenius relation with d text fe text mg 2 7 7 pm 4 2 7 times 10 7 text exp 3 20 7 pm 1 6 0 text kj text mol text rt m 2 text s d fe mg in clinopyroxene appears to be faster than diffusion involving ca species e g d ca mg while it is slower than d fe mg in other common mafic minerals spinel olivine garnet and orthopyroxene as a consequence diffusion in clinopyroxene may be the rate limiting process for the freezing of many geothermometers and compositional zoning in clinopyroxene may preserve records of a higher compared to that preserved in other coexisting mafic minerals temperature segment of the thermal history of a rock in the absence of pervasive recrystallization clinopyroxene grains will retain compositions from peak temperatures at their cores in most geological and planetary settings where peak temperatures did not exceed 1 100 c e g resetting may be expected in slowly cooled mantle rocks many plutonic mafic rocks or ultra high temperature metamorphic rocks

bedrock incision rock uplift and threshold hillslopes in the northwestern himalayas the topography of tectonically active mountain ranges reflects a poorly understood competition between bedrock uplift and erosion dating of abandoned river cut surfaces in the northwestern himalayas reveals that the indus river incises through the bedrock at extremely high rates 2 12 mm yr 1 in the surrounding mountains the average angles of hillslopes are steep and essentially independent of erosion rate suggesting control by a common threshold process in this rapidly deforming region an equilibrium is maintained between bedrock uplift and river incision with landsliding allowing hillslopes to adjust efficiently to rapid river down cutting

evolution of the santa cruz mountains california through tectonic growth and geomorphic decay i describe a strategy for modeling a specific large scale topographic feature that recognizes the discrepancy between tectonic and geomorphic scales and allows for important feedbacks between these processes in the geomorphic model the largest channel in each 4 4 km cell is explicitly treated the channel incision being driven by a stream power rule hillslopes internal to each cell respond rapidly to incision of the local channel prompting their treatment in the model as steady state forms i present analytic expressions for the local relief which is controlled by the local channel incision rate and rate constants associated with the dominant hillslope process either diffusive or landsliding initial conditions include any preexisting topography and a lightly etched a channel network the model is illustrated using the santa cruz mountains california i show that an earlier hypothesis that the santa cruz mountains are due to advection of topography past an uplift source related to a restraining bend in the san andreas fault remains viable when more specific geomorphic processes are considered the local tectonics are represented by crustal thickening necessitated by conservation of mass in the bend right lateral translation of crust with respect to the bend and flexural response to the distributed load the model produces reasonable geographic and statistical distributions of topography using rates of tectonic and geomorphic processes that are within the range of those measured locally this argues that this restraining bend is a long lived feature of the san andreas fault the slightly higher crest and shorter length of the southern santa cruz mountains may potentially be explained by slightly faster slip of the crust west of the fault with respect to the bend than of that to the east deflection is minimal beneath these relatively small mountains disallowing the interesting feedback of erosionally driven uplift the predicted pattern of exhumation could be used to guide the search for fission track sites that could aid in the testing of the model

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

topographic fingerprints of bedrock landslides bedrock landslides in mountainous regions may be triggered by either storms or earthquakes the dominant mechanism in a region affects both landscape evolution and landslide hazard we describe a simple observational test to distinguish between storm and earthquake triggers based on a probabilistic measure of hillslope morphology in areas that are dominated by storm triggered landslides steep topographic slopes are concentrated on the lowermost parts of the hillslopes storm triggers act primarily on the hillslope toes and landslides preferentially remove material from those locations giving rise to inner gorges areas where most landslides are earthquake triggered have more uniform spatial distributions of steep topographic slopes because coseismic shaking causes failures at both ridge crests and hillslope toes earthquake triggered landslides lead to planar hillslopes and rare or absent inner gorges

erosion as a driving mechanism of intracontinental mountain growth in nature mountains can grow and remain as localized tectonic features over long periods of time 10 m y by contrast according to current knowledge of lithospheric rheology and neglecting surface processes any intracontinental range with a width that exceeds that which can be supported by the strength of the lithosphere should collapse within a few tens of millions of years for example assuming a quartz dominated crustal rheology the relief of a range initially 3 km high and 300 400 km wide is reduced by half in about 15 m y as a result of lateral spreading of its crustal root we suggest that surface processes might actually prevent such a subsurface collapse removal of material from topographic heights and deposition in the foreland oppose spreading of the crustal root and could eventually drive a net influx of material toward the orogeny we performed a set of numerical experiments in order to validate this hypothesis a section of a lithosphere with a brittle elasto ductile rheology initially loaded by a mountain range is submitted to horizontal shortening and to surface processes if erosion is intense material is removed more rapidly than it can be supplied by crustal thickening below the range and the topography is rapidly smoothed for example a feature 3 km high and 300 400 km wide is halved in height in about 15 m y for an erosion coefficient k 10 3 m 2 yr the erosion rate is of the order of a few 0 1 mm yr this regime might be called erosional collapse if erosion is not active enough the crustal root spreads out laterally and subsurface collapse occurs in the third intermediate regime removal of the material by erosion is dynamically compensated by isostatic rebound and inward flow in the lower crust so that the range can grow in this mountain growth regime the range evolves toward a characteristic graded shape that primarily depends on the erosion law the erosion rate may be high e g 0 5 0 9 mm yr close to the rate of tectonic uplift e g 0 7 1 1 mm yr and few times higher than the rate of topographic uplift 0 15 0 2 mm yr these experiments show that surface processes can favor localized crustal shortening and participate in the development of an intracontinental mountain surface processes must therefore be taken into account in the interpretation and modeling of long term deformation of continental lithosphere conversely the mechanical response of the lithosphere must be accounted for when large scale topographic features are interpreted and modeled in terms of geomorphologic processes

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

denudation processes and rates in the transverse ranges southern california erosional response of a transitional landscape to external and anthropogenic forcing 1 quantification of denudation in the transverse ranges of california permits reconstruction of spatial and temporal variations in erosion that represent both the response and evolution of interacting hillslopes and channels on the southern front of the san gabriel mountains observational records of the infilling of debris basins and dams define twentieth century landscape erosion rates averaging 1 6 and 0 9 mm yr 1 respectively although all major sediment transport occurs during intense winter storms debris production on hillslopes is greatly enhanced by recurrent fires consequently in this populated region anthropogenic fires have augmented the natural erosion rates we perform a global inversion to estimate the role of precipitation intensity burned areas and local slope on catchment denudation rates after subtracting the effects of anthropogenic fires we estimate landscape denudation rates given natural fire ignition rates increased fire during the past century has augmented sediment production in debris basins by an average of 60 and individual basins show increases up to 400 to identify the dominant hillslope erosion processes the volumetric contribution of landslides was estimated using repeat aerial photographs for the same time interval over which the debris basins have been operative between 1928 and 1973 landsliding produced only 10 of the sediment in debris basins even in the long term when infrequent but volumetrically important landslides occur bedrock landslides appear to contribute a maximum of 50 to the long term landscape denudation previous mapping of soil slippage in the san dimas experimental forest within the san gabriel mountains rice et al 1969 rice and foggin 1971 indicates that shallow landsliding is likely to be the dominant modern hillslope erosion process when compared to incision rates derived from a fluvial shear stress model and to exhumation rates based on low temperature thermochronological data blythe et al 2000 modern natural erosion rates are comparable to denudation rates since the pliocene comparisons of modern erosion rates suggest that debris production on hillslopes and first order channels is directly dependent on vegetation cover and precipitation intensity for higher order channels drainage areas 2 km2 only major storms convey the sediments down valley during the past century some temporal decorrelation occurs between small scale and large scale catchments because hillslope produced sediment is stored in second order and higher channels until major storms mobilize it thus the larger fluvial network damps the episodic fire induced hillslope sediment pulses that occur within small watersheds at a longer temporal scale however uplift and denudation may have been sustained sufficiently long in much of the san gabriel mountains for the topography to reach a macroscale steady state in contrast to rapidly eroding 2 mm yr 1 ranges for which quantification of bedrock landsliding will approximate the sediment flux the san gabriel mountains occupy a niche of intermediate rates 0 1 1 0 mm yr 1 in which a broad suite of hillslope processes including shallow seated and deep seated landslides debris flows and wet and dry ravel contribute to the sediment flux

miocene rifting in the los angeles basin evidence from the puente hills half graben volcanic rocks and p wave tomography formation of the puente hills half graben in the northeastern los angeles basin and eruption of the glendora and el modeno volcanics 16 14 ma help to define the timing of extension in the basin normal faulting on the proto whittier fault ca 14 ma established the puente hills half graben in which sedimentary strata accumulated between ca 14 and 10 ma and into which diabase sills intruded north south contraction began to invert the puente hills half graben ca 7 ma leading to formation of the puente hills anticline and the whittier fault our high resolution three dimensional p wave velocity model shows two anomalous higher velocity 6 63 km s bodies at depths between 9 and 18 km which we attribute to dioritic plutons named here for whittier narrows and el modeno the stocklike whittier narrows pluton could have been a source for the glendora volcanics and the diabase sills in the puente hills half graben the sill shaped el modeno pluton was a likely source for the el modeno volcanics the northwesterly alignment of the plutons may mark the location of the northeastern los angeles basin rift boundary which is associated with the clockwise rotation of the western transverse ranges three active faults the elysian park blind thrust the puente hills blind thrust and the whittier fault converge on the whittier narrows pluton which may have played a role in their location and segmentation

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

stochastic forcing of sediment routing and storage in channel networks the stochastic field of sediment supply to the channel network of a drainage basin depends on the large scale interactions among climatically driven processes such as forest fire and rainstorms topography channel network topology and basin scale during infrequent periods of intense erosion large volumes of colluvium are concentrated in parts of a channel network particularly near tributary junctions the rivers carry bed material and wash load downstream from these storage sites at different rates the bed material travels slowly creating transient patterns of sediment transport sediment storage and channel morphology along the channel network as the concentrations of bed material migrate along the network their waveforms can undergo changes by diffusion interference at tributary junctions and loss of mass through temporary sediment storage in fans and terraces and through particle abrasion which converts bed material to wash load we investigated how these processes might influence the sediment mass balance in channels of third and higher order in a 215 km 2 drainage basin within the oregon coast range over a simulated 3000 year period with a climate typical of the late holocene we used field measurements and a simulation model to illustrate interactions between the major controls on large scale processes functioning over long periods of time in complex drainage basins

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