The sensitivity of cosmogenic radionuclide analysis to soil bulk density: implications for soil formation rates

Improving our knowledge of soil formation is critical so that we can better understand the first‐order controls on soil thickness and more effectively inform land‐management decisions. Cosmogenic radionuclide analysis has allowed soil scientists to more accurately constrain the rates at which soils form from bedrock. In such analysis, the concentration of an isotope, such as Beryllium‐10, is measured from a sample of bedrock. Because this concentration is partly governed by the lowering of the bedrock‐soil interface, a cosmogenic depth‐profile model can be fitted to infer the bedrock and surface lowering rates compatible with the measured concentrations. Given that the bedrock‐soil interface is shielded by soil, the cosmic rays responsible for the in‐situ production of the radionuclide are attenuated, with attenuation rates dependent on the density profile of this soil. Many studies have assumed that soil bulk density is either equal to that of the bedrock or constant with depth. The failure to acknowledge the variations in soil bulk density means that cosmogenically derived soil formation rates previously published may be under‐ or overestimates. Here, we deploy a new model called “CoSOILcal” to a global compilation of cosmogenic analyses of soil formation and, by making use of estimated bulk density profiles, recalculate rates of soil formation to assess the sensitivity to this important parameter. We found that where a soil mantle >0.25 m overlies the soil‐bedrock interface, accounting for the soil bulk density profile brings about a significantly slower rate of soil formation than that previously published. Moreover, the impact of using bulk density profiles on cosmogenically derived soil formation rates increases as soil thickens. These findings call into question the accuracy of our existing soil formation knowledge and we suggest that future cosmogenic radionuclide analysis must consider the bulk density profile of the overlying soil

Timing of Pleistocene glaciations in the High Atlas, Morocco: new 10Be and 36Cl exposure ages

This paper presents data from 42 new samples yielding Late Pleistocene cosmogenic 10Be and 36Cl exposure ages of moraine boulders across a series of glaciated valleys in the Toubkal Massif (4167 m a.s.l.), High Atlas, Morocco. This represents the first comprehensive Pleistocene glacial chronology in North Africa and one of the largest datasets from the Mediterranean region. The timing of these glacier advances has major implications for understanding the influence of Atlantic depressions on moisture supply to North Africa and the Mediterranean basin during the Pleistocene. The oldest and lowest moraines which span elevations from ∼1900 to 2400 m a.s.l. indicate that the maximum glacier advance occurred from MIS 5 to 3 with a combined mean 10Be and 36Cl age of 50.2 ± 19.5 ka (1 SD; n = 12, 7 outliers). The next moraine units up-valley at higher elevations (∼2200–2600 m a.s.l.) yielded exposure ages close to the global Last Glacial Maximum (LGM) with a combined mean 10Be and 36Cl age of 22.0 ± 4.9 ka (1 SD; n = 9, 7 outliers). The youngest exposure ages are from moraines that were emplaced during the Younger Dryas with a combined mean 10Be and 36Cl age of 12.3 ± 0.9 ka (1 SD; n = 7, no outliers) and are found in cirques at the highest elevations ranging from ∼2900 to 3300 m a.s.l. From moraines predating the Younger Dryas, a large number of young outliers are spread evenly between 6 and 13 ka suggesting a continuing process of exhumation or repositioning of boulders during the early to mid-Holocene. This attests to active seismic processes and possibly intense erosion during this period

Late Quaternary glaciation in the Hebrides sector of the continental shelf : cosmogenic nuclide dating of glacial events on the St Kilda archipelago

