Updated lithological ('quartz') weighing functions.
Use functions get_ages_from_server()
and get_erates_from_server()
to obtain
point-based exposure ages and erosion rates calculated by the online calculator
from sample data in your Excel spreadsheet and save the results as spreadsheet.
See updated documentation quickstart.ipynb
riversand
is a python package to calculate catchmentwide erosion rates from
cosmogenic nuclide concentrations in river sand samples. The program computes
the hypsometric statistics of the catchment area from a digital elevation model.
It uses the online erosion rate calculator by Greg Balco
(e.g. http://stoneage.hzdr.de/) to determine predicted
nuclide concentrations
The method works for in situ Be-10 and Al-26 data. It is fast (few seconds for one catchment) for all production scaling methods implemented in the online calculator (St: Lal 1991/Stone 2000; Lm: Lal/Stone with a geomagnetic correction after Nishiizumi et al. 1989; LSDn: Lifton et al. 2014) and independent of the catchment size or the resolution of the digital elevation model. It is considered robust for catchments up to approx. 600 km x 600 km; for larger catchments the effect of latitude on cosmogenic production may become significant.
The approach is described in:
Stübner, K., Balco, G., and Schmeisser, N. (2023). Riversand: a new tool for efficient computation of catchmentwide erosion rates. Radiocarbon. (link).
Definitely check out the documentation of the online calculator (e.g. here or here) and the publication Balco et al. (2008) before using this calculator.
quickstart.ipynb
step_by_step.ipynb
test_data/
: geotiffs of a 35m-resolution digital elevation model, a topographic shielding raster generated with TopoToolbox and a binary raster indicating quartz-bearing and quartz-free lithologies; shapefiles with catchment outlines; a spreadsheet with sample data.
Install latest version of riversand by running:
$ pip install riversand==1.3.1
- numpy, scipy, pandas, xarray
- rasterio, fiona, pyproj
- matplotlib