Search for: All records

We apply geologic evidence from ice-free areas in Antarctica to evaluate model simulations of ice sheet response to warm climates. This is important because such simulations are used to predict ice sheet behaviour in future warm climates, but geologic evidence of smaller-than-present past ice sheets is buried under the present ice sheet and therefore generally unavailable for model benchmarking. We leverage an alternative accessible geologic dataset for this purpose: cosmogenic-nuclide concentrations in bedrock surfaces of interior nunataks. These data produce a frequency distribution of ice thickness over multimillion-year periods, which is also simulated by ice sheet modelling. End-member transient models, parameterized with strong and weak marine ice sheet instability processes and ocean temperature forcings, simulate large and small sea-level impacts during warm periods and also predict contrasting and distinct frequency distributions of ice thickness. We identify regions of Antarctica where predicted frequency distributions reveal differences in end-member ice sheet behaviour. We then demonstrate that a single comprehensive dataset from one bedrock site in West Antarctica is sufficiently detailed to show that the data are consistent only with a weak marine ice sheet instability end-member, but other less extensive datasets are insufficient and/or ambiguous. Finally, we highlight locations where collecting additional data could constrain the amplitude of past and therefore future response to warm climates. 

Pine Island Glacier, West Antarctica, is the largest Antarctic contributor to global sea-level rise and is vulnerable to rapid retreat, yet our knowledge of its deglacial history since the Last Glacial Maximum is based largely on marine sediments that record a retreat history ending in the early Holocene. Using a suite of 10Be exposure ages from onshore glacial deposits directly adjacent to Pine Island Glacier, we show that this major glacier thinned rapidly in the early to mid-Holocene. Our results indicate that Pine Island Glacier was at least 690 m thicker than present prior to ca. 8 ka. We infer that the rapid thinning detected at the site farthest downstream records the arrival and stabilization of the retreating grounding line at that site by 8–6 ka. By combining our exposure ages and the marine record, we extend knowledge of Pine Island Glacier retreat both spatially and temporally: to 50 km from the modern grounding line and to the mid-Holocene, providing a data set that is important for future numerical ice-sheet model validation.

 

Abstract. Cosmogenic-nuclide concentrations in subglacial bedrock cores show that the West Antarctic Ice Sheet (WAIS) at a site between Thwaites and Pope glaciers was at least 35 m thinner than present in the past several thousand years and then subsequently thickened. This is important because of concern that present thinning and grounding line retreat at these and nearby glaciers in the Amundsen Sea Embayment may irreversibly lead to deglaciation of significant portions of the WAIS, with decimeter- to meter-scale sea level rise within decades to centuries. A past episode of ice sheet thinning that took place in a similar, although not identical, climate was not irreversible. We propose that the past thinning–thickening cycle was due to a glacioisostatic rebound feedback, similar to that invoked as a possible stabilizing mechanism for current grounding line retreat, in which isostatic uplift caused by Early Holocene thinning led to relative sea level fall favoring grounding line advance. 

Abstract. We collected a debris-rich ice core from a buried icemass in Ong Valley, located in the Transantarctic Mountains in Antarctica. Wemeasured cosmogenic nuclide concentrations in quartz obtained from the icecore to determine the age of the buried ice mass and infer the processesresponsible for the emplacement of the debris currently overlaying the ice.Such ice masses are valuable archives of paleoclimate proxies; however, thepreservation of ice beyond 800 kyr is rare, and therefore much effort hasbeen 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 thecosmic-ray-produced nuclides 10Be, 21Ne, and 26Al in quartzin the englacial and supraglacial debris and compare the model predictionsto measured nuclide concentrations in order to further constrain the age.Large downcore variation in measured cosmogenic nuclide concentrationssuggests that the englacial debris is sourced both from subglacially derivedmaterial and recycled paleo-surface debris that has experienced surfaceexposure prior to entrainment. We find that the upper section of the icecore is 2.95 + 0.18 / −0.22 Myr old. The average ice sublimation rate duringthis time period is 22.86 + 0.10 / −0.09 m Myr−1, and the surfaceerosion rate of the debris is 0.206 + 0.013 / −0.017 m Myr−1. Burialdating of the recycled paleo-surface debris suggests that the lower sectionof the ice core belongs to a separate, older ice mass which we estimate tobe 4.3–5.1 Myr old. The ages of these two stacked, separate ice masses canbe directly related to glacial advances of the Antarctic ice sheet andpotentially coincide with two major global glaciations during the early andlate Pliocene epoch when global temperatures and CO2 were higher thanpresent. These ancient ice masses represent new opportunities for gatheringancient climate information. 

We assess if variations in the in situ cosmogenic 26Al/10Be production ratio expected from nuclear physics are consistent with empirical data, knowledge critical for two-isotope studies. We do this using 313 samples from glacially transported boulders or scoured bedrock with presumed simple exposure histories in the Informal Cosmogenic-nuclide Exposure-age Database (ICE-D) from latitudes between 53°S to 70°N and altitudes up to 5000 m above sea level. Although there were small systematic differences in Al/Be ratios measured in different laboratories, these were not significant and are in part explained by differences in elevation distribution of samples analyzed by each laboratory. We observe a negative correlation between the 26Al/10Be production ratio and elevation (p = 0.0005), consistent with predictions based on the measured energy dependence of nuclear reaction cross-sections and the spatial variability in cosmic-ray energy spectra. We detect an increase in the production ratio with increasing latitude, but this correlation is significant only in a single variate model, and we attribute at least some of the correlation to sample elevation bias because lower latitude samples are typically from higher elevations (and vice versa). Using 6.75 as the 26Al/10Be production ratio globally will bias two-isotope results at higher elevations and perhaps higher latitudes. Data reported here support using production rate scaling that incorporates such ratio changes, such as the LSDn scheme, to minimize such biases. 

