RIVERSAND: A new tool for efficient computation of catchmentwide erosion rates

In-situ cosmogenic 10Be (and 26Al) concentrations in alluvial sediments provide a spatially averaged signal of the erosion rate of the catchment area. Catchmentwide erosion rates reflect the production rate of the entire basin, and their calculation requires knowledge of the complete production rate model. Available calculators determine production rates on a pixel-based approach and achieve computational efficiency by relying on a scaling method that ignores geomagnetic field strength variations. Here we introduce a new python-based tool that determines erosion rates based on the hypsometry of the catchment. The method relies on the fact that production rates are much more sensitive to changes in elevation than latitude. Our tool has two main advantages: (1) computation time is short (<30 seconds) and independent of the scaling method; there is no need to neglect magnetic field variations, and (2) because production rate scaling is performed by a widely used online calculator, the results are fully comparable to exposure ages or point-based erosion rates determined with the same calculator; future updates to production rate scaling are immediately effective for catchmentwide erosion rate calculation. We demonstrate in two case studies that (1) for similar scaling methods, our calculator reproduces pixel-based results within a few percent, and (2) erosion rates determined with different scaling methods may differ by >20%, differences can vary systematically with erosion rate, and using a time-constant scaling method may result in a bias in the interpretation of catchmentwide erosion rates.