An official website of the United States government
Despite considerable community effort, there is no general set of equations to model long‐term landscape evolution. In order to determine a suitable set of landscape evolution process laws for a site where postglacial erosion has incised valleys up to 50 m deep, we generate a set of alternative models and perform a multimodel analysis. The most basic model we consider includes stream power channel incision, uniform lithology, hillslope transport by linear diffusion, and surface‐water discharge proportional to drainage area. We systematically add one, two, or three elements of complexity to this model from one of four categories: hillslope processes, channel processes, surface hydrology, and representation of geologic materials. We apply methods of formal model analysis to the 37 alternative models. The global Method of Morris sensitivity analysis method is used to identify model input parameters that most and least strongly influence model outputs. Only a few parameters are identified as important, and this finding is consistent across two alternative model outputs: one based on a collection of topographic metrics and one that uses an objective function based on a topographic difference. Parameters that control channel erosion are consistently important, while hillslope diffusivity is important for only select model outputs. Uncertainty in initial and boundary conditions is associated with low sensitivity. Sensitivity analysis provides insight to model dynamics and is a critical step in using model analysis for mechanistic hypothesis testing in landscape evolution theory.
more »
« less
Inverting Topography for Landscape Evolution Model Process Representation: 2. Calibration and Validation
We present a multimodel analysis for mechanistic hypothesis testing in landscape evolution theory. The study site is a watershed with well‐constrained initial and boundary conditions in which a river network locally incised 50 m over the last 13 ka. We calibrate and validate a set of 37 landscape evolution models designed to hierarchically test elements of complexity from four categories: hillslope processes, channel processes, surface hydrology, and representation of geologic materials. Comparison of each model to a base model, which uses stream power channel incision, uniform lithology, hillslope transport by linear diffusion, and surface water discharge proportional to drainage area, serves as a formal test of which elements of complexity improve model performance. Model fit is assessed using an objective function based on a direct difference between observed and simulated modern topography. A hybrid optimization scheme identifies optimal parameters and uncertainty. Multimodel analysis determines which elements of complexity improve simulation performance. Validation tests which model improvements persist when models are applied to an independent watershed. The three most important model elements are (1) spatial variation in lithology (differentiation between shale and glacial till), (2) a fluvial erosion threshold, and (3) a nonlinear relationship between slope and hillslope sediment flux. Due to nonlinear interactions between model elements, some process representations (e.g., nonlinear hillslopes) only become important when paired with the inclusion of other processes (e.g., erosion thresholds). This emphasizes the need for caution in identifying the minimally sufficient process set. Our approach provides a general framework for hypothesis testing in landscape evolution.
more »
« less
- PAR ID:
- 10446794
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 125
- Issue:
- 7
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Models of landscape evolution provide insight into the geomorphic history 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 4 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 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 process elements: hillslope processes, surface-water hydrology, erosion by flowing water, and material properties. The 28 model programs are a systematic subset of the 2048 possible numerical models associated with 11 binary choices. Each binary choice is related to one of these four elements – for example, the use of linear or nonlinear hillslope diffusion. Terrainbento is an extensible framework: base classes that treat the elements common to all numerical models (such as input/output and boundary conditions) make it possible to create a new numerical model without reinventing these common methods. Terrainbento is built on top of the Landlab framework such that new Landlab components directly support the creation of new terrainbento model programs. Terrainbento is fully documented, has 100 % unit test coverage including numerical comparison with 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 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 numerical models.more » « less
-
Abstract Hillslope topographic change in response to climate and climate change is a key aspect of landscape evolution. The impact of short‐duration rainstorms on hillslope evolution in arid regions is persistently questioned but often not directly examined in landscape evolution studies, which are commonly based on mean climate proxies. This study focuses on hillslope surface processes responding to rainstorms in the driest regions of Earth. We present a numerical model for arid, rocky hillslopes with lithology of a softer rock layer capped by a cliff‐forming resistant layer. By representing the combined action of bedrock and clast weathering, cliff‐debris ravel, and runoff‐driven erosion, the model can reproduce commonly observed cliff‐profile morphology. Numerical experiments with a fixed base level were used to test hillslope response to cliff‐debris grain size, rainstorm intensities, and alternation between rainstorm patterns. The persistence of vertical cliffs and the pattern of sediment sorting depend on rainstorm intensities and the size of cliff debris. Numerical experiments confirm that these two variables could have driven the landscape in the Negev Desert (Israel) toward an observed spatial contrast in topographic form over the past 105–106 years. For a given total storm rain depth, short‐duration higher‐intensity rainstorms are more erosive, resulting in greater cliff retreat distances relative to longer, low‐intensity storms. Temporal alternation between rainstorm regimes produces hillslope profiles similar to those previously attributed to Quaternary oscillations in the mean climate. We suggest that arid hillslopes may undergo considerable geomorphic transitions solely by alternating intra‐storm patterns regardless of rainfall amounts.more » « less
-
Abstract Erosional perturbations from changes in climate or tectonics are recorded in the profiles of bedrock rivers, but these signals can be challenging to unravel in settings with non‐uniform lithology. In layered rocks, the surface lithology at a given location varies through time as erosion exposes different layers of rock. Recent modeling studies have used the Stream Power Model (SPM) to highlight complex variations in erosion rates that arise in bedrock rivers incising through layered rocks. However, these studies do not capture the effects of coarse sediment cover on channel evolution. We use the “Stream Power with Alluvium Conservation and Entrainment” (SPACE) model to explore how sediment cover influences landscape evolution and modulates the topographic expression of erodibility contrasts in horizontally layered rocks. We simulate river evolution through alternating layers of hard and soft rock over million‐year timescales with a constant and uniform uplift rate. Compared to the SPM, model runs with sediment cover have systematically higher channel steepness values in soft rock layers and lower channel steepness values in hard rock layers. As more sediment accumulates, the contrast in steepness between the two rock types decreases. Effective bedrock erodibilities back‐calculated assuming the SPM are strongly influenced by sediment cover. We also find that sediment cover can significantly increase total relief and timescales of adjustment toward landscape‐averaged steady‐state topography and erosion rates.more » « less
-
Abstract Features of landscape morphology—including slope, curvature, and drainage dissection—are important controls on runoff generation in upland landscapes. Over long timescales, runoff plays an essential role in shaping these same features through surface erosion. This feedback between erosion and runoff generation suggests that modeling long‐term landscape evolution together with dynamic runoff generation could provide insight into hydrological function. Here we examine the emergence of variable source area runoff generation in a new coupled hydro‐geomorphic model that accounts for water balance partitioning between surface flow, subsurface flow, and evapotranspiration as landscapes evolve over millions of years. We derive a minimal set of dimensionless numbers that provide insight into how hydrologic and geomorphic parameters together affect landscapes. Across the parameter space we investigated, model results collapsed to a single inverse relationship between the dimensionless relief and the ratio of catchment quickflow to discharge. Furthermore, we found an inverse relationship between the Hillslope number, which describes topographic relief relative to aquifer thickness, and the proportion of the landscape that was variably saturated. While the model generally produces fluvial topography visually similar to simpler landscape evolution models, certain parameter combinations produce wide valley bottom wetlands and non‐dendritic, trellis‐like drainage networks, which may reflect real conditions in some landscapes where aquifer gradients become decoupled from topography. With these results, we demonstrate the power of hydro‐geomorphic models for generating new insights into hydrological processes, and also suggest that subsurface hydrology may be integral for modeling aspects of long‐term landscape evolution.more » « less
