Catchmentwide erosion rates can now be calculated with the python package riversand
.
Input data are a digital elevation model (optional: a topographic shielding raster and a simplified map indicating the distribution of quartz-bearing rocks) in geotiff format and a shapefile with one or several catchment outlines.
The riversand
tools clip the catchment and determine the catchmentwide erosion rate based on the elevation statistics of the catchment and the cosmogenic nuclide production as predicted by the online calculator.
The calculation is fast (few seconds for one catchment) for all production scaling methods implemented in the online calculator (St, Lm, LSDn) and independent of the catchment size or the resolution of the digital elevation model.
The package is available through PyPI or github and comes with a quick tutorial and some test data.