Seasonal variability in ice velocity driven by subglacial hydrology of Drang Drung Glacier, Western Himalayas: Modeling scripts and data

Description of Files Contained in this Repository

Contact: Aleah Sommers

aleah.n.sommers@dartmouth.edu

Please note that these scripts and model output are in the format of the Ice-sheet and Sea-level System Model (ISSM). For handling, you can download and install ISSM pre-compiled binaries from this website: issm.jpl.nasa.gov/download/

Model output:

• DD_SHAKTI_winter_1yr.mat – SHAKTI subglacial hydrology spin-up simulation with zero velocity and zero meltwater inputs

• DD_coupled_winter_paper.mat – SHAKTI-ISSM coupled winter spin-up simulation

• DD_coupled_seasonal_paper.mat – SHAKTI-ISSM model output for seasonal coupled simulation

Matlab scripts to run ISSM/SHAKTI:

• runme_DD_setup.m – Set up model domain, run initial SHAKTI simulation

• runme_transientrestart.m – Continue SHAKTI simulation from end of a previous simulation, or turn on coupling to velocity and continue

• runme_DD_seasonal.m – Beginning from end of winter coupled simulation, run a seasonal simulation

Supporting data:

• DD.exp – Outline of model domain

• DD_terminus.exp – Designates terminus for boundary conditions

• DD_DEM.mat – Surface elevation (from LEGOS and OMP, 2025)

• DD_thickness.mat – Ice thickness (from Millan et al., 2022)

• DD_annual_vel.mat – Observed velocity

• Centerline_drangdrung.shp – Centerline coordinates

• DD_velocities_ZoneA_m_per_year.csv – Observed velocities, Zone A

• DD_velocities_ZoneB_m_per_year.csv – Observed velocities, Zone B

Plotting scripts:

• plot_DD_timeseries_with_obs.m

• plot_DD_seasonal_spatial.m

• plot_DD_SI_figs.m

• plot_DD_spatial_with_obs.m

• calc_DD_centerline.m

• read_plot_obs_vel_DD.m