A locally calibrated,late glacial 10Be production rate from a low-latitude,high-altitude site in the Peruvian Andes

Well-dated records of tropical glacier fluctuations are essential for developing hypotheses and testing proposed mechanisms for past climate changes. Since organic material for radiocarbon dating is typically scarce in low-latitude, high-altitude environments, surface exposure-age dating, based on the measurement of in situ produced cosmogenic nuclides, provides much of the chronologic information on tropical glacier moraines. Here, we present a locally calibrated ¹⁰Be production rate for a low-latitude, high-altitude site near Quelccaya Ice Cap (∼13.95°S, 70.89°W, 4857 m asl) in the southeastern Peruvian Andes. Using an independent age (12.35 +0.2, −0.02 ka) of the late glacial Huancané IIa moraines based on thirty-four bracketing radiocarbon ages and twelve ¹⁰Be concentrations of boulders on the moraines, we determine a local production rate of 43.28 ± 2.69 atoms gram⁻¹ year⁻¹ (at g⁻¹ yr⁻¹). Reference ¹⁰Be production rates (i.e., production rates by neutron spallation appropriate for sea-level, high-latitude sites) range from 3.97 ± 0.09 to 3.78 ± 0.09 at g⁻¹ yr⁻¹, determined using scaling after Lal (1991) and Stone (2000) and depending on our assumed boulder surface erosion rate. Since our boulder surface erosion rate estimate is a minimum value, these reference production rates are also minimum values. A secondary control site on the Huancané IIIb moraines suggests that the ¹⁰Be production rates are at least as low as, or possibly lower than, those derived from the Huancané IIa moraines. These sea-level, high-latitude production rates are at least 11–15% lower than values derived using the traditional global calibration dataset, and they are also lower than those derived from the late glacial Breque moraine in the Cordillera Blanca of Peru. However, our sea-level, high-latitude production rates agree well with recently published, locally calibrated production rates from the Arctic, New Zealand, and Patagonia. The production rates presented here should be used to calculate ¹⁰Be exposure ages in low-latitude, high-altitude locations, particularly in the tropical Andes, and should improve the ability to compare the results of studies using ¹⁰Be exposure-age dating with other chronological data.     

Discordance between cosmogenic nuclide concentrations in amalgamated sands and individual fluvial pebbles in an arid zone catchment

Based on cosmogenic ¹⁰Be and ²⁶Al analyses in 15 individual detrital quartz pebbles (16–21 mm) and cosmogenic ¹⁰Be in amalgamated medium sand (0.25–0.50 mm), all collected from the outlet of the upper Gaub River catchment in Namibia, quartz pebbles yield a substantially lower average denudation rate than those yielded by the amalgamated sand sample. ¹⁰Be and ²⁶Al concentrations in the 15 individual pebbles span nearly two orders of magnitude (0.22 ± 0.01 to 20.74 ± 0.52 × 10^{6 10}Be atoms g⁻¹ and 1.35 ± 0.09 to 72.76 ± 2.04 × 10^{6 26}Al atoms g⁻¹, respectively) and yield average denudation rates of ∼0.7 m Myr⁻¹ (¹⁰Be) and ∼0.9 m Myr⁻¹ (²⁶Al). In contrast, the amalgamated sand yields an average ¹⁰Be concentration of 0.77 ± 0.03 × 10⁶ atoms g⁻¹, and an associated mean denudation rate of 9.6 ± 1.1 m Myr⁻¹, an order of magnitude greater than the rates obtained for the amalgamated pebbles. The inconsistency between the ¹⁰Be and ²⁶Al in the pebbles and the ¹⁰Be in the amalgamated sand is likely due to the combined effect of differential sediment sourcing and longer sediment transport times for the pebbles compared to the sand-sized grains. The amalgamated sands leaving the catchment are an aggregate of grains originating from all quartz-bearing rocks in all parts of the catchment. Thus, the cosmogenic nuclide inventories of these sands record the overall average lowering rate of the landscape. The pebbles originate from quartz vein outcrops throughout the catchment, and the episodic erosion of the latter means that the pebbles will have higher nuclide inventories than the surrounding bedrock and soil, and therefore also higher than the amalgamated sand grains. The order-of-magnitude grain size bias observed in the Gaub has important implications for using cosmogenic nuclide abundances in depositional surfaces because in arid environments, akin to our study catchment, pebble-sized clasts yield substantially underestimated palaeo-denudation rates. Our results highlight the importance of carefully considering geomorphology and grain size when interpreting cosmogenic nuclide data in depositional surfaces.     

CRONUS-Earth calibration samples from the Huancané II moraines,Quelccaya Ice Cap,Peru

