Cosmogenic Radionuclides Research Articles

Rapid Uplift, yet Slow Denudation of the Suckling‐Dayman Metamorphic Core Complex in Tropical Papua New Guinea

AbstractA linear relationship between rates of physical erosion and chemical weathering is apparent in slowly eroding landscapes. Whether the relationship remains linear in rapidly eroding landscapes is less clear. Field‐based research into this relationship between erosion and weathering rates has largely been conducted in temperate climates with granitic bedrock. In tropical settings, the contribution of chemical mass loss to total denudation may approach, or even exceed, that of erosion. We report 10Be‐in‐quartz cosmogenic radionuclide and soil chemistry data from the Suckling‐Dayman Metamorphic Core Complex (SDMCC) in Papua New Guinea. Despite being exhumed at cm‐per‐year rates, its lower‐plate domed and striated morphology suggests minimal denudation, which is confirmed by our 10Be‐in‐quartz data (0.02–0.18 mm/yr). We suggest that rolling hinge‐style back‐rotation of the SDMCC's lower plate and the combination of a tropical climate and highly weatherable metabasalt bedrock have played a fundamental role in preserving the tectonic topography of this remarkable metamorphic core complex.

Cosmogenic Radionuclide Beryllium-7 Skewed Data Preprocessing for Northeast Monsoon Forecasting in Malaysia

Cosmogenic Radionuclide Beryllium-7 Skewed Data Preprocessing for Northeast Monsoon Forecasting in Malaysia

Cosmogenic radionuclides in meteorites from the Otway Massif blue ice area, Antarctica: An unusual, well‐preserved H5 chondrite strewn field

AbstractThe US Antarctic Search for Meteorites (ANSMET) discovered a dense cluster of 88 ordinary chondrites with a total mass of more than 100 kg on a blue ice area (BIA) of 1.6 × 0.3 km2 near the Otway Massif, Grosvenor Mountains, Antarctica. The larger masses (weighing up to 29 kg) were found at one end of an oval‐shaped pattern and the smaller masses (50–200 g) at the other end. We measured concentrations of the cosmogenic radionuclides 10Be (half‐life—1.36 × 106 year) and 36Cl (3.01 × 105 year) in the metal fraction of 17 H chondrites, including 14 fragments of this cluster, to verify the hypothesis that this meteorite cluster on the Otway Massif BIA represents a meteorite strewn field produced by the atmospheric breakup of a single meteoroid. The 10Be and 36Cl concentrations confirm that 10 out of 14 H chondrites from different locations within this small area are paired fragments of the same meteorite fall, while the four other H chondrites represent two additional—smaller—falls. The radionuclides suggest a pre‐atmospheric mass of 200–400 kg for the large pairing group, suggesting that 25%–50% of the meteoroid survived atmospheric entry. Based on the distribution of the paired H chondrites and evidence of their common cosmic‐ray exposure history in space, we conclude that most of the 88 meteorites within this small area represent a meteorite strewn field. The small size of the strewn field suggests that the meteoroid entered at a steep angle (>60°), while the low amount of fusion crust on most meteorite surfaces most likely indicates atmospheric break up at low altitude, while additional fragmentation of a large surviving fragment may have occurred during impact on the ice. This well‐documented strewn field provides a good opportunity to apply model simulations of the atmospheric fragmentation of this object as a function of entry angle, velocity, and meteoroid strength. Cosmogenic 14C analyses in two members of the Otway Massif pairing group yield a terrestrial age of 15.5 ± 1.5 kyr, which represents the time elapsed since this meteorite fell on Earth. The excellent preservation of an Antarctic meteorite strewn field suggests that the Otway Massif BIA represents a relatively stagnant blue ice field.

Spectral deconvolution approach for Beryllium-7 detection: Evaluating soil redistribution with NaI(Tl) detectors

