Ice Flow Velocity Research Articles

New metrics reveal the evolutionary continuum behind the morphological diversity of subglacial bedforms

Understanding the formation of subglacial bedforms is primordial to constrain ice-meltwater-bed interactions and the dynamics of past and present ice sheets. However, the difficulty to observe active subglacial bedforms below present-day ice sheets implies that their formation and evolution are essentially deduced from the inversion of morphological and sedimentological data from palaeo-subglacial bedforms. In this study, the morphological characteristics of subglacial bedforms are explored with a new approach based on the combination of three dimensionless morphometric indices i.e., circularity index, sinuosity index and elongation component ratio. We measure the spatial distribution of these indices on an unpublished database composed of ~13,500 digitized bedforms taken from two selected portions of the Irish Ice Sheet and Laurentide Ice Sheet beds considering that (1) all subglacial bedforms in each region were formed under a single ice flow configuration and (2) the bedforms may have developed either orthogonal, parallel, or oblique to ice displacement directions. Our results reveal a morphometric and spatial continuum along which ribbed bedforms of various lengths and shapes – ranging from circular to elongated flow-transverse forms – are incipient bedforms that can evolve into transitional sinuous ribbed bedforms, then into predominantly flow-parallel sinuous bedforms that progressively realign to form streamlined bedforms. Assuming the glaciological contexts of selected study areas, this continuum may provide a new geomorphic criterion to constrain spatial variations in subglacial conditions (e.g. ice flow velocity, bedrock characteristics, meltwater pressure).

An Improved Intersection Method for Positioning and Its Application to Dalk Glacier, East Antarctica

Stake measurement is a significant ground-observation method used to monitor glacier dynamics in the cryosphere and validate the reliability of remotely sensed investigations. However, some stake observation records obtained by total stations have not been efficiently utilized via the classic forward intersection method. The forward intersection method requires a pair of cognominal observation records and cannot function if stakes are only observed from one control point at a stake measurement, which commonly occurs in measurements. Therefore, we propose an improved intersection method that combines stake trajectories and sightlines, the trajectory-sightline intersection method, which requires only a single observation record. In the processing of stake identification, location and elevation from the trajectory-sightline intersection method are constrained by the velocity and surface altitude of the glacier. Then, the trajectory-sightline intersection method is applied to the stake measurement of Dalk Glacier in East Antarctica. The results show that the trajectory-sightline intersection method extracted an additional 44 intersections, thus increasing the data utilization ratio by 20.35%. Moreover, leave-p-out cross-validation indicated that the trajectory-sightline intersection method, with an average offset of 0.88±0.54 m in the horizontal position and 0.26±0.12 m at the elevation, has similar accuracy to the differential global positioning system (DGPS) applied to Livingston Island. The ice flow velocities of Dalk Glacier calculated from trajectory-sightline intersections agree with those from forward intersections (the mean velocity difference is 0.47 m·year−1). Therefore, the trajectory-sightline intersection method is an effective method for improving data utilization and obtaining stake locations and glacier velocities.

Timescales of outlet-glacier flow with negligible basal friction: theory, observations and modeling

Abstract. The timescales of the flow and retreat of Greenland's and Antarctica's outlet glaciers and their potential instabilities are arguably the largest uncertainty in future sea-level projections. Here we derive a scaling relation that allows the comparison of the timescales of observed complex ice flow fields with geometric similarity. The scaling relation is derived under the assumption of fast, laterally confined, geometrically similar outlet-glacier flow over a slippery bed, i.e., with negligible basal friction. According to the relation, the time scaling of the outlet flow is determined by the product of the inverse of (1) the fourth power of the width-to-length ratio of its confinement, (2) the third power of the confinement depth and (3) the temperature-dependent ice softness. For the outflow at the grounding line of streams with negligible basal friction, this means that the volume flux is proportional to the ice softness and the bed depth, but goes with the fourth power of the gradient of the bed and with the fifth power of the width of the stream. We show that the theoretically derived scaling relation is supported by the observed velocity scaling of outlet glaciers across Greenland as well as by idealized numerical simulations of marine ice-sheet instabilities (MISIs) as found in Antarctica. Assuming that changes in the ice-flow velocity due to ice-dynamic imbalance are proportional to the equilibrium velocity, we combine the scaling relation with a statistical analysis of the topography of 13 MISI-prone Antarctic outlets. Under these assumptions, the timescales in response to a potential destabilization are fastest for Thwaites Glacier in West Antarctica and Mellor, Ninnis and Cook Glaciers in East Antarctica; between 16 and 67 times faster than for Pine Island Glacier. While the applicability of our results is limited by several strong assumptions, the utilization and potential further development of the presented scaling approach may help to constrain timescale estimates of outlet-glacier flow, augmenting the commonly exploited and comparatively computationally expensive approach of numerical modeling.

