The implications of sampling approach and geomorphological processes for cosmogenic 10Be exposure dating of marine terraces

Abstract Concentrations of cosmogenic nuclides in rocks at the Earth’s surface are routinely measured by AMS (accelerator mass spectrometry) to obtain exposure ages. One particular application of this technique has been to derive ages for the formation of marine terraces, thus allowing constraints to be placed on rates of tectonic uplift. However, multiple rock samples from the same terrace surface have typically shown an amount of scatter in cosmogenic nuclide concentrations in excess of analytical uncertainties, potentially undermining the confidence in the resulting uplift rates. In addition, earlier works have considered the application of bedrock samples for exposure age dating marine terraces, but little attention has been directed towards the validity of sampling clasts exposed on terrace surfaces. Here, we use 10Be and 26Al measurements from a flight of marine terraces in northern Chile to compare the results from pebbles with those from bedrock samples and to investigate assumptions that sample pre-exposure and burial are negligible. We examine the influence that processes of surface deflation and pebble erosion will have on cosmogenic nuclide concentrations using a Monte Carlo model, the results of which are compared to our 10Be measurements. We find good agreement of 10Be concentrations between bedrock and pebble samples. At the sampled location, sample burial is rare and 10Be inheritance in pebble samples is low relative to the age of the terraces. Our modelling suggests that in arid environments, such as northern Chile, pebble erosion will be a secondary effect compared to surface deflation and that multiple periods of deflation will result in multimodal distributions of 10Be concentrations. More broadly, our findings show that measuring multiple surface clasts from a single marine terrace reveals the geomorphological processes influencing exposure ages, and may also be used to help identify the operation of past surface processes such as episodic deflation.