Beryllium-7 is a cosmogenic radionuclide produced in the upper atmosphere through the spallation of cosmic rays and falls to Earth mainly by wet deposition. Given its short half-life (53 days) and strong adsorption to the clay present in the soil, 7Be inventory can be evaluated to obtain information about soil redistribution over periods of up to 6 months. By performing these measurements, one can infer the efficiency of soil management and conservation procedures. The isotope presence is measured using gamma-ray spectrometry to detect its decay emission, specifically a photon with an energy of 477.6 keV. Since other isotopes found in soil emit photons with similar energies, it is common to use High-Purity Germanium (HPGe) detectors due to their high energy resolution. However, these detectors are expensive and require cryogenic cooling during operation. An alternative is the Thallium-Activated Sodium Iodide (NaI(Tl)) detector, which is less costly and does not require cryogenic cooling, but has poorer resolution and cannot separate the emissions from 228Ac (463.0 keV), 7Be, and 208Tl (510.7 keV). This work aims to apply spectral deconvolution to quantify 7Be in soil samples collected from an agricultural mega-parcel situated on a hill slope and evaluate the soil redistribution experienced within. The activities of the interfering isotopes were quantified using other emissions, and the results were used to calculate the corresponding area in the desired region of the spectra. From this, 7Be activities were calculated and compared with the HPGe results to validate the deconvolution technique. The calculated activities for the three radioisotopes showed similar correspondence between the two detectors, indicating that the applied method is adequate for evaluating soil redistribution at lower costs. The calculated net soil deposition was 36.4 ± 6.4 kgm−2 for the NaI(Tl) detector and 26.4 ± 7.9 kgm−2 for the HPGe detector during the evaluated period. Furthermore, these results suggest that a similar approach can be applied to quantify other isotopes using NaI(Tl) detectors as an alternative to HPGe detectors.

Holocene solar activity inferred from global and hemispherical cosmic-ray proxy records

Variations in solar activity have been proposed to play an important role in recent and past climate change. To study this link on longer timescales, it is essential to know how the Sun has varied over the past millennia. Direct observations of solar variability based on sunspot numbers are limited to the past 400 years, and beyond this we rely on records of cosmogenic radionuclides, such as 14C and 10Be in tree rings and ice cores. Their atmospheric production rates depend on the flux of incoming galactic cosmic rays, which is modulated by Earth’s and the Sun’s magnetic fields, the latter being linked to solar variability. Here we show that accounting for differences in hemispherical production rates, related to geomagnetic field asymmetries, helps resolve so far unexplained differences in Holocene solar activity reconstructions. We find no compelling evidence for long-term variations in solar activity and show that variations in cosmogenic radionuclide production rates on millennial timescales and longer, including the 2,400-year Hallstatt cycle, are explained by variations in the geomagnetic field. Our results also suggest an on-average stronger dipole moment during the Holocene, associated with higher field intensities in the Southern Hemisphere.

Variations in the Solar Modulation Parameter Over the Last 9.5 Thousand Years and the Tilt of the Geomagnetic Dipole

Background: Calculations of the solar modulation parameter (Φ) over the past millennia typically use the relationship between the production rate of cosmogenic isotopes, the earth's dipole moment, and the magnitude of Φ. The cosmogenic isotopes 14C and 10Be are typically used in these studies. When studying solar modulation, the cyclic change in dipole tilt is usually not taken into account, which affects estimates of past solar activity. Methods: Tree rings are a reliable basis for obtaining a radiocarbon time scale (IntCal13). However, determining the concentration of 14C in tree rings is a difficult and controversial task. The time scale derived from the 10Be production rate simulation (GICC05) is less reliable. Nevertheless, there is a way to combine the accuracy of the radiocarbon time scale with the reliability of estimates of the 10Be production rate. This method is the synchronization of the radiocarbon and beryllium-10 series. We have selected the most relevant methods for calculating the solar modulation parameter Φ for the Holocene. When calculating Φ, 10Be data synchronized with 14C data were used. The latest data on the earth's dipole moment were considered. Empirical Mode Decomposition (EMD) was used in the analysis of Φ. Results: It has been shown that the first two decomposition modes are oscillating components with periods of 710 and 208 years, the amplitudes of which increase with time, reaching a maximum of 2500 BP. From contemplation, it follows that the 710-year oscillations are apparently caused by fluctuations in the tilt of the earth's dipole. After excluding the EMD component associated with the 710-year cyclicity, a corrected series was obtained for the solar modulation parameter, free from the influence of changes in the tilt of the magnetic dipole. Conclusion: The rate of formation of cosmogenic radionuclides depends on the intensity of penetration of Galactic Cosmic Rays (GCRs) into the earth's atmosphere. Before reaching earth, GCRs must cross the heliosphere, where they are exposed to solar modulation. Adequate consideration of solar modulation parameters is important for the correct interpretation of the rate of production of cosmogenic isotopes and solar activity.

