Neon diffusion kinetics and implications for cosmogenic neon paleothermometry in feldspars

Observations of cosmogenic neon concentrations in feldspars can potentially be used to constrain the surface exposure duration or surface temperature history of geologic samples. The applicability of cosmogenic neon to either application depends on the temperature-dependent diffusivity of neon isotopes. In this work, we investigate the kinetics of neon diffusion in feldspars of different compositions and geologic origins through stepwise degassing experiments on single, proton-irradiated crystals. To understand the potential causes of complex diffusion behavior that is sometimes manifest as nonlinearity in Arrhenius plots, we compare our results to argon stepwise degassing experiments previously conducted on the same feldspars. Many of the feldspars we studied exhibit linear Arrhenius behavior for neon whereas argon degassing from the same feldspars did not. This suggests that nonlinear behavior in argon experiments is an artifact of structural changes during laboratory heating. However, other feldspars that we examined exhibit nonlinear Arrhenius behavior for neon diffusion at temperatures far below any known structural changes, which suggests that some preexisting material property is responsible for the complex behavior. In general, neon diffusion kinetics vary widely across the different feldspars studied, with estimated activation energies (Ea) ranging from 83.3 to 110.7 kJ/mol and apparent pre-exponential factors (D0) spanning three orders of magnitude from 2.4 × 10−3 to 8.9 × 10−1 cm2 s−1. As a consequence of this variability, the ability to reconstruct temperatures or exposure durations from cosmogenic neon abundances will depend on both the specific feldspar and the surface temperature conditions at the geologic site of interest

Early-to-mid Miocene erosion rates inferred from pre-Dead Sea rift Hazeva River fluvial chert pebbles using cosmogenic 21Ne

In this work, we utilize a novel application of cosmogenic 21Ne measurements in chert to compare exposure times measured in eroding surfaces in the central Jordanian Plateau with exposure times from chert pebbles transported by the Miocene Hazeva River. The Miocene Hazeva River was a large fluvial system (estimated catchment size > 100 000 km2) that drained the Arabian Plateau and Sinai Peninsula into the Mediterranean Sea during the early-to-mid Miocene. It was established after the rifting of the Red Sea uplifted the Arabian Plateau during the Oligocene. Following late-Miocene-to-early-Pliocene subsidence along the Dead Sea rift, the Hazeva drainage system was abandoned and dissected, resulting in new drainage divides on either side of the rift. We find modern erosion rates derived from cosmogenic 21Ne, 26Al, and 10Be in exposed in situ chert nodules to be extremely slow (between 2-4 mm kyr-1). Comparison between modern and paleo-erosion rates, measured in chert pebbles, is not straightforward, as cosmogenic 21Ne was acquired partly during bedrock erosion and partly during transport of these pebbles in the Hazeva River. However, 21Ne exposure times calculated in Miocene cherts are generally shorter (ranging between 0C59-0 and 242±113 kyr) compared to exposure times calculated in the currently eroding chert nodules presented here (269±49 and 378±76 kyr) and other chert surfaces currently eroding in hyperarid environments. Miocene exposure times are shorter even when considering that they account for bedrock erosion in addition to maintained transport along this large river. Shorter exposure times in Miocene cherts correspond to faster paleo-erosion rates, which we attribute to a combination of continuous surface uplift and significantly wetter climatic conditions during the early-to-mid Miocene

Pleistocene Relative Sea Levels in the Chesapeake Bay Region and Their Implications for the Next Century

Today, relative sea-level rise (3.4 mm/yr) is faster in the Chesapeake Bay region than any other location on the Atlantic coast of North America, and twice the global average eustatic rate (1.7 mm/yr). Dated interglacial deposits suggest that relative sea levels in the Chesapeake Bay region deviate from global trends over a range of timescales. Glacio-isostatic adjustment of the land surface from loading and unloading of continental ice is likely responsible for these deviations, but our understanding of the scale and timeframe over which isostatic response operates in this region remains incomplete because dated sea-level proxies are mostly limited to the Holocene and to deposits 80 ka or older. To better understand glacio-isostatic control over past and present relative sea level, we applied a suite of dating methods to the stratigraphy of the Blackwater National Wildlife Refuge, one of the most rapidly subsiding and lowest-elevation surfaces bordering Chesapeake Bay. Data indicate that the region was submerged at least for portions of marine isotope stage (MIS) 3 (ca. 60–30 ka), although multiple proxies suggest that global sea level was 40–80 m lower than present. Today MIS 3 deposits are above sea level because they were raised by the Last Glacial Maximum forebulge, but decay of that same forebulge is causing ongoing subsidence. These results suggest that glacio-isostasy controlled relative sea level in the mid-Atlantic region for tens of thousands of years following retreat of the Laurentide Ice Sheet and continues to influence relative sea level in the region. Thus, isostatically driven subsidence of the Chesapeake Bay region will continue for millennia, exacerbating the effects of global sea-level rise and impacting the region’s large population centers and valuable coastal natural resources

