RiverProAnalysis, an open-source set of Matlab scripts for river longitudinal profile analysis

Based on the analytical solutions to the steady-state and linear transient-state forms of the model equation, researchers developed a set of open-source Matlab scripts, RiverProAnalysis, to integrate the functions of drainage basin concavity calibration, channel steepness estimate, drainage divide stability determination, paleo- projection, and tectonic uplift history inversion.

The products of the tool set are in the format of image, text, and vector files, which not only can be used for visual analysis, but also be read by geographic information system software.

To illustrate the ability of the tool set, the authors took examples of two transient drainage catchments in the northern margin of the Taiyuan Basin. They analyzed the long profiles of the trunk streams, identified two generations of knickpoints, and estimated the minimum depth of river incision.

They combined all the trunk and tributary channels together to model the catchment-wide uplift history and found moderate increases in the uplift rates since the Middle Pliocene and rapid acceleration since the late Quaternary.

The results were consistent with the sedimentary records in the adjacent basin. By comparing the chi value, slope, and topographic relief of both sides of the catchment divide, they concluded the stability of the drainage divide.

The left panel shows the chi map of the two sampled catchments. The right upper panel is the Statistical histogram of the river head slope on both sides of the catchment boundary. The right low panel shows the statistical histogram of the local relief of river heads. Credit: Science China Press

The flowchart tells the functions of the program and the input data and output results of each function. Credit: Science China Press

(a) The Xilingjing catchment: Elevation misfit between the forwarded and actual channel chi-plots, under different Chi value intervals. (b) The chi-plots of all the rivers in the Xilingjing catchment and the forwarded channel chi-plot are based on chi interval of 1.0 m (black thick line). (c) Dimensionless uplift history. (d) The modeled uplift history based on different erodibility values. (e)‒(h) Uplift history inversion on the Nitun catchment. The meaning of each figure is the same from (a) to (d). Credit: Science China Press