skip to main content


Title: r.sim.terrain 1.0: a landscape evolution model with dynamic hydrology
Abstract. While there are numerical landscape evolution modelsthat simulate how steady-state flows of water and sedimentreshape topography over long periods of time,r.sim.terrain is the first tosimulate short-term topographic changefor both steady-state and dynamic flow regimesacross a range of spatial scales.This free and open-sourceGeographic Information Systems (GIS)-based topographic evolution modeluses empirical models for soil erosionand a physics-based modelfor shallow overland water flow and soil erosionto compute short-term topographic change.This model uses either a steady-stateor unsteady representation of overland flowto simulate how overland sediment mass flows reshape topographyfor a range of hydrologic soil erosion regimesbased on topographic, land cover, soil, and rainfall parameters.As demonstrated by a case studyfor the Patterson Branch subwatershedon the Fort Bragg military installation in North Carolina,r.sim.terrain simulates the development offine-scale morphological features includingephemeral gullies, rills, and hillslopes.Applications include land management, erosion control,landscape planning, and landscape restoration.  more » « less
Award ID(s):
1737563
NSF-PAR ID:
10126153
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
12
Issue:
7
ISSN:
1991-9603
Page Range / eLocation ID:
2837 to 2854
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The hydrologic dynamics and geomorphic evolution of watersheds are intimately coupled—runoff generation and water storage are controlled by topography and properties of the surface and subsurface, while also affecting the evolution of those properties over geologic time. However, the large disparity between their timescales has made it difficult to examine interdependent controls on emergent hydrogeomorphic properties, such as hillslope length, drainage density, and extent of surface saturation. In this study, we develop a new model coupling hydrology and landscape evolution to explore how runoff generation affects long‐term catchment evolution, and analyze numerical results using a nondimensional scaling framework. We focus on hydrologic processes dominating in humid climates where storm runoff primarily arises from shallow subsurface flow and from precipitation on saturated areas. The model solves hydraulic groundwater equations to predict the water‐table elevation given prescribed, constant groundwater recharge. Water in excess of the subsurface capacity for transport becomes overland flow, which generates shear stress on the surface and may detach and transport sediment. This affects the landscape form that in turn affects runoff generation. We show that (a) four dimensionless parameters describe the possible steady state landscapes that coevolve under steady recharge; (b) hillslope length increases with increasing transmissivity relative to the recharge rate; (c) three topographic metrics—steepness index, Laplacian curvature, and topographic index—together provide a basis for interpreting landscapes that have coevolved with runoff generated via shallow subsurface flow. Finally we discuss the possibilities and limitations for quantitative comparisons between the model results and real landscapes.

     
    more » « less
  2. Abstract

    Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the soil‐air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re‐orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes towards an ‘ideal hill’. In this steady state, reaction‐front valves partition water volumes into shallow and deep flowpaths. These flows dissolve low‐ and high‐solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.

     
    more » « less
  3. Abstract

    Steep landscapes evolve largely by debris flows, in addition to fluvial and hillslope processes. Abundant field observations document that debris flows incise valley bottoms and transport substantial sediment volumes, yet their contributions to steepland morphology remain uncertain. This has, in turn, limited the development of debris‐flow incision rate formulations that produce morphology consistent with natural landscapes. In many landscapes, including the San Gabriel Mountains (SGM), California, steady‐state fluvial channel longitudinal profiles are concave‐up and exhibit a power‐law relationship between channel slope and drainage area. At low drainage areas, however, valley slopes become nearly constant. These topographic forms result in a characteristically curved slope‐area signature in log‐log space. Here, we use a one‐dimensional landform evolution model that incorporates debris‐flow erosion to reproduce the relationship between this curved slope‐area signature and erosion rate in the SGM. Topographic analysis indicates that the drainage area at which steepland valleys transition to fluvial channels correlates with measured erosion rates in the SGM, and our model results reproduce these relationships. Further, the model only produces realistic valley profiles when parameters that dictate the relationship between debris‐flow erosion, valley‐bottom slope, and debris‐flow depth are within a narrow range. This result helps place constraints on the mathematical form of a debris‐flow incision law. Finally, modeled fluvial incision outpaces debris‐flow erosion at drainage areas less than those at which valleys morphologically transition from near‐invariant slopes to concave profiles. This result emphasizes the critical role of debris‐flow incision for setting steepland form, even as fluvial incision becomes the dominant incisional process.

     
    more » « less
  4. Abstract The quantification of rates for the competing forces of tectonic uplift and erosion has important implications for understanding topographic evolution. Here, we quantify the complex interplay between tectonic uplift, topographic development, and erosion recorded in the hanging walls of several active reverse faults in the Ventura basin, southern California, USA. We use cosmogenic 26Al/10Be isochron burial dating and 10Be surface exposure dating to construct a basin-wide geochronology, which includes burial dating of the Saugus Formation: an important, but poorly dated, regional Quaternary strain marker. Our ages for the top of the exposed Saugus Formation range from 0.36 +0.18/-0.22 Ma to 1.06 +0.23/-0.26 Ma, and our burial ages near the base of shallow marine deposits, which underlie the Saugus Formation, increase eastward from 0.60 +0.05/-0.06 Ma to 3.30 +0.30/-0.41 Ma. Our geochronology is used to calculate rapid long-term reverse fault slip rates of 8.6–12.6 mm yr–1 since ca. 1.0 Ma for the San Cayetano fault and 1.3–3.0 mm yr–1 since ca. 1.0 Ma for the Oak Ridge fault, which are both broadly consistent with contemporary reverse slip rates derived from mechanical models driven by global positioning system (GPS) data. We also calculate terrestrial cosmogenic nuclide (TCN)-derived, catchment-averaged erosion rates that range from 0.05–1.14 mm yr–1 and discuss the applicability of TCN-derived, catchment-averaged erosion rates in rapidly uplifting, landslide-prone landscapes. We compare patterns in erosion rates and tectonic rates to fluvial response times and geomorphic landscape parameters to show that in young, rapidly uplifting mountain belts, catchments may attain a quasi-steady-state on timescales of <105 years even if catchment-averaged erosion rates are still adjusting to tectonic forcing. 
    more » « less
  5. Abstract

    Desertification processes pose a global environmental threat, impacting 61 × 106 km2of the terrestrial land area. Changes in overland flow patterns and consequent rainwater redistribution in drylands present a potential pathway to desertification, because vegetation often relies on water inputs from runoff to sustain growth under insufficient rainfall conditions. Of particular importance are the very shallow overland flows that redistribute water, nutrients, and biological matter within arid landscapes. However, characterizing overland flow patterns remains challenging, due to their very shallow depths, their distributed nature, and the poor understanding of how these flows interact with the underlying rough soil surface. This paper describes how coupling thermal images of shallow overland flows with light detection and ranging (LiDAR) scanning of the underlying soil surface in 1 m2experimental trays allows spatial patterns of shallow overland flow to be quantified. Laboratory experiments were used to explore the behaviors of shallow overland flow as mean slope gradients and soil roughness were varied. The results show that these imaging techniques are able to capture differences in flow patterns arising across soil surfaces with varying slope, roughness, and spatial variation in infiltration properties. Several spatial indices characterizing overland flow patterns were found to correlate with runoff volume. The presence of high permeability soil patches substantially regulated overland flow. A next logical step would be to apply the thermal and LiDAR measurement techniques to the hillslope scale.

     
    more » « less