Explore the vast range of titles available.

Velocities within an icefall are typically the fastest within a glacier system and experience complex flow. The combination of convergent and fast flow, and steep slope generate a quickly changing and intensely fractured surface. This complicates velocity extraction from repeat satellite images, especially when common pattern matching procedures are used. In this study, we exploit the high temporal revisit of medium-resolution satellite images using a novel image matching technique, ensemble matching, making it possible to generate a high-resolution (30 m) velocity field from high-repeat image sequences despite challenging image conditions. We demonstrate this technique for the first time in the glaciology domain using repeat Sentinel-2 optical data over the famous Khumbu icefall, situated on the southern slopes of Mount Everest. Estimates of velocity go just over 1 m d−1, which is slower than summer velocities from noisy single pair image matching. This icefall is frequently crossed by high-altitude mountaineers who use a route confined by fixed ropes and ladders set out every season. The mountain climbers typically record their trajectory on their personal satellite navigation device. We use such volunteered geographic information to verify our velocity estimates, confirming our underestimation with ensemble matching. Besides unprecedented remotely sensed surface velocities over the icefall, we also note that the generated velocity field can aid with the planning of a safe passage through this icefall.

The PlanetScope constellation consists of ∼150 optical cubesats that are evenly distributed like strings of pearls on two orbital planes, scanning the Earth's land surface once per day with an approximate spatial image resolution of 3 m. Subsequent cubesats on each of the orbital planes image the Earth surface with a nominal time lag of approximately 90 s between them, which produces near-simultaneous image pairs over the across-track overlaps of the cubesat swaths. We exploit this short time lag between subsequent Planet cubesat images to track river ice floes on northern rivers as indicators of water surface velocities. The method is demonstrated for a 60 km long reach of the Amur River in Siberia, and a 200 km long reach of the Yukon River in Alaska. The accuracy of the estimated horizontal surface velocities is of the order of ±0.01 m s−1. The application of our approach is complicated by cloud cover and low sun angles at high latitudes during the periods where rivers typically carry ice floes, and by the fact that the near-simultaneous swath overlaps, by design, do not cover the complete Earth surface. Still, the approach enables direct remote sensing of river surface velocities for numerous cold-region rivers at a number of locations and occasionally several times per year – which is much more frequent and over much larger areas than currently feasible. We find that freeze-up conditions seem to offer ice floes that are generally more suitable for tracking, and over longer time periods, compared with typical ice break-up conditions. The coverage of river velocities obtained could be particularly useful in combination with satellite measurements of river area, and river surface height and slope.

Satellite measurements of coseismic displacements are typically based on synthetic aperture radar (SAR) interferometry or amplitude tracking, or based on optical data such as from Landsat, Sentinel-2, SPOT, ASTER, very high-resolution satellites, or air photos. Here, we evaluate a new class of optical satellite images for this purpose – data from cubesats. More specific, we investigate the PlanetScope cubesat constellation for horizontal surface displacements by the 14 November 2016 Mw 7.8 Kaikoura, New Zealand, earthquake. Single PlanetScope scenes are 2–4 m-resolution visible and near-infrared frame images of approximately 20–30 km  ×  9–15 km in size, acquired in continuous sequence along an orbit of approximately 375–475 km height. From single scenes or mosaics from before and after the earthquake, we observe surface displacements of up to almost 10 m and estimate matching accuracies from PlanetScope data between ±0.25 and ±0.7 pixels (∼ ±0.75 to ±2.0 m), depending on time interval and image product type. Thereby, the most optimistic accuracy estimate of ±0.25 pixels might actually be typical for the final, sun-synchronous, and near-polar-orbit PlanetScope constellation when unrectified data are used for matching. This accuracy, the daily revisit anticipated for the PlanetScope constellation for the entire land surface of Earth, and a number of other features, together offer new possibilities for investigating coseismic and other Earth surface displacements and managing related hazards and disasters, and complement existing SAR and optical methods. For comparison and for a better regional overview we also match the coseismic displacements by the 2016 Kaikoura earthquake using Landsat 8 and Sentinel-2 data.