The Huancané II moraines deposited by the Quelccaya Ice Cap in southern Peru were selected by the CRONUS-Earth Project as a primary site for evaluating cosmogenic-nuclide scaling methods and for calibrating production rates. The CRONUS-Earth Project is an effort to improve the state of the art for applications of cosmogenic nuclides to earth-surface chronology and processes. The Huancané II moraines are situated in the southern Peruvian Andes at about 4850 m and ∼13.9°S, 70.9°W. They are favorable for cosmogenic-nuclide calibration because of their low-latitude and high-elevation setting, because their age is very well constrained to 12.3 ± 0.1 ka by 34 radiocarbon ages on peat bracketing the moraines, and because boulder coverage by snow or soil is thought to be very unlikely. However, boulder-surface erosion by granular disintegration is observed and a ∼4% correction was applied to measured concentrations to compensate. Samples from 10 boulders were analyzed for ¹⁰Be, ²⁶Al, and ³⁶Cl. Interlaboratory bias at the ∼5% level was the largest contributor to variability of the ¹⁰Be samples, which were prepared by three laboratories (the other two nuclides were only prepared by one laboratory). Other than this issue, variability for all three nuclides was very low, with standard deviations of the analyses only slightly larger than the analytical uncertainties. The site production rates (corrected for topographic shielding, erosion, and radionuclide decay) at the mean site elevation of 4857 m were 45.5 ± 1.6 atoms ¹⁰Be (g quartz)⁻¹ yr⁻¹, 303 ± 15 atoms ²⁶Al (g quartz)⁻¹ yr⁻¹, and 1690 ± 100 atoms ³⁶Cl (g K)⁻¹ yr⁻¹. The nuclide data from this site, along with data from other primary sites, were used to calibrate the production rates of these three nuclides using seven global scaling methods. The traditional Lal formulation and the new Lifton-Sato-Dunai calibrations yield average ages for the Huancané samples that are in excellent-to-good agreement with the radiocarbon age control (within 0.7% for ¹⁰Be and ³⁶Cl and 6% for ²⁶Al). However, all of the neutron-monitor-based methods yielded ages that were too young by about 20%. The nuclide production ratios at this site are 6.74 ± 0.34 for ²⁶Al/¹⁰Be in quartz and 37.8 ± 2.3 (atoms ³⁶Cl (g K)⁻¹) (atom ¹⁰Be (g SiO₂)⁻¹)⁻¹ for ³⁶Cl/¹⁰Be, in sanidine and quartz, respectively.     

Toward the feldspar alternative for cosmogenic 10Be applications

The possibility of quantifying surface processes in mafic or volcanic environment using the potentialities offered by the in situ-produced cosmogenic nuclides, and more specifically by the in situ-produced ¹⁰Be, is often hampered by the rarity of quartz minerals in the available lithologies. As an alternative to overcome this difficulty, we explore in this work the possibility of relying on feldspar minerals rather that on quartz to perform in situ-produced ¹⁰Be measurements in such environments. Our strategy was to cross-calibrate the total production rate of ¹⁰Be in feldspar (P_10fsp) against the total production rate of ³He in pyroxene (P_3px) by measuring ³He and ¹⁰Be in cogenetic pyroxene (³He_px) and feldspar (¹⁰Be_fsp). The samples were collected from eight ignimbritic boulders, exposed from ca 120 to 600 ka at elevations ranging from 800 to 2500 m, along the preserved rock-avalanche deposits of the giant Caquilluco landslide (18°S, 70°W), Southern Peru. Along with data recently published by Blard et al. (2013a) at a close latitude (22°S) but higher elevation (ca. 4000 m), the samples yield a remarkably tight cluster of ³He_px - ¹⁰Be_fsp total production ratios whose weighted-mean is 35.6 ± 0.5 (1σ). The obtained weighted-mean ³He_px - ¹⁰Be_fsp total production ratio combined with the local ³He_py total production rate in the high tropical Andes published by Martin et al. (2017) allows to establish a total SLHL ¹⁰Be in situ-production rate in feldspar mineral (P_10fsp) of 3.57 ± 0.21 at.g⁻¹.yr⁻¹ (scaled for the LSD scaling scheme, the ERA40 atm model and the VDM of Lifton, 2016).Despite the large elevation range covered by the whole dataset (800–4300 m), no significant variation of the ³He_px - ¹⁰Be_fsp total production ratios in pyroxene and feldspar was evidenced. As an attempt to investigate the effect of the chemical composition of feldspar on the total ¹⁰Be production rate, major and trace element concentrations of the studied feldspar samples were analyzed. Unfortunately, giving the low compositional variability of our dataset, this issue is still pending.     

Anomalous cosmogenic 3He production and elevation scaling in the high Himalaya

