Glacier Flow Research Articles

Multidimensional Small Baseline Subset (Msbas) Software for Constrained and Unconstrained Deformation Analysis of Partially Coherent Dinsar and Speckle Offset Data

Differential Interferometric Synthetic Aperture Radar (DInSAR) and speckle offset data are used to measure ground deformation. Computing deformation time series from these data can be challenging due to varying radar line-of-sight acquisition geometry, rapid coherence loss, and the large datasets involved requiring computationally intensive processing. The Multidimensional Small Baseline Subset (MSBAS) Software is designed to compute deformation time series from these data, and its parallelized version is optimized to handle large datasets from modern SAR satellites. The software can produce 1D, 2D, and 3D Surface-Parallel and Aspect-Parallel Flow constrained and 3D unconstrained or 4D deformation time series. It runs on workstations and clusters utilizing OpenMP and MPI/OpenMP. This study presents novel results demonstrating MSBAS capabilities, such as landslides in Northwest Territories, Argentina, and Colombia, as well as tectonic deformation and subsidence in Turkey and Yemen. A 4D deformation time series of glacier flow and surge at the Malaspina Glacier in Alaska is showcased. Equations for solving multidimensional problems are included along with processing parameters and supplemental data used to create constraints. The manuscript serves as a comprehensive manual for the MSBAS Software Version 10. The source code is available.

Quantifying the buttressing contribution of landfast sea ice and melange to Crane Glacier, Antarctic Peninsula

Abstract. The January 2022 disintegration of multi-year landfast sea ice in the Larsen B embayment, Antarctic Peninsula, was closely followed by a significant acceleration of ice flow and ice-front retreat of numerous outlet glaciers. Crane Glacier was a notable example of this, with 6 km of its floating ice shelf lost to calving in the first month following the disintegration and a 3.4 % increase in terminus flow speeds over the same time period. In this study we quantify for the first time the buttressing stresses that were transmitted to Crane by the ice melange at the glacier outlet using the ice-flow model Úa. We constrained our model with high-resolution surface elevation profiles of the glacier and ambient melange and reconstructed the observed flow velocities by optimising the rheology rate factor throughout our model domain. This allowed us to quantify the stress regime across both the glacier and ice melange. Results showed that resistive backstresses were imparted to Crane by the ice melange with a mean buttressing ratio of ΘN=0.68 calculated at the glacier terminus (ΘN=1 implies no buttressing). In addition, diagnostic modelling showed an expected 19.2 kPa mean increase in extensional stress at the ice front following the disintegration of the ice melange. This perturbation in stress likely triggered the observed rapid calving over the near-terminus region, leading to the periodic loss of sections of Crane's buttressing ice shelf and thus further acceleration of ice flow in the subsequent months.

Impacts of GCP Distributions on UAV-PPK Photogrammetry at Sermeq Avannarleq Glacier, Greenland

Real-Time/Post-Processing Kinematic (RTK/PPK) technology has been widely applied in Unmanned Aerial Vehicle (UAV) photogrammetry in glaciological research. Considering that ground control points (GCPs) cannot be set on glaciers, evaluating the impacts of one-sided distribution is essential. In this study, 8571 images were captured at Sermeq Avannarleq glacier in western Greenland from 4 August 2021 to the 6th, covering approximately 85 km2, with the furthest distance being 13.22 km away from the coastline. Benefited by the meandering coastline, 11 roving stations roughly uniformly distributed on bare rock were surveyed with the RTK technique. PPK-geotagged images were processed in Agisoft Metashape Professional to derive the DSMs, utilizing eight different configurations of GCP distributions that gradually extended longitudinally (along the glacier flow direction) to the upper part of the glacier. The accuracy of DSMs was evaluated by referring to the validation points (VPs) that were not employed in the Bundle Block Adjustment (BBA). The results indicated that the RMSE values of the easting, northing, and height of the reconstruction model georeferenced by only PPK geotagging (no GCPs applied) were 0.038 m, 0.031 m, and 0.146 m, respectively. Applying four GCPs located at one side of the region but with both longitudinal and lateral distribution improved the RMSE values in easting, northing, and vertical to 0.037 m, 0.031 m, and 0.081 m, respectively, and these values were stable even when distributing four GCPs evenly or when increasing the number of GCPs to eleven. Moreover, the cross-validation with ICESat-2 and ArcticDEM performed only at an off-glacier region also suggested that vertical accuracy shows significant improvements for every configuration of GCPs compared to the reconstruction model optimized only by PPK, but such improvements were not obvious if the number of GCPs exceeded four. Moreover, no elevation ramps appeared in the UAV DSM, even for the GCP configuration with only two GCPs distributed at the terminus. Therefore, combining PPK with only a few GCPs but distributing in both directions of the surveying region can offer a viable solution for obtaining glacier DSMs at the coastline with decimeter-level accuracy.

