Supplementary material for "Numerical issues in modeling ice sheet instabilities such as binge-purge type cyclic ice stream surging"

• Main Author: Kevin Hank
• Format: Other/Unknown Material
• Language: unknown
• Published: 2023
• Subjects: Hudson Bay; Hudson; Hudson Strait; Ice Sheet
• Online Access: https://zenodo.org/record/7529646; https://doi.org/10.5281/zenodo.7529646

Running the models: GSM: A full description of the GSM with the associated code archive and complete documentation will soon be submitted and then linked here. We are happy to provide the code upon individual request. The input files used specifically for "Numerical issues in modeling ice sheet instabilities such as binge-purge type cyclic ice stream surging" are described here. The GSM uses bed masks to set the sediment cover of each grid cell. The bed masks for all four resolutions (25 km, 12.5 km, 6.25 km, and 3.125 km) with an abrupt transition between soft sediment and hard bedrock are: inputGSM01.bedMskLISsqHighRes25kmSymmetric.dat inputGSM02.bedMskLISsqHighRes12p5kmSymmetric.dat inputGSM03.bedMskLISsqHighRes6p25kmSymmetric.dat inputGSM04.bedMskLISsqHighRes3p125kmSymmetric.dat The GSM bed masks with a small (3.125 km) and wide (25 km) transition zone between soft sediment and hard bedrock are: inputGSM05.bedMskLISsqHighRes3p125kmSymmetricSmallSmoothTransition.dat inputGSM06.bedMskLISsqHighRes3p125kmSymmetricSmoothTransition.dat Similar to the sediment cover, the GSM uses a topography mask to set the bed elevation. The topography files with a small (3.125 km) and wide (25 km) transition zone between the pseudo-Hudson Strait and Hudson Bay at 200 m below sea level, the surrounding at sea level, and the ocean grid cells at 500 m below sea level are: inputGSM07.topo3p125kmSymmetricSmallSmooth.dat inputGSM08.topo3p125kmSymmetricSmooth.dat The GSM setup used here requires 9 input parameters. The parameters are read as parameter vectors from a parameter file. The GSM parameter vectors used in this paper are: 1.00 0.30 0.81 1.20 1.30 2.60 1.44 1.20 0.30 (parameter vector 00) 1.00 2.00 0.80 1.20 0.75 1.07 1.00 0.50 0.30 (parameter vector 01) 0.59 1.20 1.00 1.20 0.52 2.87 1.27 0.82 0.20 (parameter vector 02) 1.00 2.00 1.00 0.54 0.27 1.52 1.17 0.95 0.20 (parameter vector 03) 1.00 1.50 1.00 1.20 0.60 1.52 0.97 0.60 0.10 (parameter vector 04) C_rmu, C_fslid, lapsr, PDDmelt, hpre, PrecRef, rTnorth, n_b PISM: ...