{ "cells": [ { "cell_type": "markdown", "id": "d57527fa-3460-4b16-9df9-c37f72bf4fb1", "metadata": {}, "source": [ "# Runoff comparison Indus:" ] }, { "cell_type": "markdown", "id": "696ab51b-ae30-4cf3-9796-97cbde8e82ce", "metadata": {}, "source": [ "Here we compare the runoff estimated from Rounce et al., (2023) to the runoff estimated from OGGM (simulations done by Erik Holmgren using OGGM v. XXX), and to the runoff estimated from OGGM v1.6.\n", "\n", "--> we compare different methods for single glaciers and a single gcm here: www_lschuster/provide/runoff_example_data/runoff_method_comparison_oggm_pygem.ipynb" ] }, { "cell_type": "code", "execution_count": 1, "id": "557132b4-2b7f-4acd-a524-ee910eb7a37c", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import seaborn as sns" ] }, { "cell_type": "code", "execution_count": 2, "id": "861fcc42-6f94-48f3-b8e9-60434554a13f", "metadata": {}, "outputs": [], "source": [ "basin= 'INDUS'" ] }, { "cell_type": "code", "execution_count": 3, "id": "fe1cea5b-d005-467d-8b11-dfe899dbace0", "metadata": {}, "outputs": [], "source": [ "\n", "import geopandas as gpd\n", "pd_basin_num = gpd.read_file('/home/www/fmaussion/misc/magicc/basins_shape/glacier_basins.shp')\n", "basin_idx = pd_basin_num[pd_basin_num['RIVER_BASI'] == basin]['MRBID'].values[0]\n", "import json\n", "f = open('/home/www/fmaussion/misc/magicc/rgi_ids_per_basin.json')\n", "rgis_basin = json.load(f)[str(basin_idx)]\n", "# they are all either from RGI region 13 or 14\n", "\n", "path_oggm_erik = f'/home/www/eholmgren/data/basins/gha_basins/{basin_idx}/'\n", "gcms = ['BCC-CSM2-MR', 'CESM2', 'CESM2-WACCM', 'EC-Earth3', 'EC-Earth3-Veg',\n", " 'FGOALS-f3-L', 'GFDL-ESM4', 'INM-CM4-8', 'INM-CM5-0', 'MPI-ESM1-2-HR',\n", " 'MRI-ESM2-0', 'NorESM2-MM']" ] }, { "cell_type": "markdown", "id": "e9a744ec-5f0b-46f9-b78f-73055744b7c5", "metadata": {}, "source": [ "### Extract information of terminus location of the glaciers " ] }, { "cell_type": "code", "execution_count": 5, "id": "0fc7a554-3754-4403-a909-7436983926e4", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/tmp/ipykernel_3985733/666086889.py:3: DtypeWarning: Columns (16,17) have mixed types. Specify dtype option on import or set low_memory=False.\n", " pd_stats_13 = pd.read_csv(path_stats_13)\n" ] }, { "data": { "text/html": [ "
| \n", " | terminus_lon | \n", "terminus_lat | \n", "
|---|---|---|
| rgi_id | \n", "\n", " | \n", " |
| RGI60-13.00002 | \n", "77.951262 | \n", "35.541153 | \n", "
| RGI60-13.00188 | \n", "78.662054 | \n", "34.858492 | \n", "
| RGI60-13.00224 | \n", "79.169854 | \n", "34.436261 | \n", "
| RGI60-13.00228 | \n", "79.235076 | \n", "34.422087 | \n", "
| RGI60-13.00229 | \n", "79.407439 | \n", "34.419129 | \n", "
| ... | \n", "... | \n", "... | \n", "
| RGI60-14.27984 | \n", "78.974081 | \n", "33.337359 | \n", "
| RGI60-14.27985 | \n", "79.036696 | \n", "33.293083 | \n", "
| RGI60-14.27986 | \n", "79.038796 | \n", "33.287543 | \n", "
| RGI60-14.27987 | \n", "79.098103 | \n", "33.212602 | \n", "
| RGI60-14.27988 | \n", "79.258021 | \n", "33.217487 | \n", "
23424 rows × 2 columns
\n", "<xarray.DataArray 'glac_runoff_fixed_monthly' (model: 12, glacier: 1000,\n",
" time: 1212)>\n",
"[14544000 values with dtype=float64]\n",
"Coordinates:\n",
" RGIId (glacier) object 'RGI60-01.25001' ... 'RGI60-01.26000'\n",
" Climate_Model (model) object 'BCC-CSM2-MR' 'CESM2' ... 'NorESM2-MM'\n",
" lon (glacier) float64 -148.6 -148.4 -148.6 ... -138.1 -138.0\n",
" lat (glacier) float64 61.06 61.0 60.93 ... 59.27 59.31 59.33\n",
" * time (time) datetime64[ns] 2000-01-01 2000-02-01 ... 2100-12-01\n",
"Dimensions without coordinates: model, glacier\n",
"Attributes:\n",
" long_name: Glacier-wide runoff\n",
" unit: m3\n",
" temporal_resolution: monthly\n",
" comment: runoff from a fixed-gauge at glacier terminus that ...\n",
" grid_mapping: crs<xarray.DataArray 'area' (gcm: 12, rgi_id: 23424, year: 102)>\n",
"[28670976 values with dtype=float64]\n",
"Coordinates:\n",
" * year (year) int64 2000 2001 2002 2003 2004 ... 2097 2098 2099 2100 2101\n",
" * rgi_id (rgi_id) object 'RGI60-14.26870' ... 'RGI60-14.16511'\n",
" * gcm (gcm) object 'BCC-CSM2-MR' 'CESM2' ... 'MRI-ESM2-0' 'NorESM2-MM'\n",
" ssp <U6 'ssp126'\n",
"Attributes:\n",
" long_name: Glacier area\n",
" temporal_resolution: annual\n",
" units: m2\n",
" comment: area at start of the year\n",
" grid_mapping: crs\n",
" description:: glacier area evolution\n",
" unit:: m2\n",
" data source:: Rounce, D.R., R. Hock, and F. Maussion. 2022. Globa...\n",
" coordinates: lon lat Climate_Model RGIId<xarray.DataArray 'area' (gcm: 12, rgi_id: 1209)>\n",
"array([[ 862655.67642581, 11649.61450761, 43873.26024531, ...,\n",
" 3534114.15883217, 237472.98506712, 16855.62211932],\n",
" [1261587.68789969, 11531.82940832, 43498.57830669, ...,\n",
" 2031131.61805209, 257678.3889862 , 20311.29401576],\n",
" [1366226.00742894, 11768.24818609, 43318.10858468, ...,\n",
" 2035933.29573755, 277005.77875329, 20311.29401576],\n",
" ...,\n",
" [1358318.38756374, 11738.061399 , 67742.23828215, ...,\n",
" 6903747.58379304, 512115.07236148, 20311.29401576],\n",
" [2356095.56813629, 11485.35955636, 42604.41865651, ...,\n",
" 2937876.70362377, 790010.16224168, 20311.29401576],\n",
" [1704805.8030851 , 11732.2258039 , 66985.52942439, ...,\n",
" 3453157.22099302, 429550.01053918, 20311.29401576]])\n",
"Coordinates:\n",
" lon (rgi_id) float64 74.11 72.19 77.9 77.06 ... 75.26 76.11 74.26 76.19\n",
" lat (rgi_id) float64 36.38 35.77 31.7 32.73 ... 36.88 33.44 36.89 35.44\n",
" year int64 2101\n",
" * rgi_id (rgi_id) object 'RGI60-14.03563' ... 'RGI60-14.07503'\n",
" * gcm (gcm) object 'BCC-CSM2-MR' 'CESM2' ... 'MRI-ESM2-0' 'NorESM2-MM'\n",
" ssp <U6 'ssp585'\n",
"Attributes:\n",
" long_name: Glacier area\n",
" temporal_resolution: annual\n",
" units: m2\n",
" comment: area at start of the year\n",
" grid_mapping: crs\n",
" description:: glacier area evolution\n",
" unit:: m2\n",
" data source:: Rounce, D.R., R. Hock, and F. Maussion. 2022. Globa...<xarray.DataArray 'area' (time: 81)>\n",
"array([34357.95 , 34280.99 , 34215.387 , 34162.484 , 34178.492 ,\n",
" 34082.918 , 34020.58 , 34058.508 , 33891.47 , 33839.02 ,\n",
" 33715.01 , 33744.797 , 33618.258 , 33482.81 , 33372.523 ,\n",
" 33286.41 , 33145.11 , 33070.875 , 33215.348 , 32869.895 ,\n",
" 32705.7 , 32592.934 , 32383.156 , 32157.914 , 31994.863 ,\n",
" 31765.959 , 31908.54 , 31379.496 , 31041.059 , 30854.506 ,\n",
" 30638.01 , 30371.191 , 30047.613 , 29675.621 , 29315.379 ,\n",
" 29019.947 , 28447.281 , 28031.45 , 27606.912 , 26852.445 ,\n",
" 26377.83 , 25822.45 , 25129.336 , 24372.604 , 23721.61 ,\n",
" 23257.338 , 22656.05 , 21781.678 , 21159.922 , 20467.223 ,\n",
" 19974.305 , 19342.895 , 18817.152 , 17679.863 , 17270.387 ,\n",
" 16431.309 , 15825.784 , 15097.484 , 14615.607 , 14350.192 ,\n",
" 13731.659 , 13225.485 , 12927.309 , 12267.748 , 11778.892 ,\n",
" 11384.041 , 11758.571 , 11026.897 , 10290.394 , 9950.551 ,\n",
" 9603.716 , 9435.831 , 9043.37 , 8772.6455, 8393.718 ,\n",
" 8294.797 , 7961.6606, 8083.894 , 7520.426 , 7161.3066,\n",
" 7035.2793], dtype=float32)\n",
"Coordinates:\n",
" * time (time) float64 2.02e+03 2.021e+03 ... 2.099e+03 2.1e+03\n",
" hydro_year (time) int64 2020 2021 2022 2023 ... 2097 2098 2099 2100\n",
" hydro_month (time) int64 1 1 1 1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1 1 1 1 1\n",
" calendar_year (time) int64 2019 2020 2021 2022 ... 2096 2097 2098 2099\n",
" calendar_month (time) int64 10 10 10 10 10 10 10 ... 10 10 10 10 10 10 10<xarray.Dataset>\n",
"Dimensions: (gcm: 12, ssp: 4, hydro_yr: 81, calendar_month_2d: 12,\n",
" rgi_id: 23424)\n",
"Coordinates:\n",
" * gcm (gcm) <U13 'BCC-CSM2-MR' 'CESM2' ... 'NorESM2-MM'\n",
" * rgi_id (rgi_id) object 'RGI60-13.00002' ... 'RGI60-14.27988'\n",
" * calendar_month_2d (calendar_month_2d) int64 10 11 12 1 2 3 4 5 6 7 8 9\n",
" * hydro_yr (hydro_yr) float64 2.02e+03 2.021e+03 ... 2.1e+03\n",
" * ssp (ssp) <U6 'ssp126' 'ssp245' 'ssp370' 'ssp585'\n",
"Data variables:\n",
" runoff (ssp, gcm, hydro_yr, calendar_month_2d, rgi_id) float32 ...\n",
" calendar_year (ssp, gcm, hydro_yr) int64 2019 2020 2021 ... 2098 2099<xarray.DataArray 'hydro_yr' (hydro_yr: 81)>\n",
"array([2020., 2021., 2022., 2023., 2024., 2025., 2026., 2027., 2028., 2029.,\n",
" 2030., 2031., 2032., 2033., 2034., 2035., 2036., 2037., 2038., 2039.,\n",
" 2040., 2041., 2042., 2043., 2044., 2045., 2046., 2047., 2048., 2049.,\n",
" 2050., 2051., 2052., 2053., 2054., 2055., 2056., 2057., 2058., 2059.,\n",
" 2060., 2061., 2062., 2063., 2064., 2065., 2066., 2067., 2068., 2069.,\n",
" 2070., 2071., 2072., 2073., 2074., 2075., 2076., 2077., 2078., 2079.,\n",
" 2080., 2081., 2082., 2083., 2084., 2085., 2086., 2087., 2088., 2089.,\n",
" 2090., 2091., 2092., 2093., 2094., 2095., 2096., 2097., 2098., 2099.,\n",
" 2100.])\n",
"Coordinates:\n",
" * hydro_yr (hydro_yr) float64 2.02e+03 2.021e+03 ... 2.099e+03 2.1e+03"
],
"text/plain": [
"<xarray.DataArray 'glac_runoff_fixed_monthly' (model: 12, glacier: 1000,\n",
" year: 101)>\n",
"array([[[ 421163.32271247, 428350.40173807, 236180.93655957, ...,\n",
" 164672.80833187, 210343.38690611, 204119.75813437],\n",
" [ 419109.9438986 , 401803.846151 , 227305.2308056 , ...,\n",
" 225963.495927 , 288746.00352455, 280003.57513188],\n",
" [1200690.21356446, 1207128.29095253, 790434.24852661, ...,\n",
" 473402.22337652, 549029.18131243, 541426.94817699],\n",
" ...,\n",
" [ 121973.20055788, 120205.55944412, 81002.05416204, ...,\n",
" 55278.35690715, 70767.02924289, 68353.39559877],\n",
" [ 82400.64872411, 84388.69332114, 38323.83851968, ...,\n",
" 32573.08283819, 40373.74997177, 39766.66193446],\n",
" [ 67066.16812958, 67548.62282137, 36004.93468778, ...,\n",
" 26318.50433423, 39312.09686865, 31933.47962748]],\n",
"\n",
" [[ 213474.1875566 , 214997.28457408, 231778.00170529, ...,\n",
" 265009.59288821, 260529.3312875 , 215212.1645575 ],\n",
" [ 193730.087781 , 214297.80972054, 209145.56561322, ...,\n",
" 361543.00789529, 358550.187143 , 292168.01165931],\n",
" [ 725963.51284073, 674283.88542866, 799749.54182649, ...,\n",
" 664774.87261924, 651967.81120448, 566542.11503985],\n",
"...\n",
" [ 78493.35010314, 89572.02531429, 90249.46571689, ...,\n",
" 81771.61663878, 79838.20242421, 104513.47581776],\n",
" [ 52961.3326751 , 47532.8882202 , 51224.48943128, ...,\n",
" 48771.53509895, 48084.24070128, 62045.04366603],\n",
" [ 41497.63177282, 38071.0126381 , 42226.37690294, ...,\n",
" 37999.50941189, 37055.06549567, 48234.05358955]],\n",
"\n",
" [[ 341047.03013004, 224381.56446911, 354442.55627237, ...,\n",
" 219126.70681384, 205052.48392779, 240545.31214457],\n",
" [ 324855.09983937, 207129.25432627, 338909.65640149, ...,\n",
" 306119.02967775, 286513.31353742, 334019.27335634],\n",
" [1013776.03184228, 749739.00684619, 1036305.16864216, ...,\n",
" 574357.88998442, 547773.83981725, 605663.16537171],\n",
" ...,\n",
" [ 104880.11344552, 75876.54268797, 104281.58117388, ...,\n",
" 77793.34227983, 72341.02392792, 75983.74595444],\n",
" [ 61314.32660297, 35710.38808346, 69694.06590394, ...,\n",
" 44416.68380562, 41258.23942087, 43469.92784485],\n",
" [ 52424.6549027 , 33725.50944504, 54220.19118793, ...,\n",
" 37468.24987823, 34848.12606489, 36624.65171225]]])\n",
"Coordinates:\n",
" RGIId (glacier) object 'RGI60-13.08001' ... 'RGI60-13.09000'\n",
" Climate_Model (model) object 'BCC-CSM2-MR' 'CESM2' ... 'NorESM2-MM'\n",
" lon (glacier) float64 77.53 77.52 77.49 ... 77.08 77.07 77.06\n",
" lat (glacier) float64 41.81 41.81 41.92 ... 41.84 41.85 41.84\n",
" * year (year) int64 2000 2001 2002 2003 2004 ... 2097 2098 2099 2100\n",
"Dimensions without coordinates: model, glacier<xarray.DataArray 'glac_runoff_fixed_monthly' (model: 12, glacier: 1000,\n",
" time: 1212)>\n",
"array([[[0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.],\n",
" ...,\n",
" [0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.]],\n",
"\n",
" [[0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.],\n",
" ...,\n",
" [0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.]],\n",
"\n",
" ...,\n",
"\n",
" [[0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.],\n",
" ...,\n",
" [0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.]],\n",
"\n",
" [[0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.],\n",
" ...,\n",
" [0., 0., ..., 0., 0.],\n",
" [0., 0., ..., 0., 0.]]])\n",
"Coordinates:\n",
" RGIId (glacier) object 'RGI60-13.08001' ... 'RGI60-13.09000'\n",
" Climate_Model (model) object 'BCC-CSM2-MR' 'CESM2' ... 'NorESM2-MM'\n",
" lon (glacier) float64 77.53 77.52 77.49 ... 77.08 77.07 77.06\n",
" lat (glacier) float64 41.81 41.81 41.92 ... 41.84 41.85 41.84\n",
" * time (time) datetime64[ns] 2000-01-01 2000-02-01 ... 2100-12-01\n",
"Dimensions without coordinates: model, glacier\n",
"Attributes:\n",
" long_name: Glacier-wide runoff\n",
" unit: m3\n",
" temporal_resolution: monthly\n",
" comment: runoff from a fixed-gauge at glacier terminus that ...\n",
" grid_mapping: crs<xarray.DataArray 'runoff' (time: 100, rgi_id: 23390)>\n",
"array([[8.3875504e+07, 7.9531110e+08, 2.0042262e+08, ..., 2.8898712e+07,\n",
" 6.0809712e+07, 3.8973804e+07],\n",
" [9.6166120e+07, 8.8865837e+08, 1.7704848e+08, ..., 2.8444002e+07,\n",
" 6.9275264e+07, 3.6920380e+07],\n",
" [1.8680954e+08, 9.2112640e+08, 2.2733728e+08, ..., 3.5258640e+07,\n",
" 8.4173264e+07, 4.8226992e+07],\n",
" ...,\n",
" [3.2251334e+08, 1.2122401e+09, 2.4285963e+08, ..., 2.3531648e+07,\n",
" 6.1522008e+07, 2.8962562e+07],\n",
" [3.3573936e+08, 2.0365238e+09, 4.3380963e+08, ..., 6.3017200e+07,\n",
" 1.6475448e+08, 7.7561056e+07],\n",
" [2.2699155e+08, 1.0146522e+09, 2.0557677e+08, ..., 2.1397358e+07,\n",
" 5.5942040e+07, 2.6335692e+07]], dtype=float32)\n",
"Coordinates:\n",
" * time (time) float64 2e+03 2.001e+03 ... 2.098e+03 2.099e+03\n",
" * rgi_id (rgi_id) object 'RGI60-13.00002' ... 'RGI60-14.27988'\n",
" hydro_year (time) int64 2000 2001 2002 2003 ... 2096 2097 2098 2099\n",
" hydro_month (time) int64 4 4 4 4 4 4 4 4 4 4 4 ... 4 4 4 4 4 4 4 4 4 4 4\n",
" calendar_year (time) int64 2000 2001 2002 2003 ... 2096 2097 2098 2099\n",
" calendar_month (time) int64 1 1 1 1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1 1 1 1 1<xarray.DataArray 'glac_runoff_fixed_monthly' (gcm: 12, year: 100)>\n",
"array([[ nan, nan, 99.09062826, ..., 72.96119266,\n",
" 68.96954099, 67.13733133],\n",
" [ nan, nan, 66.94200875, ..., 73.40398869,\n",
" 63.21606232, 63.04132331],\n",
" [ nan, nan, 93.13540184, ..., 76.23517047,\n",
" 80.52481986, 74.00299028],\n",
" ...,\n",
" [ nan, nan, 82.23918759, ..., 79.42487775,\n",
" 74.59565253, 78.97827646],\n",
" [ nan, nan, 97.7674347 , ..., 56.78967033,\n",
" 61.22004234, 60.59642815],\n",
" [ nan, nan, 96.74293681, ..., 80.33039862,\n",
" 76.02856694, 70.22836774]])\n",
"Coordinates:\n",
" * gcm (gcm) object 'BCC-CSM2-MR' 'CESM2' ... 'MRI-ESM2-0' 'NorESM2-MM'\n",
" ssp <U6 'ssp126'\n",
" * year (year) int64 2000 2001 2002 2003 2004 ... 2095 2096 2097 2098 2099<xarray.DataArray 'glac_runoff_fixed_monthly' (gcm: 12, rgi_id: 23422,\n",
" time: 1212)>\n",
"array([[[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
"...\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]]])\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 2000-01-01 2000-02-01 ... 2100-12-01\n",
" * rgi_id (rgi_id) object 'RGI60-13.52348' ... 'RGI60-14.00116'\n",
" * gcm (gcm) object 'BCC-CSM2-MR' 'CESM2' ... 'MRI-ESM2-0' 'NorESM2-MM'\n",
" ssp <U6 'ssp126'<xarray.DataArray 'glac_runoff_fixed_monthly' (gcm: 12, rgi_id: 0, time: 1212)>\n",
"array([], shape=(12, 0, 1212), dtype=float64)\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 2000-01-01 2000-02-01 ... 2100-12-01\n",
" * rgi_id (rgi_id) object \n",
" * gcm (gcm) object 'BCC-CSM2-MR' 'CESM2' ... 'MRI-ESM2-0' 'NorESM2-MM'\n",
" ssp <U6 'ssp585'<xarray.DataArray 'glac_runoff_fixed_monthly' (gcm: 12, rgi_id: 23422,\n",
" time: 1212)>\n",
"array([[[ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" ...,\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan]],\n",
"\n",
" [[ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
"...\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan]],\n",
"\n",
" [[ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" ...,\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan],\n",
" [ nan, nan, nan, ...,\n",
" nan, nan, nan]]])\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 2000-01-01 2000-02-01 ... 2100-12-01\n",
" * rgi_id (rgi_id) object 'RGI60-13.52348' ... 'RGI60-14.00116'\n",
" * gcm (gcm) object 'BCC-CSM2-MR' 'CESM2' ... 'MRI-ESM2-0' 'NorESM2-MM'\n",
" ssp <U6 'ssp585'"
],
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (time: 972)\n",
"Coordinates:\n",
" * time (time) object 2020-01-16 12:00:00 ... 2100-12-1...\n",
" month (time) int64 1 2 3 4 5 6 7 8 ... 6 7 8 9 10 11 12\n",
" year (time) int64 2020 2020 2020 ... 2100 2100 2100\n",
"Data variables:\n",
" prcp (time) float64 50.99 148.0 46.39 ... 1.931 80.36\n",
" PET (time) float64 0.3812 0.2695 ... 5.302 0.3345\n",
" glacier_runoff (time) float32 0.0002651 0.008911 ... 0.0 0.0\n",
" glacier_runoff_nofac (time) float32 0.0002613 0.006933 ... 0.0 0.0\n",
" glacier_runoff_adj (time) float64 8.064e-06 0.000271 ... 0.0 0.0\n",
" glacier_runoff_adj_nofac (time) float64 7.947e-06 0.0002109 ... 0.0 0.0\n",
" prcp_adj (time) float64 49.44 143.5 44.98 ... 1.872 77.92\n",
" D (time) float64 50.61 147.7 43.33 ... -3.371 80.03\n",
" D_adj (time) float64 49.06 143.2 41.93 ... -3.43 77.58\n",
" D_adj_nofac (time) float64 49.06 143.2 41.93 ... -3.43 77.58\n",
"Attributes:\n",
" gcm: BCC-CSM2-MR_ssp126\n",
" max glaciated area: 26291694000.0<xarray.Dataset>\n",
"Dimensions: (time: 81, rgi_id: 23424, month_2d: 12)\n",
"Coordinates:\n",
" * time (time) float64 2.02e+03 2.021e+03 ... 2.1e+03\n",
" * rgi_id (rgi_id) object 'RGI60-13.00002' ... 'RGI60...\n",
" hydro_year (time) int64 2020 2021 2022 ... 2098 2099 2100\n",
" hydro_month (time) int64 1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1\n",
" calendar_year (time) int64 2019 2020 2021 ... 2097 2098 2099\n",
" calendar_month (time) int64 10 10 10 10 10 ... 10 10 10 10 10\n",
" * month_2d (month_2d) int64 1 2 3 4 5 6 7 8 9 10 11 12\n",
" calendar_month_2d (month_2d) int64 10 11 12 1 2 3 4 5 6 7 8 9\n",
"Data variables: (12/24)\n",
" volume (time, rgi_id) float32 ...\n",
" volume_bsl (time, rgi_id) float32 ...\n",
" volume_bwl (time, rgi_id) float32 ...\n",
" area (time, rgi_id) float32 ...\n",
" length (time, rgi_id) float32 ...\n",
" calving (time, rgi_id) float32 ...\n",
" ... ...\n",
" liq_prcp_on_glacier_monthly (time, month_2d, rgi_id) float32 ...\n",
" snowfall_off_glacier_monthly (time, month_2d, rgi_id) float32 ...\n",
" snowfall_on_glacier_monthly (time, month_2d, rgi_id) float32 ...\n",
" water_level (rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n",
" glen_a (rgi_id) float32 2.789e-24 ... 2.254e-24\n",
" fs (rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n",
"Attributes:\n",
" description: OGGM model output\n",
" oggm_version: 1.5.2\n",
" calendar: 365-day no leap\n",
" creation_date: 2021-10-19 16:57:59<xarray.DataArray 'glacier_runoff' (year: 81)>\n",
"array([1419.7529, 1521.861 , 1539.1827, 1519.1991, 1646.0447, 1357.7666,\n",
" 1529.0192, 1644.2972, 1473.3912, 1645.3276, 1507.8579, 1444.3009,\n",
" 1284.3569, 1450.0841, 1486.779 , 1427.2823, 1504.4945, 1438.1185,\n",
" 1308.4275, 1607.2449, 1337.0426, 1522.8254, 1575.3405, 1374.9585,\n",
" 1537.5654, 1559.8712, 1390.764 , 1674.3567, 1420.5938, 1560.2672,\n",
" 1827.844 , 1372.9092, 1577.5342, 1490.4916, 1620.1671, 1503.2883,\n",
" 1324.624 , 1445.3851, 1609.6399, 1502.6619, 1399.5876, 1422.5944,\n",
" 1346.9691, 1282.4486, 1388.2607, 1435.6373, 1336.7617, 1342.1437,\n",
" 1248.9712, 1505.0847, 1248.3435, 1469.8102, 1114.7633, 1558.7113,\n",
" 1400.9601, 1360.7726, 1312.5115, 1533.4857, 1669.1787, 1173.8094,\n",
" 1386.818 , 1479.4722, 1148.2114, 1259.6659, 1342.435 , 1175.8828,\n",
" 1307.1123, 1237.6127, 1432.2024, 1123.2223, 1376.1476, 1450.1809,\n",
" 1239.796 , 1193.6906, 1146.5719, 1198.6334, 1374.1276, 1270.9683,\n",
" 1302.4999, 1405.1603, 0. ], dtype=float32)\n",
"Coordinates:\n",
" * year (year) int64 2020 2021 2022 2023 2024 ... 2096 2097 2098 2099 2100