A Geomagnetic Estimate of Heliospheric Modulation Potential over the Last 175 Years

Galactic cosmic rays (GCRs) interact with the Earth’s atmosphere to produce energetic neutrons and cosmogenic radionuclides, such as 14C. The atmosphere is partially shielded from the interstellar GCR spectrum by both the geomagnetic and solar magnetic fields. Solar shielding is often expressed as the heliospheric modulation potential ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document}, which consolidates information about the strength and structure of the solar magnetic field into a single parameter. For the period 1951 to today, ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} can be estimated from ground-based neutron monitor data. Prior to 1950, 14C in tree rings can be used to estimate ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} and hence the solar magnetic field, back centuries or millennia. Bridging the gap in the ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} record is therefore of vital importance for long-term solar reconstructions. One method is to model ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} using the sunspot number (SN) record. However, the SN record is only an indirect measure of the Sun’s magnetic field, introducing uncertainty, and the record suffers from calibration issues. Here we present a new reconstruction of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} based on geomagnetic data, which spans both the entire duration of the neutron monitor record and stretches back to 1845, providing a significant overlap with the 14C data. We first modify and test the existing model of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} based on a number of heliospheric parameters, namely the open solar flux FS\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$F_{S}$\\end{document}, the heliospheric current sheet tilt angle α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\alpha $\\end{document}, and the global heliospheric magnetic polarity p\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$p$\\end{document}. This modified model is applied to recently updated geomagnetic estimates of FS\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$F_{S}$\\end{document} and cyclic variations of α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\alpha $\\end{document} and p\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$p$\\end{document}. This approach is shown to produce an annual estimate of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} in excellent agreement with that obtained from neutron monitors over 1951 – 2023. It also suggests that ionisation chamber estimates of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} – which have previously been used to extend the instrumental estimate back from 1951 to 1933 – are not well calibrated. Comparison of the new geomagnetic ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} with 14C estimates of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} suggests that the long-term trend is overestimated in the most recent 14C data, possibly due to hemispheric differences in the Suess effect, related to the release of carbon by the burning of fossil fuels. We suggest that the new geomagnetic estimate of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi $\\end{document} will provide an improved basis for future calibration of long-term solar activity reconstructions.

Cosmogenic nuclide-derived downcutting rates of canyons within large limestone plateaus of southern Massif Central (France) reveal a different regional speleogenesis of karst networks

Abstract. We present 35 new burial ages (27 sites) based on 26Al / 10Be ratios of terrestrial cosmogenic radionuclides measured in clasts and sediments deep within 12 caves in the southern Massif Central, France. Our results, together with previously published burial ages, verify that cave morphogenesis has been continuously active in this region for at least the past ∼ 6 Myr. Combining sample burial ages with their associated cave elevation above the modern stream bed gives a mean regional incision rate of 88 ± 5 m Ma−1 for the Grands Causses area. South of the Cevennes Fault Zone bordering the Grands Causses, the incision rate is 43 ± 5 m Ma−1, suggesting that this difference might be accommodated by the fault zone. Sediment burial ages from caves which are not located on river valley flanks or cliff walls are surprisingly too young compared to their expected ages when calculated using this regional average river incision rate. This suggests that the classical epigenic speleogenesis model that presumes a direct correlation between cave level development and regional base level lowering does not apply for the study area. Therefore, we propose that regional speleogenesis is mainly controlled by the removal of ghost rocks by headward erosion from river canyons to central parts of the plateaus, emptying incipient primokarst passages to create cave systems. Our results suggest a continuum process from hypogene primokarst composed of passages filled with ghost rock to one of epigene karst dynamics emptying these passages and creating cave networks. We propose that these processes are the major mechanism in the southern Massif Central that initiates speleogenesis and controls the geometry of the networks. In this region, tiered karst cannot be associated with the pace of incision of the major rivers but must be explained by former ghost rock (or hypogene) processes.

Synchronizing the Western Gotland Basin (Baltic Sea) and Lake Kälksjön (central Sweden) sediment records using common cosmogenic radionuclide production variations

