Criteria and tools for determining drainage divide stability

© 2018 Elsevier B.V. Watersheds are the fundamental organizing units in landscapes and thus the controls on drainage divide location and mobility are an essential facet of landscape evolution. Additionally, many common topographic analyses fundamentally assume that river network topology and divide locations are largely static, allowing channel profile form to be interpreted in terms of spatio-temporal patterns of rock uplift rate relative to base level, climate, or rock properties. Recently however, it has been suggested that drainage divides are more mobile than previously thought and that divide mobility, and resulting changes in drainage area, could potentially confound interpretations of river profiles. Ultimately, reliable metrics are needed to diagnose the mobility of divides as part of routine landscape analyses. One such recently proposed metric is cross-divide contrasts in χ a proxy for steady-state channel elevation, but cross-divide contrasts in a number of topographic metrics show promise. Here we use a series of landscape evolution simulations in which we induce divide mobility under different conditions to test the utility of a suite of topographic metrics of divide mobility and for comparison with natural examples in the eastern Greater Caucasus Mountains, the Kars Volcanic Plateau, and the western San Bernadino Mountains. Specifically, we test cross-divide contrasts in mean gradient, mean local relief, channel bed elevation, and χ all measured at, or averaged upstream of, a reference drainage area. Our results highlight that cross-divide contrasts in χ only faithfully reflect current divide mobility when uplift, rock erodibility, climate, and catchment outlet elevation are uniform across both river networks on either side of the divide, otherwise a χ-anomaly only indicates a possible future divide instability. The other metrics appear to be more reliable representations of current divide motion, but in natural landscapes, only cross-divide contrasts in mean gradient and local relief appear to consistently provide useful information. Multiple divide metrics should be considered simultaneously and across-divide values of all metrics examined quantitatively as visual assessment is not sufficiently reliable in many cases. We provide a series of Matlab tools built using TopoToolbox to facilitate routine analysis

Amplified erosion above waterfalls and oversteepened bedrock reaches

None of the conventional bedrock erosion laws can predict incision immediately upslope of a waterfall lip where the flow is accelerating toward a freefall. Considering the expected increase in flow velocity and shear stress at the lip of a waterfall, we determine erosion amplification at a waterfall lip as [equation], where [equation] is the erosion rate at the upstream end of the flow acceleration zone above a waterfall, Fr is the Froude number at this setting, and n ranges between 0.5–1.7. This amplification expression suggests that erosion at the lip could be as much as 2–5 times higher relative to erosion at a normal setting with identical hydraulic geometry. Utilizing this erosion amplification expression in numerical simulations, we demonstrate its impact on reach-scale morphology above waterfalls. Amplified erosion at the lip of a waterfall can trigger the formation of an oversteepened reach whose length is longer than the flow acceleration zone, provided incision wave velocity (Vi) at the upstream edge of the flow acceleration zone is higher than the retreat velocity of the waterfall face. Such an oversteepened reach is expected to be more pronounced when Vi increases with increasing slope. The simulations also suggest that oversteepening can eventually lead to steady state gradients adjacent to a waterfall lip provided Vi decreases with increasing slope. Flow acceleration above waterfalls can thus account, at least partially, for prevalent oversteepened bedrock reaches above waterfalls. Using the cosmogenic isotope Cl-36, we demonstrate that incision wave velocity upstream of a waterfall at the Dead Sea western escarpment is probably high enough for freefall-induced oversteepening to be feasible

Short communication: The Topographic Analysis Kit (TAK) for TopoToolbox

Quantitative analysis of digital topographic data is an increasingly important part of many studies in the geosciences. Initially, performing these analyses was a niche endeavor, requiring detailed domain knowledge and programming skills, but increasingly broad, flexible, open-source code bases have been developed to increasingly democratize topographic analysis. However, many of these analyses still require specific computing environments and/or moderate levels of knowledge of both the relevant programming language and the correct way to take these fundamental building blocks and conduct an efficient and effective topographic analysis. To partially address this, we have written the Topographic Analysis Kit (TAK), which leverages the power of one of these open code bases, TopoToolbox, to build a series of high-level topographic analysis tools to perform a variety of common topographic analyses. These analyses include the generation of maps of normalized channel steepness, or χ, and selection and statistical analysis of populations of watersheds. No programming skills or advanced mastery of MATLAB is required for effective use of TAK. In addition – to expand the utility of TAK along with the primary functions, which like the underlying TopoToolbox functions require MATLAB and several proprietary toolboxes to run – we provide compiled versions of these functions that use the free MATLAB Runtime Environment for users who do not have institutional access to MATLAB or all of the required toolboxes.</p