Glaciological history and structural evolution of the Shackleton Ice Shelf system, East Antarctica, over the past 60 years

Abstract. The discovery of Antarctica's deepest subglacial trough beneath the Denman Glacier, combined with high rates of basal melt at the grounding line, has caused significant concern over its vulnerability to retreat. Recent attention has therefore been focusing on understanding the controls driving Denman Glacier's dynamic evolution. Here we consider the Shackleton system, comprised of the Shackleton Ice Shelf, Denman Glacier, and the adjacent Scott, Northcliff, Roscoe and Apfel glaciers, about which almost nothing is known. We widen the context of previously observed dynamic changes in the Denman Glacier to the wider region of the Shackleton system, with a multi-decadal time frame and an improved biannual temporal frequency of observations in the last 7 years (2015–2022). We integrate new satellite observations of ice structure and airborne radar data with changes in ice front position and ice flow velocities to investigate changes in the system. Over the 60-year period of observation we find significant rift propagation on the Shackleton Ice Shelf and Scott Glacier and notable structural changes in the floating shear margins between the ice shelf and the outlet glaciers, as well as features indicative of ice with elevated salt concentration and brine infiltration in regions of the system. Over the period 2017–2022 we observe a significant increase in ice flow speed (up to 50 %) on the floating part of Scott Glacier, coincident with small-scale calving and rift propagation close to the ice front. We do not observe any seasonal variation or significant change in ice flow speed across the rest of the Shackleton system. Given the potential vulnerability of the system to accelerating retreat into the overdeepened, potentially sediment-filled bedrock trough, an improved understanding of the glaciological, oceanographic and geological conditions in the Shackleton system are required to improve the certainty of numerical model predictions, and we identify a number of priorities for future research. With access to these remote coastal regions a major challenge, coordinated internationally collaborative efforts are required to quantify how much the Shackleton region is likely to contribute to sea level rise in the coming centuries.

Variations in hydraulic efficiency of the subglacial drainage landsystem control surging and streaming regimes of outlet glaciers

AbstractSurging and streaming of glaciers are modulated by meltwater availability and pressure which controls mechanical coupling at their beds. Using laboratory-scale experimental modelling and palaeoglaciological mapping, we explore how subglacial drainage landsystems control meltwater drainage efficiency and ice flow velocities for terrestrial-based ice lobes resting on flat horizontal and permeable beds. Two end-members regimes, surging and streaming, appear in our experiments. The surge regime is characterised by a rapid increase of drainage efficiency through development of tunnel valleys and their tributaries, thus reducing the duration of ice flow speed-up events by lowering water pressures and increasing ice-bed coupling. Tunnel valleys connected to ice lobe margins, submarginal thrust moraines, reduced ice lobe extensions and ephemeral shear margins are the most distinctive characteristics of this regime. The stream regime is characterised by disconnected channels of smaller dimensions unable to evacuate all the meltwater: this prolonged drainage inefficiency leads to sustained high ice flow velocity and steady shear margins. Small and rectilinear meltwater channels devoid of tributaries, often disconnected from ice lobe margins, and lineation swarms are diagnostic of this regime.

Mass balance, ice volume, and flow velocity of the Vestre Grønfjordbreen (Svalbard) from 2013/14 to 2019/20

