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.
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- NSF-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
-
more » « less
-
more » « less
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 The potentially important influence of climate change on landscape evolution and on critical zone processes is not sufficiently understood. The relative contribution of hydro‐climatic factors on hillslope erosion rates may significantly vary with topography at the watershed scale. The objective of this study is to quantify the hydro‐geomorphic behavior of two contrasting landscapes in response to different climate change scenarios in the Luquillo Critical Zone Observatory, a site of particular geomorphological interest, in terms of hillslope erosion and rainfall‐triggered landslides. We investigate the extent to which hillslope erosion and landslide occurrence remain relatively invariant with future hydro‐climatic perturbations. The adjacent Mameyes and Icacos watersheds are studied, which are underlain by contrasting lithologies. A high resolution coupled hydro‐geomorphic model based on tRIBS (Triangulated Irregular Network‐based Real‐time Integrated Basin Simulator) is used. Observations of landslide activity and hillslope erosion are used to evaluate the model performance. The process‐based model quantifies feedbacks among different hydrologic processes, landslide occurrence, and topsoil erosion and deposition. Simulations suggest that the propensity for landslide occurrence in the Luquillo Mountains is controlled by tropical storms, subsurface water flow, and by non‐climatic factors, and is expected to remain significant through 2099. The Icacos watershed, which is underlain by quartz diorite, is dominated by relatively large landslides. The relative frequency of smaller landslides is higher at the Mameyes watershed, which is underlain by volcaniclastic rock. While projections of precipitation decrease at the study site may lead to moderate decline in hillslope erosion rates, the simulated erosional potential of the two diverse landscapes likely remains significant. © 2018 John Wiley & Sons, Ltd.
-
more » « less
Abstract Landscapes after wildfire commonly experience accelerated hillslope erosion, which often contributes to the mobilization and volume of debris flows. However, quantitative studies of the erosion and its relationship to rainfall, runoff, and landscape characteristics have been limited to a narrow range of physiographic conditions. We estimated the volume and delivery rate of slurry (a water‐sediment mixture) supplied to stream channels during a post‐wildfire rainstorm that generated large debris flows in six catchments above Montecito, CA, in 2018. We mapped the distribution of rills and measured their cross‐sectional geometries to quantify the influences of runoff, lithology, and hillslope characteristics on the sediment volumes released by rill erosion, and we scaled the results up to the 19.5 km2of burned hillslopes in the source catchments. We computed the likely rate of surface runoff during the rainstorm and developed an empirical model for the evolution of a representative hillslope‐spanning rill to illustrate the magnitude and speed of the erosion process. Rilling was the dominant form of erosion across the hillslopes of the source catchments, and the rapid evacuation and mixing of water and sediment during rill formation supplied a slurry with high solids concentrations to stream channels. Colluvium on shale formations was more continuous, finer‐grained, and probably less permeable than colluvium on sandstones, and these differences affected the extent and dimensions of rills. As a result, shale hillslopes were the dominant source of slurry to the debris flows and supplied over twice as much slurry per unit burn area as sandstones.
