Cosmogenic Radionuclide Concentrations Research Articles

Geochemistry, petrology and ages of the lunar meteorites Kalahari 008 and 009: New constraints on early lunar evolution

Kalahari 008 and 009 are two lunar meteorites that were found close to each other in Botswana. Kalahari 008 is a typical lunar anorthositic breccia; Kalahari 009 a monomict breccia with basaltic composition and mineralogy. Based on minor and trace elements Kalahari 009 is classified as VLT (very-low-Ti) mare basalt with extremely low contents of incompatible elements, including the REE. The Lu–Hf data define an age of 4286 ± 95 Ma indicating that Kalahari 009 is one of the oldest known basalt samples from the Moon. It provides evidence for lunar basalt volcanism prior to 4.1 Ga (pre-Nectarian) and may represent the first sample from a cryptomare. The very radiogenic initial 176Hf/ 177Hf (εHf = +12.9 ± 4.6), the low REE, Th and Ti concentrations indicate that Kalahari 009 formed from re-melting of mantle material that had undergone strong incompatible trace element depletion early in lunar history. This unusually depleted composition points toward a hitherto unsampled basalt source region for the lunar interior that may represent a new depleted endmember source for low-Ti mare basalt volcanism. Apparently, the Moon became chemically very heterogeneous at an early stage in its history and different cumulate sources are responsible for the diverse mare basalt types. Evidence that Kalahari 008 and 009 may be paired includes the similar fayalite content of their olivine, the identical initial Hf isotope composition, the exceptionally low exposure ages of both rocks and the fact that they were found close to each other. Since cryptomaria are covered by highland ejecta, it is possible that these rocks are from the boundary area, where basalt deposits are covered by highland ejecta. The concentrations of cosmogenic radionuclides and trapped noble gases are unusually low in both rocks, although Kalahari 008 contains slightly higher concentrations. A likely reason for this difference is that Kalahari 008 is a polymict breccia containing a briefly exposed regolith, while Kalahari 009 is a monomict brecciated rock that may never have been at the surface of the Moon. Altogether, the compositions of Kalahari 008 and 009 permit new insight into early lunar evolution, as both meteorites sample lunar reservoirs hitherto unsampled by spacecraft missions. The very low Th and REE content of Kalahari 009 as well as the depletion in Sm and the lack of a KREEP-like signature in Kalahari 008 point to a possible source far from the influence of the Procellarum-KREEP Terrane, possibly the lunar farside.

Effects of physical erosion on chemical denudation rates: A numerical modeling study of soil-mantled hillslopes

Many biogeochemical and Earth surface processes depend critically on chemical weathering. The immediate products of chemical weathering are present as solutes and secondary minerals in groundwater, soils, and streams, and form the nutritional foundation for terrestrial biogeochemistry. Chemical weathering also contributes to physical erosion by weakening bedrock and producing easily erodible regolith, and as the primary long-term sink for atmospheric CO 2 it modulates Earth's long-term climate via the greenhouse effect. Long-term chemical denudation rates on soil-mantled hillslopes can be estimated from cosmogenic radionuclide (CRN) concentrations in soil-borne quartz and the enrichment of a chemically inert tracer in soil relative to its parent bedrock, a technique that inherently assumes steady physical erosion over the timescale of CRN accumulation. We present a numerical model that computes changes in soil mineralogy and CRN concentrations under time-varying physical erosion rates, and we use this model to assess the accuracy of the CRN-based technique for estimating chemical denudation rates in non-steady conditions. Our modeling results suggest that CRN-based estimates of chemical denudation rates closely resemble actual chemical denudation rates averaged over the timescale of CRN accumulation, even during large-amplitude and long-period oscillations in physical erosion rates. For example, this model predicts that when physical erosion rates fluctuate sinusoidally by 50% of their mean over any period in time, CRN-based estimates of chemical denudation rates should differ from actual chemical denudation rates by less than 15%. Our model also implies that chemical denudation rates should approach zero both when physical erosion rates approach zero (because soluble minerals become depleted in the soil) and when physical erosion rates approach the maximum soil production rate (because soil thickness approaches zero). Modeled chemical denudation rates thus reach a maximum at intermediate physical erosion rates. If this relationship holds in nature, it implies that in rapidly eroding regions, further increases in physical erosion rates (e.g., due to increases in tectonic uplift rates) may not necessarily lead to faster chemical denudation on soil-mantled hillslopes.