The production rate of cosmogenic ³He in apatite, zircon, kyanite and garnet was obtained by cross-calibration against ¹⁰Be in co-existing quartz in glacial moraine boulders from the Nepalese Himalaya. The boulders have ¹⁰Be ages between 6 and 16 kyr and span elevations from 3200 to 4800 m. In all of these minerals ³He correlates with ¹⁰Be and is dominantly cosmogenic in origin. After modest correction for non-cosmogenic components, ³He/¹⁰Be systematics imply apparent sea-level high-latitude (SLHL) apparent production rates for ³He of 226 atoms g^? 1 yr^? 1 in zircon, 254 atoms g^? 1 yr^? 1 in apatite, 177 atoms g^? 1 yr^? 1 in kyanite, and 153 atoms g^? 1 yr^? 1 in garnet. These production rates are unexpectedly high compared with rates measured elsewhere in the world, and also compared with proposed element-specific production rates. For apatite and zircon, the data are sufficient to conclude that the ³He/¹⁰Be ratio increases with elevation. If this reflects different altitudinal scaling between production rates for the two isotopes then the SLHL production rates estimated by our approach are overestimates. We consider several hypotheses to explain these observations, including production of ³He via thermal neutron capture on ⁶Li, altitudinal variations in the energy spectrum of cosmic-ray neutrons, and the effects of snow cover. Because all of these effects are small, we conclude that the altitudinal variations in production rates of cosmogenic ³He and ¹⁰Be are distinct from each other at least at this location over the last ～ 10 kyr. This conclusion calls into question commonly adopted geographic scaling laws for at least some cosmogenic nuclides. If confirmed, this distinction may provide a mechanism by which to obtain paleoelevation estimates.     

A reevaluation of in situ cosmogenic 3He production rates

³He is among the most commonly measured terrestrial cosmogenic nuclides, but an incomplete understanding of the ³He production rate has limited robust interpretation of cosmogenic ³He concentrations. We use new measurements of cosmogenic ³He in olivine from a well-dated lava flow at Tabernacle Hill, Utah, USA, to calibrate the local ³He production rate. The new ³He measurements (n = 8) show excellent internal consistency and yield a sea level high latitude (SLHL) production rate of 123 ± 4 at g^?1 yr^?1 following the Lal (1991)/Stone (2000) scaling model [Lal, D., 1991. Cosmic ray labeling of erosion surfaces: in situ nuclide production rates and erosion models. Earth and Planetary Science Letters, 104, 424–439.; Stone, J.O., 2000. Air pressure and cosmogenic isotope production. Journal of Geophysical Research, 105, 23753–23759.]. We incorporate the new measurements from Tabernacle Hill in a compilation of all published production rate determinations, characterizing the mean global SLHL production rates (e.g. 120 ± 9.4 at g^?1 yr^?1 with Lal (1991)/Stone (2000)). The internal consistency of the global ³He production rate dataset is as good as the other commonly used cosmogenic nuclides. Additionally, ³He production rates in olivine and pyroxene agree within experimental error. The ³He production rates are implemented in an age and erosion rate calculator, forming a new module of the CRONUS-Earth web-based calculator, a simple platform for cosmogenic nuclide data interpretation [Balco, G., Stone, J., Lifton, N.A., and Dunai, T.J., 2008. A complete and easily accessible means of calculating surface exposure ages or erosion rates from ¹⁰Be and ²⁶Al measurements. Quaternary Geochronology, 3, 174–195.]. The ³He calculator is available online at http://www.cronuscalculators.nmt.edu/.     

Clean quartz matters for cosmogenic nuclide analyses: An exploration of the importance of sample purity using the CRONUS-N reference material

Reference materials are key for assessing inter-laboratory variability and measurement quality, and for placing analytical uncertainty bounds on sample analyses. Here, we investigate four years of data resulting from repeated processing of the CRONUS-N reference material for cosmogenic ¹⁰Be and ²⁶Al analyses. At University of Vermont, we prepared a CRONUS-N aliquot with most of our sample batches from 2013 to 2017; these reference material samples were then distributed to four different accelerator mass spectrometry facilities, yielding 73 ¹⁰Be analyses and 58 ²⁶Al analyses. We determine CRONUS-N ¹⁰Be concentrations of (2.26 ± 0.14) x 10⁵ atoms g⁻¹ (n = 73, mean, 1 SD) and ²⁶Al concentrations of (1.00 ± 0.08) x 10⁶ atoms g⁻¹ (n = 58, mean, 1 SD). We find a reproducibility of 6.3% for ¹⁰Be and 7.7% for ²⁶Al (relative standard deviations). We also document highly variable ²⁷Al and Mg concentrations and a ¹⁰Be dispersion twice as large as the mean AMS analytic uncertainty. Analyses of the CRONUS-N material with and without density separation demonstrate that non-quartz minerals are present in the material and have a large impact on measured concentrations of ²⁷Al, ¹⁰Be, and impurities; these non-quartz minerals represent only a very small portion of the total mass (0.6–0.8%) but have a disproportionally large effect on the resulting data. Our results highlight the importance of completely removing all non-quartz mineral phases from samples prior to Be/Al extraction for the determination of in situ cosmogenic ¹⁰Be and ²⁶Al concentrations.     

Grain size-dependent 10Be concentrations in alluvial stream sediment of the Huasco Valley,a semi-arid Andes region