Glacier Surges Controlled by the Close Interplay Between Subglacial Friction and Drainage

AbstractThe flow of glaciers and ice sheets is largely influenced by friction at the ice‐bed interface that can trigger rapid changes in glacier motion ranging from seasonal velocity variations to cyclic surge instabilities or even devastating glacier collapse. This wide range of transient glacier dynamics is currently not captured by models, and its implications for long‐term glacier evolution are uncertain. This highlights the need of developing improved descriptions for processes that occur at the glacier bed. Here, we present a model that describes the evolution of basal friction inspired by a “rate and state” approach, coupled to models of subglacial drainage and glacier flow, and investigate how these couplings affect the dynamics of glaciers. We show that a wide range of sliding behavior results from a feedback loop between subglacial drainage efficiency and friction which depends on the evolution of a frictional state that can be interpreted as the degree of cavitation or till porosity for hard and soft beds, respectively. In our simulations, we find that glaciers are susceptible to surging if they exhibit a transition to velocity weakening friction associated with a poor sensitivity of the drainage capacity to the frictional state. This potential materializes if the local topography and mass balance create the conditions for high water pressure to build up in an area sufficiently large to exceed a critical length. We advocate accounting for feedback loops between friction and drainage as a promising avenue for better understanding dynamical instabilities of glaciers and ice sheets.

Glacier dynamic study of the Terre de Lambert glacier in Greenland based on Landsat 8 and Sentinel-1 data

With global temperatures on the rise, the trend of glacial melting is becoming increasingly evident. This study focuses on the northeastern region of Greenland as the research subject, utilizing remote sensing technology and geographic information technology to study glacial dynamics. Initially, the research processes Landsat 8 data using QGIS to classify the widespread snow and ice cover, facilitating the observation and analysis of its temporal variations. The study then narrows its focus to Terre de Lambert for further analysis. Subsequently, using satellite data from Sentinel-1 and employing SNAP, the research conducts an in-depth analysis of glacier flow speed in specific regions. Through intricate computational models and algorithms, precise glacier flow direction and speed maps are generated. By combining historical data, the study compares and analyzes the glacier flow speed between 2018 and from June 20th to July 14th, 2023. Results indicate that the glacier flow speed in 2023 during the same time frame was faster. Throughout the research process, limitations in current data acquisition and processing methods were identified. Thus, the paper suggests improvements to data acquisition and processing tools to enhance the efficiency of glacier flow speed studies. Theres a fervent hope for the emergence of more convenient, streamlined, and customizable temporal dimension tools for glacier flow speed data in the future. The research findings aim to draw global citizens attention to the issue of glacial melting and inspire collective actions to protect our planet.

Spatial Distribution and Variation in Debris Cover and Flow Velocities of Glaciers during 1989–2022 in Tomur Peak Region, Tianshan Mountains

In this study, we utilized a feature optimization method combining texture and topographical factors with the random forest (RF) approach to identify changes in the extent of the debris cover around the Tianshan Tomur Peak between 1989 and 2022. Based on Sentinel-1 image data, we extracted glacier flow velocities using an offset tracking method and conducted a long-term analysis of flow velocities in combination with existing datasets. The debris identification results for 2022 showed that the debris-covered area in the study region was 409.2 km2, constituting 22.8% of the total glacier area. Over 34 years, the area of debris cover expanded by 69.4 km2, reflecting a growth rate of 20.0%. Analysis revealed that glaciers in the Tomur Peak area have been decelerating at an overall rate of −4.0% per decade, with the complexity of the glacier bed environment and the instability of the glacier’s internal structure contributing to significant seasonal and interannual variability in the movement speeds of individual glaciers.

Multi-scale variations of subglacial hydro-mechanical conditions at Kongsvegen glacier, Svalbard

