README v2 –  code and data for obtaining MODIS runoff limits and comparison to and between MAR and RACMO / IMAU-FDM

Differences to v1 of this repository

This version v2 of the repository includes the follwing changes compared to v1:

- contains new and updated versions of the code used in Machguth et al. (2026). Changes to the code were done in response to reviewers' requests for changes to the figures. updated code can be found mainly in the projects "Greenland_RCM_analysis" and "Notebooks". Elsewhere changes are minor or absent.

- provides also the MAR, RACMO and IMAU-FDM data. In the case of the MAR data, most variables that were not used in the study Machguth et al. (2026) were removed. This was done to reduce the data volume of this repository. Note that for easier up and download to and from this repository, the Greenland-wide MAR data were split  into four ZIP files. All included .nc files can be merged again into one folder after download.

- contains two tables (_depth_of_ice_slab_KAN_U.xlsx and _KAN_U_10m_firn_temperatures.xlsx) that were missing from v1 of this repository.

- clarifications, corrections and additional information in the README, as well as a correctly rendered Table 1.

Basic workflow

Table 1 lists the projects, the pieces of code in each project and the key output of each piece of code. This serves at illustrating the basic workflow. More detailed information on required input is provided below. The README files of the various projects provide details on how to use the code.

Table 1: General project overview. The projects and code are shown in the order they are intended to be used. In italic are parts of code that were not used for Machguth et al. (2026) but instead for Machguth et al. (2022). They might not be fully compatible anymore with runoff limits Y_r calculated along flowline-polygons.

Input and output data

All required input data is provided in this repository, with some limitations explained in the following. 

Note 1: Basically all data are provided to (i) run the runoff (Y_obs) limit detection, and (ii) to repeat all the RCM analysis and comparison to MODIS Y_obs. However, the MODIS data for all of Greenland are too volumous and thus only the MODIS data for the so called “test” data set for runoff limit detection are provided. These files include sat_modis_proc_test_l1.zip, sat_modis_proc_test_l2.zip, sat_modis_proc_test_l3.zip and sat_modis_proc_test_l2_NDWI.zip. The complete MODIS data (MOD10A1 and MOD09GA), however, can be downloaded from the data sources mentioned in Section 2.1 of Machguth et al. (2026).

Note 2: While we only provide the MODIS test data for the runoff limit detection, we here provide the full output of the Greenland-wide runoff limit detection.

Note 3: the PDF README provided refers to v1, please disregard.

Input to flowlines/crop_gdalwarp.py

arcticdem_mosaic_100m_v30_greenland_icesheet_geoidCorr.tif

Input to flowlines/flowline_seedpoints.py

seedline_v2.1.shp

Input to flowlines/flowline.py

In the case of using the test data, no seed file exists. Do not provide one, seed points will be calculated along a north-south line as specified.

greenland_vel_mosaic200_2015-2018_vx_v02-composite-crop.tif

greenland_vel_mosaic200_2015-2018_vy_v02-composite-crop.tif

dem_test_gapfilled.tif

In case of calculating flowlines Greenland-wide:

greenland_vel_mosaic500_2015-2018_vx_v02-composite-crop.tif

greenland_vel_mosaic500_2015-2018_vy_v02-composite-crop.tif

arcticdem_mosaic_500m_v30_greenland_icesheet_geoidCorr_GapFilled.tif

seedpoints_v3.4.shp

Input to MODIS_Greenland/MODIS_array_filter.py

