Also known as
Model type Tool
Model part of larger framework
Note on status model
Date note status model
Incorporated models or components:
Spatial dimensions 1D
Spatial extent Landscape-Scale
Model domain Terrestrial
One-line model description The module calculates crevasse splay morphology and water discharge outflow of a crevasse splay.
Extended model description The module is designed to calculate morphological changes and water discharge outflow of a crevasse splay that is triggered by a preset flood event and evolves afterwards. The inputs for "mainCS.m" should be daily water discharge and sediment flux series of the trunk channel upstream the crevasse splay. The outputs will be daily series for the cross-sectional parameters of the crevasse splay, and daily water discharge series of the trunk channel downstream the crevasse splay. One limitation of the present version is it only calculates the expanding and healing of a crevasse splay, while ignores the possible morphological change (demise or revival) of the trunk channel downstream the crevasse splay. Another limitation is the codes are originally written for the Lower Yellow River(a suspended load dominated river) for the purpose of calculating sediment budget in the Lower Yellow over a long timescale, say as long as hundreds years, so the present module can not be applied to other alluvial rivers without modifying those lines related to channel geometry, bankfull discharge and bank erosion(deposition).

coastal evolution,

Name Yunzhen Chen
Type of contact Model developer
Institute / Organization Nanjing University
Postal address 1 22 Hankou Road
Postal address 2
Town / City Nanjing
Postal code 210093
Country China
Email address
Phone +8613851402809

Supported platforms
Unix, Linux, 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 2009
Does model development still take place? Yes
If above answer is no, provide end year model development
Code development status
When did you indicate the 'code development status'?
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 seconds

Describe input parameters daily water discharge series;daily sediment flux series;

averaged channel cross-sectional depth, averaged channel cross-sectional width; floodplain width; manning coefficients of the channel and floodplain; longitudinal channel slope; Channel bed's super-elevation above the floodplain where sedimentation rate is close to 0; M-coefficient for erosion rate for the bottom of crevasse splay; M-coefficient for erosion rate for the two side slopes of crevasse splay; critical velocity for erosion; critical velocity for deposition; width of dike at the root; cross valley slope; settling velocity of suspended load in the channel.

Input format ASCII
Other input format
Describe output parameters averaged daily water discharge to lower reach of crevasse splay;

averaged daily crevasse splay depth; averaged daily crevasse splay width

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? No
If above answer is yes
Other visualization software

Describe processes represented by the model As a crevasse splay evolves, the slope of its outflow should be no less than the slope of lower channel; and the bottom elevation of a crevasse splay should be no lower than the elevation of lowest point of channel bed, so the bottom elevation of the lowest point that a crevasse splay is able to cut down is max(hs, Zcsb).

A ratio of Q above the bottom of crevasse splay can be distributed to outflow of crevasse splay. After flow parameters for the outflow of crevasse splay are calculated, the erosion (deposition) rate of crevasse splay can be calculated, thus the morphology of crevasse splay can be updated.

When the crevasse splay has not yet cut down to the lowest point max(hs,Zcsb), it can be both widened and deepened. When the crevasse splay has cut down to the lowest point max(hs,Zcsb), it can only be widened or silted vertically.

Describe key physical parameters and equations Key physical parameters:
  • Q: water discharge upstream crevasse splay;
  • Qcs: outflow discharge of crevasse splay;
  • Qabove: the water discharge above the bottom of crevasse splay;
  • rq: the discharge ratio of Qcs and Qabove;
  • hs: channel belt's super-elevation (the elevation of lowest point of channel bed);
  • Zcs: bottom elevation of crevasse splay;
  • Bcs: width of crevasse splay;
  • Hcs: flow depth of crevasse splay;
  • Vcs: flow velocity of crevasse splay;
  • jcs: slope of the outflow of crevasse splay;
  • Zcsb: bottom elevation of a crevasse splay whose flow slope is equal to the channel slope j;

Key physical equations:

  • Zcs<=max(hs,Zcsb);
  • rq=(1.55-1.45*Fi)*Bcs/wc+0.16*(1-2*Fi), in which Fi is the Fraud number for flow in the trunk channel, wc is width of the trunk channel;
  • Hcs=(nc*Qcs/sqrt(jcs)/Bcs)^(3/5);
  • Vcs=Qcs/Hcs/Bcs;
  • dE=M*(Vcs^2-ucre^2)/ucre^2*dt, where M is M-coefficient for erosion rate for crevasse slpay, ucre is critical velocity for erosion, dt is time step;
  • dD=Sv*(1-Vcs^2/ucrd^2)*ws/0.6*dt, where Sv is volume sediment concentration, ucrd is critical velocity for deposition, ws is settling velocity of suspended load, dt is time step.
Describe length scale and resolution constraints length scale: meters

resolution constraint: centimeter

Describe time scale and resolution constraints time scale: decades

resolution constraints: day

Describe any numerical limitations and issues

Describe available calibration data sets
Upload calibration data sets if available:
Describe available test data sets Q(m^3/s): daily water discharge series, 365 by 1 matrix

Qs(kg/s): daily sediment flux series, 365 by 1 matrix The trigger flood event is set at the 232th day

Upload test data sets if available: Media:Testdata.rar
Describe ideal data for testing As a test, you can use the attached test data Q and Qs, and run


Do you have current or future plans for collaborating with other researchers? I'm now writing some papers coauthored with James Syvitski, Irina Overeem and Albert Kettner, on historical channel avulsions and source-to-sink sediment budget models of the Yellow River basin.
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
DOI model 10.1594/IEDA/100105
For model version 0.0.1
Year version submitted 2011
Link to file
Can be coupled with:
Model info
Yunzhen Chen

Nr. of publications: 2
Total citations: 0
h-index: --"--" is not a number.
m-quotient: 0
Qrcode CrevasseFlow.png
Link to this page




Nr. of publications: 2
Total citations: 0
h-index: --"--" is not a number.
m-quotient: 0

Featured publication(s)YearModel describedType of ReferenceCitations
Chen, Y.; Overeem, I.; Syvitski, J.P.M.; Gao, S.; Kettner, A.J.; 2011. Controls of levee breaches on the Lower Yellow River during the years 1550-1855.. RCEM Publ., , 617-633.
(View/edit entry)
2011 CrevasseFlow

Model overview

Chen, Y.; 2011. CrevasseFlow, version 0.0.1.. , , 10.1594/IEDA/100105
(View/edit entry)
2011 CrevasseFlow

Source code ref.

See more publications of CrevasseFlow



Input Files

Output Files