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 Hamiltonian1 model of rock slope erosion and retreat2. 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 scheme3 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 mechanics4: 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
-
Stark, C.P., & Stark, G.J., 2022. The direction of landscape erosion. Earth Surface Dynamics, 10: 383-419. ↩
-
Howard, A.D., & Selby, M.J., 2009. Rock Slopes. In: Parsons, A.J., Abrahams, A.D. (eds). Geomorphology of Desert Environments. Springer, Dordrecht. ↩
-
Osher, S., & Fedkiw, R., 2003. Level Set Methods and Dynamic Implicit Surfaces. Springer-Verlag New York, Inc. See page 50. ↩
-
Holm, D.D., 2011. Geometric Mechanics. Part I: Dynamics and Symmetry (2nd Edition) ↩
1 maintainer
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Filter files by name, interpreter, ABI, and platform.
If you're not sure about the file name format, learn more about wheel file names.
Copy a direct link to the current filters
File details
Details for the file erosionfront-0.1.23.tar.gz.
File metadata
- Download URL: erosionfront-0.1.23.tar.gz
- Upload date:
- Size: 59.3 kB
- Tags: Source
- Uploaded using Trusted Publishing? Yes
- Uploaded via: twine/6.1.0 CPython/3.13.12
File hashes
| Algorithm | Hash digest | |
|---|---|---|
| SHA256 |
bc25cbf794cdd65f6ff945da3920c177269013015f4df36341a5ff7effc2a4e9
|
|
| MD5 |
c11e5ac848319c3ad8e3f93c61072f63
|
|
| BLAKE2b-256 |
452a5cbe6ecc87dc8c380777127ccbe23d16779482ffa6c37e69c0be4fff90f3
|
Provenance
The following attestation bundles were made for erosionfront-0.1.23.tar.gz:
Publisher:
pypi-publish.yml on cstarkjp/ErosionFront
-
Statement:
-
Statement type:
https://in-toto.io/Statement/v1 -
Predicate type:
https://docs.pypi.org/attestations/publish/v1 -
Subject name:
erosionfront-0.1.23.tar.gz -
Subject digest:
bc25cbf794cdd65f6ff945da3920c177269013015f4df36341a5ff7effc2a4e9 - Sigstore transparency entry: 1822478773
- Sigstore integration time:
-
Permalink:
cstarkjp/ErosionFront@df9c9ea74eb1bb23d4995d3ad81c57b7d27cd899 -
Branch / Tag:
refs/heads/PyPI - Owner: https://github.com/cstarkjp
-
Access:
public
-
Token Issuer:
https://token.actions.githubusercontent.com -
Runner Environment:
github-hosted -
Publication workflow:
pypi-publish.yml@df9c9ea74eb1bb23d4995d3ad81c57b7d27cd899 -
Trigger Event:
push
-
Statement type:
File details
Details for the file erosionfront-0.1.23-py3-none-any.whl.
File metadata
- Download URL: erosionfront-0.1.23-py3-none-any.whl
- Upload date:
- Size: 71.3 kB
- Tags: Python 3
- Uploaded using Trusted Publishing? Yes
- Uploaded via: twine/6.1.0 CPython/3.13.12
File hashes
| Algorithm | Hash digest | |
|---|---|---|
| SHA256 |
a98c9f07ea0efa96ada87c822a13ce7f2a2d32b23f46ccacbeb106cd29700bc6
|
|
| MD5 |
0a60a7b5b42f6651f932d841aaecbc8b
|
|
| BLAKE2b-256 |
55cf06842f6c4d00704611a74a598addb07953537bd83df561b076a324ab99e1
|
Provenance
The following attestation bundles were made for erosionfront-0.1.23-py3-none-any.whl:
Publisher:
pypi-publish.yml on cstarkjp/ErosionFront
-
Statement:
-
Statement type:
https://in-toto.io/Statement/v1 -
Predicate type:
https://docs.pypi.org/attestations/publish/v1 -
Subject name:
erosionfront-0.1.23-py3-none-any.whl -
Subject digest:
a98c9f07ea0efa96ada87c822a13ce7f2a2d32b23f46ccacbeb106cd29700bc6 - Sigstore transparency entry: 1822478777
- Sigstore integration time:
-
Permalink:
cstarkjp/ErosionFront@df9c9ea74eb1bb23d4995d3ad81c57b7d27cd899 -
Branch / Tag:
refs/heads/PyPI - Owner: https://github.com/cstarkjp
-
Access:
public
-
Token Issuer:
https://token.actions.githubusercontent.com -
Runner Environment:
github-hosted -
Publication workflow:
pypi-publish.yml@df9c9ea74eb1bb23d4995d3ad81c57b7d27cd899 -
Trigger Event:
push
-
Statement type:
