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 streampower) 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.

Gravel-bedded streams are thought to have a characteristic meandering form bordered by a self-formed, fine-grained floodplain. This ideal guides a multibillion-dollar stream restoration industry. We have mapped and dated many of the deposits along mid-Atlantic streams that formed the basis for this widely accepted model. These data, as well as historical maps and records, show instead that before European settlement, the streams were small anabranching channels within extensive vegetated wetlands that accumulated little sediment but stored substantial organic carbon. Subsequently, 1 to 5 meters of slackwater sedimentation, behind tens of thousands of 17th- to 19th-century milldams, buried the presettlement wetlands with fine sediment. These findings show that most floodplains along mid-Atlantic streams are actually fill terraces, and historically incised channels are not natural archetypes for meandering streams.

No abstract

[1] Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second-order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well-constrained, probabilistic model of storm and flood occurrence, has first-order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment-limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high-frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.INDEX TERMS: 1815 Hydrology: Erosion and sedimentation; 1824 Hydrology: Geomorphology (1625); 1869 Hydrology: Stochastic processes; 8107 Tectonophysics: Continental neotectonics; KEYWORDS: tectonic geomorphology, erosion, bedrock channels, relief, thresholds, stochastic processes Citation: Snyder, N. P., K. X Whipple, G. E. Tucker, and D. J. Merritts, Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem,

Along three rivers at the Mendocino triple junction, northern California, strath, cut, and fill terraces have formed in response to tectonic and eustatic processes. Detailed surveying and radiometric dating at multiple sites indicate that lower reaches of the rivers are dominated by the effects of oscillating sea level, primarily aggradation and formation of fill terraces during sea level high stands, alternating with deep incision during low stands. A eustasy‐driven depositional wedge extends tens of kilometers upstream on all rivers (tapering to zero thickness). This distance is greater than expected from studies of the effects of check dams on much smaller streams elsewhere, due in part to the large size of these rivers. However, the change in gradient is nearly identical to other base level rise studies: the depositional gradient is about half that of the original channel. Middle to upper reaches of each river are dominated by the effects of long‐term uplift, primarily lateral and Vertical erosion and formation of steep, unpaired strath terraces exposed only upstream of the depositional wedge. Vertical incision at a rate similar to that of uplift has occurred even during the present sea level high stand along rivers with highest uplift rates. Strath terraces have steeper gradients than the modern channel bed and do not merge with marine terraces at the river mouth; consequently, they cannot be used to determine altitudes of sea level high stands. Strath formation is a continuous process of response to long‐term uplift, and its occurrence varies spatially along a river depending on stream power, and hence position, upstream. Strath terraces are found only along certain parts of a coastal stream: upstream of the aggradational effects of oscillating sea level, and far enough downstream that stream power is in excess of that needed to transport the prevailing sediment load. For a given size river, the greater the uplift rate, the greater the rate of vertical incision and, consequently, the less the likelihood of strath terrace formation and preservation.