Model code supporting "Evidence for a more extensive Greenland Ice Sheet in southwestern Greenland during the Last Glacial Maximum"

This item contains the MATLAB code and outputs (Excel files) supporting the analysis described in the abstract. All contributions provided by Chris Sbarra, December 2020. The maximum extent and elevation of the Greenland Ice Sheet (GrIS) in southwestern Greenland during the Last Glacial Maximum (LGM...

Full description

Bibliographic Details
Main Authors: Sbarra, Christopher M., Briner, Jason P., Graham, Brandon L., Poinar, Kristin, Thomas, Elizabeth K., Young, Nicolás E.
Format: Software
Language:English
Published: 2021
Subjects:
Online Access:http://hdl.handle.net/10477/82607
Description
Summary:This item contains the MATLAB code and outputs (Excel files) supporting the analysis described in the abstract. All contributions provided by Chris Sbarra, December 2020. The maximum extent and elevation of the Greenland Ice Sheet (GrIS) in southwestern Greenland during the Last Glacial Maximum (LGM, 26-19.5 ka) is poorly constrained. Yet the size of the GrIS during the LGM helps inform estimates of past ice sheet sensitivity to climate change and provides benchmarks for ice sheet modeling. Reconstructions of LGM ice extents vary between an inner-continental-shelf minimum, a mid-shelf position, and a maximum extent at the shelf break. We use three approaches to resolve LGM ice extent in the Sisimiut sector of southwestern Greenland. First, we explore the likelihood of minimum vs. maximum GrIS reconstructions using existing relative sea level data. We use an empirical relationship between marine limit elevation and distance to LGM terminus established from other Northern Hemisphere Pleistocene ice sheets as context for interpreting marine limit data in southwestern Greenland. Our analysis supports a maximum regional GrIS extent to the shelf break during the LGM. Second, we apply a simple 1D crustal rebound model to simulate relative sea level curves for contrasting ice sheet sizes and compare these simulated curves with existing relative sea level data. The only realistic ice sheet configuration resulting in relative sea level model-data fit suggests that the GrIS terminated at the shelf break during the LGM. Lastly, we constrain the LGM ice sheet thickness using cosmogenic 10Be, 26Al and 14C exposure dating from two summit areas, one at 381 m above sea level at the coast, and another at 798 m asl 32 km inland. Twenty-four cosmogenic radionuclide measurements, combined with results of our first two approaches, reveal that our targeted summits were likely ice-covered during the LGM and became deglaciated at ~11.6 ka. Inventories of in situ 14C in bedrock at one summit point to a small degree of inherited 14C and ...