Subglacial hydrology within the Amery Ice Shelf catchment

These datasets correspond to the results from the manuscript:

Characterizing subglacial hydrology within the Amery Ice Shelf catchment using numerical modelling and satellite altimetry (Wearing et al.,)

Martin G. Wearing^1,2, Christine F. Dow³, Daniel N. Goldberg², Noel Gourmelen², Anna E. Hogg⁴, Livia Jakob^5
1ESA Centre for Earth Observation, ESRIN, Frascati, Rome, Italy
²School of Geosciences, University of Edinburgh, Edinburgh, UK
³Department of Geography and Environmental Management, University of Waterloo, Ontario, Canada
⁴School of Earth and Environment, University of Leeds, Leeds, UK
⁵Earthwave Ltd, Edinburgh, UK

The datasets are:

1) Subglacial melt rates. These are calculated using:

- basal dissipation from numerical ice sheet inversion with the ice-flow model STREAMICE (https://mitgcm.readthedocs.io/en/latest/phys_pkgs/streamice.html) run at 2km and 5km resolution.

- geothermal heat flux estimates from Shen et al., (2020) and Martos et al., (2017)

- englacial temperature from Van Liefferinge & Pattyn (2013)

Variables include in each file are:

Position coordinates (X (m) & Y (m)) in polar stereographic projection WGS 1984

Geothermal heat flux (GHF (W/m^2))

Basal friction dissipation (Basal_disp (W/m^2))

Vertical conduction (Vert_Cond (W/m^2))

Subglacial melt rate (Melt_rate (m/yr))

2) Subglacial hydrology. This is determined using the GlaDS model (Werder et al., 2013) with high and low channel conductivity and subglacial melt rate calculated using Shen GHF and 2km resolution basal dissipation.

Model results are provided in netcdf format for high and low conductivity:

Amery_subglacial_hydrology_GlaDS_high_conductivity.nc

Amery_subglacial_hydrology_GlaDS_low_conductivity.nc

Variables includes in each file are:

Node X coordinate (m)

Node Y coordinate (m)

Edge X coordinate (m)

Edge Y coordinate (m)

Channelized discharge (m3/s)

Effective pressure (Pa)

Subglacial water pressure (Pa)

Sheet discharge (x-direction) (m3/s)

Sheet discharge (y-direction) (m3/s)

sheet thickness (m)

Channel cross-sectional area (m2)

3) Ice-shelf basal melting is calculated from CryoSat-2 interferometric-swath radar altimetry acquired from 2010 to 2020, using the mass conservation approach (Gourmelen et al., 2017). Melt rates and errors are provided in the file: Amery_ice_shelf_basal_melt_rate.nc

4) The outlines of active subglacial lakes detected from CryoSat-2 interferometric-swath radar altimetry acquired from 2010 to 2020, using the method of Malczyk et al. (2020) are provided as geojson files in the zipfile: Amery_subglacial_lakes.zip