The emergent geometry of rock slopes

Summary: Simulation tools in support of a geomorphic, non-convex Hamiltonian theory of rock ramp-cliff retreat and the emergent geometry of Richter-type slopes.

Abstract

An iconic image of the American West is the desert mesa: a steep cliff, rising above a ramp-like rock slope, capped by a flat bench. This famous landform has long been assumed to develop where strong rock overlies weak, and where rockfall debris suppresses ramp erosion.

Such an explanation cannot be true in general, however, because the archetypal geometry can arise even in uniform bedrock with no talus armouring. Here we argue instead that the ramp-cliff shape is an emergent property.

Theoretical evidence comes from a simple model of scarp retreat whose combined rates of weathering and surface-normal erosion are written as a slowly varying function of gradient.

Model analysis and simulation, using geometric mechanics and level sets, reveal the sharp break in slope to form automatically as a transient shock solution of a non-convex Hamilton-Jacobi equation. Strong erodibility contrasts are not needed to explain this behaviour, but when present they lock the landform into its classic shape and allow it to persist long-term.

Comparison of differential cliff recession in geologically homogeneous versus heterogeneous bedrock at sites in Utah and Colorado confirms our hypothesis.

Level-set solution

The purpose of the Python code presented here is to derive, analyze, and numerically solve a geomorphic Hamiltonian¹ model of rock slope erosion and retreat². The code is provided as a Python library package and associated Jupyter notebooks (e.g., here and here).

Numerical solution of the model Hamilton-Jacobi equation is achieved with a level-set scheme³ that employs Lax-Friedrichs finite differencing to obtain stable viscosity solutions for a non-convex Hamiltonian. The level-set code is custom implemented in Python.

Model analysis is performed using some tools from geometric mechanics⁴: having converted the rock-slope erosion model into geomorphic Hamiltonian $\mathcal{H}(\mathbf{r}, \mathbf{p})$ form, this Hamiltonian is then used to derive Hamilton's ray tracing equations $(\partial_{\mathbf{p}}\mathcal{H}, -\partial_{\mathbf{r}}\mathcal{H})$ and the co-metric of rock slope erosion tensor $g^{ij} = \partial_{ij}\mathcal{H}$; these properties are then probed to understand model stability, notably to place bounds on the non-convexity of $\mathcal{H}$ and to identify critical angles.

References