ABSTRACT The first seven years (2013/14–2019/20) of annual and seasonal mass-balance monitoring on the glacier Vestre Grønfjordbreen (16.4 km2), located south of the town of Barentsburg on Spitsbergen, Svalbard, are presented. This part of the archipelago is one of the least glaciated on Svalbard and at the same time it experiences the most prominent glacier retreat within the last few decades. The annual mass balance of Vestre Grønfjordbreen is negative, ranging from −0.60 ± 0.18 to −2.01 ± 0.26 m w.e. The results of direct observations are compared with the geodetic mass balance for the same period (July 2015 through end of summer 2019) to identify systematic bias in the record. As the mismatch between cumulative mass balances, defined by the glaciological method (−5.66 ± 0.47 m w.e.) and computed from geodetic differencing (−5.52 ± 0.40 m w.e.), lies within the uncertainty limits, no calibration of the mass-balance series is needed. From the results of a ground-penetrating radar (GPR) survey (spring 2019), which confirmed the polythermal glacier structure, a total glacier volume of 1.987 ± 0.139 km3 was found, meaning that the cumulative mass loss during the reported seven-year period equals 8 ± 1% of the total glacier mass. Observed annual ice-flow velocities, varying from 0.50 ± 0.10 to 4.50 ± 0.10 m year−1, are consistent with low mean bed and surface slopes (5° and 8°, respectively). Correlations of mass-balance values with meteorological observations at the Barentsburg weather station are mediocre, possibly due to anomalous values recorded for 2015/16: the negative mass-balance peak reported for the other land-terminating Svalbard glaciers was not observed at Vestre Grønfjordbreen.

Mapping Ice Flow Velocity of Tidewater Glaciers in Hornsund Fiord Area with the Use of Autonomous Repeat Image Feature Tracking (2018–2022)

Dynamic climate changes are particularly apparent in polar regions. Glaciers are retreatng towards the land at a very fast pace. This study demonstrates the application of the feature tracking method in the analysis of ice flow velocity in the region of the Hornsund fiord, southern Spitsbergen, in the years 2018–2022. The calculations were based on the Geogrid and autoRIFT environments and on the Sentinel 1 images. The study also employed external data, such as a numerical terrain model and reference velocity values. The input data, e.g., the chip size and the search limit, were prepared accounting for the specific character of the investigated area. The velocities were calculated for nine biggest glaciers which terminated in the fiord. The accuracy of the results was identified by calculating the median absolute deviation (MAD) of the obtained displacement velocity values from the reference value for areas identified as stable. The study also attempted a causal analysis of the influence of weather factors on the dynamics of ice mass displacement. A systematic year-to-year decrease of the velocity was observed for the entire fiord. In the case of several glaciers, changes related to the ablation season (summer) are also clearly visible. The research results are promising and fill a research gap related to the absence of permanent monitoring and analysis of the dynamics of ice flow in polar regions. It is the first complex and precise study of glacier surface velocity changes, performed on the basis of satellite radar images for the entire Hornsund fiord.

A hierarchical network densification approach for reconstruction of historical ice velocity fields in East Antarctica

AbstractAccurate ice flow velocity data are essential for studying the mass balance of the Antarctic ice sheet. However, there is a lack of ice velocity maps of 1960s–80s in basin-wide regions or the entire ice sheet. In this study, an enhanced hierarchical network densification approach is developed for basin-wide Antarctic velocity mapping using historical ARGON and Landsat images. The produced multiple historical velocity maps from 1963 to 1989 in the region of the Fimbul and Jelbart ice shelves, East Antarctica, achieved an accuracy better than 29 m a−1. They revealed that the ice flow velocity had no significant changes over the period. Combining the surface mass balance estimate with the ice discharge estimated from our historical velocity maps and recently published velocity maps, we estimated a positive mass balance of 8.6 ± 3.9 Gt a−1 in the study area from 1963 and 2015. Our results indicate that the region's positive mass balance, as estimated in recently published studies, has been maintained since the 1960s. It is also in concordance with the low level of mass balance from 1992 to 2017 in East Antarctica. This suggests that the study area has been stable since the 1960s.

Long-Term Monitoring and Change Analysis of Pine Island Ice Shelf Based on Multi-Source Satellite Observations during 1973–2020

