A novel approach to constraining ice sheet models with glacial geomorphology

Lead Research Organisation: University of Sheffield
Department Name: Geography

Abstract

The ice sheets in Greenland and Antarctica hold enough water to raise global sea level by ~65 m, and are currently losing mass at an accelerating rate. Predictions of future mass loss are difficult to make, meaning that the contribution of ice sheets is a large source of uncertainty for estimates of future sea level rise. The numerical models which make these predictions need observations of ice sheet behaviour for validation and improvement, but the record of extant ice sheets is short. A longer record of ice sheet behaviour can be obtained by studying the regions occupied by ice sheets during the last glacial period (~120,000 to ~10,000 years ago). The ice sheets which existed during this time left behind abundant evidence for their dynamic behaviour in the form of glacial landforms. In fields such as weather forecasting, model improvement has been achieved by hindcasting - the act of simulating past conditions and comparing them to historical data. However, simulations of palaeo-ice sheets rarely utilise the glacial landform record to its full extent. This project is a deliberate attempt to bring together the rather disparate disciplines of ice sheet modelling and glacial geomorphology to improve hindcasting procedures. With this motivation, my aim is to develop and apply tools for validating ice sheet models with glacial geomorphology.
During the last glacial, the Greenland and Antarctic ice sheets were larger and much of North America, Northwest Europe and Patagonia was glaciated. Information regarding the behaviour of these major ice sheets is currently disparately stored. I will collate information regarding the ice sheets that covered these five regions into a single standardised database. From this, I will produce the first global scale empirical reconstruction of ice sheet extent and behaviour during the last glacial. I will also use a state of the art ice sheet model to simulate the behaviour of these ice sheets. The ice sheets will be simulated hundreds of times, each simulation with a separate set of input values, to capture the full range of likely ice configurations given the current limitations of our knowledge regarding how ice sheets flow and past climate conditions.
To evaluate my model simulations, I will develop a new set of tools which will quantify how well the model predicts the properties of palaeo-ice sheets inferred from the dataset of glacial geomorphology. To filter out extremely unrealistic model simulations, I will first compare model output to my reconstructed ice sheet extent. An overall score of model performance will then be calculated from the ability to reproduce the timing, flow direction, changes in flow pattern, margin position, erosive capabilities and basal pressure gradients that are recorded by glacial geomorphology.
I will then combine the simulations which best conform to the geomorphological data to produce a high resolution simulation of the behaviour of the major ice sheets which existed during the last glacial cycle. The simulated ice sheets will serve as important analogues for existing ice sheets, demonstrating how controls on ice sheet behaviour such as topography, sea level and climate conditioned these ice sheets to retreat, and how their changing geometries contributed to past rapid sea level rise, such as meltwater pulses and major iceberg discharge events. My model-data comparison procedure will also identify regions where the model struggles to replicate the palaeo-ice sheet behaviour. This will highlight deficiencies in the model, providing targets for model development. Furthermore, models of the future Antarctic and Greenland ice sheets rely upon accurately simulating their past extent and behaviour. The procedures and model experiments I conduct will improve how ice sheet models adhere to this history. Therefore, my project will bring about a step-change in how ice sheet models are tested and calibrated and offer a new framework for other researchers to utilise.

Planned Impact

There are 4 groups who will benefit from this project:
1. Scientists interested in the past, present and future of ice sheets.
2. Policymakers
3. Teachers and school Children
4. The General Public and the Media
How will they benefit?
1. Scientists will benefit from improved methods for validating, constraining and testing ice sheet models. How the improved understanding of past, present and future ice sheets will benefit scientists is outlined in the 'Academic Beneficiaries' section of this proforma.
2. Rising sea levels are a globally significant concern, and ice sheets are set to be the main contributors to this trend. Policymakers will benefit from the reduced uncertainty in ice sheet models, and thus improved accuracy in predictions of sea level rise, brought about by this project. The Intergovernmental Panel on Climate Change (IPCC) provides policymakers with information from the scientific literature, which I will contribute to. I will also attend relevant NERC workshops to engage with policy makers.
3. Glacial geomorphology and ice sheets are now part of the GCSE and A-Level curriculum. Schools will benefit from access to the geomorphological data generated by this project and the model simulations of past-ice sheets produced. Schools will be reached through the Geographical Association, who provide school teachers with teaching materials and professional advice.
4. Ice sheets have been the focus of much media attention, highlighting the interest of the general public in this topic. Indeed, some of my previous work (Kingslake et al., 2017; Nature) received widespread international media coverage (The Washington Post, Buzzfeed, The Independent). As palaeo-ice sheets used to cover populous areas, I envisage that the improved picture of ice sheet behaviour in the past developed in this project will further engage the media and the public. This will provide them with an improved understanding of the role of ice sheets in sculpting the landscapes surrounding them, and the importance of contemporary ice sheets. The general public will be directly targeted by the development of website materials, leaflets and phone-based applications which explain the glaciation of the Yorkshire Dales National Park. The public will also be engaged through NERC public engagement events.

Publications

10 25 50

