Skip to content

Commit f716373

Browse files
weiji14weiji14
committed
📝 Refresh Whillans IX figures on main README.md up to 20201224
Update the main README.md file with new figures over subglacial lake Whillans IX up to 20201224. Looks like the ice surface over the active subglacial lake has returned back to the same elevation pre-filling (and may have sank down a little bit more actually)! Also made some tiny updates to the atlxi_lake and atlxi_xover scripts
1 parent 38949ad commit f716373

File tree

6 files changed

+35
-36
lines changed

6 files changed

+35
-36
lines changed

‎README.md‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -10,16 +10,16 @@ in Antarctica using remote sensing and machine learning.
1010
[![Dependabot Status](https://api.dependabot.com/badges/status?host=github&repo=weiji14/deepicedrain)](https://dependabot.com)
1111
![License](https://img.shields.io/github/license/weiji14/deepicedrain)
1212

13-
| Ice Surface Elevation trends over Antactica | Active Subglacial Lake filling event |
13+
| Ice Surface Elevation trends over Antactica | Active Subglacial Lake fill-drain event |
1414
|---|---|
15-
| ![ICESat-2 ATL11 rate of height change over time in Antarctica 2019-03-29 to 2020-12-24](https://user-images.githubusercontent.com/23487320/123902132-65cfd680-d9c0-11eb-88d6-4e0e8c5abc47.png) | ![dsm_whillans_ix_cycles_3-9.gif](https://user-images.githubusercontent.com/23487320/110536564-7b599000-8186-11eb-9ae2-aca8d76f7313.gif) |
15+
| ![ICESat-2 ATL11 rate of height change over time in Antarctica 2019-03-29 to 2020-12-24](https://user-images.githubusercontent.com/23487320/123902132-65cfd680-d9c0-11eb-88d6-4e0e8c5abc47.png) | ![dsm_whillans_ix_cycles_3-9.gif](https://user-images.githubusercontent.com/23487320/124219379-5ed7ce00-db50-11eb-95d0-f1f660d4d688.gif) |
1616

1717
![DeepIceDrain Pipeline Part 1 Exploratory Data Analysis](https://yuml.me/diagram/scruffy;dir:LR/class/[Land-Ice-Elevation|atl06_play.ipynb]->[Convert|atl06_to_atl11.ipynb],[Convert]->[Land-Ice-Height-time-series|atl11_play.ipynb])
1818
![DeepIceDrain Pipeline Part 2 Subglacial Lake Analysis](https://yuml.me/diagram/scruffy;dir:LR/class/[Height-Change-over-Time-(dhdt)|atlxi_dhdt.ipynb],[Height-Change-over-Time-(dhdt)]->[Subglacial-Lake-Finder|atlxi_lake.ipynb],[Subglacial-Lake-Finder]->[Crossover-Analysis|atlxi_xover.ipynb])
1919

2020
| Along track view of an ATL11 Ground Track | Elevation time-series at Crossover Points |
2121
|---|---|
22-
| ![alongtrack_whillans_ix_1080_pt3](https://user-images.githubusercontent.com/23487320/110536370-41888980-8186-11eb-96e6-1ce92aa9966b.png) | ![crossover_anomaly_whillans_ix_2018-10-14_2020-11-11](https://user-images.githubusercontent.com/23487320/110536098-efdfff00-8185-11eb-97d9-065dd59b5727.png) |
22+
| ![alongtrack_whillans_ix_1080_pt3](https://user-images.githubusercontent.com/23487320/124219416-744cf800-db50-11eb-83a1-45e1e1159ba6.png) | ![crossover_anomaly_whillans_ix_2019-03-29_2020-12-24](https://user-images.githubusercontent.com/23487320/124219432-7a42d900-db50-11eb-92b4-c83728b8dc1c.png) |
2323

2424

2525

‎atlxi_lake.ipynb‎

Lines changed: 12 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -873,7 +873,7 @@
873873
},
874874
"outputs": [],
875875
"source": [
876-
"# Save or load dhdt data from Parquet file\n",
876+
"# Load dhdt data from Parquet file\n",
877877
"placename: str = \"siple_coast\" # \"slessor_downstream\" # \"Recovery\" # \"Whillans\"\n",
878878
"df_dhdt: cudf.DataFrame = cudf.read_parquet(\n",
879879
" f\"ATLXI/df_dhdt_{placename.lower()}.parquet\"\n",
@@ -890,26 +890,26 @@
890890
"output_type": "stream",
891891
"text": [
892892
"index 0\n",
893-
"geometry POLYGON ((-444773.5763656513 -545238.396390805...\n",
893+
"geometry POLYGON ((-444731.6953220846 -545129.683759524...\n",
894894
"basin_name Whillans\n",
895-
"refgtracks 74|135|196|266|327|388|516|577|638|769|830|891...\n",
896-
"num_points 4276\n",
897-
"maxabsdhdt 7.26322\n",
898-
"inner_dhdt 1.567\n",
899-
"mean_dhdt 1.91527\n",
900-
"outer_dhdt 0.350537\n",
901-
"outer_std 0.160555\n",
902-
"outer_mad 0.0827204\n",
895+
"refgtracks 74|135|196|266|327|388|577|638|769|830|1019|10...\n",
896+
"num_points 3422\n",
897+
"maxabsdhdt 6.731061\n",
898+
"inner_dhdt 1.152791\n",
899+
"mean_dhdt 1.365484\n",
900+
"outer_dhdt 0.338404\n",
901+
"outer_std 0.151085\n",
902+
"outer_mad 0.081393\n",
903903
"Name: 0, dtype: object\n"
904904
]
905905
},
906906
{
907907
"data": {
908908
"image/svg+xml": [
909-
"<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"300\" height=\"300\" viewBox=\"-456964.4459335574 -545977.9425100257 19967.745218940254 17925.926309279283\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-1074029.958710772)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"133.1183014596017\" opacity=\"0.6\" d=\"M -444773.57636565133,-545238.3963908057 L -445426.23271775024,-545062.6933081353 L -447894.70407518453,-544175.6198389902 L -450167.5050020933,-543148.723574959 L -451936.7037078315,-541987.963628023 L -454210.39452298917,-540418.8503920811 L -455612.023089771,-538987.7429525732 L -456224.8998143374,-537787.353574816 L -454509.0218480702,-532882.2389798673 L -452726.7733748321,-531671.9759466153 L -450107.561159974,-530560.6250711954 L -448142.8086421694,-529956.804608957 L -442105.3774746713,-528791.5623199664 L -441605.42441635387,-528869.279102371 L -439571.86208080617,-530579.3740726819 L -439131.22288676386,-531351.7919350173 L -438656.8873741209,-534025.102285925 L -438422.0274304641,-535974.3104041764 L -438272.83840486367,-537429.4367084408 L -437736.2468338372,-542726.6039209291 L -437776.2069919053,-542836.0370197135 L -444681.03551864496,-545210.4538980428 L -444773.57636565133,-545238.3963908057 z\" /></g></svg>"
909+
"<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"300\" height=\"300\" viewBox=\"-456630.37936665915 -545831.6143030439 18952.12467501691 17100.473795304657\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-1074562.754810783)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"126.34749783344606\" opacity=\"0.6\" d=\"M -444731.6953220846,-545129.6837595248 L -445011.25885498157,-545106.2968520958 L -445376.93950812094,-545031.6591196847 L -447409.3526933879,-544382.2805507262 L -450352.5856256729,-542856.9347974786 L -451613.65405471186,-542088.3567467475 L -453168.7111508558,-541005.8354607922 L -454541.14624465484,-540030.3540272998 L -455129.4508626142,-539161.6296595022 L -455928.44882314,-537936.6574710397 L -455889.53593075665,-537286.6599730677 L -454243.5090539303,-535519.7292632827 L -453242.79123223276,-534521.465268024 L -449104.8354354054,-531478.2867125623 L -448620.22966874443,-531145.8302439799 L -443146.6610664141,-529433.0710512584 L -441727.20288979204,-529844.5068711725 L -439855.5902259358,-530889.0817186871 L -439736.6727377982,-531015.2193078512 L -439630.14850606647,-531141.0021021602 L -438656.6297759172,-534025.5166154493 L -438380.1852351614,-535865.5798414922 L -439494.1268296778,-541257.5824335428 L -441276.22981477395,-542729.7436845967 L -444731.6953220846,-545129.6837595248 z\" /></g></svg>"
910910
],
911911
"text/plain": [
912-
"<shapely.geometry.polygon.Polygon at 0x7fd0bc5b05e0>"
912+
"<shapely.geometry.polygon.Polygon at 0x7fbc201cd670>"
913913
]
914914
},
915915
"execution_count": 12,

‎atlxi_lake.py‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -351,7 +351,7 @@
351351
# # Select a subglacial lake to examine
352352

353353
# %%
354-
# Save or load dhdt data from Parquet file
354+
# Load dhdt data from Parquet file
355355
placename: str = "siple_coast" # "slessor_downstream" # "Recovery" # "Whillans"
356356
df_dhdt: cudf.DataFrame = cudf.read_parquet(
357357
f"ATLXI/df_dhdt_{placename.lower()}.parquet"

‎atlxi_xover.ipynb‎

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -57,7 +57,7 @@
5757
"tag: str = \"X2SYS\"\n",
5858
"os.environ[\"X2SYS_HOME\"] = os.path.abspath(tag)\n",
5959
"client = dask.distributed.Client(\n",
60-
" n_workers=8, threads_per_worker=1, env={\"X2SYS_HOME\": os.environ[\"X2SYS_HOME\"]}\n",
60+
" n_workers=4, threads_per_worker=1, env={\"X2SYS_HOME\": os.environ[\"X2SYS_HOME\"]}\n",
6161
")\n",
6262
"client"
6363
]
@@ -68,7 +68,7 @@
6868
"metadata": {},
6969
"outputs": [],
7070
"source": [
71-
"min_date, max_date = (\"2018-10-14\", \"2020-09-30\")"
71+
"min_date, max_date = (\"2019-03-29\", \"2020-12-24\")"
7272
]
7373
},
7474
{
@@ -105,7 +105,7 @@
105105
"outputs": [],
106106
"source": [
107107
"# Save or load dhdt data from Parquet file\n",
108-
"placename: str = \"siple_coast\" # \"slessor_downstream\" # \"Recovery\" # \"Whillans\"\n",
108+
"placename: str = \"whillans_upstream\" # \"slessor_downstream\"\n",
109109
"df_dhdt: pd.DataFrame = pd.read_parquet(f\"ATLXI/df_dhdt_{placename.lower()}.parquet\")"
110110
]
111111
},
@@ -116,12 +116,12 @@
116116
"outputs": [],
117117
"source": [
118118
"# Choose one Antarctic active subglacial lake polygon with EPSG:3031 coordinates\n",
119-
"lake_name: str = \"Lake 12\"\n",
119+
"lake_name: str = \"Whillans IX\"\n",
120120
"lake_catalog = deepicedrain.catalog.subglacial_lakes()\n",
121-
"lake_ids: list = (\n",
121+
"lake_ids, transect_id = (\n",
122122
" pd.json_normalize(lake_catalog.metadata[\"lakedict\"])\n",
123-
" .query(\"lakename == @lake_name\")\n",
124-
" .ids.iloc[0]\n",
123+
" .query(\"lakename == @lake_name\")[[\"ids\", \"transect\"]]\n",
124+
" .iloc[0]\n",
125125
")\n",
126126
"lake = (\n",
127127
" lake_catalog.read()\n",
@@ -147,7 +147,7 @@
147147
"source": [
148148
"# Subset data to lake of interest\n",
149149
"placename: str = region.name.lower().replace(\" \", \"_\")\n",
150-
"df_lake: cudf.DataFrame = region.subset(data=df_dhdt) # bbox subset\n",
150+
"df_lake: pd.DataFrame = region.subset(data=df_dhdt) # bbox subset\n",
151151
"gdf_lake = gpd.GeoDataFrame(\n",
152152
" df_lake, geometry=gpd.points_from_xy(x=df_lake.x, y=df_lake.y, crs=3031)\n",
153153
")\n",
@@ -362,7 +362,7 @@
362362
" text=f\"Track {track1} and {track2} crossover\",\n",
363363
" position=\"TC\",\n",
364364
" offset=\"jTC0c/0.2c\",\n",
365-
" V=\"q\",\n",
365+
" verbose=\"q\",\n",
366366
")\n",
367367
"# Plot data points\n",
368368
"fig.plot(x=df_max.t, y=df_max.h, style=\"c0.15c\", color=\"darkblue\", pen=\"thin\")\n",

‎atlxi_xover.py‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -7,7 +7,7 @@
77
# extension: .py
88
# format_name: hydrogen
99
# format_version: '1.3'
10-
# jupytext_version: 1.7.1
10+
# jupytext_version: 1.11.3
1111
# kernelspec:
1212
# display_name: deepicedrain
1313
# language: python
@@ -53,13 +53,13 @@
5353
tag: str = "X2SYS"
5454
os.environ["X2SYS_HOME"] = os.path.abspath(tag)
5555
client = dask.distributed.Client(
56-
n_workers=8, threads_per_worker=1, env={"X2SYS_HOME": os.environ["X2SYS_HOME"]}
56+
n_workers=4, threads_per_worker=1, env={"X2SYS_HOME": os.environ["X2SYS_HOME"]}
5757
)
5858
client
5959

6060

6161
# %%
62-
min_date, max_date = ("2018-10-14", "2020-09-30")
62+
min_date, max_date = ("2019-03-29", "2020-12-24")
6363

6464
# %%
6565
# Initialize X2SYS database in the X2SYS/ICESAT2 folder
@@ -78,18 +78,18 @@
7878

7979
# %%
8080
# Save or load dhdt data from Parquet file
81-
placename: str = "siple_coast" # "slessor_downstream" # "Recovery" # "Whillans"
81+
placename: str = "whillans_upstream" # "slessor_downstream"
8282
df_dhdt: pd.DataFrame = pd.read_parquet(f"ATLXI/df_dhdt_{placename.lower()}.parquet")
8383

8484

8585
# %%
8686
# Choose one Antarctic active subglacial lake polygon with EPSG:3031 coordinates
87-
lake_name: str = "Lake 12"
87+
lake_name: str = "Whillans IX"
8888
lake_catalog = deepicedrain.catalog.subglacial_lakes()
89-
lake_ids: list = (
89+
lake_ids, transect_id = (
9090
pd.json_normalize(lake_catalog.metadata["lakedict"])
91-
.query("lakename == @lake_name")
92-
.ids.iloc[0]
91+
.query("lakename == @lake_name")[["ids", "transect"]]
92+
.iloc[0]
9393
)
9494
lake = (
9595
lake_catalog.read()
@@ -107,7 +107,7 @@
107107
# %%
108108
# Subset data to lake of interest
109109
placename: str = region.name.lower().replace(" ", "_")
110-
df_lake: cudf.DataFrame = region.subset(data=df_dhdt) # bbox subset
110+
df_lake: pd.DataFrame = region.subset(data=df_dhdt) # bbox subset
111111
gdf_lake = gpd.GeoDataFrame(
112112
df_lake, geometry=gpd.points_from_xy(x=df_lake.x, y=df_lake.y, crs=3031)
113113
)
@@ -270,7 +270,7 @@
270270
text=f"Track {track1} and {track2} crossover",
271271
position="TC",
272272
offset="jTC0c/0.2c",
273-
V="q",
273+
verbose="q",
274274
)
275275
# Plot data points
276276
fig.plot(x=df_max.t, y=df_max.h, style="c0.15c", color="darkblue", pen="thin")

‎deepicedrain/features/subglacial_lakes.feature‎

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -33,7 +33,6 @@ Feature: Mapping Antarctic subglacial lakes
3333
# | Kamb 8 | whillans_upstream | 3-9 | 157.5 | 45 |
3434
# | Kamb 1 | whillans_upstream | 3-9 | 157.5 | 45 |
3535
| Kamb 34 | whillans_upstream | 4-9 | 157.5 | 45 |
36-
# | Kamb 12 | siple_coast | 3-8 | 157.5 | 45 |
3736
# | MacAyeal 1 | siple_coast | 3-8 | 157.5 | 60 |
3837
# | Slessor 45 | slessor_downstream | 3-9 | 202.5 | 60 |
3938
# | Slessor 23 | slessor_downstream | 3-9 | 202.5 | 60 |

0 commit comments

Comments
 (0)