Glacial erosion and relief production in the Eastern Sierra

Abstract The proposal that climate change can drive the uplift of mountain summits hinges on the requirement that glacial erosion significantly enhances the relief of a previously fluvially sculpted mountain range. We have tested this hypothesis through a systematic investigation of neighbouring glaciated and nonglaciated drainage basins on the eastern side of the Sierra Nevada, CA. We present a simple, objective method for investigating the relief structure of a drainage basin, which shows noticeable differences in the spatial distribution of relief between nonglaciated and glaciated basins. Glaciated basins on the eastern side of the Sierra Nevada have only f80 m greater mean geophysical relief than nonglaciated basins. This &apos;&apos;extra&apos;&apos; relief, though, is attributable principally to the larger size of the glaciated basins, as geophysical relief generally increases with basin size. The glaciers on this side of the range were only responsible for relief production if they substantially increased headward erosion rates into low relief topography, such as an elevated plateau, and thus enlarged previously fluvial basins. We carried out a preliminary morphometric analysis to elucidate the importance of this effect and found that the glaciers of the eastern Sierra Nevada may have eroded headward at considerably faster rates than rivers, but only when they were not obstructed from doing so by either competing larger glaciers in adjacent valleys or transfluent ice at the head of the basin. Our results also suggest that, in temperate regions, alpine glaciers are capable of eroding downward at faster rates than rivers above the equilibrium line altitude (ELA). Although we can rule out significant peak uplift in response to local relief production, in the special case of the Sierra Nevada the concentration of mass removal above the ELA could have contributed to flexural uplift at the edge of a tilting block.

Can springs cut canyons into rock?

Amphitheater-headed valleys on Earth and Mars are often assumed to result from erosion by emerging spring water (i.e., seepage erosion or groundwater sapping) rather than by surface runoff. The origin of such valleys has implications for landscape evolution on Earth and the hydrologic cycle and associated potential for life on other planets. In this paper we explore the evidence for seepage erosion in bedrock to address whether valley morphology can be used as a diagnostic indicator of seepage erosion. Seepage erosion is an important process in loose sediment where hydraulic forces cause grain detachment, often resulting in amphitheater-headed valleys. However, the extension of these processes to resistant rock is uncertain. In sedimentary rocks, groundwater might control the shape and rate of valley formation. It is possible, however, that seepage plays only a secondary role to runoff processes. This seems likely in basaltic valleys on Earth, where little evidence exists for seepage erosion. Since the ability of seepage to erode bedrock valleys remains unclear and because many amphitheater-headed valleys were probably carved by other processes, seepage erosion should not be inferred based solely on valley form

Constraints on the tectonic and landscape evolution of the Bhutan Himalaya from thermochronometry

The observed geomorphology and calculated thermal histories of the Bhutan Himalaya provide an excellent platform to test ideas regarding the influence of tectonics and climate on the evolution of a convergentmountain range. However, little consensus has been reached regarding the late Cenozoic history of the Bhutan Himalaya. Some researchers have argued that observed geologic relationships show slowing deformation rates, such that the range is decaying from a geomorphic perspective, while others see the range as growing and steepening. We suggest that a better understanding is possible through the integrated interpretation of geomorphic and thermochronometric data from the comparison of predictions from models of landscape evolution and thermal-kinematic models of orogenic systems. New thermochronometric data throughout Bhutan aremost consistent with a significant decrease in erosion rates, from2 to 3 km/Ma down to 0.1–0.3 km/Ma, around 6–4Ma. We interpret this pattern as a decrease in rock uplift rates due to the activation of contractional structures of the Shillong Plateau, an uplifted region approximately 100 km south of Bhutan. However, low-relief, fluvial landscapes throughout the Bhutanese hinterland record a late pulse of surface uplift likely due to a recent increase in rock uplift rates. Constraints from our youngest thermochronometers suggest that this increase in rock uplift and surface uplift occurred within the last 1.75Ma. These results imply that the dynamics of the Bhutan Himalaya and Shillong Plateau have been linked during the late Cenozoic, with structural elements of both regions active in variable ways and times over that interval

