skip to main content

Title: HyLands 1.0: a hybrid landscape evolution model to simulate the impact of landslides and landslide-derived sediment on landscape evolution
Abstract. Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) more » model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics. « less
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Geoscientific Model Development
Page Range or eLocation-ID:
3863 to 3886
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract In most landscape evolution models, extreme rainfall enhances river incision. In steep landscapes, however, these events trigger landslides that can buffer incision via increased sediment delivery and aggradation. We quantify landslide sediment aggradation and erosional buffering with a natural experiment in southern Taiwan where a northward gradient in tectonic activity drives increasing landscape steepness. We find that landscape response to extreme rainfall during the 2009 typhoon Morakot varied along this gradient, where steep areas experienced widespread channel sediment aggradation of >10 m and less steep areas did not noticeably aggrade. We model sediment export to estimate a sediment removal timeline and find that steep, tectonically active areas with the most aggradation may take centuries to resume bedrock incision. Expected sediment cover duration reflects tectonic uplift. We find that despite high stream power, sediment cover may keep steep channels from eroding bedrock for up to half of a given time period. This work highlights the importance of dynamic sediment cover in landscape evolution and suggests a mechanism by which erosional efficiency in tectonically active landscapes may decrease as landscape steepness increases.
  2. Interactions between deep Earth geodynamics and Earth surface processes are well documented at various scales, but many challenges remain in how inversion of a fluvially incised landscape should be interpreted in terms of long-term geodynamics or how deep Earth dynamics impact natural hazards. Here, we present results from geomorphic stream channel metrics and modeling of long profiles of streams draining the Tyrrhenian (northern) flank of Sicily (Italy), to assess the inferred, rapid, west-to-east horizontal translation of the Calabrian forearc. A detachment-limited stream power model-based determination of landscape response time and knickpoint migration provides an independent prediction for transient base level fall associated with the sweeping forearc over the past ~4 Ma. The model shows that two pulses of time-transgressive, west-to-east propagating base level fall occurred in the drainages of northern Sicily, where parallel north-flowing streams are arranged across the migrating path of the forearc. The long profile analysis indicates that the paired uplift pulses last ~1 Ma and are separated in time by ~1.5 Ma, consistent with the west-to-east passage first of the forearc high, followed by dynamic uplift in its wake due to sub-lithospheric mantle flow, as proposed in other plate boundary settings. The ongoing surficial response to thesemore »dynamics is represented by river incision, knickpoint migration, and drainage divide migration. Furthermore, these processes steepened the landscape, leading to an increase in active landsliding and contributing to the natural hazards in this region.« less
  3. Valla, Pierre (Ed.)
    Abstract Over the past few decades, tectonic geomorphology has been widely implemented to constrain spatial and temporal patterns of fault slip, especially where existing geologic or geodetic data are poor. We apply this practice along the eastern margin of Bull Mountain, Southwest Montana, where 15 transient channels are eroding into the flat, upstream relict landscape in response to an ongoing period of increased base level fall along the Western North Boulder fault. We aim to improve constraints on the spatial and temporal slip rates across the Western North Boulder fault zone by applying channel morphometrics, cosmogenic erosion rates, bedrock characteristics, and calibrated reproductions of the modern river profiles using a 1-dimensional stream power incision model that undergoes a change in the rate of base level fall. We perform over 104 base level fall simulations to explore a wide range of fault slip dynamics and stream power parameters. Our best fit simulations suggest that the Western North Boulder fault started as individual fault segments along the middle to southern regions of Bull Mountain that nucleated around 6.2 to 2.5 Ma, respectively. This was followed by the nucleation of fault segments in the northern region around 1.5 to 0.4 Ma. We recreate themore »evolution of the Western North Boulder fault to show that through time, these individual segments propagate at the fault tips and link together to span over 40 km, with a maximum slip of 462 m in the central portion of the fault. Fault slip rates range from 0.02 to 0.45 mm/yr along strike and are consistent with estimates for other active faults in the region. We find that the timing of fault initiation coincides well with the migration of the Yellowstone hotspot across the nearby Idaho-Montana border and thus attribute the initiation of extension to the crustal bulge from the migrating hotspot. Overall, we provide the first quantitative constraints on fault initiation and evolution of the Western North Boulder fault, perhaps the farthest north basin in the Northern Basin and Range province that such constraints exist. We show that river profiles are powerful tools for documenting the spatial and temporal patterns of normal fault evolution, especially where other geologic/geodetic methods are limited, proving to be a vital tool for accurate tectonic hazard assessments.« less
  4. Bedrock topography is a key boundary condition for ice sheet modeling, and determining changes in subglacial topography through time can provide insight into the timing of ice sheet development, the magnitude of glacial erosion, and the co-development of glaciers and glacial topography. West Antarctica hosts an unusually high geothermal gradient supported by hot, low-viscosity mantle which likely enhanced the lithospheric response to West Antarctic Ice Sheet (WAIS) cycles of growth and increased the sensitivity of thermochronometers to landscape evolution on million-year timescales. Thus, a valuable record of glacial landscape change might be recovered from apatite fission track [AFT 80-130°C range] and (U-Th)/He [AHe; 50-90°C] dating, provided that landscape evolution can be distinguished from tectonic signals, including the effects of faults. This study utilizes AFT-AHe thermochronology and thermo-kinematic Pecube modeling to investigate interactions between the hot geotherm, glacial erosion, and inferred crustal structures in the Ford Ranges and the DeVicq Glacier trough in western and central Marie Byrd Land (MBL), respectively. The Ford Ranges host glacial troughs (up to 3km relief) dissecting a low-relief erosional surface. Previous work suggests a majority of bedrock exhumation and cooling occurred at/by 80 Ma. However, new data hint at renewed exhumation linked to glacial incisionmore »since WAIS formation at 34 or 20 Ma. Prior (U-Th)/He zircon dates from exposures of crystalline bedrock span 90 – 67 Ma. New AHe bedrock dates are 41 to 26 Ma, while two glacial erratics (presumed to be eroded from walls or floor of glacial troughs) yielded AHe dates of 37 Ma and 16 Ma. Initial modeling results suggest a tectonic boundary between the Ford Ranges and Edward VII Peninsula separating regions with distinct exhumation histories. The boundary may cause differential WAIS incision at 34 or 20 Ma, a possibility being investigated with new models. The DeVicq Glacier trough (>3.5km relief) coincides with a prominent crustal lineament but lacks temperature-time information compared to other regions. The crustal structure may have accommodated motion between elevated central MBL and the subdued crust of the Ford Ranges. Here, owing to the lack of onshore non-volcanic bedrock exposure, we have employed AHe and AFT dating of glacial sediment marine core samples offshore of the DeVicq Glacier to investigate the timing and rates of exhumation of the bedrock carved by the DeVicq trough, with initial results revealing detrital AHe ages as young as 24 Ma. Our new Pecube models test a series of thermal, tectonic, and landscape evolution scenarios against a suite of thermochronologic data, allowing us to assess the timing of glacial incision and WAIS initiation in the Ford Ranges, and to seek evidence of an inferred tectonic boundary at DeVicq Trough. Modeling efforts will be aided by new AHe and AFT analyses from ongoing work. These models combine topographic, tectonic, thermal, and key thermochronologic datasets to produce new insight into the unique cryosphere-lithosphere interactions affecting landscape change in West Antarctica.« less
  5. Abstract. River avulsions are an important mechanism by which sediment is routed andemplaced in foreland basins. However, because avulsions occur infrequently,we lack observational data that might inform where, when, and why avulsionsoccur and these issues are instead often investigated by rule-basednumerical models. These models have historically simplified or neglected the effects of abandoned channels on avulsion dynamics, even though fluvialmegafans in foreland basins are characteristically covered in abandonedchannels. Here, we investigate the pervasiveness of abandoned channels onmodern fluvial megafan surfaces. Then, we present a physically basedcellular model that parameterizes interactions between a single avulsingriver and abandoned channels in a foreland basin setting. We investigate how abandoned channels affect avulsion setup, pathfinding, and landscapeevolution. We demonstrate and discuss how the processes of abandoned channel inheritance and transient knickpoint propagation post-avulsion serve to shortcut the time necessary to set up successive avulsions. Then, we address the idea that abandoned channels can both repel and attract future pathfinding flows under different conditions. By measuring the distance between the mountain front and each avulsion over long (106 to 107 years) timescales, we show that increasing abandoned channel repulsion serves to push avulsions farther from the mountain front, while increasing attraction pulls avulsions proximally. Abandoned channels domore »not persist forever, and we test possible channel healing scenarios (deposition-only, erosion-only, and far-field-directed) and show that only the final scenario achieves dynamic equilibrium without completely filling accommodation space. We also observe megafan growth occurring via ∼100 000-year cycles of lobe switching but only in our runs that employ deposition-only or erosion-only healing modes. Finally, we highlight opportunities for future field work and remote sensing efforts to inform our understanding of the role that floodplain topography, including abandoned channels, plays on avulsion dynamics.« less