/sat_modis_proc_GR_l1test/*

Input to MODIS_Greenland/MODIS_stddev_spatial.py

/sat_modis_proc_test_l2/*

Input to MODIS_Greenland/MODIS_mean_stddev.py

/sat_modis_proc_test_l3/*

dem_test_gapfilled.tif

Input to MODIS_Greenland/ MODIS_NDWI_tiff_to_nc.py

Input data were computed directly on Google Earth Engine, they have not been preserved, only output data exist.

Input to MODIS_Greenland/MODIS_find_slush_limit.py

/sat_modis_proc_test_l3/*

/sat_modis_proc_test_l2_NDWI/*

mask_greenland_icesheet/dem_test_gapfilled.tif

Ys_polygons__test_W20km.shp

flowlines__test_W20km.shp

test_MOD10A1.l3.v4_yrs2000-2021_doy126-136_stddev_median.tif

Input to Greenland_RCM_analysis_prep/modis_vs_mar.py

This code is a direct copy of notebooks that were used to compare Greenland-wide max Y_obs to max Y_rcm MAR. It has not been tested, not all input data might be included in the repository.

/flash/tedstona/MARv.HorstRCMStudy_20240624/*

flowline Polygons

/flash/tedstona/_list_PolyIDs.xlsx

Input to Greenland_RCM_analysis_prep/modis_vs_racmo.py

This code is a direct copy of notebooks that were used to compare Greenland-wide max Y_obs to max Y_rcm MAR. It has not been tested, not all input data might be included in the repository.

/flash/tedstona/RACMO/1km/runoff/*

arcticdem_mosaic_500m_v30_greenland_icesheet_GeoidCorr_GapFilled_RACMO1km.tif

racmo_polys.nc

Input to MODIS_Greenland_analysis/MODIS_A_SLmax.py

_test_slush-limit_output_table.xlsx

or the Greenland-wide output data:

_GR_slush-limit_output_table.xlsx

Input to MODIS_Greenland_analysis/MODIS_A_SLanalyis.py

KAN_U_hourly_v3_fewer_columns.xlsx

KAN_M_hourly_v3_fewer_columns.xlsx

(selected_SL_years_and_stripes_20km.xlsx: not found and was also only used in Machguth et al., 2022)

__test_table_complete_annual_max_SL.xlsx

__test_slush-limit_output_table_OnlyValidEntries.xlsx

or the Greenland-wide output data:

__GR_table_complete_annual_max_SL.xlsx

__GR_test_slush-limit_output_table_OnlyValidEntries.xlsx

Input to MODIS_Greenland_analysis/MODIS_comp_RCM_maxYs.py

It has not been tested whether this file also works with the “test” data.

RACMO2.3p2_ERA5_3h_FGRN055.1km-rlim-RUa1mm.xlsx

MAR-v.20240624-rlim-RUa1mm.xlsx

_stripes_with_aquifers.xlsx

__test_table_simple_annual_max_SL.xlsx

__test_table_complete_annual_max_SL.xlsx

or the Greenland-wide output data:

__GR_table_simple_annual_max_SL.xlsx

__GR_table_complete_annual_max_SL.xlsx

Input to MODIS_Greenland_analysis/MODIS_comp_RCM_Ys.py

It has not been tested whether this file also works with the “test” data.

__GR_slush-limit_output_table_OnlyValidEntries.xlsx

flowlines_daily_rlims_RACMO_1mmEvents_10mmAnnual_2000_2021.xlsx

flowlines_daily_rlims_MAR-v.20240624_1mmEvents_10mmAnnual_2000_2021.xlsx

Input to Greenland_RCM_analysis/RCM_analysis_comp_Ktransect.py

It has not been tested whether this file also works with the “test” data.

/MARv3.14_Greenland-wide_continuous_reduced/*

/FDM_Greenland-wide/* 

/RACMO-FDM_K-Transect_updated/* 

/RACMO_Greenland-wide/* 

/MARv3.14_K-transect_continuous/*

FGRN055_Masks.nc

__GR_slush-limit_output_table_OnlyValidEntries.xlsx

__GR_table_simple_annual_max_SL.xlsx

__GR_table_complete_annual_max_SL.xlsx

Input to _notebooks\K-transect_MAR_Greenland-wide_to_K-transect.ipynb

This notebooks applies the Bresenham algorithm to extract those RCM grid cells that are located closest along the K-Transect (or any other Transect). The output of this notebook is already provided in this repository as /MARv3.14_K-transect_continuous/* . (The input to _notebooks\K-transect_MAR_Greenland-wide_to_K-transect.ipynb would be all netCDF files of /MAR_v3.14_reduced_*/).

While this notebook can be applied to any other Greenland-wide dataset, note that (i) the IMAU-FDM data in /RACMO-FDM_K-Transect_updated/ were obtained directly from IMAU (due to the imense data volume of Greenland-wide IMAU-FDM output for all firn layers) and (ii) that the RACMO grid cells along the K-Transect are extracted on the fly in RCM_analysis_comp_Ktransect.py, also using the Bresenham algorithm.

Input to _notebooks\K-transect_RACMO-MAR_depth-data-plot.ipynb

_depth_of_ice_slab_KAN_U.xlsx

_KAN_U_10m_firn_temperatures.xlsx

/RACMO-FDM_K-Transect_updated/*

/MARv3.14_K-transect_continuous/*

Input to _notebooks\scatter_plots_production.pynb

table_influence_difference_maxysrcm_incl-albedo_v7.4.xlsx which is being written to output by Greenland_RCM_analysis\RCM_analysis_comp_Ktransect_tables.py

Output data

The basic output files of the runoff limit detection for all of Greenland are:

_GR_slush-limit_output_table.xlsx (written by MODIS_Greenland/MODIS_find_slush_limit.py)

__GR_slush-limit_output_table_OnlyValidEntries.xlsx (by MODIS_Greenland_analysis/MODIS_A_SLmax.py)

__GR_table_simple_annual_max_SL.xlsx  (by MODIS_Greenland_analysis/MODIS_A_SLmax.py)

__GR_table_complete_annual_max_SL.xlsx  (by MODIS_Greenland_analysis/MODIS_A_SLmax.py)

The basic output files of the runoff limit detection for the “test” data set are:

_test_slush-limit_output_table.xlsx (written by MODIS_Greenland/MODIS_find_slush_limit.py)

__test_slush-limit_output_table_OnlyValidEntries.xlsx (by MODIS_Greenland_analysis/MODIS_A_SLmax.py)

__test_table_simple_annual_max_SL.xlsx  (by MODIS_Greenland_analysis/MODIS_A_SLmax.py)

__test_table_complete_annual_max_SL.xlsx  (by MODIS_Greenland_analysis/MODIS_A_SLmax.py)