📈 Cascade plots of a few subglacial lakes

Not to be confused with waterfall plots, nor anything really. This is still a fairly basic graph showing an inferred cascade of Whillans 7 -> Whillans IX -> Subglacial Lake Whillans -> Whillans 12. The numbers don't quite sum up perfectly on the chart, so I might need to look at a few more adjacent lakes. This should do for now though!

Author: weiji14

atlxi_xover.ipynb

@@ -430,6 +430,60 @@
     "fig.show()"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## Combined ice volume displacement plot\n",
+    "\n",
+    "Showing how subglacial water cascades down a drainage basin!"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "lines_to_next_cell": 0
+   },
+   "outputs": [],
+   "source": [
+    "fig = pygmt.Figure()\n",
+    "fig.basemap(\n",
+    "    region=f\"2019-02-28/2020-09-30/-0.3/0.5\",\n",
+    "    frame=[\"wSnE\", \"xaf\", 'yaf+l\"Ice Volume Displacement (km@+3@+)\"'],\n",
+    ")\n",
+    "pygmt.makecpt(cmap=\"davosS\", color_model=\"+c\", series=(-2, 4, 0.5))\n",
+    "for i, (_placename, linestyle) in enumerate(\n",
+    "    iterable=zip(\n",
+    "        [\"whillans_ix\", \"subglacial_lake_whillans\", \"whillans_12\", \"whillans_7\"],\n",
+    "        [\"\", \".-\", \"-\", \"..-\"],\n",
+    "    )\n",
+    "):\n",
+    "    fig.plot(\n",
+    "        data=f\"figures/{_placename}/ice_dvol_dt_{_placename}.txt\",\n",
+    "        cmap=True,\n",
+    "        pen=f\"thick,{linestyle}\",\n",
+    "        zvalue=i,\n",
+    "        label=_placename,\n",
+    "        columns=\"0,3\",  # time column (0), ice_dvol column (3)\n",
+    "    )\n",
+    "fig.text(\n",
+    "    position=\"TL\",\n",
+    "    offset=\"j0.2c\",\n",
+    "    text=\"Whillans Ice Stream Central Catchment active subglacial lakes\",\n",
+    ")\n",
+    "fig.legend(position=\"jML+jML+o0.2c\", box=\"+gwhite\")\n",
+    "fig.savefig(\"figures/cascade_whillans_ice_stream.png\")\n",
+    "fig.show()"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,

atlxi_xover.py

@@ -113,6 +113,7 @@
 )
 df_lake: pd.DataFrame = df_lake.loc[gdf_lake.within(lake.geometry)]  # polygon subset
 
+
 # %%
 # Run crossover analysis on all tracks
 track_dict: dict = deepicedrain.split_tracks(df=df_lake)
@@ -222,7 +223,6 @@
 fig.savefig(f"figures/{placename}/crossover_area_{placename}_{min_date}_{max_date}.png")
 fig.show()
 
-
 # %% [markdown]
 # ### Plot Crossover Elevation time-series
 #
@@ -320,3 +320,40 @@
 fig.show()
 
 # %%
+
+# %% [markdown]
+# ## Combined ice volume displacement plot
+#
+# Showing how subglacial water cascades down a drainage basin!
+
+# %%
+fig = pygmt.Figure()
+fig.basemap(
+    region=f"2019-02-28/2020-09-30/-0.3/0.5",
+    frame=["wSnE", "xaf", 'yaf+l"Ice Volume Displacement (km@+3@+)"'],
+)
+pygmt.makecpt(cmap="davosS", color_model="+c", series=(-2, 4, 0.5))
+for i, (_placename, linestyle) in enumerate(
+    iterable=zip(
+        ["whillans_ix", "subglacial_lake_whillans", "whillans_12", "whillans_7"],
+        ["", ".-", "-", "..-"],
+    )
+):
+    fig.plot(
+        data=f"figures/{_placename}/ice_dvol_dt_{_placename}.txt",
+        cmap=True,
+        pen=f"thick,{linestyle}",
+        zvalue=i,
+        label=_placename,
+        columns="0,3",  # time column (0), ice_dvol column (3)
+    )
+fig.text(
+    position="TL",
+    offset="j0.2c",
+    text="Whillans Ice Stream Central Catchment active subglacial lakes",
+)
+fig.legend(position="jML+jML+o0.2c", box="+gwhite")
+fig.savefig("figures/cascade_whillans_ice_stream.png")
+fig.show()
+
+# %%

deepicedrain/features/subglacial_lakes.feature

@@ -61,4 +61,6 @@ Feature: Mapping Antarctic subglacial lakes
     Examples:
     | lake_name                | location            |
     | Whillans 7               | whillans_upstream   |
+    # | Whillans IX              | whillans_upstream   |
+    # | Subglacial Lake Whillans | whillans_downstream |
     | Whillans 12              | whillans_downstream |