The authors thank NERC-CIAF for funding analysis of the 10Be and 36Cl exposure ages (Allocation 9116/0412), the Carnegie Trust for the Universities of Scotland for a grant towards travel expenses.The St Kilda archipelago lies ~65 km west of the Outer Hebrides and ~60 km east of the Atlantic shelf break, and represents a key site for testing the assertion that during the Last Local Glacial Maximum (LLGM; c. 27 ka) the British–Irish Ice Sheet (BIIS) extended to near the shelf edge in all sectors. Two consistent cosmogenic 36Cl exposure ages averaging (≥) 81.6±7.8 ka for perched boulders at 290 m altitude demonstrate that the last ice sheet failed to over-run high ground on the largest island, Hirta. 36Cl and 10Be exposure ages for glacially emplaced boulders on low ground indicate deposition by small, locally nourished glaciers that last occupied a north-facing valley (Gleann Mòr) at c. 30.9±3.2 ka, prior to extension of the last ice sheet to the outer shelf, and a south-facing valley (Village Bay) at c. 19.2±2.3 ka, several millennia after the LLGM. Our dating evidence is consistent with previous interpretations of lithostratigraphical, seismostratigraphical and geomorphological evidence and confirms that the last ice sheet failed to encroach on St Kilda. A simple ice-flow model demonstrates that even if thin, low-gradient ice lobes encircled the archipelago during the LLGM, the ice margin can only have reached the outermost moraine banks, ~40 km west of St Kilda, under extremely low (<2 kPa) driving stresses, implying either surge-like transient streaming behaviour at the ice-sheet margin or that the moraine banks relate to an earlier, more extensive ice sheet. The final glaciation of the Village Bay area at c. 19.2±2.3 ka was out of phase with the behaviour of the BIIS, which was undergoing net retreat during this period.PostprintPeer reviewe

Combining surface exposure dating and burial dating from paired cosmogenic depth profiles. Example of El Límite alluvial fan in Huércal-Overa basin (SE Iberia)

Cosmogenic nuclide depth-profiles are used to calculate the age of landforms, the rates at which erosion has affected them since their formation and, in case of deposits, the paleo-erosion rate in the source area. However, two difficulties are typically encountered: 1) old deposits or strongly affected by cosmogenic nuclide inheritance often appear to be saturated, and 2) a full propagation of uncertainties often yields poorly constrained ages. Here we show how to combine surface-exposure-dating and burial-dating techniques in the same profile to get more accurate age results and to constrain the extent of pre-depositional burial periods. A 10Be-26Al depth-profile measured in an alluvial fan of SE Iberia is presented as a natural example