publication icon
Barr I (2019) The dynamics of mountain erosion: Cirque growth slows as landscapes age in Earth Surface Processes and Landforms

 
Description The ice sheets we have in the world today, in Greenland and Antarctica, are losing ice. This loss of ice storage in the ice sheets, contributes to global sea level rise. To predict how these ice masses are going to respond to our warming world, we need reliable computer models of ice flow. Currently, these models are trained to limited observations of Greenland and Antarctica. To expand the database of observations, my fellowship focussed on the evidence left behind by ice sheets that existed during the last glacial period (roughly 100,000 to 10,000 years ago). As these ice sheets grew and shrunk over Canada, Northern Europe and Patagonia, they left behind an archive which allows us to infer ice behaviour in changing climates. This is stored in landforms, sediments and rocks that have been created and deposited by these past ice sheets. In my fellowship, I have developed tools for testing ice flow models against this archive, creating novel approaches toward combining the previously disparate observational evidence from models. This has led to knew insights into the behaviour of these past ice sheets, and has created a framework for calibrating predictive models of Greenland and Antarctica against past ice extent data.
Exploitation Route The tools for constraining models are being used by others. The sea level community are using my predictions of palaeo-British ice (Bradley, Ely et al., 2023 for example).
Sectors Environment

 
Title Data for: Comparing observations of subglacial linations to the output of ice sheet models: The Likelihood of Accordant Lineations Analysis (LALA) Tool 
Description NetCDF file of a model simulation of the former British Irish Ice Sheet. Supplement to: Archer, R.E., Ely, J.C., Heaton, T.J., Butcher, F.E.G., Hughes, A.L.C. and Clark, C.D. Comparing observations of subglacial linations to the output of ice sheet models: The Likelihood of Accordant Lineations Analysis (LALA) Tool. Earth Surface Processes and Landforms The folder is provided as a .tar.gz file. This can be uncompressed using a command (linux) such as: tar -xf archive.tar.gz or in Windows using 7zip. Within the folder is output from a numerical model simulation of the former British Irish Ice Sheet, modelled using PISM: https://www.pism.io/ The netcdf file contains several variables, output every 100 years: thk - ice thickness uvel - ice velocity in the x direction vvel - ice velocity in the y direction mask - a mask of ice cover This can be viewed using the ncview software and manipulated in python using packages such as netCDF4 and xarray. 
Type Of Material Database/Collection of data 
Year Produced 2023 
Provided To Others? Yes  
URL https://figshare.shef.ac.uk/articles/dataset/Data_for_Comparing_observations_of_subglacial_linations...
 
Title Data for: Exploring the ingredients required to successfully model the placement, generation, and evolution of ice streams in the British-Irish Ice Sheet 
Description Contained are files related to Gandy et al., "Exploring the ingredients required to successfully model the placement, generation, and evolution of ice streams in the British-Irish Ice Sheet". Directory "inputs" contains required required files for reproducing experiments described. Directory "advance_retreat" contains .hdf5 output files of the ADVANCE and RETREAT experiments. Directory "advance_retreat_geo" contains .hdf5 output files of the ADVANCE_GEOTHERMAL and RETREAT_GEOTHERMAL experiments. Details on viewing BISICLES output files are found here; http://davis.lbl.gov/Manuals/BISICLES-DOCS/index.html 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes  
URL https://data.mendeley.com/datasets/w5tn37pkym
 
Title Data for: Exploring the ingredients required to successfully model the placement, generation, and evolution of ice streams in the British-Irish Ice Sheet 
Description Contained are files related to Gandy et al., "Exploring the ingredients required to successfully model the placement, generation, and evolution of ice streams in the British-Irish Ice Sheet". Directory "inputs" contains required required files for reproducing experiments described. Directory "advance_retreat" contains .hdf5 output files of the ADVANCE and RETREAT experiments. Directory "advance_retreat_geo" contains .hdf5 output files of the ADVANCE_GEOTHERMAL and RETREAT_GEOTHERMAL experiments. Details on viewing BISICLES output files are found here; http://davis.lbl.gov/Manuals/BISICLES-DOCS/index.html 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes  
URL https://data.mendeley.com/datasets/w5tn37pkym/1
 
Title Dataset from: Quantifying the Uncertainty in the Eurasian Ice-Sheet Geometry at the Penultimate Glacial Maximum (Marine Isotope Stage 6) 
Description The North Sea Last Interglacial sea level is sensitive to the fingerprint of mass loss from polar ice sheets. However, the signal is complicated by the influence of glacial isostatic adjustment driven by Penultimate Glacial Period ice-sheet changes, and yet these ice-sheet geometries remain significantly uncertain. Here, we produce new reconstructions of the Eurasian ice sheet during the Penultimate Glacial Maximum by employing large ensemble experiments from a simple ice-sheet model that depends solely on basal sheer stress, ice extent, and topography. To explore the range of uncertainty in possible ice geometries, we use a parameterised shear-stress map as input that has been developed to incorporate bedrock characteristics and the influence of ice-sheet basal processes. We perform Bayesian uncertainty quantification, utilising Gaussian Process emulation, to calibrate against global ice-sheet reconstructions of the last deglaciation and rule out combinations of input parameters that produce unrealistic ice sheets. The refined parameter space is then applied to the Penultimate Glacial Maximum to create an ensemble of constrained 3D Eurasian ice-sheet geometries. Our reconstructed Penultimate Glacial Maximum Eurasian ice-sheet volume is 48 ± 8 m sle. We find that the Barents-Kara Sea region displays both the largest mean volume and volume uncertainty of 24 ± 8 m sle while the British-Irish sector's volume of 1.7 ± 0.2 m sle is smallest. Our new workflow may be applied to other locations and periods where ice-sheet histories have limited empirical data. 
Type Of Material Database/Collection of data 
Year Produced 2023 
Provided To Others? Yes  
URL https://archive.researchdata.leeds.ac.uk/1183/