Abstract. The flow of glaciers is largely controlled by changes at the ice–bed interface, where basal slip and sediment deformation drive basal glacier motion. Determining subglacial conditions and their responses to hydraulic forcing remains challenging due to the difficulty of accessing the glacier bed. Here, we monitor the interplay between surface runoff and hydro-mechanical conditions at the base of the Kongsvegen glacier in Svalbard. From July 2021 to August 2022, we measured both subglacial water pressure and till strength. Additionally, we derived median values of subglacial hydraulic gradient and radius of channelized subglacial drainage system from seismic power, recorded at the glacier surface. To characterize the variations in the subglacial conditions caused by changes in surface runoff, we investigate the variations of the following hydro-mechanical properties: measured water pressure, measured sediment ploughing forces, and derived hydraulic gradient and radius, over seasonal, multi-day, and diurnal timescales. We discuss our results in light of existing theories of subglacial hydrology and till mechanics to describe subglacial conditions. We find that during the short, low-melt-rate season in 2021, the subglacial drainage system evolved at equilibrium with runoff, increasing its capacity as the melt season progressed. In contrast, during the long and high-melt-rate season in 2022, the subglacial drainage system evolved transiently to respond to the abrupt and large water supply. We suggest that in the latter configuration, the drainage capacity of the preferential drainage axis was exceeded, promoting the expansion of hydraulically connected regions and local weakening of ice–bed coupling and, hence, enhanced sliding.

Glacial lakes control ice flow: new insights from satellite SAR PO-SBAS observations in Duiya Glacier, southern Tibetan Plateau

ABSTRACT This study investigates the effects of glacial lakes on the flow velocities of Duiya Glacier, southern Tibetan Plateau, by employing satellite Synthetic Aperture Radar (SAR) Pixel-Offset-Tracking Small-Baseline-Subset (PO-SBAS) technology. Duiya Glacier, which terminates in a lake, exhibits approximately 10 m/yr higher horizontal speed and a higher melting rate with a vertical Non-Surface-Parallel-Flow (nSPF) velocity of – 7 m/yr than those of its land-terminating counterpart – West Duosangpu Glacier. Notably, Duiya Glacier experiences significant seasonal velocity fluctuations with accelerated flow between May and October. By integrating glacier geometry, changes in glacier boundaries and the extent of proglacial lakes, the distribution of supraglacial lakes, and climatic variables, we reveal that proglacial and supraglacial lakes play a crucial role in increasing the flow velocities of Duiya Glacier. Duiya Glacier flows faster because of increased subglacial water pressure resulting from water influx from these lakes. This phenomenon becomes conspicuously evident during May – October, when increased meltwater due to increased temperatures and precipitation further elevates the subglacial water pressure. Our method highlights the potential for understanding the impact of glacial lakes on glacier movements at a large scale, leveraging the capabilities of satellite SAR PO-SBAS technology for continuous, wide-scale, high temporal resolution, and 3D velocity monitoring.

Submerged imprint of glacier dynamics in the NW sector of Lago Viedma (southern Patagonia, Argentina)

ABSTRACTLago Viedma is a proglacial lake into which the Viedma Glacier flows from the Southern Patagonian Ice Field. This glacier has rapidly lost mass in recent times, and its deglacial history is reflected in the lake's subsurface. New high‐resolution multi‐channel seismic profiles acquired in the northwestern sector of the lake have allowed us to reconstruct the bathymetry of the area and identify several small sub‐basins, which have a maximum depth of 240 m in this sector. Four seismic facies have been recognized, separated by erosional unconformities, reflecting the depositional conditions in this sector of the basin during different Quaternary sedimentation phases. There is a transitional phase from ice‐contacted deposits to subglacial deposits, probably associated with a subglacial fan, and finally a phase of lacustrine sedimentation. In addition, three depositional stages were identified within the lacustrine deposits, indicating a different sedimentary input, and the morphologies observed at the lake bottom suggest the existence of at least ten stagnations of the glacial margin. This study helps to improve the history of the retreat of the Viedma Glacier in this area during the Pleistocene/Holocene and provides a basis for further geophysical measurements aimed at mapping this remote lacustrine environment.

A modified viscous flow law for natural glacier ice: Scaling from laboratories to ice sheets

