Modelling the effects of ice transport and sediment sources on the form of detrital thermochronological age probability distributions from glacial settings

The impact of glaciers on the Quaternary evolution of mountainous landscapes remains controversial. Although in situ or bedrock low-temperature thermochronology offers insights on past rock exhumation and landscape erosion, the method also suffers from potential biases due to the difficulty of sampling bedrock buried under glaciers. Detrital thermochronology attempts to overcome this issue by sampling sediments at e.g. the catchment outlet, a component of which may originate from beneath the ice. However, detrital age distributions not only reflect the catchment exhumation, but also spatially variable patterns and rates of surface erosion and sediment transport. In this study, we use a new version of a glacial landscape evolution model, iSOSIA, to address the effect of erosion and sediment transport by ice on the form of synthetic detrital age distributions. Sediments are tracked as Lagrangian particles formed by bedrock erosion, and their transport is restricted to ice or hillslope processes, neglecting subglacial hydrology, until they are deposited. We base our model on the Tiedemann Glacier (British Columbia, Canada), which has simple morphological characteristics, such as a linear form and no connectivity to large tributary glaciers. Synthetic detrital age distributions are generated by specifying an erosion history, then sampling sediment particles at the frontal moraine of the modelled glacier. Results show that sediment sources, reflecting different processes such as glacier and hillslope erosion, can have distinct bedrock age distribution signatures, and estimating such distributions should help to identify predominant sources in the sampling site. However, discrepancies between the detrital and bedrock age distributions occur due to (i) the selective storage of a large proportion of sediments in small tributary glaciers and in lateral moraines, (ii) the large range of particle transport times due to varying transport lengths and strong variability of glacier ice velocity, (iii) the heterogeneous pattern of erosion, and (iv) the advective nature of glacier sediment transport along ice streamlines. This last factor leads to a poor lateral mixing of particle detrital signatures inside the frontal moraine, and then local sampling of the frontal moraine is likely to reflect local sources upstream. Therefore, sampling randomly across the moraine is preferred for a more representative view of the catchment age distribution. Finally, systematic comparison