This repository contains algorithms, datasets, and techniques referenced in following publication and digital object identifier:
King, J., Howell, S., Brady, M., Toose, P., Derksen, C., Haas, C., and Beckers, J.: Local-scale variability of snow density on Arctic sea ice, The Cryosphere, https://doi.org/10.5194/tc-14-4323-2020, 14(12), 4323–4339, 2020.
Local-scale variations in snow density and layering on Arctic sea ice were characterized using a combination of traditional snow pit and SnowMicroPen (SMP) measurements. In total, 14 sites were evaluated within the Canadian Arctic Archipelago and Arctic Ocean on both first (FYI) and multi-year (MYI) sea ice. Sites contained multiple snow pits with coincident SMP profiles as well as unidirectional SMP transects. An existing SMP density model was recalibrated using manual density cutter measurements (n=186) to identify best-fit parameters for the observed conditions. Cross-validation of the revised SMP model showed errors comparable to the expected baseline for manual density measurements (RMSE=34 kg m-3 or 10.9%) and strong retrieval skill (R2=0.78). The density model was then applied to SMP transect measurements to characterize variations at spatial scales of up to 100 m. A supervised classification trained on snow pit stratigraphy allowed separation of the SMP density estimates by layer-type. The resulting dataset contains 58,882 layer-classified estimates of snow density on sea ice representing 147 m of vertical variation and equivalent to more than 600 individual snow pits. An average bulk density of 310 kg m-3 was estimated with clear separation between FYI and MYI environments. Lower densities on MYI (277 kg m-3) corresponded with increased depth hoar composition (49.2%), in strong contrast to composition of the thin FYI snowpack (19.8%). Spatial auto-correlation analysis showed layered composition on FYI snowpack to persist over long distances while composition on MYI rapidly decorrelated at distances less than 16 m. Application of the SMP profiles to determine propagation bias in radar altimetry showed the potential errors of 0.5 cm when climatology is used over known snow density.
Figure 1c: Unidirectional SnowMicoPen (SMP) transects were collected at multiple sites to evaluate spatial variability of snowpack properties with sets of 10 profiles separated at distances of 0.1, 1, and 10 m, in sequence
There are several Jupyter notebooks within this repository which reproduce the various algorithms and figures used in the publication. To run the notebooks yourself, there are a couple of options:
- Download and install Miniconda
- Clone this repository:
$ git clone https://github.com/kingjml/SMP-Sea-Ice.git
- Create the
smp-sea-ice
environment using this repository's environment specification:$ conda env create -f environment.yml
- Activate the
smp-sea-ice
environment, and launch Jupyter to explore the notebooks:$ conda activate smp-sea-ice (smp-sea-ice) $ jupyter notebook index.ipynb
- (Optional) If you would like to run the brute-force measurement alignment routine in notebook #1 yourself, you will need to retrieve the SMP datasets from this link