Glacier flow modulates sea level and is governed largely by the viscous deformation of ice. Multiple molecular-scale mechanisms facilitate viscous deformation, but it remains unclear how each contributes to glacier-scale deformation. Here, we present a model of ice deformation that bridges laboratory and glacier scales, unifies existing estimates of the viscous parameters, and provides a framework for estimating the parameters from observations and incorporating flow laws derived from laboratory observations into glacier-flow models. Our results yield a map of the dominant deformation mechanisms in the Antarctic Ice Sheet, showing that, contrary to long-standing assumptions, dislocation creep, characterized by a value of the stress exponent [Formula: see text], likely dominates in all fast-flowing areas. This increase from the canonical value of [Formula: see text] dramatically alters the climate conditions under which marine ice sheets may become unstable and drive rapid rates of sea-level rise.

Dynamics of the spatiotemporal velocity of glaciers on the eastern slope of Mount Gongga, China, under climate change

The quantitative assessment of glacier flow velocity dynamics plays a pivotal role in understanding its response mechanisms concerning climate warming. This work provides a systematic quantitative assessment of the deceleration status of glaciers in this region by investigating the motion evolution of typical glaciers in Mount Gongga in recent years, thereby revealing the seasonal dynamics and inter-annual evolution over an extensive time span. We used the optical flow-small baseline subset (OF-SBAS) method to compute the time-series velocities of the Hailuogou Glacier and the Mozigou Glacier using 178 archived Sentinel-1 satellite synthetic aperture radar (SAR) images from 2014 to 2021. The findings revealed a prominent seasonal pattern in glacier motion, characterised by cyclic variations in velocity from cold to warm seasons. Moreover, we identified variations in velocities across distinct regions of the glacier surface, underscored by the lag in the peak time node of glacier flow with increasing elevation. This pattern may have been determined by a combination of internal and external factors, including mass accumulation and ablation-driven subglacial drainage, as well as the glacier geomorphological setting. Furthermore, during 2015–2021, the glaciers on the eastern slope of Mount Gongga exhibited an overarching trend of deceleration. Notably, the ablation area of the Hailuogou Glacier recorded the most substantial deceleration, exceeding 8% per year. This study underscores the efficacy of the OF-SBAS method in extracting long-term glacier velocities. This work also establishes a robust foundation for the analysis of spatiotemporal fluctuations in glacier movement within the context of climate warming.

Hydrological control of the surging behaviour of the Ghujerab River Head Glacier, Karakoram (2019–2023): Insights from high-temporal-resolution remote sensing monitoring

Study regionThe northern region of the Karakoram Range. Study focusKarakoram is a region in High Mountain Asia with many surge-type glaciers. This study employed over 200 high-temporal-resolution remote sensing images and investigated the variations in elevation and velocity of the Ghujerab River Head Glacier (GRHG) from 2019 to 2023. Furthermore, we elucidated the potential controlling mechanisms. New hydrological insights for the regionOur findings revealed that the GRHG, akin to typical surge-type glaciers in Karakoram, started to surge in the spring and finished surging in the summer, with a duration of less than two years. Throughout the surging process, the glacier transferred a mass of 0.11±0.003 km³ from the reservoir area to the receiving area, resulting in a thickening of 91.59±1.04 m at the glacier terminus and thinning of 11.78±1.04 m in the upper glacier. By analysing the mass balance and glacier surface albedo during surging, we proposed that climatic disturbances in the glacier region provided essential material inputs for the surge. Additionally, based on the seasonal evolution pattern of glacier flow velocity, we inferred a close correlation between surging and variations in subglacial hydrology. The duration of acceleration and deceleration during glacier surging, as well as a comparison with existing studies, further support our conclusion. Future research integrating multi-source remote sensing and onsite observations can support numerical simulations to quantitatively reveal the key processes occurring beneath and within glaciers during surge events.

Watershed water policy based on a glacier water service perspective in the typical glacier basins, Western China