Terrestrial cosmogenic nuclide concentrations in sediment are used to quantify mean denudation rates in catchments. This article explores the differences between the ¹⁰Be concentration in fine (sand) and in coarse (1–3 or 5–10 cm pebbles) river sediment. Sand and pebbles were sampled at four locations in the Huasco Valley, in the arid Chilean Andes. Sand has ¹⁰Be concentrations between 4.8 and 8.3·10⁵ at g⁻¹, while pebbles have smaller concentrations between 2.2 and 3.3·10⁵ at g⁻¹. It appears that the different concentrations, systematically measured between sand and pebbles, are the result of different denudation rates, linked with the geomorphologic processes that originated them. We propose that the ¹⁰Be concentrations in sand are determined by the mean denudation rate of all of the geomorphologic processes taking place in the catchment, including debris flow processes as well as slower processes such as hill slope diffusion. In contrast, the concentrations in pebbles are probably related to debris flows occurring in steep slopes. The mean denudation rates calculated in the catchment are between 30 and 50 m/Myr, while the denudation rates associated with debris flow are between 59 and 81 m/Myr. These denudation rates are consistent with those calculated using different methods, such as measuring eroded volumes.     

In situ cosmogenic 10Be production-rate calibration from the Southern Alps,New Zealand

We present a ¹⁰Be production-rate calibration derived from an early Holocene debris-flow deposit at about 1000 m above sea level in the central Southern Alps, New Zealand, in the mid-latitude Southern Hemisphere. Ten radiocarbon ages on macrofossils from a soil horizon buried by the deposit date the deposit to 9690 ± 50 calendar years before AD2008. Surface ¹⁰Be concentrations of seven large boulders partially embedded in the stable surface of the deposit are tightly distributed, yielding a standard deviation of ～2%. Conversion of the ¹⁰Be measurements to sea level/high-latitude values using each of five standard scaling methods indicates ¹⁰Be production rates of 3.84 ± 0.08, 3.87 ± 0.08, 3.83 ± 0.08, 4.15 ± 0.09, and 3.74 ± 0.08 atoms g^?1 a^?1, relative to the ‘07KNSTD’ ¹⁰Be AMS standard, and including only the local time-integrated production-rate uncertainties. When including a sea level high-latitude scaling uncertainty the overall error is ～2.5% (1σ) for each rate. To test the regional applicability of this production-rate calibration, we measured ¹⁰Be concentrations in a set of nearby moraines deposited before 18 060 ± 200 years before AD2008. The ¹⁰Be ages are only consistent with minimum-limiting ¹⁴C age data when calculated using the new production rates. This also suggests that terrestrial in situ cosmogenic-nuclide production did not change significantly from Last Glacial Maximum to Holocene time in New Zealand. Our production rates agree well with those of a recent calibration study from northeastern North America, but are 12–14% lower than other commonly adopted values. The production-rate values presented here can be used elsewhere in New Zealand for rock surfaces exposed during or since the last glacial period.     

Cosmogenic isotope burial dating of fluvial sediments from the Lower Rhine Embayment,Germany

Cosmogenic isotope burial dating, using ¹⁰Be and ²⁶Al, was applied to Plio–Pleistocene fluvial successions from the Lower Rhine Embayment, Germany. The approach consists of three principal steps: (1) measurement of cosmogenic nuclides in depth profiles, (2) modelling of hypothetical nuclide concentrations based on a first-order conceptualisation of the geological context and the principal succession of depositions and subsequent erosional and burial phases, and (3) using parameter estimation to identify values for the a priori unknown model parameters (burial age, initial nuclide concentrations, terrace erosion rates) that result in minimal disagreement between hypothetical and measured nuclide concentrations.The Late Pliocene Kieseloolite Formation was dated to 3650 ± 1490 ka and the Early Pleistocene Waalre Formation to 900 ± 280 ka. The unconformably overlying Upper Terrace Formation revealed ages of 740 ± 210 ka and 750 ± 250 ka for the two different sites. These findings are in good agreement with independent age control derived by bio-, magneto-, and litho-stratigraphy. Furthermore, identifiability and uncertainty analysis reveal a relationship between burial depth and sensitivity of isotope concentrations to burial age and erosion rate. Our results indicate that using shallower buried samples would enable a considerably more robust estimation of the burial age and the terrace erosion rate. Uncertainties arose mainly from nuclide measurements, and not from the uncertainties derived from modelling or insufficient knowledge of nuclide production and decay properties.     

Depth-dependence of the production rate of in situ 14C in quartz from the Leymon High core,Spain