Channel dynamics, sediment transport, and the slope of alluvial fans: experimental study, The

A B S T R A C T We present the results of an experimental study of alluvial fan sedimentation under conditions of constant inflow water discharge Q w , sediment supply Q so , and median grain size D. The study was designed to complement and test a recently formulated model of alluvial fan sedimentation and to emphasize the interactions between, and controls on, flow channelization and equilibrium fan slope. Flow channelization and fan sedimentation were studied under conditions of nearly steady, uniform aggradation. Steady conditions were achieved by imposing a steadily rising base level, just in balance with the average sediment aggradation rate. Experimental inflows covered a wide range of conditions, allowing examination of the effects of variations in Q w , Q so , and D in both bedload-and suspension-dominated environments. Experimental results were most consistent with an expanding-flow channel model. Key experimental findings successfully and quantitatively predicted by the expanding-flow theory include: (1) straight to slightly concave radial profiles of bedload-dominated fans; (2) distinctly convexo-concave profiles of suspension-dominated fans; (3) a strong, inverse relationship between Q w and fan slope; (4) a strong, but secondary, relationship between Q so and fan slope; and (5) near-independence of D and fan slope so long as transport stage is high and bedload transport dominant. However, potential scale effects in the experiments arose from reduced flow Reynolds numbers and incorrect geometric scaling of channel widths; no confident conclusions regarding the debate over the relative importance of &apos;&apos;sheetfloods&apos;&apos; and braided channel flows can be drawn from the experimental data. Extrapolation to field scale is best accomplished through appropriate application of the theoretical model herein confirmed against experimental data

Lattice-Boltzmann Method for Geophysical Plastic Flows

We explore possible applications of the Lattice-Boltzmann Method for the simulation of geophysical flows. This fluid solver, while successful in other fields, is still rarely used for geotechnical applications. We show how the standard method can be modified to represent free-surface realization of mudflows, debris flows, and in general any plastic flow, through the implementation of a Bingham constitutive model. The chapter is completed by an example of a full-scale simulation of a plastic fluid flowing down an inclined channel and depositing on a flat surface. An application is given, where the fluid interacts with a vertical obstacle in the channel.Comment: in W. Wu, R.I. Borja (Edts.) Recent advances in modelling landslides and debris flow, Springer Series in Geomechanics and Geoengineering (2014), ISBN 978-3-319-11052-3, pp. 131-14

Late Cenozoic evolution of the eastern margin of the Tibetan Plateau: Inferences from ^(40)Ar/^(39)Ar and (U-Th)/He thermochronology

High topography in central Asia is perhaps the most fundamental expression of the Cenozoic Indo-Asian collision, yet an understanding of the timing and rates of development of the Tibetan Plateau remains elusive. Here we investigate the Cenozoic thermal histories of rocks along the eastern margin of the plateau adjacent to the Sichuan Basin in an effort to determine when the steep topographic escarpment that characterizes this margin developed. Temperature-time paths inferred from ^(40)Ar/^(39)Ar thermochronology of biotite, multiple diffusion domain modeling of alkali feldspar ^(40)Ar release spectra, and (U-Th)/He thermochronology of zircon and apatite imply that rocks at the present-day topographic front of the plateau underwent slow cooling (30°–50°C/m.y.) coincident with exhumation from inferred depths of ∼8–10 km, at denudation rates of 1–2 mm/yr. Samples from the interior of the plateau continued to cool relatively slowly during the same time period (∼3°C/m.y.), suggesting limited exhumation (1–2 km). However, these samples record a slight increase in cooling rate (from <1 to ∼3°C/m.y.) at some time during the middle Tertiary; the tectonic significance of this change remains uncertain. Regardless, late Cenozoic denudation in this region appears to have been markedly heterogeneous, with the highest rates of exhumation focused at the topographic front of the plateau margin. We infer that the onset of rapid cooling at the plateau margin reflects the erosional response to the development of regionally significant topographic gradients between the plateau and the stable Sichuan Basin and thus marks the onset of deformation related to the development of the Tibetan Plateau in this region. The present margin of the plateau adjacent to and north of the Sichuan Basin is apparently no older than the late Miocene or early Pliocene (∼5–12 Ma)