Abstract. Diffusion properties of cosmogenic 3He in quartz at Earth surface temperatures offer the potential to directly reconstruct the evolution of pastin situ temperatures from formerly glaciated areas, which is important information for improving our understanding of glacier–climateinteractions. In this study, we apply cosmogenic 3He paleothermometry to rock surfaces gradually exposed from the Last Glacial Maximum(LGM) to the Holocene period along two deglaciation profiles in the European Alps (Mont Blanc and Aar massifs). Laboratory experiments conducted onone representative sample per site indicate significant differences in 3He diffusion kinetics between the two sites, with quasi-linearArrhenius behavior observed in quartz from the Mont Blanc site and complex Arrhenius behavior observed in quartz from the Aar site, which weinterpret to indicate the presence of multiple diffusion domains (MDD). Assuming the same diffusion kinetics apply to all quartz samples along eachprofile, forward model simulations indicate that the cosmogenic 3He abundance in all the investigated samples should be at equilibrium withpresent-day temperature conditions. However, measured cosmogenic 3He concentrations in samples exposed since before the Holocene indicate anapparent 3He thermal signal significantly colder than today. This observed 3He thermal signal cannot be explained with a realisticpost-LGM mean annual temperature evolution in the European Alps at the study sites. One hypothesis is that the diffusion kinetics and MDD modelapplied may not provide sufficiently accurate, quantitative paleo-temperature estimates in these samples; thus, while a pre-Holocene 3Hethermal signal is indeed preserved in the quartz, the helium diffusivity would be lower at Alpine surface temperatures than our diffusion modelspredict. Alternatively, if the modeled helium diffusion kinetics is accurate, the observed 3He abundances may reflect a complexgeomorphic and/or paleoclimatic evolution, with much more recent ground temperature changes associated with the degradation of alpine permafrost. 

The rapidly retreating Thwaites and Pine Island glaciers together dominate present-day ice loss from the West Antarctic Ice Sheet and are implicated in runaway deglaciation scenarios. Knowledge of whether these glaciers were substantially smaller in the mid-Holocene and subsequently recovered to their present extents is important for assessing whether current ice recession is irreversible. Here we reconstruct relative sea-level change from radiocarbon-dated raised beaches at sites immediately seawards of these glaciers, allowing us to examine the response of the earth to loading and unloading of ice in the Amundsen Sea region. We find that relative sea level fell steadily over the past 5.5 kyr without rate changes that would characterize large-scale ice re-expansion. Moreover, current bedrock uplift rates are an order of magnitude greater than the rate of long-term relative sea-level fall, suggesting a change in regional crustal unloading and implying that the present deglaciation may be unprecedented in the past ~5.5 kyr. While we cannot preclude minor grounding-line fluctuations, our data are explained most easily by early Holocene deglaciation followed by relatively stable ice positions until recent times and imply that Thwaites and Pine Island glaciers have not been substantially smaller than present during the past 5.5 kyr.

 

Abstract. Evidence for the timing and pace of past grounding lineretreat of the Thwaites Glacier system in the Amundsen Sea embayment (ASE)of Antarctica provides constraints for models that are used to predict thefuture trajectory of the West Antarctic Ice Sheet (WAIS). Existingcosmogenic nuclide surface exposure ages suggest that Pope Glacier, a formertributary of Thwaites Glacier, experienced rapid thinning in the early tomid-Holocene. There are relatively few exposure ages from the lower ice-freesections of Mt. Murphy (<300 m a.s.l.; metres above sea level) that are uncomplicated byeither nuclide inheritance or scatter due to localised topographiccomplexities; this makes the trajectory for the latter stages ofdeglaciation uncertain. This paper presents 12 new 10Be exposure agesfrom erratic cobbles collected from the western flank of Mt. Murphy, within160 m of the modern ice surface and 1 km from the present grounding line.The ages comprise two tightly clustered populations with mean deglaciationages of 7.1 ± 0.1 and 6.4 ± 0.1 ka (1 SE). Linear regressionanalysis applied to the age–elevation array of all available exposure agesfrom Mt. Murphy indicates that the median rate of thinning of Pope Glacierwas 0.27 m yr−1 between 8.1–6.3 ka, occurring 1.5 times faster thanpreviously thought. Furthermore, this analysis better constrains theuncertainty (95 % confidence interval) in the timing of deglaciation atthe base of the Mt. Murphy vertical profile (∼ 80 m above themodern ice surface), shifting it to earlier in the Holocene (from 5.2 ± 0.7 to 6.3 ± 0.4 ka). Taken together, the results presentedhere suggest that early- to mid-Holocene thinning of Pope Glacier occurredover a shorter interval than previously assumed and permit a longer durationover which subsequent late Holocene re-thickening could have occurred. 

Abstract. We use 25 new measurements of in situ produced cosmogenic 26Al and 10Bein 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 quartzextracted from central Cuban river sand range from3.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 soluteconcentrations and modeled runoff, exceed measured cosmogenic-10Be-derived erosion rates in 18 of 23 basins. This disparity mandatesthat in this environment landscape-scale mass loss is not fully representedby 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 consistentwith the presence of a deep, mixed, regolith layer providing extendedstorage times on slopes and/or burial and extended storage during fluvialtransport. 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.