To constrain the depth-dependence of in situ ¹⁴C production we measured the cosmogenic ¹⁴C concentration of quartz separates along a quartzite core from the Leymon High site in northwest Spain. A total of 16 quartz samples were measured over a depth range of 1–1545 cm (3–4017 g cm⁻²). The obtained ¹⁴C profile was modeled using a neutron production rate, exponentially decreasing with depth, and a fast and negative muon production parameterized as a function of the local muon flux as derived by (Heisinger et al., 2002a, 2002b). This model yields a total negative muon capture probability of 1.72 (+0.22/−0.56) × 10⁻² and a fast muon reaction cross section of 0 (+11.8/−0.00) μb. Rescaled to sea level high latitude using a Lal/Stone scaling scheme, these estimates yield a surface muon production rate of 3.31 (+0.43/−1.07) and 0 (+0.42/−0.00) at·g⁻¹ yr⁻¹ for negative muon capture and fast muons, respectively. This is the first muon production estimate for in situ ¹⁴C from a natural setting and is within uncertainty of the previous experimental estimates. The present contribution also provides new long-term blank and standard (PP-4, CRONUS-A & CRONUS-N) in situ ¹⁴C data from the ETH Zürich ¹⁴C extraction line.     

Understanding long-term soil processes using meteoric 10Be: A first attempt on loessic deposits

Meteoric ¹⁰Be, due to its high affinity with soil and sediment particles, is widely used in geomorphologic and environmental studies attempting to evaluate the soil production/denudation rates or soil ages up to 10⁷ years. However, the evolution of the ¹⁰Be distribution as a function of depth is poorly known in soils as revealed by recent reviews (Graly et al., 2010; Willenbring and von Blanckenburg, 2010). In this study, ¹⁰Be concentrations in the bulk and the 0–2 μm (lutum) granulometric fraction of samples along Luvisols profiles developed from loess in Northern France have been measured. The bulk ¹⁰Be concentrations are significantly higher in one of the three sites, likely reflecting differences in the inherited ¹⁰Be concentrations of the loess parent material as well as in the accumulation rates of the later. However, the bulk ¹⁰Be concentrations along all profiles are significantly correlated with the lutum (0–2 μm fraction) content, the maximum ¹⁰Be concentrations being evidenced in the Bt-horizon. Dominant adsorption of ¹⁰Be to the lutum has been furthermore corroborated by the mass-balance calculations with as much as 79.8 ± 9.0% of ¹⁰Be being associated with the lutum. Contrary to the bulk ¹⁰Be concentrations, the lutum ¹⁰Be concentrations showed several maxima coinciding with shifts in the coarse to fine silt ratio. This was interpreted as a change in the loess deposit dynamic. Finally, using numerical modeling approach based on the advection-diffusion equation, an average downward migration of ¹⁰Be by clay translocation was estimated. It ranges from 0.01 to 0.08 cm yr⁻¹. Inherited ¹⁰Be in the loess parent material represented from 64 to 71% of the total ¹⁰Be content in the simulated soils. Vertical ¹⁰Be distributions and their maximum concentrations in the Bt-horizon thus mainly result from redistribution of the inherited ¹⁰Be by clay translocation and bioturbation.     

Sidewall erosion: Insights from in situ-produced 10Be concentrations measured on supraglacial clasts (Mont Blanc massif,France)

Sidewall erosion because of rockfalls is one of the most efficient erosional processes in the highest parts of mountain ranges; it is therefore important to quantify sidewall erosion to understand the long-term evolution of mountainous topography. In this study, we analyse how the ¹⁰Be concentration of supraglacial debris can be used to quantify sidewall erosion in a glaciated catchment. We first analyse, in a glaciated catchment, the cascade of processes that move a rock from a rockwall to a supraglacial location and propose a quantitative estimate of the number of rockfalls statistically mixed in a supraglacial sand sample. This model incorporates the size of the rockwall, a power law distribution of the size of the rockfalls and the mean glacial transport velocity. In the case of the Bossons glacier catchment (Mont Blanc massif), the ¹⁰Be concentrations obtained for supraglacial samples vary from 1.97 ± 0.24 to 23.82 ± 1.68 × 10⁴ atoms g⁻¹. Our analysis suggests that part of the ¹⁰Be concentration dispersion is related to an insufficient number of amalgamated rockfalls that does not erase the stochastic nature of the sidewall erosion. In the latter case, the concentration of several collected samples is averaged to increase the number of statistically amalgamated rockfalls. Variable and robust ¹⁰Be-derived rockwall retreat rates are obtained for three distinct rockfall zones in the Bossons catchment and are 0.19 ± 0.08 mm year⁻¹, 0.54 ± 0.1 mm year⁻¹ and 1.08 ± 0.17 mm year⁻¹. The mean ¹⁰Be retreat rate for the whole catchment (ca. 0.65 mm year⁻¹) is close to the present-day erosion rate derived from other methods. © 2019 John Wiley & Sons, Ltd.     

Nature of decadal-scale sediment accumulation on the western shelf of the Mississippi River delta