Pine Island Glacier (PIG) is one of the largest contributors to sea level rise in Antarctica. Continuous thinning and frequent calving imply significant destabilization of Pine Island Glacier Ice Shelf (PIGIS). To understand the mechanism of its accelerated disintegration and its future development, we conducted a long-term monitoring and comprehensive analysis of PIGIS, including ice flow velocity, ice shelf fronts, ocean water temperature, rifts, and surface strain rates, based on multi-source satellite observations during 1973–2020. The results reveal that: (1) ice flow velocities of PIGIS increased from 2.3 km/yr in 1973 to 4.5 km/yr in 2020, with two rapid acceleration periods of 1995–2009 and 2017–2020, and its change was highly correlated to the ocean water temperature variation. (2) At least 13 calving events occurred during 1973–2020, with four unprecedented successive retreats in 2015, 2017, 2018, and 2020. (3) The acceleration of ice shelf rifting and calving may correlate to the destruction of shear margins, while this damage was likely a response to the warming of bottom seawater. The weakening southern shear margin may continue to recede, indicating that the instability of PIGIS will continue.

ANALYSIS OF OVERESTIMATION IN HISTORICAL ICE FLOW VELOCITY MAPS IN WESTERN PACIFIC OCEAN SECTOR, ANTARCTICA

Abstract. High accuracy reconstruction of historical ice flow velocity fields is crucial in mass balance research of the Antarctic Ice Sheet by utilizing the input-output approach. A historical flow velocity of the Western Pacific Ocean sector region of East Antarctica from 1963 to 1989 was mapped and then corrected for its velocity overestimation. In this study, we analyzed the spatial distribution of the corrected values, and further assessed the relationship between the corrected values and related factors including timespan of image pairs, ice flow velocity, the spatial acceleration of ice flow velocity, and surface slope. The results indicate that the corrected ice flow velocity points are mostly dispersed between a buffer 25 km upstream and a buffer 25 km downstream the grounding line, with the largest mean value emerging in the region between the grounding line and its 25 km downstream. The corrected values exhibit linear correlation with three conditions: 1) when the ice flow velocity range is 0 – 1586 m/y, 2) spatial acceleration is 0 – 69 (m/y)/km, and 3) the slope is 0 – 1.95 degrees and the R2 is higher than 0.7. However, the correlation between timespan and corrected values is not obvious. The corrected values for the floating region have a greater linear correlation with all three factors than the inland region. Ice flow velocity, spatial acceleration, and surface slope all have an influence on the distribution of the corrected values of the reconstructed historical ice flow velocity maps, and may further affect the assessment of the mass balance of the Antarctic Ice Sheet.

ICE FLOW VELOCITY MAPPING IN GREENLAND USING HISTORICAL IMAGES FROM 1960s TO 1980s: SCHEME DESIGN

Abstract. Ice flow velocity over long time series in Greenland plays a vital role in estimating and predicting the mass balance of the Greenland Ice Sheet and its contribution to global sea level rise. However, there are few Greenland ice flow velocity products with large spatial coverage available showing the Greenland ice flow velocity pattern before the 1990s. We proposed three methods, including parallax decomposition, grid-based NCC image matching, feature and gird-based image matching with constraints for estimation of surface velocity in East Antarctica based on ARGON KH-5 and LANDSAT imagery, and a systematic compilation method for the ice surface velocity in East Antarctica from the 1960s to 1980s. Based on the above methods and cartographic processes, this study designed a framework for the mapping of the historical ice flow velocity of the Greenland ice sheet. Currently, the early LANDSAT images covering several glaciers in North Greenland and Northwest Greenland have been processed and applied to velocity mapping using the cartographic process proposed in this study, and some preliminary results have been obtained.

EXTRACTION OF ICE FLOW VELOCITY BY COMBINATION OF DINSAR AND OFFSET TRACKING METHODS FOR PINE GLACIER, WEST ANTARCTICA