Cosmogenic Nuclide Systematics and the CRONUScalc Program

As cosmogenic nuclide applications continue to expand, the need for a common basis for calculation becomes increasingly important. In order to accurately compare between results from different nuclides, a single method of calculation is necessary. Calculators exist in numerous forms with none matching the needs of the CRONUS-Earth project to provide a simple and consistent method to interpret data from most commonly used cosmogenic nuclides. A new program written for this purpose, CRONUScalc, is presented here. This unified code presents a method applicable to 10Be, 26Al, 36Cl, 3He, and 14C, with 21Ne in testing. The base code predicts the concentration of a sample at a particular depth for a particular time in the past, which can be used for many applications. The multi-purpose code already includes functions for performing production rate calibrations as well as calculating erosion rates and surface exposure ages for single samples and depth profiles. The code is available under the GNU General Public License agreement and can be downloaded and modified to deal with specific atypical scenarios

Abrupt mid-Holocene ice loss in the western Weddell Sea Embayment of Antarctica

The glacial history of the westernmost Weddell Sea sector of Antarctica since the Last Glacial Maximum is virtually unknown, and yet it has been identified as critical for improving reliability of glacio-isostatic adjustment models that are required to correct satellite-derived estimates of ice sheet mass balance. Better knowledge of the glacial history of this region is also important for validating ice sheet models that are used to predict future contribution of the Antarctic ice sheet to sea level rise. Here we present a new Holocene deglacial chronology from a site on the Lassiter Coast of the Antarctic Peninsula, which is situated in the western Weddell Sea sector. Samples from 12 erratic cobbles and 18 bedrock surfaces from a series of presently-exposed ridges were analysed for cosmogenic 10Be exposure dating, and a smaller suite of 7 bedrock samples for in situ 14C dating. The resulting 10Be ages are predominantly in the range 80–690 ka, whereas bedrock yielded much younger in situ 14C ages, in the range 6.0–7.5 ka for samples collected from 138–385 m above the modern ice surface. From these we infer that the ice sheet experienced a period of abrupt thinning over a short time interval (no more than 2700 years) in the mid-Holocene, resulting in lowering of its surface by at least 250 m. Any late Holocene change in ice sheet thickness — such as re-advance, postulated by several modelling studies — must lie below the present ice sheet surface. The substantial difference in exposure ages derived from 10Be and 14C dating for the same samples additionally implies ubiquitous 10Be inheritance acquired during ice-free periods prior to the last deglaciation, an interpretation that is consistent with our glacial-geomorphological field observations for former cold-based ice cover. The results of this study provide evidence for an episode of abrupt ice sheet surface lowering in the mid-Holocene, similar in rate, timing and magnitude to at least two other locations in Antarctica

New Last Glacial Maximum Ice Thickness constraints for the Weddell Sea Embayment, Antarctica