Sediment delivered to coastal systems by rivers (15×10⁹ tons) plays a key role in the global carbon and nutrient cycles, as deltas and continental shelves are considered to be the main repositories of organic matter in marine sediments. The Mississippi River, delivering more than 60% of the total dissolved and suspended materials from the conterminous US, dominates coastal and margin processes in the northern Gulf of Mexico. Draining approximately 41% of the conterminous US, the Mississippi and Atchafalaya river system deliver approximately 2×10⁸ tons of suspended matter to the northern Gulf shelf each year. Unlike previous work, this study provides a comprehensive evaluation of sediment accumulation covering majority of the shelf (<150 m water depth) west of the Mississippi Delta from 92 cores collected throughout the last 15 years. This provides a unique and invaluable data set of the spatial and modern temporal variations of the sediment accumulation in this dynamic coastal environment.Three types of ²¹⁰Pb profiles were observed from short cores (15–45 cm) collected on the shelf. Proximal to Southwest Pass in 30–100 m water depths, non-steady-state profiles were observed indicating rapid accumulation. Sediment accumulation rates in this area are typically >2.5 cm yr⁻¹ (>1.8 g cm⁻² yr⁻¹). Kasten cores (∼200 cm in length) collected near Southwest Pass also indicate rapid deposition (>4 cm yr⁻¹; >3 g cm⁻² yr⁻¹) on a longer timescale than that captured in the box cores. Near shore (<20 m), profiles are dominated by sediments reworked by waves and currents with no accumulation (the exception is an area just south of Barataria Bay where accumulation occurs). The remainder of the shelf (distal of Southwest Pass) is dominated by steady-state accumulation beneath a ∼10-cm thick mixed layer. Sediment accumulation rates for the distal shelf are typically <0.7 cm yr⁻¹ (<0.5 g cm⁻² yr⁻¹). A preliminary sediment budget based on the distribution of ²¹⁰Pb accumulation rates indicates that 40–50% of the sediment delivered by the river is transported out of the study region. Sediment is moved to distal regions of the shelf/slope through two different mechanisms. Along-isobath sediment movement occurs by normal resuspension processes west of the delta, whereas delivery of sediments south and southwest of the delta may be also be influenced by mass movement events on varying timescales.     

10Be‐derived assessment of accelerated erosion in a glacially conditioned inner gorge,Entlebuch, Central Alps of Switzerland

Inner gorges often result from the propagation of erosional waves related to glacial/interglacial climate shifts. However, only few studies have quantified the modern erosional response to this glacial conditioning. Here, we report in situ ¹⁰Be data from the 64 km² Entlen catchment (Swiss Alps). This basin hosts a 7 km long central inner gorge with two tributaries that are >100 m‐deeply incised into thick glacial till and bedrock. The ¹⁰Be concentrations measured at the downstream end of the gorge yield a catchment‐wide erosion rate of 0.42 ± 0.04 mm yr^‐1, while erosion rates are consistently lower upstream of the inner gorge, ranging from 0.14 ± 0.01 mm yr^‐1 to 0.23 ± 0.02 mm yr^‐1. However, ¹⁰Be‐based sediment budget calculations yield rates of ~1.3 mm yr^‐1 for the inner gorge of the trunk stream. Likewise, in the two incised tributary reaches, erosion rates are ~2.0 mm yr^‐1 and ~1.9 mm yr^‐1. Moreover, at the erosional front of the gorge, we measured bedrock incision rates ranging from ~2.5 mm yr^‐1 to ~3.8 mm yr^‐1. These rates, however, are too low to infer a post‐glacial age (15–20 ka) for the gorge initiation. This would require erosion rates that are between 2 and 6 times higher than present‐day estimates. However, the downcutting into unconsolidated glacial till favored high erosion rates through knickzone propagation immediately after the retreat of the LGM glaciers, and subsequent hillslope relaxation led to a progressive decrease in erosion rates. This hypothesis of a two‐ to sixfold decrease in erosion rates does not conflict with the ¹⁰Be‐based erosion rate budgets, because the modern erosional time scale recorded by ¹⁰Be cover the past 2–3 ka only. These results point to the acceleration of Holocene erosion in response to the glacial overprint of the landscape. Copyright © 2012 John Wiley & Sons, Ltd.     

P–PINI: A cosmogenic nuclide burial dating method for landscapes undergoing non-steady erosion