Multi-archive studies of climate events and archive-specific response times require synchronous time scales. Aligning common variations in the cosmogenic radionuclide production rate via curve fitting methods provides a tool for the continuous synchronization of natural environmental archives down to decadal precision. Based on this approach, we synchronize 10Be records from Western Gotland Basin (WGB, Baltic Sea) and Lake Kälksjön (KKJ, central Sweden) sediments to the 14C production time series from the IntCal20 calibration curve during the Mid-Holocene period ~6400 to 5200 a BP. Before the synchronization, we assess and reduce non-production variability in the 10Be records by using 10Be/9Be ratios and removing common variability with the TOC record from KKJ sediments based on regression analysis. The synchronizations to the IntCal20 14C production time scale suggest decadal to multi-decadal refinements of the WGB and KKJ chronologies. These refinements reduce the previously centennial chronological uncertainties of both archives to about ± 20 (WGB) and ±40 (KKJ) years. Combining proxy time series from the synchronized archives enables us to interpret a period of ventilation in the deep central Baltic Sea basins from ~6250 to 6000 a BP as possibly caused by inter-annual cooling reducing vertical water temperature gradients allowing deep water formation during exceptionally cold winters.

Mammal-bearing Pleistocene deposits, Vranić, southwestern Pannonian Basin System

Vranić site is a sand quarry that is located on the southern foothills of Papuk Mountain. From bottom to top of the succession, three sedimentary units have been recognized as: Unit 1 containing massive sand with scattered gravel-sized clasts, marl cobbles and boulders, and abundant marine mammal and fish fossils; Unit 2 consisting of sand intercalated with silt, clay and gravel, which may be horizontally bedded; and Unit 3, which is an erosionally-based lenticular matrix to clast-supported structureless gravel. The basal part of Unit 1 shows numerous reworked skeletal remains of Miocene marine mammals. Cosmogenic radionuclides constrain the age of burial of Unit 1 to 895 ± 211 ka, while the source area of Unit 1 had a quick denudation.The massive sands deposited on the slopes of Papuk Mountain were vulnerable to erosion due to the absence of clay. Heavy rainfall or water from snow melting created flash floods that infiltrated the sands, thereby causing slope destabilization and deformation. This process led to slumps that were transformed into a sandy debris flow. This sediment was probably deposited during the interglacial marine isotope age (MIS) 21 period when the scarce vegetation and a warmer climate favored the melting of permafrost ice and consequently triggered slope movements during MIS 22. The reworked skeletal material sampled at the Vranić site comprises fossilized fish along with mammal bones and teeth. Thus, the Vranić site represents an important site for cetacean fossil remains and serves as an important data source for marine life in the Central Paratethys.

Study of cosmogenic activation in 76Ge enriched germanium detectors during fabrication and transportation above ground

Rare event search experiments using germanium detectors are operated in underground laboratories to minimize the background induced by cosmic rays. However, the cosmogenic activation in germanium crystals on the ground during fabrication and transportation generates long half-life radionuclides and contributes a considerable background. We simulated the production rates of cosmogenic radionuclides in germanium and calculated the specific activities of cosmogenic radionuclides according to the scheduled fabrication and transportation processes of 76Ge enriched germanium detectors. The impact of cosmogenic background in germanium crystals for the next generation CDEX experiment was assessed with the scheduled exposure history above ground.

Case studies of three geological archives for rare radionuclide measurements using accelerator mass spectrometry

Long-lived radionuclides in our environment provide important information on natural and anthropogenic processes. Their presence and concentration reflect the balance of production and decay. Geological archives store such information and the nuclides can be chemically extracted from the bulk sample. Accelerator mass spectrometry (AMS) represents a sensitive method to quantify those nuclides at natural levels. Three different terrestrial archives are discussed here as examples for radionuclide extraction using various chemical separation methods for subsequent AMS measurements. We focus on sample preparation for the cosmogenic radionuclides 10Be and 26Al, various anthropogenic actinide isotopes such as U, Pu, and Am as well as the astrophysically interesting nuclides 41Ca, 53Mn, and 60Fe. The processed materials cover samples with masses between a few mg and up to a few hundred kg and protocols are presented for the quantitative extraction of some 10,000 atoms of cosmogenic or interstellar origin per sample and even as low as a few hundred actinide atoms.

Simulations of 7Be and 10Be with the GEOS-Chem global model v14.0.2 using state-of-the-art production rates

