1,569 research outputs found
Recommended from our members
A lattice grain model of hillslope evolution
Thispaperdescribesandexploresanewcontinuous-timestochasticcellularautomatonmodelofhill- slope evolution. The Grain Hill model provides a computational framework with which to study slope forms that arise from stochastic disturbance and rock weathering events. The model operates on a hexagonal lattice, with cell states representing fluid, rock, and grain aggregates that are either stationary or in a state of motion in one of the six cardinal lattice directions. Cells representing near-surface soil material undergo stochastic disturbance events, in which initially stationary material is put into motion. Net downslope transport emerges from the greater likelihood for disturbed material to move downhill than to move uphill. Cells representing rock undergo stochas- tic weathering events in which the rock is converted into regolith. The model can reproduce a range of common slope forms, from fully soil mantled to rocky or partially mantled, and from convex-upward to planar shapes. An optional additional state represents large blocks that cannot be displaced upward by disturbance events. With the addition of this state, the model captures the morphology of hogbacks, scarps, and similar features. In its simplest form, the model has only three process parameters, which represent disturbance frequency, characteris- tic disturbance depth, and base-level lowering rate, respectively. Incorporating physical weathering of rock adds one additional parameter, representing the characteristic rock weathering rate. These parameters are not arbitrary but rather have a direct link with corresponding parameters in continuum theory. Comparison between observed and modeled slope forms demonstrates that the model can reproduce both the shape and scale of real hillslope profiles. Model experiments highlight the importance of regolith cover fraction in governing both the downslope mass transport rate and the rate of physical weathering. Equilibrium rocky hillslope profiles are possible even when the rate of base-level lowering exceeds the nominal bare-rock weathering rate, because increases in both slope gradient and roughness can allow for rock weathering rates that are greater than the flat-surface maximum. Examples of transient relaxation of steep, rocky slopes predict the formation of a regolith-mantled pediment that migrates headward through time while maintaining a sharp slope break
Modeling the shape and evolution of normal-fault facets
Facets formed along the footwalls of active normal-fault blocks display a variety of longitudinal profile forms, with variations in gradient, shape, degree of soil cover, and presence or absence of a slope break at the fault trace. We show that a two-dimensional, process-oriented cellular automaton model of facet profile evolution can account for the observed morphologic diversity. The model uses two dimensionless parameters to represent fault slip, progressive rock weathering, and downslope colluvial-soil transport driven by gravity and stochastic disturbance events. The parameters represent rock weathering and soil disturbance rates, respectively, scaled by fault slip rate; both can be derived from field-estimated rate coefficients. In the model's transport-limited regime, slope gradient depends on the ratio of disturbance to slip rate, with a maximum that represents the angle of repose for colluvium. In this regime, facet evolution is consistent with nonlinear diffusion models of soil-mantled hillslope evolution. Under the weathering-limited regime, bedrock becomes partly exposed but microtopography helps trap some colluvium even when facet gradient exceeds the threshold angle. Whereas the model predicts a continuous gradient from footwall to colluvial wedge under transport-limited behavior, fully weathering-limited facets tend to develop a slope break between footwall and basal colluvium as a result of reduced transport efficiency on the rocky footwall slope. To the extent that the model provides a reasonable analogy for natural facets, its behavior suggests that facet profile morphology can provide useful constraints on relative potential rates of rock weathering, soil disturbance, and fault slip
Recommended from our members
Developing and exploring a theory for the lateral erosion of bedrock channels for use in landscape evolution models
Understanding how a bedrock river erodes its banks laterally is a frontier in geomorphology. Theories for the vertical incision of bedrock channels are widely implemented in the current generation of landscape evolution models. However, in general existing models do not seek to implement the lateral migration of bedrock channel walls. This is problematic, as modeling geomorphic processes such as terrace formation and hillslopeāchannel coupling depends on the accurate simulation of valley widening. We have developed and implemented a theory for the lateral migration of bedrock channel walls in a catchment-scale landscape evolution model. Two model formulations are presented, one representing the slow process of widening a bedrock canyon and the other representing undercutting, slumping, and rapid downstream sediment transport that occurs in softer bedrock. Model experiments were run with a range of values for bedrock erodibility and tendency towards transport- or detachment-limited behavior and varying magnitudes of sediment flux and water discharge in order to determine the role that each plays in the development of wide bedrock valleys. The results show that this simple, physics-based theory for the lateral erosion of bedrock channels produces bedrock valleys that are many times wider than the grid discretization scale. This theory for the lateral erosion of bedrock channel walls and the numerical implementation of the theory in a catchment-scale landscape evolution model is a significant first step towards understanding the factors that control the rates and spatial extent of wide bedrock valleys
Recommended from our members
Statistical Classification of SelfāOrganized Snow Surfaces
Windāswept snow selfāorganizes into bedforms. These bedforms affect local and global energy fluxes but have not been incorporated into Earth system models because the conditions governing their development are not well understood. To address this difficulty, we created statistical classifiers, drawn from 736 hr of timeālapse footage in the Colorado Front Range, that predict bedform presence as a function of wind speed and time since snowfall. These classifiers provide the first quantitative predictions of bedform and sastrugi presence in varying weather conditions. We find that the likelihood that a snow surface is covered by bedforms increases with time since snowfall and with wind speed and that the likelihood that a surface is covered by sastrugi increases with time and with the highest wind speeds. Our observations will be useful to Earth system modelers and represent a new step toward understanding selfāorganized processes that ornament 8% of the surface of the planet
Recommended from our members
VariableāThreshold Behavior in Rivers Arising From HillslopeāDerived Blocks
Geomorphologists often rely on simple models of river channel incision for predicting rates of landscape evolution and channel response to perturbations, as well as extracting climatic and tectonic signals from river longitudinal profiles. Recent work has shown that large, hillslopeāderived blocks delivered to rivers may noticeably alter the form and evolution of river profiles from the behavior predicted by the most common models. Here we use a 1āD model of river reach erosion and hillslope block delivery to explore the conditions under which block delivery strongly influences channel evolution. We use global sensitivity analysis to understand which model parameters most strongly affect the channel longitudinal profile. We explore the effects of blocks on the relationship between erosion rate and channel gradient, and on the erosion rateāchannel steepness exponent Ļ, and find that block effects result in highly variable slope and Ļ over the range of erosion rates and climatic conditions (discharge mean and variability) tested. The influence of blocks on erosion rateāslope scaling may be approximated by a piecewise model: The erosion threshold imposed by blocks scales linearly with erosion rate when blocks are infrequently mobile and remains constant when blocks are frequently mobile. We explore the implications of this variableāthreshold model for the erosion rateāchannel steepness relationship and find that erosion rateādependent thresholds imposed by hillslopeāderived blocks cause significant departures from previous models but may be consistent with existing field data sets. Our work has implications for landscape evolution modeling and the inversion of channel profiles for forcing information
Recommended from our members
The Babelizer: language interoperability for model; coupling in the geosciences
The babelizer is a Python utility that generates code to import libraries from other languages into Python. Target libraries must expose a Basic Model Interface (BMI) (Hutton et al., 2020; Peckham et al., 2013) and be written in C, C++, or Fortran, although the babelizer is extendable, so other languages can be added in the future. The babelizer provides a streamlined mechanism for bringing scientific models into a common language where they can communicate with one another as components of an integrated model.</p
Recommended from our members
CellLab-CTS 2015: continuous-time stochastic cellular automaton modeling using Landlab
CellLab-CTS 2015 is a Python-language software library for creating two-dimensional, continuous-time stochastic (CTS) cellular automaton models. The model domain consists of a set of grid nodes, with each node assigned an integer state code that represents its condition or composition. Adjacent pairs of nodes may undergo transitions to different states, according to a user-defined average transition rate. A model is created by writing a Python code that defines the possible states, the transitions, and the rates of those transitions. The code instantiates, initializes, and runs one of four object classes that represent different types of CTS models. CellLab-CTS provides the option of using either square or hexagonal grid cells. The software provides the ability to treat particular grid-node states as moving particles, and to track their position over time. Grid nodes may also be assigned user-defined properties, which the user can update after each transition through the use of a callback function. As a component of the Landlab modeling framework, CellLab-CTS models take advantage of a suite of Landlab's tools and capabilities, such as support for standardized input and output
Recommended from our members
The Landlab v1.0 OverlandFlow component: a Python tool for computing shallow-water flow across watersheds
Representation of flowing water in landscape evolution models (LEMs) is often simplified compared to hydrodynamic models, as LEMs make assumptions reducing physical complexity in favor of computational efficiency. The Landlab modeling framework can be used to bridge the divide between complex runoff models and more traditional LEMs, creating a new type of framework not commonly used in the geomorphology or hydrology communities. Landlab is a Python-language library that includes tools and process components that can be used to create models of Earth-surface dynamics over a range of temporal and spatial scales. The Landlab OverlandFlow component is based on a simplified inertial approximation of the shallow water equations, following the solution of de Almeida et al. (2012). This explicit two-dimensional hydrodynamic algorithm simulates a flood wave across a model domain, where water discharge and flow depth are calculated at all locations within a structured (raster) grid. Here we illustrate how the OverlandFlow component contained within Landlab can be applied as a simplified event-based runoff model and how to couple the runoff model with an incision model operating on decadal timescales. Examples of flow routing on both real and synthetic landscapes are shown. Hydrographs from a single storm at multiple locations in the Spring Creek watershed, Colorado, USA, are illustrated, along with a map of shear stress applied on the land surface by flowing water. The OverlandFlow component can also be coupled with the Landlab DetachmentLtdErosion component to illustrate how the nonsteady flow routing regime impacts incision across a watershed. The hydrograph and incision results are compared to simulations driven by steady-state runoff. Results from the coupled runoff and incision model indicate that runoff dynamics can impact landscape relief and channel concavity, suggesting that on landscape evolution timescales, the OverlandFlow model may lead to differences in simulated topography in comparison with traditional methods. The exploratory test cases described within demonstrate how the OverlandFlow component can be used in both hydrologic and geomorphic applications
Recommended from our members
The evolution of snow bedforms in the Colorado Front Range and the processes that shape them
When wind blows over dry snow, the snow surface self-organizes into bedforms such as dunes, ripples, snow waves, and sastrugi. These bedforms govern the interaction between wind, heat, and the snowpack, but thus far they have attracted few scientific studies. We present the first time-lapse documentation of snow bedform movement and evolution, as part of a series of detailed observations of snow bedform movement in the Colorado Front Range. We show examples of the movement of snow ripples, snow waves, barchan dunes, snow steps, and sastrugi. We also introduce a previously undocumented bedform: the stealth dune. These observations show that (1) snow dunes accelerate minute-by-minute in response to gusts, (2) sastrugi and snow steps present steep edges to the wind and migrate downwind as those edges erode, (3) snow waves and dunes deposit layers of cohesive snow in their wake, and (4) bedforms evolve along complex cyclic trajectories. These observations provide the basis for new conceptual models of bedform evolution, based on the relative fluxes of snowfall, aeolian transport, erosion, and snow sintering across and into the surface. We find that many snow bedforms are generated by complex interactions between these processes. The prototypical example is the snow wave, in which deposition, sintering, and erosion occur in transverse stripes across the snowscape.</p
- ā¦