Unraveling complex exposure-burial histories of bedrock surfaces under ice sheets by integrating cosmogenic nuclide concentrations with climate proxy records

The production, accumulation, and decay of cosmogenic radionuclides in rock surfaces subjected to episodes of exposure and burial by ice results in nuclide concentrations in present day rock surfaces that can be used to address a variety of questions in glacial geomorphology and Quaternary geology. Of particular importance is the fact that these nuclide concentrations reflect both the timing of initial exposure of the rock surface and the chronology of subsequent exposure, burial, and erosion episodes. For landscapes where geomorphic evidence indicates that little/no erosion occurred, constraining the timing of initial exposure and the number of phases of exposure and burial that a rock surface has been subjected to is possible using multiple cosmogenic radionuclide concentrations combined with proxies for the timing and duration of periods of ice cover, such as ice core or marine isotope records. However, interpretations based on this approach require determination of an appropriate cutoff value to separate the proxy record into ice-free and ice-covered conditions and assessment of the sensitivity of the results to different cutoff values. We have developed a numerical model to evaluate variations in total exposure and burial durations as a function of different proxy records and cutoff values. This program is available at http://www.missouri.edu/~liyk/ClimateProxyCurve.zip. Initial results for sites in West Antarctica and northern Sweden show that the method provides a quick and robust way to derive best-fit cutoff values and chronologies of burial and exposure, and small changes in cutoff values can result in significant shifts in results. The method described here provides new insight into the interpretation and reliability of multiple nuclide samples. This approach also has the potential to provide improved constraints for ice sheet dynamics and landscape evolution, and a means to assess the sensitivity of calculated initial exposure dates to assumptions about ice sheet history.

A threshold in soil formation at Earth’s arid–hyperarid transition

The soils of the Atacama Desert in northern Chile have long been known to contain large quantities of unusual salts, yet the processes that form these soils are not yet fully understood. We examined the morphology and geochemistry of soils on post-Miocene fans and stream terraces along a south-to-north (27° to 24° S) rainfall transect that spans the arid to hyperarid transition (21 to ∼2 mm rain y −1). Landform ages are ⩾ 2 My based on cosmogenic radionuclide concentrations in surface boulders, and Ar isotopes in interbedded volcanic ash deposits near the driest site indicate a maximum age of 2.1 My. A chemical mass balance analysis that explicitly accounts for atmospheric additions was used to quantify net changes in mass and volume as a function of rainfall. In the arid (21 mm rain y −1) soil, total mass loss to weathering of silicate alluvium and dust (−1030 kg m −2) is offset by net addition of salts (+170 kg m −2). The most hyperarid soil has accumulated 830 kg m −2 of atmospheric salts (including 260 kg sulfate m −2 and 90 kg chloride m −2), resulting in unusually high volumetric expansion (120%) for a soil of this age. The composition of both airborne particles and atmospheric deposition in passive traps indicates that the geochemistry of the driest soil reflects accumulated atmospheric influxes coupled with limited in-soil chemical transformation and loss. Long-term rates of atmospheric solute addition were derived from the ion inventories in the driest soil, divided by the landform age, and compared to measured contemporary rates. With decreasing rainfall, the soil salt inventories increase, and the retained salts are both more soluble and present at shallower depths. All soils generally exhibit vertical variation in their chemistry, suggesting slow and stochastic downward water movement, and greater climate variability over the past 2 My than is reflected in recent (∼100 y) rainfall averages. The geochemistry of these soils shows that the transition from arid to hyperarid rainfall levels marks a fundamental geochemical threshold: in wetter soils, the rate and character of chemical weathering results in net mass loss and associated volumetric collapse after 10 5 to 10 6 years, while continuous accumulation of atmospheric solutes in hyperarid soils over similar timescales results in dramatic volumetric expansion. The specific geochemistry of hyperarid soils is a function of atmospheric sources, and is expected to vary accordingly at other hyperarid sites. This work identifies key processes in hyperarid soil formation that are likely to be independent of location, and suggests that analogous processes may occur on Mars.

