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

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

Antarctic ice sheet paleo-constraint database

We present a database of observational constraints on past Antarctic Ice Sheet changes during the last glacial cycle intended to consolidate the observations that represent our understanding of past Antarctic changes, for state-space estimation, and paleo-model calibrations. The database is a major expansion of the initial work of Briggs and Tarasov (2013). It includes new data types and multi-tier data quality assessment. The updated constraint database “AntICE2” consists of observations of past grounded and floating ice sheet extent, past ice thickness, past relative sea level, borehole temperature profiles, and present-day bedrock displacement rates. In addition to paleo-observations, the present-day ice sheet geometry and surface ice velocities are incorporated to constrain the present-day ice sheet configuration. The method by which the data is curated using explicitly defined criteria is detailed. Moreover, the observational uncertainties are specified. The methodology by which the constraint database can be applied to evaluate a given ice sheet reconstruction is discussed. The implementation of the “AntICE2” database for Antarctic Ice Sheet model calibrations will improve Antarctic Ice Sheet predictions during past warm and cold periods and yield more robust paleo model spin ups for forecasting future ice sheet changes

Exposure-age record of Holocene ice sheet and ice shelf change in the northeast Antarctic Peninsula

A

New 10Be exposure ages improve Holocene ice sheet thinning history near the grounding line of Pope Glacier, Antarctica

Evidence for the timing and pace of past grounding line retreat of the Thwaites Glacier system in the Amundsen Sea embayment (ASE) of Antarctica provides constraints for models that are used to predict the future trajectory of the West Antarctic Ice Sheet (WAIS). Existing cosmogenic nuclide surface exposure ages suggest that Pope Glacier, a former tributary of Thwaites Glacier, experienced rapid thinning in the early to mid-Holocene. There are relatively few exposure ages from the lower ice-free sections of Mt. Murphy (<300 m a.s.l.; metres above sea level) that are uncomplicated by either nuclide inheritance or scatter due to localised topographic complexities; this makes the trajectory for the latter stages of deglaciation uncertain. This paper presents 12 new 10Be exposure ages from erratic cobbles collected from the western flank of Mt. Murphy, within 160 m of the modern ice surface and 1 km from the present grounding line. The ages comprise two tightly clustered populations with mean deglaciation ages of 7.1 ± 0.1 and 6.4 ± 0.1 ka (1 SE). Linear regression analysis applied to the age–elevation array of all available exposure ages from Mt. Murphy indicates that the median rate of thinning of Pope Glacier was 0.27 m yr−1 between 8.1–6.3 ka, occurring 1.5 times faster than previously thought. Furthermore, this analysis better constrains the uncertainty (95 % confidence interval) in the timing of deglaciation at the base of the Mt. Murphy vertical profile (∼ 80 m above the modern ice surface), shifting it to earlier in the Holocene (from 5.2 ± 0.7 to 6.3 ± 0.4 ka). Taken together, the results presented here suggest that early- to mid-Holocene thinning of Pope Glacier occurred over a shorter interval than previously assumed and permit a longer duration over which subsequent late Holocene re-thickening could have occurred