Also known as
Model type Single
Model part of larger framework
Incorporated models or components:
Spatial dimensions 1D
Spatial extent Reach-Scale
Model domain Terrestrial, Hydrology
One-line model description This model is designed to simulate longitudinal profiles with headward advancing headcuts.
Extended model description This model is designed to simulate longitudinal profiles with headward advancing headcuts. This model simulates gully erosion on the centennial-scale given information such as average rainfall and infiltration rates. The modeler also specifies a headcut erosion rate and or a rule for headcut retreat (either discharge-dependent or height-dependent retreat).

Gully, elevation profile, Matlab, arid alluvial channel,

Name Francis Rengers
Type of contact Model developer
Institute / Organization University of Colorado
Postal address 1 UCB 399
Postal address 2
Town / City Boulder
Postal code 80309-0399
State Colorado
Country United States
Email address

Supported platforms
Mac OS, Windows
Other platform
Programming language


Other program language
Code optimized Single Processor
Multiple processors implemented
Nr of distributed processors
Nr of shared processors
Start year development
Does model development still take place? No
If above answer is no, provide end year model development 2012
Code development status As is, no updates are provided
When did you indicate the 'code development status'? 2020
Model availability As code
Source code availability
(Or provide future intension)
Through CSDMS repository
Source web address
Source csdms web address
Program license type GPL v2
Program license type other
Memory requirements
Typical run time Depends on the settings. Can be several hours.

Describe input parameters The input file is a text file and users are required to input:

Time (in model years) dT (the time step in fractions of a year) tauc (Critical Shear stress for portions of the channel that are vegetated in Pascals) taucWepp (Critical Shear stress for portions of the channel that are soil in Pascals) lenzone (the length of channel that is bare soil in Meters) Pmmphr (Rainfall to be used for erosion in Millimeters per Hour) tval (this is the number of loop iterations before a profile is saved as output) Immphr (Infiltration to be used for erosion in Millimeters per Hour)

One additional input: One must supply the input length of the channel as a matlab data array called xcell.mat

Input format ASCII
Other input format
Describe output parameters The output are elevation values that represent the gully channel profile.
Output format ASCII
Other output format
Pre-processing software needed? No
Describe pre-processing software
Post-processing software needed? No
Describe post-processing software
Visualization software needed? Yes
If above answer is yes Matlab
Other visualization software In order to plot the output data, it is probably easiest in Matlab

Describe processes represented by the model This code will erode cells according to a shear stress and also deposit sediment based on the concentration of sediment in a modeled water column. Additionally it has a headcut that migrates upstream and as the headcut erodes it deposits sediment downstream that the model must erode.
Describe key physical parameters and equations manning's n for vegetation

manning's n for soil critical shear stress for vegetation critical shear stress for soil rainfall infiltration

Describe length scale and resolution constraints This works for channels of hundreds of meters with 1 m cells
Describe time scale and resolution constraints It is typically run for times no less than 25 years and no greater than 700 years.
Describe any numerical limitations and issues Numerical instabilities occur if the time step is too large.

Describe available calibration data sets
Upload calibration data sets if available:
Describe available test data sets
Upload test data sets if available:
Describe ideal data for testing

Do you have current or future plans for collaborating with other researchers?
Is there a manual available? No
Upload manual if available:
Model website if any
Model forum / discussion board

This part will be filled out by CSDMS staff

OpenMI compliant No but possible
BMI compliant No but possible
WMT component No but possible
PyMT component
Is this a data component
Can be coupled with:
Model info
Nr. of publications: 1
Total citations: 40
h-index: 1
m-quotient: 0.11
Qrcode GullyErosionProfiler1D.png
Link to this page




Nr. of publications: 1
Total citations: 40
h-index: 1
m-quotient: 0.11

Featured publication(s)YearModel describedType of ReferenceCitations
Rengers, Francis K.; Tucker, G. E.; 2014. Analysis and modeling of gully headcut dynamics, North American high plains. Journal of Geophysical Research: Earth Surface, 119, 983–1003. 10.1002/2013JF002962
(View/edit entry)
2014 GullyErosionProfiler1D
Model overview 40
See more publications of GullyErosionProfiler1D



Input Files

Output Files