Interpreting erosion rates from cosmogenic radionuclide concentrations measured in rapidly eroding terrain

AbstractA combination of numerical analysis and 10Be concentrations measured in sediment samples from the high‐relief Torrente catchment, southern Spain, allows us to investigate the sampling requirements for determining erosion rates using cosmogenic nuclides in high‐relief, landslide‐dominated terrain. We use simple modelling to quantify the effect of particle spalling and/or landsliding on erosion rates determined using a cosmogenic in‐situ produced isotope. Analytical results show that the cosmogenic nuclide concentration of a surface experiencing regular detachment of a grain or block may be considered to be in steady state, and ‘in‐situ’ erosion rates estimated, when an appropriate number of spatially independent samples are amalgamated. We present equations that enable calculation of the number of bedrock samples that must be amalgamated for the estimation of mean erosion rates on an outcrop experiencing regular detachment of a grain or chip of thickness L every T years. Our findings confirm that mean catchment erosion rates may be reliably estimated from 10Be concentrations in fluvial sediment in high‐relief rapidly eroding terrain. These catchment‐wide integrated erosion rates can be calculated where erosion is primarily accomplished through shallow (<3 m) spalling processes; where deep‐seated (>3 m) landslides are the dominant mode of erosion only minimum erosion rates can be determined. Lastly, we present erosion rate measurements from the Torrente catchment that reveal variation of two orders of magnitude (0·03–1·6 m ka−1) quantifying the high degree of spatial variation in erosion rates expected within rapidly uplifting catchments. Copyright © 2006 John Wiley & Sons, Ltd.

Cosmogenic nuclide evidence for minimal erosion across two subglacial sliding boundaries of the late glacial Fennoscandian ice sheet

The existence of sliding and frozen bed areas under ice sheets is significant in understanding basal thermal regimes, patterns of erosion and landform development, and in constraining boundary conditions for the reconstructions of ice sheets. Recognition of subglacial boundaries between sliding and frozen-bed areas for former ice sheets is typically based on distinct morphological contrasts between areas with glacial landform assemblages and relict areas showing little alteration of pre-existing features. Some of these boundaries, especially on continental shield areas, however, are clearly visible from air photos but have minimal topographic expression. Understanding the chronology and erosional development of such boundaries is important to provide insight into the pattern and persistence of basal conditions under ice sheets. Geomorphic evidence and cosmogenic radionuclide concentrations of bedrock outcrops on either side of two sliding boundaries on Ultevis and Arvestuottar, low-relief upland plateaus in northern Sweden, are consistent with negligible erosion in relict landscape (frozen bed) areas due to the last glaciation, but also indicate insignificant erosion in the sliding areas. Such a pattern and magnitude of landscape modification indicates that sliding was short lived in these areas, likely as a transient phase during deglaciation. These sites demonstrate that short periods of sliding are in some cases sufficient to produce landscapes that are recognized as ‘glacial’ from air photos. Thus, regions of sliding identified on shield areas must be viewed as the cumulative total area that has experienced sliding at any time during a glaciation. The actual extent of sliding areas during any single ice sheet phase is presumably considerably less than this cumulative total, which has important implications for establishing appropriate basal boundary conditions for ice sheet reconstructions.

Use of cosmogenic radionuclides as a sediment tracer in the Santa Cruz littoral cell, California, United States

We demonstrate that cosmogenic radionuclides can be employed to trace sediment when the sources have sufficiently distinct concentrations, and develop a theoretical mixing model for a rocky coastline littoral system. We combine the resulting mixing model with existing cosmogenic radionuclide methods that quantify river inputs and terrace ages to determine the major components of the long-term littoral sediment budget of the Santa Cruz, California, coastline. Sediment derived from coastal basins eroding at 0.2 mm/yr has a much lower concentration than sediment derived from 60‐84 ka terrace sands atop backwearing seacliffs. The complex pattern of cosmogenic radionuclide concentrations in littoral sands along .100 km of coastline can be explained by mixing sediment derived from seacliffs backwearing at 10 cm/yr with that delivered by rivers having widely different amounts of sediment discharge.

Pleistocene relief production in Laramide mountain ranges, western United States