Abstract Taking the Yangtze River Source Basin (YRSB) and Shule River Basin (SRB) as two typical cases, the sustainability of the water resources in these two basins was evaluated using the level of water stress (LWS) from sustainable development goal 6.4.2, and the regulating effect of the glacier runoff on the LWS was quantified. From 2000 to 2030, the level of socioeconomic development in the YRSB is low, and the total water consumption is only about 0.18 × 108 m3, whereas the SRB has a relatively high level of socioeconomic development and total water consumption is about 10 × 108 m3, i.e., 50 times higher than that in the YRSB. For the aforementioned reasons, the SRB's LWS is much higher than the YRSB's, resulting in a very low sustainability of water resources. As natural assets, glaciers flow downstream in the runoff mode, so compensation at the watershed scale should be considered. In the basin, the optimal allocation of water resources is needed. At the inter-basin scale, the compensation mechanism of glacier water resources needs to be improved.

The use of terrestrial monoscopic time-lapse cameras for surveying glacier flow velocity

Glacier flow velocity is crucial information to assess the impact of global warming on glaciers since it is related to the ice thickness, its variations provide information on mass balance and, in general, on the current state of glacier “health”. Moreover, flow anomalies are often an indicator of glacier instabilities. Therefore, more attention has been dedicated to surveying glacier kinematics. Historically, flow velocity was the first quantitative variable measured on glaciers since the 19th century. In the last decades, terrestrial monoscopic time-lapse digital cameras have permitted to conduct automatic surveying at high spatial and temporal resolutions using digital image correlation. Even though terrestrial time-lapse imagery is currently a consolidated technique in glacier monitoring, possible strategies, limitations and potentialities have never been systematically reviewed. This work aims to illustrate the typical procedures required to monitor glacier flow velocity using terrestrial monoscopic time-lapse cameras. We describe possible inconveniences that can arise during the survey and provide some guidelines to minimise such issues. We present six study cases in the European Alps (Mont Blanc and Bernina massifs) that feature different monitoring equipment, site geometry and glacier morphodynamics to illustrate possible solutions for terrestrial imagery monitoring. The results of this review highlight the high benefit-to-cost ratio of terrestrial time-lapse cameras in glacier flow surveying.

Link between glacial striation morphology and induced drag

Abrasion acts to smooth glacial terrains and leaves behind linear scratch-like features (striations) on bedrock landscapes. Striations are often used as measures of glacier flow directions, but their morphology can also provide information about the subglacial stress conditions that produced the features. While striations are often abundant in the field, the processes that create them can be opaque and hard to examine in situ because they occur under thick layers of flowing ice. To alleviate that difficulty and provide information for interpretation of the populations of striations that are observed in the field, we conducted a set of laboratory experiments in which a ring of temperate debris-laden ice was slid atop a planar marble bed under various contact force conditions that led to the creation of hundreds of striations. During the experiment, numerous glaciological properties were continuously measured, including the resistive drag. Following the completion of the experiments, the marble beds were extracted, and the striations were measured for length and categorized by morphological type, and a subset was measured using a high-resolution white-light profilometer. These experiments showed that, similar to field observations, type 2 striations were initially the most abundant; however, we found that type 3 striations became the most abundant at large displacements. We found good correlation between the abundance of striations as a function of displacement and measured drag as a function of displacement. When taken together, these results suggest that, in natural settings, ice flow around small roughness elements in glacier beds can “reset” the basal debris field, causing striations to become more abundant in their wake. As roughness is linked to quarrying, abrasion rates may increase in areas of increased quarrying.

Accumulation rates over the past 260 years archived in Elbrus ice core, Caucasus

Abstract. In this study, we present a seasonally resolved accumulation record spanning from 1750 to 2009 Common Era (CE), based on a 181.8 m ice core obtained from the Elbrus Western Plateau in the Caucasus. We implemented various methods to account for uncertainties associated with glacier flow, layer thinning, and dating. Additionally, we applied a novel approach to calculate a seasonal calendar for meteorological data, enabling comparison with ice core records. The reconstructed accumulation data were compared with available meteorological data, gridded precipitation records, and paleo-reanalysis data. Reconstructed accumulation is representative for a large region south of the Eastern European plain and Black Sea region with summer precipitation being the primary driver of precipitation variability. We identified a statistically significant relationship between changes in regional precipitation and fluctuations in the North Atlantic Oscillation (NAO) index, which is, however, not stable over the entire period covered by the ice core.

Bayesian estimation of glacier surface elevation changes from DEMs

