Cosmogenic Research Articles

A tool for the ages: The Probabilistic Cosmogenic Age Analysis Tool (P-CAAT)

While revolutionary to the geomorphic community, the application of terrestrial cosmogenic nuclide (TCN) dating is complicated by geological uncertainties, which often lead to skewed or poorly clustered TCN age distributions. Although a range of statistical approaches are typically used to detect and remove outliers, few are optimized for analysis of TCN datasets. Many are mean- or median-based and therefore explicitly assume a single probability distribution (e.g., Mean Squared Weighted Deviates, Chauvenet's Criterion, etc.). Given the ubiquity of pre- and post-depositional modification of rock surfaces, which occur at different rates in different geomorphic settings, these approaches struggle with multimodal distributions which often characterize TCN datasets. In addition, most statistical approaches do not propagate measurement or production rate uncertainties, which become increasingly important as dataset size or clustering increases. Finally, most approaches provide arithmetic single solutions, irrespective of geologic context. To address these limitations, we present the Probabilistic Cosmogenic Age Analysis Tool (P-CAAT), a new approach for outlier detection and landform age analysis. This tool incorporates both sample age and geologic uncertainties and uses Monte Carlo simulations to eliminate dataset skewness by isolating component normal distributions from a cumulative probability density estimate for datasets with three or more samples. This approach allows geologic context to inform post-analysis interpretations, as researchers can assign landform ages based upon statistically distinct subpopulations, informed by the characteristics of geomorphic systems (e.g., exhumation of boulders as moraines degrade through time). To evaluate the effectiveness of P-CAAT, we analyzed a range of synthetic TCN datasets and compared the results to commonly used statistical approaches for outlier detection. Irrespective of dataset size or clustering, P-CAAT outperformed other approaches and returned accurate solutions that improve in precision as sample size increases. To enable more comprehensive utilization of our approach, P-CAAT is packaged with a GUI interface and is available for download at kgs. uky.edu/anorthite/PCAAT.

Terrestrial cosmogenic 10Be dating of the Última Esperanza ice lobe moraines (52°S, Patagonia) indicates the global Last Glacial Maximum (LGM) extent was half of the local LGM

Whether glaciers in the southern hemisphere were asynchronous to those in the north during the global Last Glacial Maximum (gLGM; 26.5–19 ka) is still debated. In Patagonia (South America), numerous ice lobes attained their maximum extents during diverse episodes of the last glacial cycle. To understand the variations in the gLGM vs local LGM record, we studied the Lago Aníbal Pinto area (52°00′ S, 72°40′ E; Chile), where several moraine belts were deposited by one of the eastward-flowing southern Patagonian Ice Sheet (PIS) outlet ice lobes; the Última Esperanza. We report eight 10Be terrestrial cosmogenic nuclides (TCN) surface ages from granitic moraine boulders. Our weighted average age obtained from the southern part of the Río Turbio moraine belt yields 50.7 ± 2.4 ka (oldest boulder age; 53.8 ± 5.3 ka) and confirms the greatest extent of the local Last Glacial Maximum (lLGM) during Marine Isotope Stage 3 (MIS 3) in the previously dated northern moraines from the same belt. Our 10Be TCN age (32.6 ± 2.2 ka) derived from the Dos Lagunas moraine, which makes up the northernmost margin of the Última Esperanza ice lobe's Arauco advance, also validates the MIS 3 timing. Following the formation of Arauco moraines, the Última Esperanza ice lobe was split into three main tributaries in the south, which formed three restricted and previously undated moraine complexes. We dated one of them, the Aníbal Pinto moraine complex. Whereas the highest moraine yielded the oldest age (28.3 ± 2.2 ka), lower moraine's surface that was later truncated as a lacustrine erosional platform, yielded younger boulder ages (weighted average age = 18.9 ± 1.0 ka; oldest boulder age = 19.0 ± 1.5 ka), indicating they were deposited under the Pinto Lake level. Our new age data allow us to propose a new chronology for the Aníbal Pinto moraine complex and consolidate previously published ages from other moraine belts. We attribute the Aníbal Pinto moraine complex to early gLGM (MIS 2), emphasising that the gLGM was half the extent of lLGM in the Última Esperanza ice lobe that underpins its interhemispheric asynchronous character.

The Chronology of Early Human Settlement in Three Gorges Region, China—Contribution of Coupled Electron Spin Resonance and Uranium-Series Dating Method

The Three Gorges region (TGR) located in the geographic center of China, is a transition zone between mountain and plain areas, and a probable migration corridor for hominins and other mammals between South and North China. Detailed chronological information of paleoanthropological evidence in this area could help us better understand the human evolution in East Asia. The OSL and U-series dating methods are two conventional dating methods generally adopted to date such sites; however, their applications were limited by the dating range—restricted to several hundred of millennia and ambiguous stratigraphic relationship between the archaeological remains and the dating target materials. Cosmogenic nuclide burial dating of quartzite stone artifacts and coupled electron spin resonance and uranium series dating (ESR/U-series) of fossil teeth have the potential to date Early–Middle Pleistocene hominin sites in Asia and were applied increasingly in China in recent years. However, the application of cosmogenic 26Al/10Be burial dating is limited in TGR because most sites are dominated by limestone, leading to the scarcity of the quartz component. In this case, the coupled ESR/U-series method plays a more important role in the establishment of the chronology of human settlement. In TGR, by using the coupled ESR/U-series method, we have dated seven important Early and Middle Pleistocene hominin settlement sites, including Longgupo, Jianshi, Yunxian, Meipu, Bailongdong, Changyang, and Yumidong sites. Based on our dating results, we propose that hominins were settled in TGR probably from the early stage of Early Pleistocene (∼2.5−2.2 Ma) at the Longgupo site to the late Middle Pleistocene to Late Pleistocene of the Yumidong site (∼274−14 ka) and very likely to spread to other parts of East Asia during this time period. In view of the potential of coupled ESR/U-series dating on fossil teeth from the hominin sites in the TGR, future work may consider the micro damage or non-destructive analysis of enamel fragment with the ESR method and laser ablation ICP-MS techniques that will make possible the direct dating of precious human fossils in China.

Cosmogenic nuclide dating of two stacked ice masses: Ong Valley, Antarctica

Abstract. We collected a debris-rich ice core from a buried ice mass in Ong Valley, located in the Transantarctic Mountains in Antarctica. We measured cosmogenic nuclide concentrations in quartz obtained from the ice core to determine the age of the buried ice mass and infer the processes responsible for the emplacement of the debris currently overlaying the ice. Such ice masses are valuable archives of paleoclimate proxies; however, the preservation of ice beyond 800 kyr is rare, and therefore much effort has been recently focused on finding ice that is older than 1 Myr. In Ong Valley, the large, buried ice mass has been previously dated at > 1.1 Ma. Here we provide a forward model that predicts the accumulation of the cosmic-ray-produced nuclides 10Be, 21Ne, and 26Al in quartz in the englacial and supraglacial debris and compare the model predictions to measured nuclide concentrations in order to further constrain the age. Large downcore variation in measured cosmogenic nuclide concentrations suggests that the englacial debris is sourced both from subglacially derived material and recycled paleo-surface debris that has experienced surface exposure prior to entrainment. We find that the upper section of the ice core is 2.95 + 0.18 / −0.22 Myr old. The average ice sublimation rate during this time period is 22.86 + 0.10 / −0.09 m Myr−1, and the surface erosion rate of the debris is 0.206 + 0.013 / −0.017 m Myr−1. Burial dating of the recycled paleo-surface debris suggests that the lower section of the ice core belongs to a separate, older ice mass which we estimate to be 4.3–5.1 Myr old. The ages of these two stacked, separate ice masses can be directly related to glacial advances of the Antarctic ice sheet and potentially coincide with two major global glaciations during the early and late Pliocene epoch when global temperatures and CO2 were higher than present. These ancient ice masses represent new opportunities for gathering ancient climate information.

Fluvial evolution in a growing thrust-fold range of the Yumu Shan, NE Tibetan Plateau

In an active thrust-fold belt, the drainage evolution, such as the alternation between transverse and longitudinal rivers, will influence the surface erosion within the range and the sedimentary sequence in the foreland basin. However, knowledge of the controls of changes in transverse or longitudinal drainage patterns in an uplifting range is limited. The Yumu Shan, on the northeastern margin of the Tibetan Plateau, is an example for a young and active thrust-fold range that shows evidence for drainage reversals since ∼3 Ma. The materials we used are the late Cenozoic depositional sequences around the Yumu Shan. Analyses of the lithology, grain size, roundness, and clast orientation of gravel deposits on the northern and southern flanks, combined with the results of cosmogenic nuclide burial dating, reveal that a transverse river, the Liyuan River, crossed the Yumu Shan range during 3.0–1.8 Ma by eroding weak Cenozoic sedimentary rocks. After 1.8 Ma, however, this transverse river was deflected and flowed along the edge of the range, and a longitudinal river developed along the back flank of the range. This process was synchronous with the exposure of resistant Paleozoic rocks by the continuous uplift of the range. After comparison with other potential controlling factors for drainage reversals, such as variations in rock uplift, climate or upstream aggradation, we find that the change in bedrock resistance was the most likely control on the change from a transverse to a longitudinal drainage pattern in this relatively young thrust-fold belt. Furthermore, through comparison with a series of active fold-thrust belts worldwide, we propose a simple conceptual model for the drainage evolution of active fold-thrust belts. In the earlier stage of the uplift, the exposed soft sedimentary rocks will cause transverse (antecedent) rivers to cross the range. After the continuous uplift of the range, resistant rocks cause the deflection of a transverse river and its change to a longitudinal river along the back flank of the range. In the final stage, after the several million years to around ten million years of uplift, with the development of a mature range, the local base level of the foreland will greatly fall, promoting a return to a transverse drainage pattern by headward erosion and capture of longitudinal rivers.

230Th dating of gypsum from lacustrine, brackish-marine and terrestrial environments

Gypsum is an important environmental archive in hyperarid settings, but its application in earth sciences has been hampered by the limited availability of suitable dating methods. Here we re-evaluate the potential of 230Th/U dating for sedimentary gypsum as well as anhydrite samples in different depositional environments. We provide a robust analytical protocol based on a simple dissolution using nitric acid followed by a two-column separation procedure using anion exchange resin. Isotope analyses were performed using MC-ICPMS and a mixed 229Th-233U-236U tracer. We applied the method to a suite of samples that likely reached secular equilibrium as well as to three suites of younger (<300 ka old) samples, two from the Atacama Desert in Chile and one from a drill core of a saline playa lake in southern Spain. We employed a selective sampling approach targeting the least visually altered and clearest parts of the samples largely devoid of detritus and analysed multiple subsamples.The results demonstrate the general applicability of the 230Th/U method to gypsum as well as anhydrite, and recrystallization of gypsum to anhydrite does not significantly affect the U-series isotope system. However, there is evidence for open-system behaviour, which needs to be carefully assessed. Our results suggest that α-recoil processes as well as uptake or loss of uranium are important processes biasing the results of gypsum dating. Furthermore, we see signs of multiple-phase contamination in our case studies involving detrital components and authigenic phases. Based on these results we propose geochemical criteria in to identify biased 230Th/U gypsum ages.The application of our protocol to a suite of selenite samples deposited by brine outflows on the lowest terrace of Cerro Soledad (Atacama Desert, Chile) provides a mean age of 212 ± 8 ka, indicating a relatively humid climate interval with local rainfall in this area that is coinciding with marine isotope stage 7C. Selenites deposited in the paleo-lake Soledad system (Atacama Desert, Chile) yield 230Th/U ages ranging from 250 to 320 ka, in agreement with previous age estimates for the drainage of this paleo-lake obtained by cosmogenic nuclide dating. Application to a suite of gypsum crystals recovered from a drill core in the Laguna de Fuente de Piedra (South Spain) shows good agreement with their stratigraphic context, however, multiple-phase contamination produces larger age uncertainties for these ages.

Cosmogenic nuclide weathering biases: corrections and potential for denudation and weathering rate measurements

Abstract. Cosmogenic radionuclides (CRNs) are the standard tool to derive centennial-to-millennial timescale denudation rates; however, it has been demonstrated that chemical weathering in some settings can bias CRNs as a proxy for landscape denudation. Currently, studies investigating CRN weathering biases have mostly focused on the largely insoluble target mineral quartz in felsic lithologies. Here, we examine the response of CRN build-up for both soluble and insoluble target minerals under different weathering scenarios. We assume a simple box model in which bedrock is converted to a well-mixed regolith at a constant rate, and denudation occurs by regolith erosion and weathering either in the regolith or along the regolith–bedrock interface, as is common in carbonate bedrock. We show that weathering along the regolith–bedrock interface increases CRN concentrations compared to a no-weathering case and how independently derived weathering rates or degrees can be used to correct for this bias. If weathering is concentrated within the regolith, insoluble target minerals will have a longer regolith residence time and higher nuclide concentration than soluble target minerals. This bias can be identified and corrected using paired-nuclide measurements of minerals with different solubility coupled with knowledge of either the bedrock or regolith mineralogy to derive denudation and long-term weathering rates. Similarly, single-nuclide measurements on soluble or insoluble minerals can be corrected to determine denudation rates if a weathering rate and compositional data are available. Our model highlights that for soluble target minerals, the relationship between nuclide accumulation and denudation is not monotonic. We use this understanding to map the conditions of regolith mass, weathering, and denudation rates at which weathering corrections for cosmogenic nuclides become large and ambiguous, as well as identify environments in which the bias is mostly negligible and CRN concentrations reliably reflect landscape denudation. We highlight how measurements of CRNs from soluble target minerals, coupled with bedrock and regolith mineralogy, can help to expand the range of landscapes for which centennial-to-millennial timescale denudation and weathering rates can be obtained.

Cosmogenic nuclide and solute flux data from central Cuban rivers emphasize the importance of both physical and chemical mass loss from tropical landscapes

Abstract. We use 25 new measurements of in situ produced cosmogenic 26Al and 10Be in river sand, paired with estimates of dissolved load flux in river water, to characterize the processes and pace of landscape change in central Cuba. Long-term erosion rates inferred from 10Be concentrations in quartz extracted from central Cuban river sand range from 3.4–189 Mg km−2 yr−1 (mean 59, median 45). Dissolved loads (10–176 Mg km−2 yr−1; mean 92, median 97), calculated from stream solute concentrations and modeled runoff, exceed measured cosmogenic-10Be-derived erosion rates in 18 of 23 basins. This disparity mandates that in this environment landscape-scale mass loss is not fully represented by the cosmogenic nuclide measurements. The 26Al / 10Be ratios are lower than expected for steady-state exposure or erosion in 16 of 24 samples. Depressed 26Al / 10Be ratios occur in many of the basins that have the greatest disparity between dissolved loads (high) and erosion rates inferred from cosmogenic nuclide concentrations (low). Depressed 26Al / 10Be ratios are consistent with the presence of a deep, mixed, regolith layer providing extended storage times on slopes and/or burial and extended storage during fluvial transport. River water chemical analyses indicate that many basins with lower 26Al / 10Be ratios and high 10Be concentrations are underlain at least in part by evaporitic rocks that rapidly dissolve. Our data show that when assessing mass loss in humid tropical landscapes, accounting for the contribution of rock dissolution at depth is particularly important. In such warm, wet climates, mineral dissolution can occur many meters below the surface, beyond the penetration depth of most cosmic rays and thus the production of most cosmogenic nuclides. Our data suggest the importance of estimating solute fluxes and measuring paired cosmogenic nuclides to better understand the processes and rates of mass transfer at a basin scale.

Late Cenozoic locally landslide-dammed lakes across the Middle Yangtze River

Understanding the drainage history of the Middle Yangtze River is critical for determining its interaction of tectonics, climate, and denudation. Previous studies have focused on interpreting the provenance of the Xigeda lacustrine sediments, however, the poorly constrained ages of these sediments mean that the evolutionary history remains elusive. Here, we present new cosmogenic 10Be, 21Ne and 26Al data from the Dadu and Anning river sediments, indicating deposition ages of paleo-lake sediments between ~1 and ~0.5 Ma. Paleo-denudation rates together with neon isotopic composition and topographic analysis suggest that there was no connection between the Dadu and Anning rivers during the Middle Pleistocene. Combined with field investigations and cosmogenic nuclides ages, we suggest that paleo-lakes formed upstream of channel blocking landslides, which in turn could have been the result of external forcings, such as earthquakes, debris flows, increased precipitation, and glaciation. Previous thermochronology data constrain the fluvial incision of the paleo-Dadu-Anning River to have enhanced between ~13–9 Ma. Subsequently, our cosmogenic results show paleo-lakes were formed (~1 Ma) after the paleo-Dadu-Anning River capture (~2 Ma). Finally, the draining of paleo-lakes occurred at ~0.5 Ma.

Regional Differences in Carbon-14 Data of the 993 CE Cosmic Ray Event

Cosmogenic nuclides such as 14C from tree rings and 10Be and 36Cl from ice cores are excellent proxies for the past extremely large solar energetic particle (SEP) events, which are dozens of times larger than the largest SEP event in the history of observation. So far, several rapid 14C increases have been discovered, which are considered to have originated from extreme SEP events (or set of successive SEP events) from verifications using multiple cosmogenic nuclide analyses in natural archives. Although these events are characterized by a rapid increase in cosmogenic nuclide concentrations, 14C data recorded worldwide do not always show similar variations, especially during the 993 CE event, where a rapid increase was recorded in either 992–993 CE or 993–994 CE in several records. We present new 14C data of the Japanese cedar sample for the 993 CE event. Although the latest data show no significant increase in 1 year, an overall increase pattern is consistent with the previously reported 14C data of the Japanese cedar, which supports that a significant 14C increase occurred from 993 to 994 CE in the Japanese sample. Given the dominant 14C production in high latitudes by SEPs, the difference in timing of increase may be a transport effect in the atmosphere. Moreover, the difference in the timing of the 14C increase can cause a 1-year age-determination error using the 993 CE radiocarbon spike. Compared with the 14C data between tree samples from high latitude and midlatitude, including Japan, high-latitude data can capture 14C changes originating from SEP events more quickly and clearly and may be more suitable for a SEP event exploration in the past.

A question of time: historical overview and modern thought on Quaternary dating methods to produce fluvial chronologies

Since its early stage one century ago, Quaternary geochronology has witnessed major progress with the development of a vast array of relative, correlation and numerical dating methods. Fluvial geomorphologists have simultaneously striven to quantitatively assess processes on very different timescales (101-106 years). Dating fluvial/alluvial landforms and/or associated material/archives rapidly turned out to be one of the main challenges accordingly. Here, a bipartite reflection is proposed on how Quaternary dating methods have allowed the production of fluvial chronologies over the last century, including at the present time. Firstly, a progressive methodological shift has been witnessed. Whilst age controls were mostly attained via relative and correlation dating until the late xxth century, numerical dating methods have increasingly taken over the prevailing role since the last two or three decades. Two explanatory causes are then proposed. The bipartite evolution of relative dating methods over the last decades is explored: whilst some methods were further developed in the fluvial context (e.g., amino-acid racemization), some others were vehemently questioned (e.g., lichenometry). The late flourishing of radiometric dating methods during the second half of the xxth century is also considered. Secondly, the current status of Quaternary geochronology is doubly questioned: (i) Is the die cast for “disused” relative dating methods? and (ii) Can radiometric dating methods solve all chronological issues related to fluvial environments? Surface exposure dating of alluvial surfaces via cosmogenic nuclides is then thoroughly discussed in this respect, focusing on some methodological issues. A rehabilitation of some methods fallen into disuse, such as weathering indicators, is then proposed to support numerical age outputs. Finally, geochronological approaches based on multiple dating tools, either by crosschecking the outputs of different methods or coupling numerical dating methods, are fostered.

Modeling the Transport and Deposition of 10Be Produced by the Strongest Solar Proton Event During the Holocene

AbstractProminent excursions in the number of cosmogenic nuclides (e.g., 10Be) around 774 CE/775 document the most severe solar proton event (SPE) throughout the Holocene. Its manifestation in ice cores is valuable for geochronology, but also for solar‐terrestrial physics and climate modeling. Using the ECHAM/MESSy Atmospheric Chemistry (EMAC) climate model in combination with the Warning System for Aviation Exposure to SEP (WASAVIES), we investigate the transport, mixing, and deposition of the cosmogenic nuclide 10Be produced by the 774 CE/775 SPE. By comparing the model results to the reconstructed 10Be time series from four ice core records, we study the atmospheric pathways of 10Be from its stratospheric source to its sink at Earth's surface. The reconstructed post‐SPE evolution of the 10Be surface fluxes at the ice core sites is well captured by the model. The downward transport of the 10Be atoms is controlled by the Brewer‐Dobson circulation in the stratosphere and cross‐tropopause transport via tropopause folds or large‐scale sinking. Clear hemispheric differences in the transport and deposition processes are identified. In both polar regions the 10Be surface fluxes peak in summertime, with a larger influence of wet deposition on the seasonal 10Be surface flux in Greenland than in Antarctica. Differences in the peak 10Be surface flux following the 774 CE/775 SPE at the drilling sites are explained by specific meteorological conditions depending on the geographic locations of the sites.

Cosmogenic nuclide dating of Australopithecus at Sterkfontein, South Africa

Sterkfontein is the most prolific single source of Australopithecus fossils, the vast majority of which were recovered from Member 4, a cave breccia now exposed by erosion and weathering at the landscape surface. A few other Australopithecus fossils, including the StW 573 skeleton, come from subterranean deposits [T. C. Partridge etal., Science 300, 607-612 (2003); R. J. Clarke, K. Kuman, J. Hum. Evol. 134, 102634 (2019)]. Here, we report a cosmogenic nuclide isochron burial date of 3.41 ± 0.11 million years (My) within the lower middle part of Member 4, and simple burial dates of 3.49 ± 0.19 My in the upper middle part of Member 4 and 3.61 ± 0.09 My in Jacovec Cavern. Together with a previously published isochron burial date of 3.67 ± 0.16 My for StW 573 [D. E. Granger etal., Nature 522, 85-88 (2015)], these results place nearly the entire Australopithecus assemblage at Sterkfontein in the mid-Pliocene, contemporaneous with Australopithecus afarensis in East Africa. Our ages for the fossil-bearing breccia in Member 4 are considerably older than the previous ages of ca. 2.1 to 2.6 My interpreted from flowstones associated with the same deposit. We show that these previously dated flowstones are stratigraphically intrusive within Member 4 and that they therefore underestimate the true age of the fossils.

Long-term denudation rate of karstic North Dalmatian Plain (Croatia) calculated from 36Cl cosmogenic nuclides

The development of a karst landscape results from complex interactions between lithology, climate, hydrology, soil, vegetation and tectonics. Weathering and erosion of carbonate rocks leads to denudation of karst landscapes. As dissolution of carbonate rocks is often considered to be the main process governing carbonate weathering, other processes are often overlooked. Here we present research done in the North Dalmatian Plain, a carbonate erosive surface located in the Dinaric karst region. Although the study site is composed of two different carbonate lithologies having different weathering style, there is no evident lithological impact on the topography of the erosive surface. Analyses of 36Cl concentration were performed in ten proximal bedrock samples from both lithologies and resulted in long-term denudation rates ranging from 14.7 to 22.8 m/Ma. Since no statistical significance was found between samples from different lithologies (all samples belong to a single normal population) and they have the same geomorphological context and climate features, variable denudation rates are attributed to local (sample specific) differences. In the study site there are no large outstanding rock residuals or patches with deep soil profiles. Thus, to maintain the levelled erosive surface, local differential denudation rates have to vary with time. We hypothesize that lichens and pedogenic carbonates have a significant role in modulating local differences in denudation rates. Our study shows that even at outcrop scale those differences can be significant, and the study of sample populations is preferred to single or limited number of analyses. Thus, the long-term denudation rate of the North Dalmatian Plain, including its local variability is 18.91 ± 0.81 m/Ma. Despite classical studies on karst terrains assume that dissolution is the main process responsible for development of these landscapes, our research highlights the importance of physical weathering in combination with dissolution of carbonates as main controls on the denudation of karst landscapes.

Alpine rock glacier activity over Holocene to modern timescales (western French Alps)

Abstract. Active rock glaciers are some of the most frequent cryospheric landforms in midlatitude high-elevation mountain ranges. Their activity strongly influences the hydrology and geomorphology of alpine environments over short (years to decades) and long (centuries to millennia) timescales. Being conspicuous expressions of mountain permafrost and important water reserves in the form of ground ice, rock glaciers are seen as increasingly important actors in the geomorphological and hydrological evolution of mountain systems, especially in the context of current climate change. Over geological timescales, rock glaciers both reflect paleoclimate conditions and transport rock boulders produced by headwall erosion, and they therefore participate in shaping high mountain slopes. However, the dynamics of rock glaciers and their evolution over different timescales remain under-constrained. In this study, we adopt a multi-method approach, including field observations, remote sensing, and geochronology, to investigate the rock glacier system of the Vallon de la Route (Combeynot Massif, western French Alps). Remotely sensed images and correlation techniques are used to document the displacement field of the rock glacier over timescales ranging from days to decades. Additionally, to estimate displacement over periods from centuries to millennia, we employ terrestrial cosmogenic nuclide (quartz 10Be) surface-exposure dating on rock boulder surfaces located along the central flow line of the rock glacier, targeting different longitudinal positions from the headwall to the rock glacier terminus. The remote sensing analysis demonstrates that between 1960 and 2018 the two lower units of the rock glacier were motionless, the transitional unit presented an integrated surface velocity of 0.03±0.02 m a−1, and the two upper active units above 2600 m a.s.l. showed a velocity between 0.14±0.08 and 0.15±0.05 m a−1. Our results show 10Be surface-exposure ages ranging from 13.10±0.51 to 1.88±0.14 ka. The spatial distribution of dated rock glacier boulders reveals a first-order inverse correlation between 10Be surface-exposure age and elevation and a positive correlation with horizontal distance to the headwall. These observations support the hypothesis of rock boulders falling from the headwall and remaining on the glacier surface as they are transported down valley, and they may therefore be used to estimate rock glacier surface velocity over geological timescales. Our results also suggest that the rock glacier is characterized by two major phases of activity. The first phase, starting around 12 ka, displays a 10Be age gradient with a rock glacier surface velocity of about 0.45 m a−1, following a quiescent period between ca. 6.2 and 3.4 ka before the emplacement of the present-day upper two active units. Climatic conditions have favored an integrated rock glacier motion of around 0.18 m a−1 between 3.4 ka and present day. These results allow us to quantify back-wearing rates of the headwall of between 1.0 and 2.5 mm a−1, higher than catchment-integrated denudation rates estimated over millennial timescales. This suggests that the rock glacier system promotes the maintenance of high rock wall erosion by acting as debris conveyor and allowing freshly exposed bedrock surfaces to be affected by erosion processes.

The impact of different atmospheric phenomena to cosmogenic 22Na/7Be ratio

The isotopic ratio of two cosmogenic isotopes, 22Na and 7Be, has been recognized as a possible tracer for atmospheric dynamics. We studied how different atmospheric phenomena affect the 22Na/7Be ratio that are known to affect the ground level 7Be concentrations. The 22Na/7Be ratio was studied in the following context: estimation of aerosol residence time, in the intrusions of stratospheric air after Sudden Stratospheric Warming (SSW) events, seasonal variations due to different atmospheric mixing conditions and finally the interannual variations caused by the large-scale atmospheric circulation such as North Atlantic Oscillation (NAO), Arctic Oscillation (AO), Scandinavian Pattern (SCAND) and Quasi-Biennial Oscillation (QBO). This study utilized airborne radioactivity monitoring data from three Finnish stations located in Rovaniemi (66.51 °N, 25.73 °E), Kajaani (64.08 °N, 27.71 °E) and Kotka (60.48 °N, 26.92 °E). The theoretical calculations were performed to see how the 22Na/7Be ratio changes as a function of altitude during when produced by the galactic cosmic rays. The calculations showed that the 22Na/7Be ratio is higher in the stratosphere and lower in the troposphere. At mid-stratosphere level the ratio was 0.001, at tropopause level approximately 0.0007 and at ground level approximately 0.00045. Due to the different half-lives of 22Na and 7Be the mean atmospheric residence time can be estimefrom the radioactive decay law by assuming an initial activity concentration ratio of 0.0007 and using the measured ground level 22Na/7Be ratios. This resulted in the median atmospheric residence time as 91–95 days. The atmospheric impact on 22Na/7Be ratio was studied at different time scales starting from short-term impact from SSW, the known seasonal cycle and variations at interannual timescale. SSW events caused rapid increases in the ground level 7Be and 22Na concentrations while the impact on 22Na/7Be ratio was relatively small. A clear seasonal cycle in the 22Na/7Be ratio was observed where the 22Na/7Be ratio was at low value during November–February, increased during March–May, stayed at stable high value during June–July and steadily decreased during August–October. The long-term average for 22Na/7Be ratio was established to be 1.1 × 10−4 by using the Kotka data from 1993 to 2020. During 1995–2009, the 22Na/7Be ratio was below the long-term average, and during 2010–2020, generally above. The interannual variations in the 22Na/7Be ratio were studied by using wavelet transform analysis of the Kotka data. The analysis showed 4-6-year and 11-year periodicities. The wavelet coherence analyses with AO, NAO, SCAND, QBO and Southern Oscillation (SO) indices showed that the shorter periodicity was caused by the NAO/AO effects the atmospheric circulation in Scandinavia and hence the transport of the isotopes from upper atmosphere to ground level. The observed coherence with the Southern Oscillation Index (SOI) starting from 2012 onwards appeared at an approximately 3–4year scale while the AO and NAO had a coherence at an approximately 6-year scale. The observed 11-year periodicity was caused by the Solar cycle which modulates the 7Be and 22Na production in the upper atmosphere.

Technical note: Accelerator mass spectrometry of &amp;lt;sup&amp;gt;10&amp;lt;/sup&amp;gt;Be and &amp;lt;sup&amp;gt;26&amp;lt;/sup&amp;gt;Al at low nuclide concentrations

Abstract. Accelerator mass spectrometry (AMS) is currently the standard technique to measure cosmogenic 10Be and 26Al concentrations, but the challenge with measuring low nuclide concentrations is to combine high AMS measurement efficiency with low backgrounds. The current standard measurement setup at ANSTO uses the 3+ charge state with Ar stripper gas at 6 MV for Be and 4 MV for Al, achieving ion transmission through the accelerator for 10Be3+ and 26Al3+ of around 35 % and 40 %, respectively. Traditionally, 26Al measurement uncertainties are larger than those for 10Be. Here, however, we show that 26Al can be measured to similar precision as 10Be even for samples with 26Al / 27Al ratios in the range of 10−15, provided that measurement times are sufficiently long. For example, we can achieve uncertainties of 5 % for 26Al / 27Al ratios around 1×10-14, typical for samples of late Holocene age or samples with long burial histories. We also provide empirical functions between the isotope ratio and achievable measurement precision, which allow predictive capabilities for future projects and serve as a benchmark for inter-laboratory comparisons. For the smallest signals, not only is understanding the source of 10Be or 26Al background events required to select the most appropriate blank correction method but also the impact of the data reduction algorithms on the obtained nuclide concentration becomes pronounced. Here we discuss approaches to background correction and recommend quality assurance practices that guide the most appropriate background correction method. Our sensitivity analysis demonstrates a 30 % difference between different background correction methods for samples with 26Al / 27Al ratios below 10−14. Finally, we show that when the measured signal is small and the number of rare isotope counts is also low, differing 26Al or 10Be concentrations may be obtained from the same data if alternate data reduction algorithms are used. Differences in the resulting isotope concentration can be 50 % or more if only very few (≲ 10) counts were recorded or about 30 % if single measurement is shorter than 10 min. Our study presents a comprehensive method for analysis of cosmogenic 10Be and 26Al samples down to isotope concentrations of a few thousand atoms per gram of sample, which opens the door to new and more varied applications of cosmogenic nuclide analysis.

Late Pleistocene glacial advances, equilibrium-line altitude changes and paleoclimate in the Jakupica Mts (North Macedonia)

In the Jakupica Mts a plateau glacier was reconstructed (max. area ∼45 km2, max thickness: ∼260 m). The study area comprises six formerly glaciated valleys, five of which were fed by the plateau glacier and one had an independent cirque when local glaciation reached its maximum ice extent (MIE). The equilibrium line altitude (ELA) of the most extended glacial phase was at 2075+37/-25 m asl. The 10Be cosmic ray exposure (CRE) age of this phase was estimated at 19.3+1.7/-1.3 ka, conformable with the Last Glacial Maximum (LGM). CRE ages from the next moraine generation placed the first phase of deglaciation to 18.2+1.0/-3.0 ka. The samples from the moraine of the penultimate deglaciation phase provided CRE ages with large scatter and biased towards old ages, which is probably the result of inherited cosmogenic nuclide concentrations within the rock.Glacio-climatological modelling was performed for the MIE, which has a well-established LGM age. The degree-day model was used to calculate the amount of accumulation required to sustain the glaciological equilibrium assuming a certain temperature drop at the ELA. The degree-day model constrained by the pollen-based July paleo-temperature reconstructions yielded an annual total melt at the LGM ELA comparable to or slightly higher than the current mean annual precipitation at the same elevation. These wetter LGM conditions inferred from the paleo-glaciological evidence in Jakupica Mts suggest an enhanced moisture advection in the region.

Reconciling the apparent absence of a Last Glacial Maximum alpine glacial advance, Yukon Territory, Canada, through cosmogenic beryllium-10 and carbon-14 measurements

Abstract. We present a new in situ produced cosmogenic beryllium-10 and carbon-14 nuclide chronology from two sets (outer and inner) of alpine glacier moraines from the Grey Hunter massif of southern Yukon Territory, Canada. The chronology of moraines deposited by alpine glaciers outside the limits of the Last Glacial Maximum (LGM) ice sheets potentially provides a less-ambiguous archive of mass balance, and hence climate, than can be inferred from the extents of ice sheets themselves. Results for both nuclides are inconclusive for the outer moraines, with evidence for pre-LGM deposition (beryllium-10) and Holocene deposition (carbon-14). Beryllium-10 results from the inner moraine are suggestive of canonical LGM deposition but with relatively high scatter. Conversely, in situ carbon-14 results from the inner moraines are tightly clustered and suggestive of terminal Younger Dryas deposition. We explore plausible scenarios leading to the observed differences between nuclides and find that the most parsimonious explanation for the outer moraines is that of pre-LGM deposition, but many of the sampled boulder surfaces were not exhumed from within the moraine until the Holocene. Our results thus imply that the inner and outer moraines sampled pre- and post-date the canonical LGM and that moraines dating to the LGM are lacking likely due to overriding by the subsequent Late Glacial/earliest Holocene advance.

Large-energy single hits at JUNO from atmospheric neutrinos and dark matter

Large liquid scintillator detectors, such as JUNO, present a new opportunity to study neutral current events from the low-energy end of the atmospheric neutrinos, and possible new physics signals due to light dark matter. We carefully study the possibility of detecting ``large-energy singles'' (LES), i.e., events with visible scintillation energy $&gt;15\text{ }\text{ }\mathrm{MeV}$, but no other associated tags. For an effective exposure of $20\text{ }\text{ }\mathrm{kton}\ensuremath{-}\mathrm{yr}$ and considering only Standard Model physics, we expect the LES sample to contain $\ensuremath{\sim}40$ events from scattering on free protons and $\ensuremath{\sim}108$ events from interaction with carbon, from neutral-current interactions of atmospheric neutrinos. Backgrounds, largely due to $\ensuremath{\beta}$ decays of cosmogenic isotopes, are shown to be significant only below 15 MeV visible energy. The LES sample at JUNO can competitively probe a variety of new physics scenarios, such as boosted dark matter and annihilation of galactic dark matter to sterile neutrinos.