Gently sloped summits and ridges (collectively referred to as summit flats) are abundant in many Laramide ranges in the western United States. The erosion rate of summit flats is ~10 m/m.y., on the basis of the concentrations of cosmogenic radionuclides. Because erosion rates in valleys between summit flats are an order of magnitude faster, relief within these ranges is currently increasing by about 100 m/m.y. If summit-flat erosion is slower than rock uplift driven by the isostatic response to valley erosion, then this relief production could result in increased summit elevations. The mean depth of material eroded from a smooth surface fit to existing summit flats varies from 280 to 340 m in four Laramide ranges, based on geographic information system (GIS) analyses of digital elevation models. This erosion would result in a maximum of 250‐300 m of rock uplift, assuming Airy isostasy. However, because the Laramide ranges examined here are narrow relative to the flexural wavelength of the lithosphere, erosionally driven rock uplift is limited to ~ 50‐100 m. Over the past several million years, summit erosion would approximately offset this rock uplift. Therefore, we conclude that summit elevations have remained essentially constant even though several hundred meters of relief has been produced. On the basis of valley and summit erosion rates and the average depth of erosion, we estimate that relief production in Laramide ranges began at ca. 3 Ma. We hypothesize that this relief production was climatically driven and was associated with the onset or enhancement of alpine glaciation in these ranges.

Age of Allan Hills 82102, a meteorite found inside the ice

THE recovery of thousands of meteorites in the Antarctic since 1969 has not only greatly increased knowledge of the meteorites themselves but has also provided a new tool for glaciology. Most Antarctic meteorites are found on blue ice areas where old ice is continuously ablated. A measurement of the age of this ice helps us understand meteorite accumulation mechanisms and the dynamics of ice movement. The terrestrial age of a meteorite, the time period since the date of meteorite fall, can be determined from the reduction in concentration of cosmogenic radionuclides during the time the meteorite has been shielded by the Earth's atmosphere. Here we report the terrestrial age of a meteorite that was recovered from below the surface of the ice and argue that this represents a measurement of the age of the ice itself. We measured the cosmogenic radionuclides 10Be (half-life = 1.5 Myr), 14C (5,730 yr), 26A1 (0.71 Myr), 36C1 (0.30 Myr) and 53Mn (3.7 Myr) in the meteorite and 10Be and 36C1 in the ice. We obtained a terrestrial age of 11,000 years for the meteorite which suggests that the snow accumulation area where it fell was only a few tens of kilometres away.

Cosmic-ray-produced 53Mn in thirty-one meteorites

Cosmic-ray-produced 53Mn (t ½ = 3.7 × 10 6years) has been determined by neutron activation analysis in twenty-two chondrites including three Antarctic meteorites: Yamato-7301 (j), -7305 (k) and -7304 (m). 53Mn was also measured in four mesosiderites, three iron meteorites, Bencubbin (unique) and Udei Station (iron with silicate inclusions). In addition, preliminary results for 10Be (1.6 × 10 6 years) were obtained in the Yamato meteorites using a low-background needle GM counter. Based on published values of rare gas ages, corrections were made for undersaturation; the average specific saturation activities of 53Mn were found in the range 450 ± 90dpm 53Mn/kg Fe in most of the chondrites and 490 ± 75dpm 53Mn/kg Fe in the mesosiderites. Two meteorites had extremely low contents of 53Mn: 102 ± 6 in Yamato-7301 and 48 ± 3dpm 53Mn/kg Fe in Bondoc. The Bondoc mesosiderite was already known to have a low concentration of cosmogenic radionuclides due to its large pre-atmospheric size. Several possible mechanisms are discussed to explain the low 53Mn activity in Yamato-7301: (1) long terrestrial age of about 7 m.y.; (2) low production rate of 53Mn due to heavy pre-atmospheric shielding (>70cm); (3) multi-stage irradiation history resulting in an undersaturation of 53Mn; and (4) a mechanism in which two or three of the above factors are combined. The ratio of 53Mn production rate in Ni to that in Fe has been estimated to be 1/3, based on the measurements of 53Mn in the metallic and silicate phases of St. Séverin meteorite, as well as on published results of some high-energy bombardment experiments.