{ "cells": [ { "cell_type": "markdown", "id": "9a981ba7-9fbd-494d-a2cf-25d0259c9473", "metadata": {}, "source": [ "# Figure 7\n", "\n", "This notebook produces Figure 7 and other auxiliary figures.\n", "\n", "## Get data\n", "\n", "(1) from the Ice Thickness Models Intercomparison eXperiment (ITMIX) data set (Farinotti et al., 2017), available at https://doi.org/10.5905/ethz-1007-92. Redistributed under the CC BY-NC-SA 4.0 License (see the previous link for details).\n", "\n", "Reference: Farinotti, D., Brinkerhoff, D. J., Clarke, G. K. C., Fürst, J. J., Frey, H., Gantayat, P., Gillet-Chaulet, F., Girard, C., Huss, M., Leclercq, P. W., Linsbauer, A., Machguth, H., Martin, C., Maussion, F., Morlighem, M., Mosbeux, C., Pandit, A., Portmann, A., Rabatel, A., …Andreassen, L. M. (2017). How accurate are estimates of glacier ice thickness? Results from ITMIX, the Ice Thickness Models Intercomparison eXperiment. Cryosphere, 11(2), 949–970. https://doi.org/10.5194/tc-11-949-2017" ] }, { "cell_type": "code", "execution_count": 1, "id": "7a91ca36-c996-4547-b6b2-d82a60888518", "metadata": {}, "outputs": [], "source": [ "%%bash\n", "\n", "if [ ! -d ../data/ITMIX/ ]; then\n", " wget -P ../data/ITMIX/ --content-disposition --no-verbose 'http://www.geo.cornell.edu/eas/gstudent/wz278/paperdraft2021/ITMIX_dataset.zip' \n", " unzip ../data/ITMIX/ITMIX_dataset.zip -d ../data/ITMIX/\n", " unzip ../data/ITMIX/01_ITMIX_input_data.zip -d ../data/ITMIX/01_ITMIX_input_data/\n", " unzip ../data/ITMIX/02_ITMIX_results.zip -d ../data/ITMIX/02_ITMIX_results/\n", " unzip ../data/ITMIX/03_ITMIX_measured_thickness.zip -d ../data/ITMIX/03_ITMIX_measured_thickness/\n", " rm -rf ../data/ITMIX/*.zip\n", "fi" ] }, { "cell_type": "markdown", "id": "a7467f7c-2d64-4ad1-a865-02f586473eb3", "metadata": {}, "source": [ "(2) From ITS_LIVE: A NASA MEaSUREs project to provide automated, low latency, global glacier flow and elevation change datasets. Original data available at https://its-live.jpl.nasa.gov/. We use the 2018 speed mosaics.\n", "\n", "References: \n", "\n", "1. Gardner, A. S., M. A. Fahnestock, and T. A. Scambos (2019) [Accessed October 21, 2021]. ITS_LIVE Regional Glacier and Ice Sheet Surface Velocities. Data archived at National Snow and Ice Data Center. https://doi.org/10.5067/6II6VW8LLWJ7\n", "2. Gardner, A. S., G. Moholdt, T. Scambos, M. Fahnstock, S. Ligtenberg, M. van den Broeke, and J. Nilsson (2018). Increased West Antarctic and unchanged East Antarctic ice discharge over the last 7 years, Cryosphere, 12(2), 521–547. https://doi.org/10.5194/tc-12-521-2018" ] }, { "cell_type": "code", "execution_count": 2, "id": "7f811b8b-16ec-467a-9d40-17fda1bb7bda", "metadata": {}, "outputs": [], "source": [ "%%bash\n", "\n", "if [ ! -f ../data/SRA_G0240_2018_v_EPSG32633.tif ]; then\n", " wget http://www.geo.cornell.edu/eas/gstudent/wz278/paperdraft2021/SRA_G0240_2018_v_EPSG32633.tif -P ../data/ --no-verbose\n", "fi" ] }, { "cell_type": "markdown", "id": "8d871ef2-4275-428a-a7c7-13d7baa277ad", "metadata": {}, "source": [ "(3) The flowline shapefile available at `../data/flowline_statpts4_clean_EPSG32633.shp` for flowline vertices, which is included in the data folder.\n", "\n", "## Analysis" ] }, { "cell_type": "code", "execution_count": 3, "id": "7bc13070-52ea-4007-a569-c2f083297a81", "metadata": {}, "outputs": [], "source": [ "import geopandas as gpd\n", "import rasterio as rio\n", "import numpy as np\n", "import pejzero\n", "import matplotlib.pyplot as plt\n", "from matplotlib.colors import LinearSegmentedColormap" ] }, { "cell_type": "markdown", "id": "20350657-a8c9-4537-8cce-0d7b44ba2fe2", "metadata": {}, "source": [ "Here are all file locations:" ] }, { "cell_type": "code", "execution_count": 4, "id": "43bb9963-1b85-46db-9ccc-ca94093a9916", "metadata": {}, "outputs": [], "source": [ "flowline_file = '../data/flowline_statpts4_clean_EPSG32633.shp'\n", "bed_file = '../data/ITMIX/02_ITMIX_results/Farinotti_Austfonna_bedrock.asc'\n", "surface_file = '../data/ITMIX/01_ITMIX_input_data/Austfonna/02_surface_Austfonna_2007_UTM33.asc'\n", "speed_file = '../data/ITMIX/01_ITMIX_input_data/Austfonna/06_speed_Austfonna_1995-1996_UTM33.asc'\n", "speed2018_file = '../data/SRA_G0240_2018_v_EPSG32633.tif'" ] }, { "cell_type": "markdown", "id": "ed521616-43e4-46db-a126-838f041ac2c1", "metadata": {}, "source": [ "### (1) Data ingest\n", "\n", "Sample all the necessary measurements at the vertex locations:" ] }, { "cell_type": "code", "execution_count": 5, "id": "34712b68-54c7-4a9e-ae3c-c149402def59", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
| \n", " | glacier_id | \n", "geometry | \n", "d | \n", "s | \n", "b | \n", "u | \n", "udiff | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "11 | \n", "LINESTRING (661744.501 8812864.548, 661772.209... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[124.80580139160156, 125.9896011352539, 126.51... | \n", "[6.389999866485596, 6.630000114440918, 6.76000... | \n", "[9.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0... | \n", "[2.063974380493164, 1.063974380493164, 1.06397... | \n", "
| 1 | \n", "12 | \n", "LINESTRING (664882.963 8810146.886, 664950.156... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[105.45079803466797, 107.12020111083984, 107.8... | \n", "[-66.0, -75.98999786376953, -78.91999816894531... | \n", "[12.0, 13.0, 13.0, 13.0, 13.0, 14.0, 14.0, 14.... | \n", "[0.3671684265136719, 0.06629180908203125, 0.06... | \n", "
| 2 | \n", "13 | \n", "LINESTRING (669704.623 8810585.218, 669770.209... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[130.61459350585938, 132.28590393066406, 132.4... | \n", "[-107.91999816894531, -113.83999633789062, -11... | \n", "[18.0, 17.0, 17.0, 17.0, 17.0, 17.0, 17.0, 17.... | \n", "[-6.172802925109863, -5.395359039306641, -5.77... | \n", "
| 3 | \n", "14 | \n", "LINESTRING (673035.951 8808972.154, 673095.381... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[127.32119750976562, 127.9738998413086, 128.56... | \n", "[-41.369998931884766, -45.77000045776367, -47.... | \n", "[14.0, 14.0, 14.0, 14.0, 14.0, 14.0, 14.0, 14.... | \n", "[-0.8997259140014648, -0.8997259140014648, -0.... | \n", "
| 4 | \n", "15 | \n", "LINESTRING (676753.012 8808446.155, 676815.893... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[144.21249389648438, 146.13330078125, 147.1508... | \n", "[-46.41999816894531, -50.880001068115234, -53.... | \n", "[9.0, 8.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, ... | \n", "[5.6007080078125, 6.6007080078125, 7.600708007... | \n", "
| 5 | \n", "16 | \n", "LINESTRING (680680.473 8808481.221, 680771.006... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[142.48159790039062, 142.54119873046875, 142.3... | \n", "[-84.48999786376953, -88.70999908447266, -89.3... | \n", "[8.0, 8.0, 8.0, 8.0, 8.0, 8.0, 8.0, 8.0, 8.0, ... | \n", "[27.496936798095703, 27.496936798095703, 27.49... | \n", "
| 6 | \n", "31 | \n", "LINESTRING (709733.163 8827112.114, 709698.322... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "
| 7 | \n", "32 | \n", "LINESTRING (712538.492 8829608.857, 712479.960... | \n", "[0.0, 50.0, 100.0, 150.0, 200.0, 250.0, 300.0,... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "[nan, nan, nan, nan, nan, nan, nan, nan, nan, ... | \n", "