Abstract. West Antarctica is the main contributor to global sea level rise at present and in the coming decades, since it occupies 80% of mass loss in Antarctica. In this paper, a combination of Differential Interferometric Synthetic Aperture Radar (DInSAR) and offset tracking technology is used to extract the ice flow velocity of the Pine Island Glacier (PIG), a typical glacier in West Antarctica. Due to the large deformation gradients in PIG, DInSAR technology is used to extract ice flow velocity in bare rock and mountains, then offset tracking technology is used to extract ice flow velocity in areas where glaciers collapse frequently. Finally, the above two results are mosaiced into a new image of the interannual ice flow velocity of PIG in 2017. Through qualitative and quantitative evaluation, it is found that the ice flow velocity extracted by the combination has high accuracy in both high and low velocity areas. In summary, we concluded that the combination of DInSAR and Offset tracking can obtain reliable ice velocity products in glaciers that change rapidly. This combination is of scientific significance for monitoring the movement and evolution of glaciers in the West Antarctica.

Subglacial meltwater routes of the Fennoscandian Ice Sheet

ABSTRACT Subglacial drainage systems are crucial elements of glaciers and ice sheets because they modulate ice flow velocity. However, logistical challenges of measuring subglacial processes beneath contemporary ice and natural limitations in long-term monitoring hinder our understanding about their spatio-temporal evolution. Subglacial meltwater landforms created by palaeo-ice sheets are records of past subglacial drainage systems and offer the potential to study their large-scale development throughout deglaciation. Although collectively recording subglacial drainage, individual meltwater landforms such as eskers, meltwater channels and meltwater corridors, which comprise tunnel valleys and meltwater tracks (assemblages of landforms in broad, elongated paths with irregular surface texture), have mostly been investigated as separate entities. Using high-resolution (1–2 m) digital elevation models, we map integrated networks of subglacial meltwater landforms, herein called subglacial meltwater routes, on an ice-sheet scale in Fennoscandia. Our map provides a basis for future research on the long-term evolution of subglacial drainage networks and its effect on ice dynamics of the Fennoscandian Ice Sheet.

Overestimation and adjustment of Antarctic ice flow velocity fields reconstructed from historical satellite imagery

Abstract. Antarctic ice velocity maps describe the ice flow dynamics of the ice sheet and are one of the primary components used to estimate the Antarctic mass balance and contribution to global sea level changes. In comparison to velocity maps derived from recent satellite images of monthly to weekly time spans, historical maps, from before the 1990s, generally cover longer time spans, e.g., over 10 years, due to the scarce spatial and temporal coverage of earlier satellite image data. We found velocity overestimations (OEs) in such long-span maps that can be mainly attributed to velocity gradients and time span of the images used. In general, they are less significant in slow-flowing grounded regions with low spatial accelerations. Instead, they take effect in places of high ice dynamics, for example, near grounding lines and often in ice shelf fronts. Velocities in these areas are important for estimating ice sheet mass balance and analyzing ice shelf instability. We propose an innovative Lagrangian velocity-based method for OE correction without the use of field observations or additional image data. The method is validated by using a set of ground truth velocity maps for the Totten Glacier and Pine Island Glacier which are produced from high-quality Landsat 8 images from 2013 to 2020. Subsequently, the validated method is applied to a historical velocity map of the David Glacier region from images from 1972–1989 acquired during Landsat 1, 4, and 5 satellite missions. It is demonstrated that velocity overestimations of up to 39 m a−1 for David Glacier and 195 m a−1 for Pine Island Glacier can be effectively corrected. Furthermore, temporal acceleration information, e.g., on basal melting and calving activities, is preserved in the corrected velocity maps and can be used for long-term ice flow dynamics analysis. Our experiment results in the Pine Island Glacier (PIG) show that OEs of a 15-year span can reach up to 1300 m a−1 along the grounding line and cause an overestimated grounding line (GL) flux of 11.5 Gt a−1 if not corrected. The magnitudes of the OEs contained in both velocity and mass balance estimates are significant. When used alongside recent velocity maps of 1990s–2010s, they may lead to underestimated long-term changes for assessment and forecast modeling of the global climate change impact on the Antarctic ice sheet. Therefore, the OEs in the long-span historical maps must be seriously examined and corrected. We recommend that overestimations of more than the velocity mapping uncertainty (1σ) be corrected. This velocity overestimation correction method can be applied to the production of regional and ice-sheet-wide historical velocity maps from long-term satellite images.