Accurate estimates of glacier surface elevation changes are paramount for various aspects of the study of the cryosphere, from glacier flow and thickness estimates to hydrological forecasts and projections of sea-level-rise. We present a novel probabilistic framework to filter outliers and estimate uncertainties in glacier surface elevation changes computed from the subtraction of digital elevation models (DEM). Our methodology frames outlier filtering as a Bayesian inference problem, thus characterizing the state of knowledge on glacier surface elevation changes through the posterior distribution as the combination of glacier volume variation observations and prior knowledge arising from previously collected data and/or modeled results. We validate this technique with experiments using Gaussian random fields to generate artificial noise in glacier surface elevation variation observations and show that the model satisfactorily culls the simulated outliers. Surface elevation change estimates are consistent with results computed from widely-used outlier filtering and uncertainty estimation techniques. The Bayesian framework allows unifying DEM error models with physical considerations on glacier surface elevation changes within a simple, statistically coherent model preventing temporal correlation and additional biases in other techniques. On the basis of these results, we discuss the implications of DEM uncertainty and offer suggestions for the glaciological community.

Application of a structured decision‐making process in cryospheric hazard planning: Case study of Bering Glacier surges on local state planning in Alaska

AbstractSurging glaciers are glaciers that experience rapidly accelerated glacier flow over a comparatively short period of time. Though relatively rare worldwide, Alaska is home to the largest number of surge‐type glaciers globally. However, their impact on the broader socioecological system in the state is both poorly understood and under‐researched, which poses a challenge in developing appropriate sustainability decisions in Alaska. We investigated how the surge patterns of the Bering Glacier in Alaska have potentially devastating effects on the local ecological biodiversity of its watershed via a structured decision‐making analysis of the different possible consequences. Specifically, this analysis was conducted to explore the various outcomes of a Bering Glacier surge particularly if humans have an increased presence near the glacier due to the area potentially becoming a state park. This work explored the benefits of applying a risk and decision analytical framework in a cryosphere context, to better understand the socioeconomic impact of glacier surges. This is a novel approach in which a decision analysis tool was used to better understand an environmental sustainability challenge, offering an innovative method to support the achievement of the United Nations Sustainability Development Goals in Alaska. We therefore emphasise the need for integrated biophysical and socioeconomic analyses when it comes to understanding glacier hazards. Our research highlights the importance of understanding and researching biophysical changes as well as using a structured decision‐making process for complicated hazard planning scenarios, exemplified via glaciated regions in Alaska, in order to create adaptation strategies that are sustainable and encompass the range of possible outcomes.

Melting of glacier ice enhanced by bursting air bubbles

Feedbacks between ice melt, glacier flow and ocean circulation can rapidly accelerate ice loss at tidewater glaciers and alter projections of sea-level rise. At the core of these projections is a model for ice melt that neglects the fact that glacier ice contains pressurized bubbles of air due to its formation from compressed snow. Current model estimates can underpredict glacier melt at termini outside the region influenced by the subglacial discharge plume by a factor of 10–100 compared with observations. Here we use laboratory-scale experiments and theoretical arguments to show that the bursting of pressurized bubbles from glacier ice could be a source of this discrepancy. These bubbles eject air into the seawater, delivering additional buoyancy and impulses of turbulent kinetic energy to the boundary layer, accelerating ice melt. We show that real glacier ice melts 2.25 times faster than clear bubble-free ice when driven by natural convection in a laboratory setting. We extend these results to the geophysical scale to show how bubble dynamics contribute to ice melt from tidewater glaciers. Consequently, these results could increase the accuracy of modelled predictions of ice loss to better constrain sea-level rise projections globally.

GLOBAL GEODYNAMIC MODEL OF THE EARTH AND ITS APPLICATION FOR ANTARCTICA

A geodynamic model of the modern Earth is constructed based on the SMEAN 2 global seismic tomography model. On the basis of mantle temperature anomalies, numerical modeling of spherical mantle convection is performed for the parameters of the real Earth. For Antarctica, an explanation was obtained for the existence of the West Antarctic Rift System, which contains one of the largest and least known volcanic provinces on Earth, which is consistent with the increased surface heat flow and modern volcanism of the West Antarctic Rift System and part of the Antarctic Peninsula. The increased heat flow and volcanoes in this region are causing instability and accelerating the flow of glaciers from the West Antarctic Ice Sheet into the ocean, which could lead to a significant rise in global sea level.