Search

Constraining Earth's sediment mass balance over geologic time requires a quantitative understanding of how landscapes respond to transient tectonic perturbations. However, the mechanisms by which bedrock lithology governs landscape response remain poorly understood. Rock type influences the size of sediment delivered to river channels, which controls how efficiently rivers respond to tectonic forcing. The Mendocino triple junction region of northern California, USA, is one landscape in which large boulders, delivered by hillslope failures to channels, may alter the pace of landscape response to a pulse of rock uplift. Boulders frequently delivered by earthflows in one lithology, the Franciscan mélange, have been hypothesized to steepen channels and slow river response to rock uplift, helping to preserve high-elevation, low-relief topography. Channels in other units (the Coastal Belt and the Franciscan schist) may experience little or no erosion inhibition due to boulder delivery. Here we investigate spatial patterns in channel steepness, an indicator of erosion resistance, and how it varies between mélange and non-mélange channels. We then ask whether lithologically controlled boulder delivery to rivers is a possible cause of steepness variations. We find that mélange channels are steeper than Coastal Belt channels but not steeper than schist channels. Though channels in all units steepen with increasing proximity to mapped hillslope failures, absolute steepness values near failures are much higher (~2x) in the mélange and schist than in Coastal Belt units. This could reflect reduced rock erodibility or increased erosion rates in the mélange and schist, or disproportionate steepening due to enhanced boulder delivery by hillslope failures in those units. To investigate the possible influence of lithology-dependent boulder delivery, we map boulders at failure toes in the three units. We find that boulder size, frequency, and concentration are greatest in mélange channels and that Coastal Belt channels have the lowest concentrations. Using our field data to parameterize a mathematical model for channel slope response to boulder delivery, we find that the modeled influence of boulders in the mélange could be strong enough to account for some observed differences in channel steepness between lithologies. At the landscape scale, we lack the data to fully disentangle boulder-induced steepening from that due to spatially varying erosion rates and in situ rock erodibility. However, our boulder mapping and modeling results suggest that lithology-dependent boulder delivery to channels could retard landscape adjustment to tectonic forcing in the mélange and potentially also in the schist. Boulder delivery may modulate landscape response to tectonics and help preserve high-elevation, low-relief topography at the Mendocino triple junction and elsewhere. [ABSTRACT FROM AUTHOR]

Geologists frequently debate the origin of iconic river canyons, as well as the extent to which they record climatic and tectonic signals. Fluvial and hillslope processes work in concert to control canyon evolution; rivers both set the boundary conditions for adjoining hillslopes and respond to delivery of hillslope-derived sediment. But, what happens when canyon walls deliver boulders that are too large for a river to carry? River canyons commonly host large blocks of rock derived from resistant hillslope lithologies. Blocks have recently been shown to control the shapes of hillslopes and channels by inhibiting sediment transport and bedrock erosion. Here, we developed the first process-based model for canyon evolution that incorporates the roles of blocks in both hillslope and channel processes. Our model reveals that two-way negative channel-hillslope feedbacks driven by block delivery to the river result in characteristic plan-view and cross-sectional river canyon forms. Internal negative feedbacks strongly reduce the rate at which erosional signals pass through landscapes, leading to persistent local unsteadiness even under steady tectonic and climatic forcing. Surprisingly, while the presence of blocks in the channel initially slows incision rates, the subsequent removal of blocks from the oversteepened channel substantially increases incision rates. This interplay between channel and hillslope dynamics results in highly variable long-term erosion rates. These autogenic channel-hillslope dynamics can mask external signals, such as changes in rock uplift rate, complicating the interpretation of landscape morphology and erosion histories. [ABSTRACT FROM AUTHOR]

Models of landscape evolution provide insight into the development of specific field areas, create testable predictions of landform development, demonstrate the consequences of current geomorphic process theory, and spark imagination through hypothetical scenarios. While the last four decades have brought the proliferation of many alternative formulations for the redistribution of mass by Earth surface processes, relatively few studies have systematically compared and tested these alternative equations. We present a new Python modeling package, terrainbento 1.0, that enables multi-model comparison, sensitivity analysis, and calibration of Earth surface process models. terrainbento provides a set of 28 model programs that implement alternative transport laws related to four model elements: hillslope processes, surface-water hydrology, erosion by flowing water, and material properties. The 28 model programs stem from 13 binary choices related to one of these four elements -- for example, the use of linear or non-linear hillslope diffusion. terrainbento is an extensible framework: model base classes that treat the elements common to all models (such as input/output and boundary conditions) make it possible to create a new model without re-inventing these common methods. terrainbento is built on top of the Landlab framework, such that new Landlab components directly support the creation of new terrainbento models. terrainbento is fully documented, has 100% unit test coverage including numerical comparison with % all available analytical solutions for process models, and continuous integration testing. We support future users and developers with introductory Jupyter notebooks and a template for creating new terrainbento model programs. In this paper, we describe the package structure, process model theory, and software implementation of terrainbento. Finally, we illustrate the utility of terrainbento with a benchmark example highlighting the differences in steady state topography between five different process models. [ABSTRACT FROM AUTHOR]

Quantifying off-fault deformation (OFD) rates on geomorphic time scales (10²-105 yr) along strike-slip faults is critical for resolving discrepancies between geologic and geodetic slip-rate estimates, improving knowledge of seismic hazard, and understanding the influence of tectonic motion on landscapes. Quantifying OFD over these time scales is challenging without displacement markers such as offset terraces or geologic contacts. We present a landscape evolution model coupled with distributed lateral tectonic shear to show how drainage basins sheared by lateral tectonic motion can reveal OFD rates. The model shows that OFD rate can control the orientation of drainage basin topography: the faster the OFD rate, the greater the deflection of drainage basins toward a fault-parallel orientation. We apply the model to the southern San Andreas fault near the Mecca Hills (southern California, USA), where drainage basins change in orientation with proximity to the fault. Comparison of observed and modeled topography suggests that the OFD rate in the Mecca Hills follows an exponential-like spatial pattern with a maximum rate nearest the fault of 3.5 ± 1.5 mm/yr, which decays to ~0 at ~600 m distance from the fault. This rate is applicable since the initiation of differential rock uplift in the Mecca Hills ca. 760 ka. Our results suggest that OFD in this 800 m study area may be as high as 10% of total plate motion. This example demonstrates that curved drainage basins may be used to estimate OFD rates along strike-slip faults. [ABSTRACT FROM AUTHOR]