With its temporal resolution of 10 days (five days with two satellites, and significantly more at high latitudes), its swath width of 290 km, and its 10 m and 20 m spatial resolution bands from the visible to the shortwave infrared, the European Sentinel-2 satellites have significant potential for glacier remote sensing, in particular mapping of glacier outlines and facies, and velocity measurements. Testing Level 1C commissioning and ramp-up phase data for initial sensor quality experiences, we find a high radiometric performance, but with slight striping effects under certain conditions. Through co-registration of repeat Sentinal-2 data we also find lateral offset patterns and noise on the order of a few metres. Neither of these issues will complicate most typical glaciological applications. Absolute geo-location of the data investigated was on the order of one pixel at the time of writing. The most severe geometric problem stems from vertical errors of the DEM used for ortho-rectifying Sentinel-2 data. These errors propagate into locally varying lateral offsets in the images, up to several pixels with respect to other georeferenced data, or between Sentinel-2 data from different orbits. Finally, we characterize the potential and limitations of tracking glacier flow from repeat Sentinel-2 data using a set of typical glaciers in different environments: Aletsch Glacier, Swiss Alps; Fox Glacier, New Zealand; Jakobshavn Isbree, Greenland; Antarctic Peninsula at the Larsen C ice shelf.

Optical satellite products are available at different processing levels. Of these products, terrain corrected (i.e., orthorectified) products are the ones mostly used for glacier displacement estimation. For terrain correction, a digital elevation model (DEM) is used that typically stems from various data sources with variable qualities, from dispersed time instances, or with different spatial resolutions. Consequently, terrain representation used for orthorectifying satellite images is often in disagreement with reality at image acquisition. Normally, the lateral orthoprojection offsets resulting from vertical DEM errors are taken into account in the geolocation error budget of the corrected images, or may even be neglected. The largest offsets of this type are often found over glaciers, as these may show strong elevation changes over time and thus large elevation errors in the reference DEM with respect to image acquisition. The detection and correction of such orthorectification offsets is further complicated by ice flow which adds a second offset component to the displacement vectors between orthorectified data. Vice versa, measurement of glacier flow is complicated by the inherent superposition of ice movement vectors and orthorectification offset vectors. In this study, we try to estimate these orthorectification offsets in the presence of terrain movement and translate them to elevation biases in the reference surface. We demonstrate our method using three different sites which include very dynamic glaciers. For the Oriental Glacier, an outlet of the Southern Patagonian icefield, Landsat 7 and 8 data from different orbits enabled the identification of trends related to elevation change. For the Aletsch Glacier, Swiss Alps, we assess the terrain offsets of both Landsat 8 and Sentinel-2A: a superior DEM appears to be used for Landsat in comparison to Sentinel-2, however a systematic bias is observed in the snow covered areas. Lastly, we demonstrate our methodology in a pipeline structure; displacement estimates for the Helheim-glacier, in Greenland, are mapped and corrected for orthorectification offsets between data from different orbits, which enables a twice as dense a temporal resolution of velocity data, as compared to the standard method of measuring velocities from repeat-orbit data only. In addition, we introduce and implement a novel matching method which uses image triplets. By formulating the three image displacements as a convolution, a geometric constraint can be exploited. Such a constraint enhances the reliability of the displacement estimations. Furthermore the implementation is simple and computationally swift.

Negribreen, a tidewater glacier located in central eastern Svalbard, began actively surging after it experienced an initial collapse in summer 2016. The surge resulted in horizontal surface velocities of more than 25 m d−1, making it one of the fastest-flowing glaciers in the archipelago. The last surge of Negribreen likely occurred in the 1930s, but due to a long quiescent phase, investigations of this glacier have been limited. As Negribreen is part of the Negribreen Glacier System, one of the largest glacier systems in Svalbard, investigating its current surge event provides important information on surge behaviour among tidewater glaciers within the region. Here, we demonstrate the surge development and discuss triggering mechanisms using time series of digital elevation models (1969–2018), surface velocities (1995–2018), crevasse patterns and glacier extents from various data sources. We find that the active surge results from a four-stage process. Stage 1 (quiescent phase) involves a long-term, gradual geometry change due to high subglacial friction towards the terminus. These changes allow the onset of Stage 2, an accelerating frontal destabilization, which ultimately results in the collapse (Stage 3) and active surge (Stage 4).

