whigg/fjords_grl_2022

Code for Slater 2022, GRL, Characteristic depths, fluxes and timescales for Greenland's tidewater glacier fjords from subglacial discharge-driven upwelling during summer

fjords_GRL_2022

This directory contains code and data required to reproduce the results in the paper:

Characteristic depths, fluxes and timescales for Greenland's tidewater glacier fjords from subglacial discharge-driven upwelling during summer
D. A. Slater, D. Carroll, H. Oliver, M. J. Hopwood, F. Straneo, M. Wood, J. K. Willis, M. Morlighem
Geophysical Research Letters, 2022

Code

It contains the following code files (all MATLAB):

process_calvingfronts.m - defines calving fronts
process_fjords.m - defines fjords
process_runoff.m - gets subglacial runoff from Mankoff et al. 2020
process_glaciers.m - defines glaciers
link_data.m - brings datasets together (calving fronts, fjords, glaciers, CTDs)
process_plume.m - runs plume model for all glaciers
process_seasonality.m - investigates plume seasonality

rho.m - plume model equation of state
run_plume.m - plume model top level code

makeplots.m - makes plots used in the manuscript
maketable.m - makes the csv data table in the SI
fjord_plots.m - makes fjord-by-fjord map plots in the SI

polarstereo_inv.m - converts polar stereographic to lat/lon
latlon2utm.m - converts lat/lon to polar stereographic
customlegend.m - simple code for custom plot legend
errplot.m - a simple custom plotting function called in makeplots.m
plotbox.m - another simple custom plotting function called in makeplots.m
saveplot.m - function for saving plots as png
saveplot_pdf.m - function for saving plots as pdf

Required data

It contains the following data files (all .mat):

ice_ocean_sectors.mat - sector polygons derived in Slater et al. 2019
all_OMG_CTD_from_csv_ver_1.3.mat - OMG CTD data
carroll_2016_CTDs.mat - CTD data from Carroll et al. 2016
mortensen_meire.mat - CTD data from Mortensen & Meire 2020
sarqardleq.mat - CTD data from Straneo 2019

The following files are used but are not included in this directory:

BedMachine Greenland v4
Plotting colormaps

Created outputs

The code will create the following outputs

calvingfronts.mat - these are the calving fronts created by process_calvingfronts.m
fjords.mat - this is the fjord mask created by process_fjords.m
runoff.mat - subglacial runoff resulting from process_runoff.m
glaciers.mat - this is the glaciers structure created by process_glaciers.m
twglaciers.mat - this is the main output, containing the 136 tidewater glaciers and their associated plume properties etc

Running order

To reproduce the results, run the scripts in the order process_calvingfronts.m
process_fjords.m
process_runoff.m
process_glaciers.m
link_data.m
process_plume.m
process_seasonality.m

To make the plots in the paper and SI, run makeplots.m and/or fjord_plots.m. To make the data table run maketable.m.