Abstract. The cosmogenic radionuclides 7Be and 10Be are useful tracers for atmospheric transport studies. Combining 7Be and 10Be measurements with an atmospheric transport model can not only improve our understanding of the radionuclide transport and deposition processes but also provide an evaluation of the transport process in the model. To simulate these aerosol tracers, it is critical to evaluate the influence of radionuclide production uncertainties on simulations. Here we use the GEOS-Chem chemical transport model driven by the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) reanalysis to simulate 7Be and 10Be with the state-of-the-art production rate from the CRAC:Be (Cosmic Ray Atmospheric Cascade: Beryllium) model considering realistic spatial geomagnetic cutoff rigidities (denoted as P16spa). We also perform two sensitivity simulations: one with the default production rate in GEOS-Chem based on an empirical approach (denoted as LP67) and the other with the production rate from the CRAC:Be but considering only geomagnetic cutoff rigidities for a geocentric axial dipole (denoted as P16). The model results are comprehensively evaluated with a large number of measurements including surface air concentrations and deposition fluxes. The simulation with the P16spa production can reproduce the absolute values and temporal variability of 7Be and 10Be surface concentrations and deposition fluxes on annual and sub-annual scales, as well as the vertical profiles of air concentrations. The simulation with the LP67 production tends to overestimate the absolute values of 7Be and 10Be concentrations. The P16 simulation suggests less than 10 % differences compared to P16spa but a significant positive bias (∼18 %) in the 7Be deposition fluxes over East Asia. We find that the deposition fluxes are more sensitive to the production in the troposphere and downward transport from the stratosphere. Independent of the production models, surface air concentrations and deposition fluxes from all simulations show similar seasonal variations, suggesting a dominant meteorological influence. The model can also reasonably simulate the stratosphere–troposphere exchange process of 7Be and 10Be by producing stratospheric contribution and 10Be/7Be ratio values that agree with measurements. Finally, we illustrate the importance of including the time-varying solar modulations in the production calculation, which significantly improve the agreement between model results and measurements, especially at mid-latitudes and high latitudes. Reduced uncertainties in the production rates, as demonstrated in this study, improve the utility of 7Be and 10Be as aerosol tracers for evaluating and testing transport and scavenging processes in global models. For future GEOS-Chem simulations of 7Be and 10Be, we recommend using the P16spa (versus default LP67) production rate.

Cosmogenic Nuclide Tracking of Sediment Recycling From a Frontal Siwalik Range in the Northwestern Himalaya

AbstractThe Himalayan orogen exports millions of tons of sediment annually to the Indo‐Gangetic foreland basin, derived from both hinterland and foreland fold‐thrust belts (FTB). Although erosion rates in the hinterland are well‐constrained, erosion rates in the foreland FTB and, by extension, the sediment flux have remained poorly constrained. Here, we quantified erosion rates and sediment flux from the Mohand Range in the northwestern Himalaya by modeling and measuring the cosmogenic radionuclide (CRN) 10Be and 26Al concentrations in modern fluvial sediments. Our model uses local geological and geophysical constraints and accounts for CRN inheritance and sediment recycling, which enables us to determine the relative contributions of the hinterland and foreland FTB sources to the CRN budget of the proximal foreland deposits. Our model predictions closely match measured concentrations for a crustal shortening rate across the Mohand Range of 8.0 ± 0.5 mm yr−1 (i.e., approximately 50% of the total shortening across the Himalaya at this longitude) since Ma. This shortening implies a spatial gradient in erosion rates ranging from 0.42 ± 0.03 to 4.92 ± 0.34 mm yr−1, controlled by the geometry of the underlying structure. This erosion pattern corresponds to an annual sediment recycling of ∼2.0 megatons from the Mohand Range to the downstream Yamuna foreland. Converted to sediment fluxes per unit width along the Himalaya, the foreland FTB accounts for ∼21% ± 5% of the total flux entering the foreland. Because these sediments have lower 10Be concentrations than hinterland‐derived sediment, they would lead to ∼14% overestimation of 10Be‐derived erosion rates, based on Yamuna sediments in the proximal foreland.

36Cl, a new tool to assess soil carbon dynamics

Soil organic carbon is one of the largest surface pools of carbon that humans can manage in order to partially mitigate annual anthropogenic CO2 emissions. A significant element to assess soil sequestration potential is the carbon age, which is evaluated by modelling or experimentally using carbon isotopes. Results, however, are not consistent. The 14C derived approach seems to overestimate by a factor of 6–10 the average carbon age in soils estimated by modeling and 13C approaches and thus the sequestration potential. A fully independent method is needed. The cosmogenic chlorine nuclide, 36Cl, is a potential alternative. 36Cl is a naturally occurring cosmogenic radionuclide with a production that increased by three orders of magnitude during nuclear bomb tests. Part of this production is retained by soil organic matter in organochloride form and hence acts as a tracer of the fate of soil organic carbon. We here quantify the fraction and the duration of 36Cl retained in the soil and we show that retention time increases with depth from 20 to 322 years, in agreement with both modelling and 13C-derived estimates. This work demonstrates that 36Cl retention duration can be a proxy for the age of soil organic carbon.