Controls on Larsen C Ice Shelf Retreat From a 60‐Year Satellite Data Record

AbstractRapid retreat of the Larsen A and B ice shelves has provided important clues about the ice shelf destabilization processes. The Larsen C Ice Shelf, the largest remaining ice shelf on the Antarctic Peninsula, may also be vulnerable to future collapse in a warming climate. Here, we utilize multisource satellite images collected over 1963–2020 to derive multidecadal time series of ice front, flow velocities, and critical rift features over Larsen C, with the aim of understanding the controls on its retreat. We complement these observations with modeling experiments using the Ice‐sheet and Sea‐level System Model to examine how front geometry conditions and mechanical weakening due to rifts affect ice shelf dynamics. Over the past six decades, Larsen C lost over 20% of its area, dominated by rift‐induced tabular iceberg calving. The Bawden Ice Rise and Gipps Ice Rise are critical areas for rift formation, through their impact on the longitudinal deviatoric stress field. Mechanical weakening around Gipps Ice Rise is found to be an important control on localized flow acceleration and the propagation of two rifts that caused a major calving event in 2017. Capturing the time‐varying effects of rifts on ice rigidity in ice shelf models is essential for making realistic predictions of ice shelf flow dynamics and instability. In the context of the Larsen A and Larsen B collapses, we infer a chronology of destabilization processes for embayment‐confined ice shelves, which provides a useful framework for understanding the historical and future destabilization of Antarctic ice shelves.

Retrieval of Svalbard ice flow velocities using Sentinel 1A/1B three-pass Differential SAR Interferometry

Glacier velocity is an important parameter to characterize glacier dynamics and to derive ice thickness and mass balance. The 2-pass/3-pass Differential Synthetic Aperture Radar Interferometry (DInSAR) techniques are advantageous in estimating glacier movements. However, the 2-pass DInSAR requires an external Digital Elevation Model (DEM), whereas the 3-pass DInSAR does not. The 3-pass DInSAR technique is adopted with Sentinel-1A/1B to map the Svalbard glacier flow velocities. Furthermore, the effect of glacier surface elevation change on the 2-pass DInSAR results were revealed by comparing glacier velocity with the two time period topographic information. The coherence indicates that Sentinel-1A/1B has a high potential to infer operational glacier velocity using the 3-pass DInSAR method with an average atmospheric uncertainty of 0.24 cm/day over the studied region. The precision of the 3-pass DInSAR is analyzed by comparing 2-pass derived line-of-sight (LOS) velocity on the same dates. The 3-pass DInSAR achieves high-resolution and detailed information.

Validation for Ice Flow Velocity Variations of Shirase Glacier Derived From PALSAR-2 Offset Tracking

We have shown that the flow velocity of Shirase Glacier can be estimated by applying the offset tracking method to the amplitude images derived from the phased array type L-band synthetic aperture radar type-2 (PALSAR-2) onboard the Advanced Land Observing Satellite-2. Although the offset tracking is widely used, the method is sometimes not applicable effectively and the reason is not always clear. Here, we consider a possibility of ionospheric contamination and show the validation results of the estimated flow velocity applied the offset tracking method to the PALSAR-2 data, based on <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in-situ</i> measurement of the flow velocity from the GNSS receiver on the autonomous phase-sensitive radio echo sounder simultaneously with the PALSAR-2 observation. As a result, the ionospheric contamination can yield a spurious error of ±0.2 km a <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">−1</sup> . After statistically removing the erroneous pairs, the RMSE was 0.049 km a <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">−1</sup> derived from the flow velocity error between the <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in-situ</i> measurement results and estimated from satellite, and the estimated flow velocity obtained using PALSAR-2 data is proved to be in good agreement with the ground truth. Obtained spatial and temporal variations reveal signature of glacier dynamics, which prove the efficacy of the offset tracking method. The obtained ice-flow velocity increases rapidly from the upstream region to the coast, but its velocity is roughly constant over a region, 10-km long about the grounding line (GL), then gradually tends to increase again downstream from the GL. This trend has continued largely unchanged over24 years since 1996.

