Impacts of bed topography resolution on sea-level rise projections from coupled subglacial hydrology and ice dynamics for Thwaites Glacier, Antarctica

Ice sheet models require explicit knowledge of the underlying bed. However, much remains unknown regarding the subglacial environment due to difficulties associated with measuring it. Extensive radar surveys have been conducted across Antarctica, but the requirement of full-coverage bed topography for models necessitates interpolation over gaps between existing observations, which often span kilometers or more. Advances in modeling capabilities now allow for the application of dynamic coupling between subglacial hydrology and ice dynamics in models of Antarctica. While a bed resolution of ~1 km is recommended for modeling Antarctic ice dynamics, it has been suggested that finer spatial resolutions are necessary to resolve subglacial water flow. We use a coupled model configuration to generate projections of glacier evolution, including the subglacial hydrologic system, for Thwaites Glacier, West Antarctica, initiated with several different bed topographies. We find that the specific bed..., This dataset includes the complete output from 9 simulations conducted using the Ice-sheet and Sea-level System Model (ISSM). The simulations explore variations in the representation of subglacial hydrology and bed topography resolution. Alongside the model outputs, we include all necessary files and code used to produce the figures presented in the manuscript.
Code/Software
This work was done using the Ice-sheet and Sea-level System Model (ISSM) and the Glacier Drainage System Model (GlaDS). GlaDS is implemented within ISSM, which is open access and available for download at https://issm.jpl.nasa.gov/. These scripts were run using ISSM version 4.24 and MATLAB version R2021b.
, # Impacts of bed topography resolution on sea-level rise projections from coupled subglacial hydrology and ice dynamics for Thwaites Glacier, Antarctica

Dataset DOI: [10.5061/dryad.xd2547dvb](10.5061/dryad.xd2547dvb)

### Description of the data and file structure

The data provided here is organized into 2 folders and 1 .mat file: Results and Figures. The file TG_Initialize_Transient_500mBed_NGlaDS.mat contains the model state just prior to running the model relaxation for the simulation using effective pressure from GlaDS and the 500-meter resolution bed topography. This file contains all fields and initial conditions necessary to launch the transient simulation corresponding to this bed and effective pressure combination. It includes results from the inversion and thermal model steps, as well as all required parameter values.

### Model_Forcing: DATA_repo.zip

#### Results

This folder contains all relatent model output used in our analysis organized into 6 subfolders: Bed, Geometry...,