26Al and 10Be concentrations from alluvial drill cores across the Indo-Gangetic plain reveal multimillion-year sediment-transport lag times

Sediment burial can lead to significant lag times during sediment transport and is an important component of the sediment cycle. The Indo-Gangetic Plain (IGP) is a large sediment conveyer belt and understanding sediment lag times on various time scales is important for short- and long-term sediment transport and denudation-rate estimations. In this study, we have quantified the transport lag time in terms of burial duration with 14 different alluvial drill cores from 4 major Himalayan river basins, including the Ghaggar, Paleo-Yamuna, Ganges, and Kosi rivers. We rely on 56 samples using the paired cosmogenic radionuclides, 26Al and 10Be. The coring locations of paleo-Yamuna and Kosi, along with one of the Ghaggar cores are within ∼100 km from the mountain front, while the Ganges cores are ∼350 km away. The coring locations cover an E-W distance of ∼1200 km on the IGP and the depositional ages range from ∼120 kyr BP to the present day based on previous age estimations. Out of 56 samples, 39 samples show a burial signal ranging from 0.61±0.33 Myr to 3.75±1.9 Myr, with outlier values as large as ∼6 Myr.We show that the Kosi samples due to their proximity to the mountain front show no burial or a lower mean burial duration of ∼1.01 Myr, compared to the samples from the coring locations well within the IGP, such as the Ghaggar (∼2.48 Myr) and the Ganges (∼2.36 Myr). However, despite its proximity to the mountain front, the paleo-Yamuna samples have a mean burial duration of ∼2.20 Myr. Along with the floodplains and megafan surface of the IGP, we identify two potential major sediment-lag sources in the Himalayan catchment: The high and low elevation broad valleys and the Siwaliks. Increased storage spaces in the valleys north of the main Himalayan crest, including drainage-captured catchments are likely to provide long burial durations such as in the case of the Kosi and Ghaggar basins. The upper and middle Siwaliks with burial duration in the range of up to ∼6 Myr constitute a second source. We attribute the no-burial signals of 7 samples of the Kosi basin to rapid sediment generation and transportation, driven by dynamic hillslope erosion and fluvial dynamics.Importantly, our data suggest that even samples close to the mountain front exhibit million years of burial duration. These findings suggest that denudation rates derived from modern river sands or in paleo settings will need to account for burial duration and associated processes.

Effects of Global Geomagnetic Field Variations Over the Past 100,000 Years on Cosmogenic Radionuclide Production Rates in the Earth's Atmosphere

AbstractThe production rates of cosmogenic radionuclides, such as 10Be, 14C, and 36Cl, in the Earth's atmosphere vary with the geomagnetic field and solar activity. For the first time, the production rates of several cosmogenic nuclides are estimated for the past 100 ka based on global, time‐dependent geomagnetic field models and a moderate solar‐activity level. In particular, the production rates were high with no notable latitudinal dependence during the Laschamps geomagnetic excursion (41 ka BP). The mean global production of 10Be over the Laschamps excursion was more than two times greater than the present‐day one, whereas the increase was 1.9 times for the Norwegian‐Greenland Sea excursion (∼65 ka), and only 1.3 times for the Mono Lake/Auckland excursion (∼34 ka). All analyzed geomagnetic field models covering the past 100 ka, including the modern and Holocene epochs, lead to hemispheric asymmetry in the production rates, persistent overall time ranges, and reflected in the time‐averaged nuclide production rates. Production rates predicted by the geomagnetic field models are in good agreement with actual measurements from ice cores and sediment records. These global, long‐term production rates are important for a wide range of studies that employ cosmogenic nuclides as a proxy/tracer of different Earth system processes.

Synchronizing ice-core and U ∕ Th timescales in the Last Glacial Maximum using Hulu Cave 14C and new 10Be measurements from Greenland and Antarctica