A composite 10Be, IR-50 and 14C chronology of the pre-LGM full ice extent of the western Patagonian Ice Sheet in the Isla de Chiloé, south Chile (42ºS

Unanswered questions about the glacier and climate history preceding the global Last Glacial Maximum (LGM) in the southern temperate latitudes remain. The Marine Isotope Stage (MIS) 3 is normally understood as a global interstadial period; nonetheless its climate was punctuated by conspicuous variability, and its signature has not been resolved beyond the polar realms. In this paper, we compile a 10Be depth profile, single grain infrared (IR) stimulated luminescence dating and 14C samples to derive a new glacier record for the principal outwash plain complex, deposited by the western Patagonian Ice Sheet (PIS) during the last glacial period (Llanquihue Glaciation) on the Isla de Chiloé, southern Chile (42∘ S). In this region, the Golfo de Corcovado Ice Lobe left a distinct geomorphic and stratigraphic imprint, suitable for reconstructing former ice dynamics and timing of past climate change. Our data indicate that maximum glaciation occurred by 57.8±4.7 ka without reaching the Pacific Ocean coast. Ice readavanced and buttressed against the eastern side of the Cordillera de la Costa again by 26.0±2.9 ka. Our data further support the notion of a large ice extent during parts of the MIS 3 in Patagonia and New Zealand but appear to contradict near contemporaneous interstadial evidence in the southern midlatitudes, including Chiloé. We propose that the PIS expanded to its full-glacial Llanquihue moraines, recording a rapid response of southern mountain glaciers to the millennial-scale climate stadials that punctuated the MIS 3 at the poles and elsewhere.</p

Controls on the erosion of the continental margin of southeast Brazil from cosmogenic 10Be in river sediments

The Atlantic Ocean coast region of southeast Brazil contains two coast-parallel mountain ranges (the Serra do Mar and Serra da Mantiqueira) generated by tectonic activity pulses tens of millions years after the main continental rift event occurred around 120 Ma. Although the short-term erosion rates for the region are established, the relative importance of the factors controlling erosion is poorly constrained. We combine new and published catchment-averaged erosion rates (n = 48) using in situ-produced 10Be concentrations in quartz from river sediments to establish the regional erosion pattern. The river catchments are (i) escarpment topography, (ii) high-altitude low-relief and (iii) mixed topography, which record how escarpment fronts are migrating inland. Ocean-facing coastal escarpment catchments of the Serra do Mar (ε = 18–53 m/Ma) can be eroded approximately twice as fast as continent-facing escarpment catchments in the Serra do Mar and Serra da Mantiqueira (ε = 7–24 m/Ma). The correlation between the normalized channel steepness index (ksn) and slope angle indicates that river incision and hillslope erosion processes combine to maintain the high relief. The Serra do Mar catchments define a mean slope angle threshold indicating that landslides are the dominant erosional process when slope angles in excess of ~30°. Tectonic activity is low and plays no significant role in driving erosion. A first-order relationship between erosion rate and precipitation-temperature across the region implies that climate plays a key role in soil production, river incision and in triggering erosional processes. Although the high topographic relief is a pre-condition for the occurrence of significant erosion, the climatic condition is the outlining factor of the regional variation in erosion rates

Regional mid-Pleistocene glaciation in central Patagonia

Southern South America contains a glacial geomorphological record that spans the past million years and has the potential to provide palaeoclimate information for several glacial periods in Earth's history. In central Patagonia, two major outlet glaciers of the former Patagonian Ice Sheet carved deep basins ∼50 km wide and extending over 100 km into the Andean plain east of the mountain front. A succession of nested glacial moraines offers the possibility of determining when the ice lobes advanced and whether such advances occurred synchronously. The existing chronology, which was obtained using different methods in each valley, indicates the penultimate moraines differ in age by a full glacial cycle. Here, we test this hypothesis further using a uniform methodology that combines cosmogenic nuclide ages from moraine boulders, moraine cobbles and outwash cobbles. 10Be concentrations in eighteen outwash cobbles from the Moreno outwash terrace in the Lago Buenos Aires valley yield surface exposure ages of 169–269 ka. We find 10Be inheritance is low and therefore use the oldest surface cobbles to date the deposit at 260–270 ka, which is indistinguishable from the age obtained in the neighbouring Lago Pueyrredón valley. This suggests a regionally significant glaciation during Marine Isotope Stage 8, and broad interhemispheric synchrony of glacial maxima during the mid to late Pleistocene. Finally, we find the dated outwash terrace is 70–100 ka older than the associated moraines. On the basis of geomorphological observations, we suggest this difference can be explained by exhumation of moraine boulders

Arable soil formation and erosion: a hillslope-based cosmogenic nuclide study in the United Kingdom

Arable soils are critical resources that support multiple ecosystem services. They are frequently threatened, however, by accelerated erosion. Subsequently, policy to ensure their long-term security is an urgent societal priority. Although their long-term security relies upon a balance between the rates of soil loss and formation, there have been few investigations of the formation rates of soils supporting arable agriculture. This paper addresses this knowledge gap by presenting the first isotopically constrained soil formation rates for an arable (Nottinghamshire, UK) and coniferous woodland hillslope (Shropshire, UK). Rates ranged from 0.026 to 0.096 mm yr−1 across the two sites. These rates fall within the range of previously published rates for soils in temperate climates and on sandstone lithologies but significantly differed from those measured in the only other UK-based study. We suggest this is due to the parent material at our sites being more susceptible to weathering. Furthermore, soil formation rates were found to be greatest for aeolian-derived sandstone when compared with fluvially derived lithology raising questions about the extent to which the petrographic composition of the parent material governs rates of soil formation. On the hillslope currently supporting arable agriculture, we utilized cosmogenically derived rates of soil formation and erosion in a first-order lifespan model and found, in a worst-case scenario, that the backslope A horizon could be eroded in 138 years with bedrock exposure occurring in 212 years under the current management regime. These findings represent the first quantitative estimate of cultivated soil lifespans in the UK

Moraine crest or slope: An analysis of the effects of boulder position on cosmogenic exposure age

Terrestrial cosmogenic nuclide dating of ice-marginal moraines can provide unique insights into Quaternary glacial history. However, pre- and post-depositional exposure histories of moraine boulders can introduce geologic uncertainty to numerical landform ages. To avoid geologic outliers, boulders are typically selected based on their depositional context and individual characteristics but while these criteria have good qualitative reasoning, many have not been tested quantitatively. Of these, boulder location is critical, as boulders located on moraine crests are prioritised, while those on moraine slopes are typically rejected. This study provides the first quantitative assessment of the relative utility of moraine crest and moraine slope sampling using new and published 10Be and 36Cl ages (n = 19) and Schmidt hammer sampling (SH; n = 635 moraine boulders, ∼19,050 SH R-values) in the northern and southern Pyrenees. These data show that for many of the studied moraines, the spatial distribution of “good” boulders is effectively random, with no consistent clustering on moraine crests, ice-proximal or -distal slopes. In turn, and in contrast to prior work, there is no clear penalty to either moraine crest or moraine slope sampling. Instead, we argue that landform stability exerts a greater influence on exposure age distributions than the characteristics of individual boulders. For the studied landforms, post-depositional stability is strongly influenced by sedimentology, with prolonged degradation of matrix-rich unconsolidated moraines while boulder-rich, matrix-poor moraines stabilised rapidly after deposition. While this pattern is unlikely to hold true in all settings, these data indicate that differences between landforms can be more significant than differences at the intra-landform scale. As ad hoc assessment of landform stability is extremely challenging based on geomorphological evidence alone, preliminary SH sampling, as utilised here, is a useful method to assess the temporal distribution of boulder exposure ages and to prioritise individual boulders for subsequent analysis

Extensive MIS 3 glaciation in southernmost Patagonia revealed by cosmogenic nuclide dating of outwash sediments

The timing and extent of former glacial advances can demonstrate leads and lags during periods of climatic change and their forcing, but this requires robust glacial chronologies. In parts of southernmost Patagonia, dating pre-global Last Glacial Maximum (gLGM) ice limits has proven difficult due to post-deposition processes affecting the build-up of cosmogenic nuclides in moraine boulders. Here we provide ages for the Río Cullen and San Sebastián glacial limits of the former Bahía Inútil–San Sebastián (BI-SSb) ice lobe on Tierra del Fuego (53–54°S), previously hypothesised to represent advances during Marine Isotope Stages (MIS) 12 and 10, respectively. Our approach uses cosmogenic 10Be and 26Al exposure dating, but targets glacial outwash associated with these limits and uses depth-profiles and surface cobble samples, thereby accounting for surface deflation and inheritance. The data reveal that the limits formed more recently than previously thought, giving ages of 45.6 ka for the Río Cullen, and 30.1 ka for the San Sebastián limits. These dates indicate extensive glaciation in southern Patagonia during MIS 3, prior to the well-constrained, but much less extensive MIS 2 (gLGM) limit. This suggests the pattern of ice advances in the region was different to northern Patagonia, with the terrestrial limits relating to the last glacial cycle, rather than progressively less extensive glaciations over hundreds of thousands of years. However, the dates are consistent with MIS 3 glaciation elsewhere in the southern mid-latitudes, and the combination of cooler summers and warmer winters with increased precipitation, may have caused extensive glaciation prior to the gLGM