skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

This content will become publicly available on January 1, 2024

Title: GLAcier Feature Tracking testkit (GLAFT): a statistically and physically based framework for evaluating glacier velocity products derived from optical satellite image feature tracking
Abstract. Glacier velocity measurements are essential to understand ice flow mechanics, monitor natural hazards, and make accurate projections of future sea-level rise. Despite these important applications, the method most commonly used to derive glacier velocity maps, feature tracking, relies on empirical parameter choices that rarely account for glacier physics or uncertainty. Here we test two statistics- and physics-based metrics to evaluate velocity maps derived from optical satellite images of Kaskawulsh Glacier, Yukon, Canada, using a range of existing feature-tracking workflows. Based on inter-comparisons with ground truth data, velocity maps with metrics falling within our recommended ranges contain fewer erroneous measurements and more spatially correlated noise than velocity maps with metrics that deviate from those ranges. Thus, these metric ranges are suitable for refining feature-tracking workflows and evaluating the resulting velocity products. We have released an open-source software package for computing and visualizing these metrics, the GLAcier Feature Tracking testkit (GLAFT).  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Page Range / eLocation ID:
4063 to 4078
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. We present Glacier Image Velocimetry (GIV), an open-source and easy-to-use software toolkit for rapidly calculating high-spatial-resolutionglacier velocity fields. Glacier ice velocity fields reveal flow dynamics, ice-flux changes, and (with additional data and modelling) icethickness. Obtaining glacier velocity measurements over wide areas with field techniques is labour intensive and often associated with safetyrisks. The recent increased availability of high-resolution, short-repeat-time optical imagery allows us to obtain ice displacement fields using“feature tracking” based on matching persistent irregularities on the ice surface between images and hence, surface velocity over time. GIV isfully parallelized and automatically detects, filters, and extracts velocities from large datasets of images. Through this coupled toolchain and aneasy-to-use GUI, GIV can rapidly analyse hundreds to thousands of image pairs on a laptop or desktop computer. We present four example applicationsof the GIV toolkit in which we complement a glaciology field campaign (Glaciar Perito Moreno, Argentina) and calculate the velocity fields of smallmid-latitude (Glacier d'Argentière, France) and tropical glaciers (Volcán Chimborazo, Ecuador), as well as very large glaciers (Vavilov Ice Cap,Russia). Fully commented MATLAB code and a stand-alone app for GIV are available from GitHub and Zenodo (see, Van Wyk de Vries, 2021a). 
    more » « less
  2. Abstract. The Thwaites Eastern Ice Shelf buttresses a significant portion of Thwaites Glacier through contact with a pinning point 40 km offshore of the present grounding line. Predicting future rates of Thwaites Glacier’s contribution to sea-level rise depends on the evolution of this pinning point and the resultant change in the ice-shelf stress field since the breakup of the Thwaites Western Glacier Tongue in 2009. Here we use Landsat-8 feature tracking of ice velocity in combination with ice-sheet model perturbation experiments to show how past changes in flow velocity have been governed in large part by changes in lateral shear and pinning point interactions with the Thwaites Western Glacier Tongue. We then use recent satellite altimetry data from ICESat-2 to show that Thwaites Glacier’s grounding line has continued to retreat rapidly; in particular, the grounded area of the pinning point is greatly reduced from earlier mappings in 2014, and grounded ice elevations are continuing to decrease. This loss has created two pinned areas with ice flow now funneled between them. If current rates of surface lowering persist, the Thwaites Eastern Ice Shelf will unpin from the seafloor in less than a decade, despite our finding from airborne radar data that the seafloor underneath the pinning point is about 200 m shallower than previously reported. Advection of relatively thin and mechanically damaged ice onto the remaining portions of the pinning point and feedback mechanisms involving basal melting may further accelerate the unpinning. As a result, ice discharge will likely increase up to 10 % along a 45 km stretch of the grounding line that is currently buttressed by the Thwaites Eastern Ice Shelf. 
    more » « less
  3. Abstract. Rock glaciers are a prominent component of many alpine landscapes andconstitute a significant water resource in some arid mountainenvironments. Here, we employ satellite-based interferometric syntheticaperture radar (InSAR) between 2016 and 2019 to identify and monitor activeand transitional rock glaciers in the Uinta Mountains (Utah, USA), an area of∼3000 km2. We used mean velocity maps to generate aninventory for the Uinta Mountains containing 205 active and transitional rockglaciers. These rock glaciers are 11.9 ha in area on average andlocated at a mean elevation of 3308 m, where mean annual airtemperature is −0.25 ∘C. The mean downslope velocity for theinventory is 1.94 cm yr−1, but individual rock glaciers have velocities ranging from0.35 to 6.04 cm yr−1. To search for relationships with climaticdrivers, we investigated the time-dependent motion of three rock glaciers. Wefound that rock glacier motion has a significant seasonal component, withrates that are more than 5 times faster during the late summer compared to therest of the year. Rock glacier velocities also appear to be correlated withthe snow water equivalent of the previous winter's snowpack. Our resultsdemonstrate the ability to use satellite InSAR to monitor rock glaciers overlarge areas and provide insight into the environmental factors that controltheir kinematics. 
    more » « less
  4. Abstract

    Conservation planning and decision‐making rely on evaluations of biodiversity status and threats that are based upon species' distribution estimates. However, gaps exist regarding automated tools to delineate species' current ranges from distribution estimates and use those estimates to calculate both species‐ and community‐level biodiversity metrics. Here, we introduce changeRangeR, an R package that facilitates workflows to reproducibly transform estimates of species' distributions into metrics relevant for conservation. For example, by combining predictions from species distribution models (SDMs) with other maps of environmental data (e.g., suitable forest cover), researchers can characterize the proportion of a species' range that is under protection, metrics used under the IUCN Criteria A and B guidelines (Area of Occupancy and Extent of Occurrence), and other more general metrics such as taxonomic and phylogenetic diversity and endemism. Further, changeRangeR facilitates temporal comparisons among biodiversity metrics to inform efforts toward complementarity and consideration of future scenarios in conservation decisions. changeRangeR also provides tools to determine the effects of modeling decisions through sensitivity tests. Transparent and repeatable workflows for calculating biodiversity change metrics from SDMs such as those provided by changeRangeR are essential to inform conservation decision‐making efforts and represent key extensions for SDM methodology and associated metadata documentation.

    more » « less
  5. Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBMwas first published in 1989; here, we describe the first fully open-sourceWBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of themodel functionality, utility, and evaluation of simulated global riverdischarge and irrigation water use. This new version adds a novel suite ofwater source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions. 
    more » « less