Existing methods of cosmogenic nuclide burial dating perform well provided that sediment sources undergo steady rates of erosion and the samples experience continuous exposure to cosmic rays. These premises exert important limitations on the applicability of the methods. And yet, high mountain sediment sources are rife with transient processes, such as non-steady erosion by glacial quarrying and/or landsliding, or temporary cosmic-ray shielding beneath glaciers and/or sediment. As well as breaching the premises of existing burial dating methods, such processes yield samples with low nuclide abundances and variable ²⁶Al/¹⁰Be ratios that may foil both isochron and simple burial-age solutions. P–PINI (Particle-Pathway Inversion of Nuclide Inventories) is a new dating tool designed for dating the burial of sediments sourced from landscapes characterized by abrupt, non-steady erosion, discontinuous exposure, and catchments with elevation-dependent ²⁶Al/¹⁰Be production ratios. P–PINI merges a Monte Carlo simulator with established cosmogenic nuclide production equations to simulate millions of samples (¹⁰Be–²⁶Al inventories). The simulated samples are compared statistically with ¹⁰Be–²⁶Al measured in field samples to define the most probable burial age. Here, we target three published ¹⁰Be–²⁶Al datasets to demonstrate the versatility of the P–PINI model for dating fluvial and glacial sediments. (1) The first case serves as a robust validation of P–PINI. For the Pulu fluvial gravels (China), we obtain a burial age of 1.27 ± 0.10 Ma (1σ), which accords with the isochron burial age and two independent chronometers reported in Zhao et al. (2016) Quaternary Geochronology 34, 75–80. The second and third cases, however, reveal marked divergence between P–PINI and isochron-derived ages. (2) For the fluvial Nenana Gravel (USA), we obtain a minimum-limiting burial age of 4.5 ± 0.7 Ma (1σ), which is compatible with unroofing of the Alaska Range starting ∼ 6 Ma, while calling into question the Early Pleistocene isochron burial age presented in Sortor et al. (2021) Geology 49, 1473–1477. (3) For the Bünten Till (Switzerland), we obtain a limiting burial age of <204 ka (95th percentile range), which conforms with the classical notion of the most extensive glaciation in the northern Alpine Foreland assigned to the Riss glaciation (sensu marine isotope stage 6) contrary to the isochron burial age presented in Dieleman et al. (2022) Geosciences, 12, 39. Discrepancies between P–PINI and the isochron ages are rooted in the challenges posed by the diverse pre-burial ²⁶Al/¹⁰Be ratios produced under conditions characteristic of high mountain landscapes; i.e., non-steady erosion, discontinuous cosmic-ray exposure, and elevation-dependent ²⁶Al/¹⁰Be production ratios in the source region, which are incompatible with the isochron method, but easily accommodated by the stochastic design of P–PINI.     

Meteoric 10Be speciation in subglacial sediments of East Antarctica

A sequential chemical extraction procedure was developed and tested to investigate the utility of meteoric ¹⁰Be as a tracer for authigenic mineral formation beneath the East Antarctic Ice Sheet. Subglacial meltwater is widely available under the Antarctic Ice Sheet and dissolved gases within it have the potential to drive chemical weathering processes in the subglacial environment. Meteoric ¹⁰Be is a cosmogenic nuclide with a half-life of 1.39·10⁶ years that is incorporated into glacier ice, therefore its abundance in the subglacial environment in Antarctica is meltwater dependent. It is known to adsorb to fine-grained particles in aqueous solution, precipitate with amorphous oxides/hydroxides, and/or be incorporated into authigenic clay minerals during chemical weathering. The presence of ¹⁰Be in chemical weathering products derived from beneath the ice therefore indicates chemical weathering processes in the subglacial environment. Freshly emerging subglacial sediments from the Mt. Achernar blue ice moraine were subject to chemical extractions where these weathering phases were isolated and ¹⁰Be concentrations therein quantified. Optimization of the phase isolation was developed by examining the effects of each extraction on the sample mineralogy and chemical composition. Experiments on ¹⁰Be desorption revealed that pH 3.2–3.5 was optimal for the extraction of adsorbed ¹⁰Be. Vigorous disaggregation of the samples before grain size separations and acid extractions is crucial due to the incorporation of the nuclide in clay minerals and its preferential absorption to clay-sized particles. ¹⁰Be concentrations of 2–22·10⁷ atoms·g⁻¹ measured in oxides and clay minerals in freshly emerging sediments strongly indicate subglacial chemical weathering in the catchment of the Mt. Achernar moraine. Based on total ¹⁰Be sample concentrations, local basal melt rates, and ¹⁰Be ice concentrations, sediment-meltwater contact in the subglacial environment is on the order of thousands of years per gram of underlying fine sediment. Strong correlation (R = 0.97) between ¹⁰Be and smectite abundance in the sediments supports authigenic clay formation in the subglacial environment. This suggests meteoric ¹⁰Be is a useful tool to characterize subglacial geochemical weathering processes under the Antarctic Ice Sheet.     

Paleoseismology of the 2016 MW 6.1 Petermann earthquake source: Implications for intraplate earthquake behaviour and the geomorphic longevity of bedrock fault scarps in a low strain-rate cratonic region