In this study, we combine remote sensing, in situ and model-derived datasets from 1966 to 2014 to calculate the mass-balance components of Kronebreen, a fast-flowing tidewater glacier in Svalbard. For the well-surveyed period 2009–2014, we are able to close the glacier mass budget within the prescribed errors. During these 5 years, the glacier geodetic mass balance was −0.69 ± 0.12 m w.e. a−1, while the mass budget method led to a total mass balance of −0.92 ± 0.16 m w.e. a−1, as a consequence of a strong frontal ablation (−0.78 ± 0.11 m w.e. a−1), and a slightly negative climatic mass balance (−0.14 ± 0.11 m w.e. a−1). The trend towards more negative climatic mass balance between 1966–1990 (+0.20 ± 0.05 m w.e. a−1) and 2009–2014 is not reflected in the geodetic mass balance trend. Therefore, we suspect a reduction in ice-discharge in the most recent period. Yet, these multidecadal changes in ice-discharge cannot be measured from the available observations and thus are only estimated with relatively large errors as a residual of the mass continuity equation. Our study presents the multidecadal evolution of the dynamics and mass balance of a tidewater glacier and illustrates the errors introduced by inferring one unmeasured mass-balance component from the others.

Greenland basal topographic data show a segmented valley extending from Petermann Fjord into the centre of Greenland; however, the locations of radar scan lines, used to create the bedrock topography data, indicate that valley segmentation is due to data interpolation. Therefore, as a thought experiment, simulations where the valley is opened are used to investigate its effects on basal water movement and distribution. The simulations indicate that the opening of this valley can result in an uninterrupted water pathway from the interior to Petermann Fjord. Along its length, the path of the valley progresses gradually down an ice surface slope, causing a lowering of ice overburden pressure that could enable water flow along its path. The fact that the valley base appears to be relatively flat and follows a path near the interior ice divide that roughly intersects the east and west basal hydrological basins is presented as evidence that its present day form may have developed in conjunction with an overlying ice sheet. Experiments where basal melting is increased solely within the deep interior near the known large area of basal melting result in an increase in the flux of water northwards along the entire valley. The results are consistent with a long subglacial river; however, considerable uncertainty remains over aspects such as whether adequate water is available at the bed, whether water escapes from the valley or is refrozen, and over what form a hydrological conduit could take along the valley base.

For studying the flow of glaciers and their response to climate change it is important to detect glacier surges. Here, we compute within Google Earth Engine the normalized differences between winter maxima of Sentinel-1 C-band radar backscatter image stacks over subsequent years. We arrive at a global map of annual backscatter changes, which are for glaciers in most cases related to changed crevassing associated with surge-type activity. For our demonstration period 2018–2019 we detected 69 surging glaciers, with many of them not classified so far as surge type. Comparison with glacier surface velocities shows that we reliably find known surge activities. Our method can support operational monitoring of glacier surges and some other special events such as large rock and snow avalanches.

Satellite measurements of coseismic displacements are typically based on synthetic aperture radar (SAR) interferometry or amplitude tracking, or based on optical data such as from Landsat, Sentinel-2, SPOT, ASTER, very high-resolution satellites, or air photos. Here, we evaluate a new class of optical satellite images for this purpose – data from cubesats. More specific, we investigate the PlanetScope cubesat constellation for horizontal surface displacements by the 14 November 2016 Mw 7.8 Kaikoura, New Zealand, earthquake. Single PlanetScope scenes are 2–4 m-resolution visible and near-infrared frame images of approximately 20–30 km  ×  9–15 km in size, acquired in continuous sequence along an orbit of approximately 375–475 km height. From single scenes or mosaics from before and after the earthquake, we observe surface displacements of up to almost 10 m and estimate matching accuracies from PlanetScope data between ±0.25 and ±0.7 pixels (∼ ±0.75 to ±2.0 m), depending on time interval and image product type. Thereby, the most optimistic accuracy estimate of ±0.25 pixels might actually be typical for the final, sun-synchronous, and near-polar-orbit PlanetScope constellation when unrectified data are used for matching. This accuracy, the daily revisit anticipated for the PlanetScope constellation for the entire land surface of Earth, and a number of other features, together offer new possibilities for investigating coseismic and other Earth surface displacements and managing related hazards and disasters, and complement existing SAR and optical methods. For comparison and for a better regional overview we also match the coseismic displacements by the 2016 Kaikoura earthquake using Landsat 8 and Sentinel-2 data.