We describe new Last Glacial Maximum (LGM) ice thickness constraints for three locations spanning the Weddell Sea Embayment (WSE) of Antarctica. Samples collected from the Shackleton Range, Pensacola Mountains, and the Lassiter Coast constrain the LGM thickness of the Slessor Glacier, Foundation Ice Stream, and grounded ice proximal to the modern Ronne Ice Shelf edge on the Antarctic Peninsula, respectively. Previous attempts to reconstruct LGM-to-present ice thickness changes around the WSE used measurements of long-lived cosmogenic nuclides, primarily Be-10. An absence of post-LGM apparent exposure ages at many sites led to LGM thickness reconstructions that were spatially highly variable and inconsistent with flow line modelling. Estimates for the contribution of the ice sheet occupying the WSE at the LGM to global sea level since deglaciation vary by an order of magnitude, from 1.4 to 14.1m of sea level equivalent. Here we use a short-lived cosmogenic nuclide, in situ-produced C-14, which is less susceptible to inheritance problems than Be-10 and other long-lived nuclides. We use in situ C-14 to evaluate the possibility that sites with no post-LGM exposure ages are biased by cosmogenic nuclide inheritance due to surface preservation by cold-based ice and non-deposition of LGM-aged drift. Our measurements show that the Slessor Glacier was between 310 and up to 655m thicker than present at the LGM. The Foundation Ice Stream was at least 800m thicker, and ice on the Lassiter Coast was at least 385m thicker than present at the LGM. With evidence for LGM thickening at all of our study sites, our in situ C-14 measurements indicate that the long-lived nuclide measurements of previous studies were influenced by cosmogenic nuclide inheritance. Our inferred LGM configuration, which is primarily based on minimum ice thickness constraints and thus does not constrain an upper limit, indicates a relatively modest contribution to sea level rise since the LGM of < 4.6 m, and possibly as little as < 1.5 m

Validation of earthquake ground-motion models in southern California, USA, using precariously balanced rocks

Accurate estimates of earthquake ground shaking rely on uncertain ground-motion models derived from limited instrumental recordings of historical earthquakes. A critical issue is that there is currently no method to empirically validate the resultant ground-motion estimates of these models at the timescale of rare, large earthquakes; this lack of validation causes great uncertainty in ground-motion estimates. Here, we address this issue and validate ground-motion estimates for southern California utilizing the unexceeded ground motions recorded by 20 precariously balanced rocks. We used cosmogenic 10Be exposure dating to model the age of the precariously balanced rocks, which ranged from ca. 1 ka to ca. 50 ka, and calculated their probability of toppling at different ground-motion levels. With this rock data, we then validated the earthquake ground motions estimated by the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) seismic-source characterization and the Next Generation Attenuation (NGA)-West2 ground-motion models. We found that no ground-motion model estimated levels of earthquake ground shaking consistent with the observed continued existence of all 20 precariously balanced rocks. The ground-motion model I14 estimated ground-motion levels that were inconsistent with the most rocks; therefore, I14 was invalidated and removed. At a 2475 year mean return period, the removal of this invalid ground-motion model resulted in a 2−7% reduction in the mean and a 10−36% reduction in the 5th−95th fractile uncertainty of the ground-motion estimates. Our findings demonstrate the value of empirical data from precariously balanced rocks as a validation tool for removing invalid ground-motion models and, in turn, reducing the uncertainty in earthquake ground-motion estimates

Glacial geology of the Hudson Mountains, Amundsen Sea sector, West Antarctica

The Hudson Mountains are situated in the eastern Amundsen Sea sector of the West Antarctic Ice Sheet, adjacent to Pine Island Glacier. They form a volcanic field of 17 stratovolcanoes and parasitic vents, preserved as nunataks. Two former tributaries of Pine Island Glacier (Larter and Lucchitta glaciers) flow through the mountains. Here we present a detailed study of the glacial geology of the area. We describe field observations and measurements of geomorphological features from 15 of the nunataks, meltwater ponds found on the surface of three nunataks and supraglacial features (ice dolines) from two sites near the present grounding line. Together these provide constraints on the past ice sheet extent, flow pathways and thermal regime, and enhance our understanding of the present hydrological regime – all of which are important as context for the observed modern ice sheet behaviour. We find evidence suggesting that all nunataks in the Hudson Mountains were covered by ice during the Last Glacial Maximum (defined here as 26.5-19 ka; Clark et al., 2009) and have since deglaciated. Faceted and polished erratic cobbles and boulders of exotic lithologies (syenites, alkali granites, granites, granodiorites, tonalites and gabbros) are numerous and perched on nunatak surfaces. A marked difference between the dominant erratic lithologies on nunataks adjacent to Pine Island Glacier (granite) and Lucchitta Glacier (granodiorite-tonalite) indicates that the ice sheet was transporting clasts from at least two distinct upstream source regions. The similarity in degree of weathering suggests, however, that all the erratics were transported by one phase of (warm-based) glaciation; their presence on or close to the summits of all except one nunatak indicates that the ice sheet during that time was at least 700 m thicker than present. These results are consistent with ice sheet model simulations which suggest that all nunataks in the Hudson Mountains were completely submerged by the Last Glacial Maximum ice sheet