The 20 May 2016 M_W 6.1 Petermann earthquake in central Australia generated a 21 km surface rupture with 0.1 to 1 m vertical displacements across a low-relief landscape. No paleo-scarps or potentially analogous topographic features are evident in pre-earthquake Worldview-1 and Worldview-2 satellite data. Two excavations across the surface rupture expose near-surface fault geometry and mixed aeolian-sheetwash sediment faulted only in the 2016 earthquake. A 10.6 ± 0.4 ka optically stimulated luminescence (OSL) age of sheetwash sediment provides a minimum estimate for the period of quiescence prior to 2016 rupture. Seven cosmogenic beryllium-10 (¹⁰Be) bedrock erosion rates are derived for samples < 5 km distance from the surface rupture on the hanging-wall and foot-wall, and three from samples 19 to 50 km from the surface rupture. No distinction is found between fault proximal rates (1.3 ± 0.1 to 2.6 ± 0.2 m Myr⁻¹) and distal samples (1.4 ± 0.1 to 2.3 ± 0.2 m Myr⁻¹). The thickness of rock fragments (2–5 cm) coseismically displaced in the Petermann earthquake perturbs the steady-state bedrock erosion rate by only 1 to 3%, less than the erosion rate uncertainty estimated for each sample (7–12%). Using ¹⁰Be erosion rates and scarp height measurements we estimate approximately 0.5 to 1 Myr of differential erosion is required to return to pre-earthquake topography. By inference any pre-2016 fault-related topography likely required a similar time for removal. We conclude that the Petermann earthquake was the first on this fault in the last ca. 0.5–1 Myr. Extrapolating single nuclide erosion rates across this timescale introduces large uncertainties, and we cannot resolve whether 2016 represents the first ever surface rupture on this fault, or a > 1 Myr interseismic period. Either option reinforces the importance of including distributed earthquake sources in fault displacement and seismic hazard analyses.     

Volume and mass changes of the Greenland ice sheet inferred from ICESat and GRACE

This study examines the recent evolution of the Greenland ice sheet and its six major drainage basins. Based on laser altimetry data acquired by the Ice, Cloud and Land Elevation Satellite (ICESat), covering the period September–November 2003 to February–March 2008, ice surface height changes and their temporal variations were inferred. Our refined repeat track analysis is solely based on ICESat data and is independent of external elevation models, since it accounts for both ice height changes and the local topography. From the high resolution ice height change pattern we infer an overall mean surface height trend of −0.12 ± 0.006 m yr⁻¹. Furthermore, the largest changes could be identified at coastal margins of the ice sheet, exhibiting rates of more than −2 m yr⁻¹. The total ice volume change of the entire ice sheet amounts to −205.4 ± 10.6 km³ yr⁻¹. In addition, we assessed mass changes from 78 monthly Gravity Recovery and Climate Experiment (GRACE) solutions. The Release-04 gravity field solutions of GeoForschungsZentrum Potsdam cover the period between August 2002 and June 2009. We applied an adjusted regional integration approach in order to minimize the leakage effects. Attention was paid to an optimized filtering which reduces error effects from different sources. The overall error assessment accounts for GRACE errors as well as for errors due to imperfect model reductions. In particular, errors caused by uncertainties in the glacial isostatic adjustment models could be identified as the largest source of errors. Finally, we determined both seasonal and long-term mass change rates. The latter amounts to an overall ice mass change of −191.2 ± 20.9 Gt yr⁻¹ corresponding to 0.53 ± 0.06 mm yr⁻¹ equivalent eustatic sea level rise. From the combination of the volume and mass change estimates we determined a mean density of the lost mass to be 930 ± 11 kg m⁻³. This value supports our applied density assumption 900 ± 30 kg m⁻³ which was used to perform the volume–mass-conversion of our ICESat results. Hence, mass change estimates from two independent observation techniques were inferred and are generally in good agreement.     

An inter-laboratory comparison of cosmogenic 3He and radiogenic 4He in the CRONUS-P pyroxene standard

This study reports an inter-laboratory comparison of the ³He and ⁴He concentrations measured in the pyroxene material CRONUS-P. This forms part of the CRONUS-Earth and CRONUS-EU programs, which also produced a series of natural reference materials for in situ produced ²⁶Al, ¹⁰Be, ¹⁴C, ²¹Ne and ³⁶Cl.Six laboratories (GFZ Potsdam, Caltech Pasadena, CRPG Nancy, SUERC Glasgow, BGC Berkeley, Lamont New York) participated in this intercomparison experiment, analyzing between 5 and 22 aliquots each. Intra-laboratory results yield ³He concentrations that are consistent with the reported analytical uncertainties, which suggests that ³He is homogeneous within CRONUS-P. The inter-laboratory dataset (66 determinations from the 6 different labs) is characterized by a global weighted mean of (5.02 ± 0.12) × 10⁹ at g⁻¹ with an overdispersion of 5.6% (2σ). ⁴He is characterized by a larger variability than ³He, and by an inter-lab global weighted mean of (3.60 ± 0.18) × 10¹³ at g⁻¹ (2σ) with an overdispersion of 10.4% (2σ).There are, however, some systematic differences between the six laboratories. More precisely, 2 laboratories obtained mean ³He concentrations that are about 6% higher than the clustered other 4 laboratories. This systematic bias is larger than the analytical uncertainty and not related to the CRONUS-P material (see Schaefer et al., 2015). Reasons for these inter-laboratory offsets are difficult to identify but are discussed below. To improve the precision of cosmogenic ³He dating, we suggest that future studies presenting cosmogenic ³He results also report the ³He concentration measured in the CRONUS-P material in the lab(s) used in a given study.