Glacial influence on late Pleistocene 10Be-derived paleo-erosion rates in the north-western Himalaya, India

Terrestrial cosmogenic nuclide concentrations in fluvial deposits allow estimation of paleo-erosion rates and reconstruction of the response of landscapes to climatic perturbations. In partly ice-covered landscapes, however, incorporation of subglacially-derived sediments that were shielded by ice from cosmic can lead to erroneous erosion rate calculations. Here, we combine in situ-produced Be-10-derived erosion rates, based on sand and pebbles from a fluvial fill terrace and the modern riverbed in the upper Yamuna catchment, with numerical ice flow modelling to quantify this bias. New luminescence and surface exposure ages suggest that aggradation of the exposed deposits occurred between 29.9 +/- 2.5 ka and 14.8 +/- 2.8 ka. During most of the deposition, glaciers probably covered similar to 19% of the catchment. Be-10 concentrations of terrace sand samples differ from those of pebble samples. We obtained the lowest erosion rates from quartzite pebbles, which stem from low elevations, and the highest erosion rates from crystalline pebbles, which stem from high elevations in the Yamuna catchment. We explain these different erosion rates by differences in the steepness of the source areas, an effect that prevails throughout the entire aggradation period despite significant former ice-cover. Sand samples, which are thought to be derived from all elevation parts of the catchment, however show lower Be-10 concentrations during the aggradation compared to present-day. We argue that this difference is due to a substantial subglacial origin of the sand during the aggradation period, and not necessarily related to enhanced erosion. We conclude that aggradation of the valley fill in the Yamuna catchment is most likely due to reduced discharge, and only marginally related to higher erosion rates during the late Pleistocene. (C) 2020 Elsevier B.V. All rights reserved.