TY - JOUR
T1 - Keeping the edge
T2 - A numerical method that avoids knickpoint smearing when solving the stream power law
AU - Campforts, Benjamin
AU - Govers, Gerard
PY - 2015/7/1
Y1 - 2015/7/1
N2 - ©2015. American Geophysical Union. All Rights Reserved.The stream power equation is commonly used to model river incision into bedrock. Although specific conditions allow an analytical approach, finite difference methods (FDMs) are most frequently used to solve this equation. FDMs inevitably suffer from numerical smearing which may affect their suitability for transient river incision modeling. We propose the use of a finite volume method (FVM) which is total variation diminishing (TVD) to simulate river incision in a more accurate way. The TVD-FVM is designed to simulate sharp discontinuities, making it very suitable to simulate river incision pulses. We show that the TVD-FVM is much better capable of preserving propagating knickpoints than FDMs, using Niagara Falls as an example. Comparison of numerical results obtained using the TVD-FVM with analytical solutions shows a very good agreement. Furthermore, the uncertainty associated with parameter calibration is dramatically reduced when the TVD-FVM is applied. The high accuracy of the TVD-FDM allows correct simulation of transient incision waves as a consequence of older uplift pulses. This implies that the TVD-FVM is much more suitable than FDMs to reconstruct regional uplift histories from current river profile morphology and to simulate river incision processes in general.
AB - ©2015. American Geophysical Union. All Rights Reserved.The stream power equation is commonly used to model river incision into bedrock. Although specific conditions allow an analytical approach, finite difference methods (FDMs) are most frequently used to solve this equation. FDMs inevitably suffer from numerical smearing which may affect their suitability for transient river incision modeling. We propose the use of a finite volume method (FVM) which is total variation diminishing (TVD) to simulate river incision in a more accurate way. The TVD-FVM is designed to simulate sharp discontinuities, making it very suitable to simulate river incision pulses. We show that the TVD-FVM is much better capable of preserving propagating knickpoints than FDMs, using Niagara Falls as an example. Comparison of numerical results obtained using the TVD-FVM with analytical solutions shows a very good agreement. Furthermore, the uncertainty associated with parameter calibration is dramatically reduced when the TVD-FVM is applied. The high accuracy of the TVD-FDM allows correct simulation of transient incision waves as a consequence of older uplift pulses. This implies that the TVD-FVM is much more suitable than FDMs to reconstruct regional uplift histories from current river profile morphology and to simulate river incision processes in general.
UR - http://www.scopus.com/inward/record.url?scp=84939264393&partnerID=8YFLogxK
U2 - 10.1002/2014JF003376
DO - 10.1002/2014JF003376
M3 - Article
SN - 2169-9003
VL - 120
SP - 1189
EP - 1205
JO - Journal of Geophysical Research: Earth Surface
JF - Journal of Geophysical Research: Earth Surface
IS - 7
ER -