An official website of the United States government
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
Inverting Topography for Landscape Evolution Model Process Representation: 1. Conceptualization and Sensitivity Analysis
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
- PAR ID:
- 10446798
- 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
-
-
Temporal and spatial variations of tectonic rock uplift are generally thought to be the main controls on long-term erosion rates in various landscapes. However, rivers continuously lengthen and capture drainages in strike-slip fault systems due to ongoing motion across the fault, which can induce changes in landscape forms, drainage networks, and local erosion rates. Located along the restraining bend of the San Andreas Fault, the San Bernardino Mountains provide a suitable location for assessing the influence of topographic disequilibrium from perturbations by tectonic forcing and channel reorganization on measured erosion rates. In this study, we measured 17 new basin-averaged erosion rates using cosmogenic 10Be in river sands (hereafter, 10Be-derived erosion rates) and compiled 31 10Be-derived erosion rates from previous work. We quantify the degree of topographic disequilibrium using topographic analysis by examining hillslope and channel decoupling, the areal extent of pre-uplift surface, and drainage divide asymmetry across various landscapes. Similar to previous work, we find that erosion rates generally increase from north to south across the San Bernardino Mountains, reflecting a southward increase in tectonic activity. However, a comparison between 10Be-derived erosion rates and various topographic metrics in the southern San Bernardino Mountains suggests that the presence of transient landscape features such as relict topography and drainage-divide migration may explain local variations in 10Be-derived erosion rates. Our work shows that coupled analysis of erosion rates and topographic metrics provides tools for assessing the influence of tectonic uplift and channel reorganization on landscape evolution and 10Be-derived erosion rates in an evolving strike-slip restraining bend.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
-
Abstract Channel networks are important landscape features that transport water, sediment, and nutrients. Their emergence and evolution are controlled by the competition between hillslope and fluvial processes on landscapes. Investigating the geomorphic and topologic properties of these networks is crucial for quantifying the roles of processes in creating distinct patterns of channel networks and developing models to predict the network dynamics under changing environment. Here, we study the response of landscapes to changing climatic forcing via numerical‐modeling and the topographic analysis of natural landscapes. We use a physically‐based numerical landscape evolution model to investigate the channel network structure for varying hillslope and fluvial processes represented by different magnitudes of soil transport () and fluvial incision () coefficients. We show that landscapes with the same Péclet number (defined as the ratio of the timescales of advective (fluvial) to diffusive (hillslope) processes) and thus the same characteristic length scale may exhibit different geomorphic and topologic characteristics. Specifically, increasingDandK(mimicking humid conditions) or decreasingDandK(mimicking dry conditions), while keeping the same Péclet number, results in distinct branching structures. These changes lead to an exponential decrease in relief under humid conditions and an increase under dry conditions. For smaller and combinations, higher number of branching channels is observed, whereas for larger and combinations, higher number of side‐branching channels is obtained. These results align with topographic analysis of natural landscapes, suggesting that varying climatic conditions imprint distinct signatures on the branching structure of channel networks.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