Mid-Term Monitoring of Glacier’s Variations with UAVs: The Example of the Belvedere Glacier

Recently, Unmanned Aerial Vehicles (UAV) have opened up unparalleled opportunities for alpine glacier monitoring, as they allow for reconstructing extensive and high-resolution 3D models. In order to evaluate annual ice flow velocities and volume variations, six yearly measurements were carried out between 2015 and 2020 on the debris-covered Belvedere Glacier (Anzasca Valley, Italian Alps) with low-cost fixed-wing UAVs and quadcopters. Every year, ground control points and check points were measured with GNSS. Images acquired from UAV were processed with Structure-from-Motion and Multi-View Stereo algorithms to build photogrammetric models, orthophotos and digital surface models, with decimetric accuracy. Annual glacier velocities were derived by combining manually-tracked features on orthophotos with GNSS measurements. Velocities ranging between 17 m y−1 and 22 m y−1 were found in the central part of the glacier, whereas values between 2 m y−1 and 7 m y−1 were found in the accumulation area and at the glacier terminus. Between 2 × 106 m3 and 3.5 × 106 m3 of ice volume were lost every year. A pair of intra-year measurements (October 2017–July 2018) highlighted that winter and spring volume reduction was ∼1/4 of the average annual ice loss. The Belvedere monitoring activity proved that decimetric-accurate glacier models can be derived with low-cost UAVs and photogrammetry, limiting in-situ operations. Moreover, UAVs require minimal data acquisition costs and allow for great surveying flexibility, compared to traditional techniques. Information about annual flow velocities and ice volume variations of the Belvedere Glacier may have great value for further understanding glacier dynamics, compute mass balances, or it might be used as input for glacier flow modelling.

Modeling the Dynamics of Supraglacial Rivers and Distributed Meltwater Flow With the Subaerial Drainage System (SaDS) Model

AbstractMeltwater produced at the surface of glaciers and ice sheets has important implications for basal sliding rates and therefore ice flow velocities. In order to determine the role of supraglacial water in ice dynamics and predict future changes, we first need to understand and be able to accurately predict moulin input rates. To this end, we present the Subaerial Drainage System (SaDS) model. SaDS is a dynamic model that couples supraglacial runoff in the bare‐ice ablation zone in a distributed sheet with flow in discrete channels. Flow in the distributed sheet drives melt through potential energy dissipation, allowing a channel network to form naturally with no prior assumptions about channel locations. We apply the model to a synthetic ice sheet margin and carry out a suite of sensitivity tests. Modeled moulin inputs show expected behaviors including large diurnal variability, multi‐hour lags following peak surface melt, and demonstrate complex and diverse seasonal dynamics. The sensitivity tests illustrate the range of possible model behaviors and constrain the parameter values for which the model predicts physically realistic moulin inputs. We also apply the model to a ∼20 × 27 km2 catchment on the southwestern Greenland Ice Sheet using RACMO melt forcing and previously mapped moulin locations. Modeled supraglacial lake and stream locations match those mapped from Landsat 8 images, and moulin inputs show varied daily and seasonal dynamics. These results demonstrate that the model is a promising tool to provide moulin inputs for subglacial and ice dynamic studies.

Complementary Approaches Towards a Universal Model of Glacier Surges

Although many convincing, diverse, and sometimes competing models of glacier surging have been proposed, the observed behavior of surging glaciers does not fit into distinct categories, and suggests the presence of a universal mechanism driving all surges. On the one hand, recent simulations of oscillatory flow behavior through the description of transient basal drag hint at a fundamental underlying process. On the other hand, the proposition of a unified model of oscillatory flow through the concept of enthalpy adopts a systems based view, in an attempt to rather unify different mechanisms through a single universal measure. While these two general approaches differ in perspective, they are not mutually exclusive, and seem likely to complement each other. A framework incorporating both approaches would see the mechanics of basal drag describing ice flow velocity and surge propagation as a function of forcing by conditions at the glacier bed, in turn modulated through the unified measure of enthalpy.