Moulins across the Greenland ice sheet

This repository provides the datasets and codes corresponding to the manuscript "Satellite mapping of the moulin distribution across the entire Greenland Ice Sheet". Detailed discriptions are provided below. 01_Moulin_Map_across_the_Greenland_Ice_Sheet.shp: Moulins (and terminal lakes) across the GrIS mapped from Sentinel-2 imagery in the summer of 2019.02_GrIS_Basin_with_Attributes.shp: The six basins of GrIS with attributes used in the analysis, including daily runoff, runoff area and ice velocity.03_GrIS_Basin_with_Bin_Area.shp: The areas of elevation bins in the six basins of GrIS.04_Supraglacial_lakes_across_the_GrIS_Dunmire2021.shp: Convex hulls of supraglacial lakes provided by Dunmire et al. (2021).05_Moulin_Samples.shp: Moulin samples in statistically significant hotspot and coldspot regions.06_Non-Moulin_Samples.shp: Non-moulin samples across the GrIS.07_Moulin_Statistics_across_Zwally_Basins.shp: Moulin numbers and densities, as well as ice sheet attributes, analyzed across the GrIS drainage systems delineated by Zwally et al. (2012). Code_1_Analysis of GrIS-wide moulin distribution.ipynb: A Jupyter Notebook to analyze moulin distribution across the GrIS, reproducing tables and figures presented in the manuscript with the provided data.Code_2_Analysis of potential controls on moulin distribution.ipynb: A Jupyter Notebook to analyze potential controls on moulin distribution, reproducing tables and figures presented in the manuscript with the provided data. Data Preparation: This folder includes data and codes supporting moulin identification. calculate NDWI from Sentinel-2 imagery.py: A Python script to calculate the normalized difference water index (NDWI) from Sentinel-2 satellite imagery. calculate Drainage Networks from ArcticDEM.py: A Python script to model supraglacial drainage networks based on ArcticDEM. calculate Strain Rate from XY ice velocity.py: A Python script to calculate strain rates from ice velocity components in the x and y directions. GrIS_Flow_regime_4.tif: Flow regime 4 used to mask regions where moulin identification was not conducted.