Abstract. Between 15 and 27 kyr b2k (thousands of years before 2000 CE) during the last glacial, Greenland experienced a prolonged cold stadial phase, interrupted by two short-lived warm interstadials. Greenland ice-core calcium data show two periods, preceding the interstadials, of anomalously high atmospheric dust loading, the origin of which is not well understood. At approximately the same time as the Greenland dust peaks, the Chinese Hulu Cave speleothems exhibit a climatic signal suggested to be a response to Heinrich Event 2, a period of enhanced ice-rafted debris deposition in the North Atlantic. In the climatic signal of Antarctic ice cores, moreover, a relative warming occurs between 23 and 24.5 kyr b2k that is generally interpreted as a counterpart to a cool climate phase in the Northern Hemisphere. Proposed centennial-scale offsets between the polar ice-core timescales and the speleothem timescale hamper the precise reconstruction of the global sequence of these climatic events. Here, we examine two new 10Be datasets from Greenland and Antarctic ice cores to test the agreement between different timescales, by taking advantage of the globally synchronous cosmogenic radionuclide production rates. Evidence of an event similar to the Maunder Solar Minimum is found in the new 10Be datasets, supported by lower-resolution radionuclide data from Greenland and 14C in the Hulu Cave speleothem, representing a good synchronization candidate at around 22 kyr b2k. By matching the respective 10Be data, we determine the offset between the Greenland ice-core chronology, GICC05, and the Antarctic chronology for the West Antarctic Ice Sheet Divide ice core (WDC), WD2014, to be 125 ± 40 years. Furthermore, via radionuclide wiggle-matching, we determine the offset between the Hulu speleothem and ice-core timescales to be 375 years for GICC05 (75–625 years at 68 % confidence) and 225 years for WD2014 (−25–425 years at 68 % confidence). The rather wide uncertainties are intrinsic to the wiggle-matching algorithm and the limitations set by data resolution. The undercounting of annual layers in GICC05 inferred from the offset is hypothesized to have been caused by a combination of underdetected annual layers, especially during periods with low winter precipitation, and misinterpreted unusual patterns in the annual signal during the extremely cold period often referred to as Heinrich Stadial 1.

Precipitation‐Driven Gamma Radiation Enhancement Over the Atlantic Ocean

AbstractGamma radiation over the Atlantic Ocean was measured continuously from January to May 2020 by a NaI(Tl) detector installed on board the Portuguese navy's ship NRP Sagres. Enhancements in the gamma radiation values are identified automatically by an algorithm for detection of anomalies in mean and variance as well as by visual inspection. The anomalies are typically +50% above the background level and relatively rare events (∼<10% of the days). All the detected anomalies are associated with simultaneous precipitation events, consistent with the wet deposition of scavenged radionuclides. The enhancements are detected in the open ocean even at large distances (+500 km) from the nearest coastline. Back trajectories reveal that half of these events are associated with air masses experiencing continental land influences, but the other half do not display evidence of recent land contact. The enhancements in gamma radiation very far from land and with no evidence of continental fetch from back trajectories are difficult to explain as resulting only from radionuclides with a terrestrial source such as radon and its progeny. Further investigation and additional measurements are needed to improve understanding on the sources of ambient radioactivity in the open ocean and assess whether gamma radiation in the marine environment is influenced not only by radionuclides of terrestrial origin, but also cosmogenic radionuclides, like Beryllium‐7, formed in the upper atmosphere but with the ability to be transported downward and serve as a tracer of the aerosols to which it attaches.

Cosmogenic Radionuclide-Beryllium 7 (7Be) for Monsoon Rainfall Forecasting in Malaysia: A Systematic Literature Review

Within the region of India, monsoon prediction was reported to use the concentration of 7Be cosmogenic radionuclide that contains crucial information such as atmospheric behaviour change. This review aims to identify and analyse the existing research on the usage of 7Be Cosmogenic radionuclide in assisting weather forecasts, especially rain onset and withdrawal in the monsoon season and the correlations of rainfall and monsoonal variation in Malaysia. A comprehensive search of major scientific databases was conducted, and relevant studies were selected based on inclusion and exclusion criteria between January 2017 and August 2022. The review found that the method of utilising 7Be as a parameter in monsoon forecasting for India is newly developed and only be assessed using the trans-equatorial method. The review also highlighted the expectation that the 7Be concentration analysis throughout the years could provide new insight into understanding monsoon behaviour and can be further utilised in forecasting the rainfall intensity during a wet season in Malaysia.