Glacier Ice Research Articles

Simulation of flexural-gravity wave propagation for elastic plates in shallow water using an energy-stable finite difference method with weakly enforced boundary and interface conditions

We introduce an energy stable, high-order-accurate finite difference approximation of the dynamic, pure bending Kirchhoff plate equations for complex geometries and spatially variable properties. We utilize the summation-by-parts (SBP) framework to discretize the biharmonic operator with variable coefficients, with attention given to free and clamped boundary conditions and corner conditions. Energy conservation is established by combining SBP boundary closures with weak enforcement of the boundary and interface conditions using a penalty (simultaneous approximation term, SAT) technique. Then we couple the plate equations to the shallow water equations to study flexural-gravity wave propagation, and prove that the semi-discrete system of equations is self-adjoint. We demonstrate the stability and accuracy properties of the method on curvilinear multiblock grids using the method of manufactured solutions. The method, which we provide in an open-source code, is then used to model ocean wave interactions with the Thwaites Glacier and Pine Island Ice Shelf in the Amundsen Sea off the coast of West Antarctica.

The Twannberg iron meteorite strewn field in the Swiss Jura mountains: insights for Quaternary environmental conditions

The ~ 10 km2 strewn field of the Twannberg type IIG iron meteorite is located in the Swiss Jura Mountains, 30 km northwest of Bern. The strewn field has been mapped by a group of citizen scientists since 2006, yielding more than 2000 meteorite fragments with a total mass of 152.7 kg until the end of 2022. With a terrestrial age of 176 ± 19 ka and a minimum pre-atmospheric mass of ~ 250 t, the Twannberg meteorite is a local time marker in an area with a poorly-known paleoenvironmental history. The Twannberg strewn field is located just outside of the maximum extent of ice during the Last Glacial Maximum (LGM). On the Mont Sujet, meteorites are size-sorted in a 6-km long section of the primary strewn field (altitude 945–1370 m a.s.l.), indicating a fall direction from east-northeast to west-southwest (azimuth approximately 250°). On the Twannberg plateau and in the Twannbach gorge, meteorites are not size-sorted and occur in a ~ 5.7-km long area associated with till and recent stream sediments (altitude 430–1075 m a.s.l.). The mass distribution of meteorites on the Twannberg plateau demonstrate that these meteorites were not found where they fell but that they must have been transported up to several km by glacier ice flow after the fall. The distribution of meteorites and of glacially transported Alpine clasts on the Mont Sujet and on the Chasseral chain indicates the presence of local ice caps and of an approximately 200-m higher Alpine ice surface with respect to the LGM at the time of fall. This high ice level during MIS 6 (Marine Isotopic Stage 6, 191–130 ka) indicated by the meteorite distribution is consistent with surface exposure ages of 50–144 ka from nearby resting erratic boulders at altitudes of up to 1290 m a.s.l., including the newly dated Jobert boulder (63 ka). These boulders indicate an ice level ~ 400 m higher than during LGM at a time not later than MIS 6. Post-LGM luminescence ages of loess-containing meteorites on the Mont Sujet and 14C ages of materials associated with meteorite finds indicate relatively young pedoturbation and increased oxidation of meteorites since ~ 7300 cal BP, possibly correlated with deforestation and enhanced erosion resulting from increased human activities since the Neolithic. This study shows that Twannberg meteorites in their palaeoenvironmental context provide valuable information about ice levels and transport directions during MIS 6 and about their interaction with the post-LGM environmental conditions. The unique Twannberg strewn field has the potential to reveal more valuable information.

Heterogeneity in glacio-hydrological processes and estimation of different components in streamflow from central Himalayan glaciers

Study regionThe study region includes two glaciers from Alaknanda and Bhagirathi river basins in Central Himalaya. Study focusThe study focuses on the analysis of high-resolution isotopic data sets of different components of the streamflow with ground-based meteorological observations from Automatic Weather Stations (AWSs) at two glaciers (∼4000 m asl). New hydrological insights for the regionThe glaciers in the Himalaya are difficult to access due to their topography and climate and require complicated logistics to work in the region. Hence, the understanding of hydrological processes in this region is limited. The debris cover (sublimation of surface ice), orientation and microclimatic conditions (temperature, wind and rain) of the two glaciers control the isotope signatures of the glacier surface ice in the Himalayan region, indicating heterogeneity and complexity in the isotopic compositions. The studies estimating the contribution of different components to the streamflow downstream using generalized values of stable isotopes (glacier ice, snow) are complicated, as several glaciers contribute to the total runoff in large basins. The stable isotopes of streamflow indicate the contribution of snow and ice melt during early ablation (May-Jun.); rainfall and ice melt during the ISM (Jul.-Aug.) and ice melt during late ablation (Sep.-Oct.). The contribution of snow-glacier melt and rainfall for the ablation season (Jun.-Oct.) was 89% and 11%, respectively. The separation of the hydrograph is complex, site and time-specific, which needs attention.

Estimation of glacier depth and ice volume of Kabul Basin, Afghanistan

Estimation of glacier depth and ice volume of Kabul Basin, Afghanistan

Glacier dynamics assessment in Eastern Dhauliganga basin (1994–2018), Kumaun Himalaya, India

AbstractThe glacier database provides an opportunity to track its dynamics and address issues linked with water resources, hazards and sustainability. To better understand the dynamics of Eastern Dhauliganga basin, we mapped the clean and debris‐covered glaciers by using satellite images from Landsat‐series and Sentinel 2. The band‐ratio based methods—NDSI (Normalised Difference Snow Index) and NDVI (Normalised Difference Vegetation Index) were used to map the clean ice glaciers. The debris‐covered glaciers were delineated using topographic parameters and thermal band of Landsat TM. The glacier extent and length change were monitored for 25 glaciers in the basin, of which 13 are clean ice and 12 are debris‐covered glaciers. During the time period from 1994 to 2018, higher area loss was observed for clean ice glaciers (±20.6%) compared with debris‐covered glaciers (±10.9%). The clean‐ice glaciers retreated at a rate of 4.8–33.7 m/year, whereas debris‐covered glaciers retreated at a rate of 2.9–28.8 m/year during 1994–2018.

Spatial differences of ice volume across High Mountain Asia

Advanced knowledge of glacier ice volume is vital for water resource assessment. Previous studies have focused on the estimation of ice volume, but the quantitative understanding of the spatial variability of ice volume across High Mountain regions is currently lacking. Here, we used global-scale ice thickness, debris cover and equilibrium line data to analyse ice-volume differences at various scales across High Mountain Asia (HMA). The results showed that 6.3% of the HMA glaciers are covered by debris, with debris area and volume accounting for 9% and 13.8% of the total glacier area and volume, respectively. An average debris-cover volume ratio of 13% was observed. The spatial distribution of ice volume across the HMA varies considerably from region to region. The ice volume is predominately distributed on north-facing slopes and accounts for approximately 38% of the total. It is very common in Altay and Sayan, East Tian Shan, West Kunlun, East Kunlun and Qilian Shan. Meanwhile, ice volumes in the Himalayas and Hengduan Shan are mainly distributed on the southeast aspect. Relative weight functions showed that glacier area, maximum length and average thickness are closely related to ice volume, with average relative weights of 63.7%, 22.5% and 9.8%, respectively. This study is important for the evolution of glacier volume and water resource assessment.

Basal melt rates and ocean circulation under the Ryder Glacier ice tongue and their response to climate warming: a high-resolution modelling study

Abstract. The oceanic forcing of basal melt under floating ice shelves in Greenland and Antarctica is one of the major sources of uncertainty in climate ice sheet modelling. We use a high-resolution, nonhydrostatic configuration of the Massachusetts Institute of Technology general circulation model (MITgcm) to investigate basal melt rates and melt-driven circulation in the Sherard Osborn Fjord under the floating tongue of Ryder Glacier, northwestern Greenland. The control model configuration, based on the first-ever observational survey by Ryder 2019 Expedition, yielded melt rates consistent with independent satellite estimates. A protocol of model sensitivity experiments quantified the response to oceanic thermal forcing due to warming Atlantic Water and to the buoyancy input from the subglacial discharge of surface fresh water. We found that the average basal melt rates show a nonlinear response to oceanic forcing in the lower range of ocean temperatures, while the response becomes indistinguishable from linear for higher ocean temperatures, which unifies the results from previous modelling studies of other marine-terminating glaciers. The melt rate response to subglacial discharge is sublinear, consistent with other studies. The melt rates and circulation below the ice tongue exhibit a spatial pattern that is determined by the ambient density stratification.

Glacier–rock glacier interactions in the eastern Hindu Kush, Nuristan, Afghanistan [35.92,71.13] in the period 1976–2019

ABSTRACT Landsystem relationships between glaciers and rock debris supply in a mountain landscape domain, (), are described. Decimal latitude-longitude [dLL] geolocations are used to identify features and transects in an information landscape. Geo-located features are coded, enabling transects between a 1976 expedition and 2019 Google Earth imagery to be compared. Rock debris is progressively added to 1-3 km long glaciers which become debris-covered. Cirque glaciers eventually assume rock glacier (RG) forms when supraglacial debris loads are high. Some rock glacier snouts reach main valley floors and still advance over meadows. This behaviour is attributed to high geomorphic activity producing rock detritus and transport to glaciers in the early Little Ice Age. The advances of rock glacier snouts are a consequence of thinning; low-angle glaciers still moving beneath debris-covered glaciers (GLd) covers. Persistent melt pools continue to develop within the surface debris cover of glaciers and rock glaciers and expose glacier ice. All the rock glaciers are below the regional snowline and permafrost can be discounted for rock glacier formation. Scree slope (SS) development may ultimately be sufficient to cover bare glacier ice moving from a glacier (GL) to debris-covered glacier (GLd) to rock glacier (RG). Reverse slopes in the debris at the foot of screes mark the mass continuum of glacier flow below the debris cover, not the ‘rooting zone’ of a permafrost-derived RG. Scree slopes themselves show no evidence of rock glacier-like flow. A simple glacier ice-debris transport continuum model is sufficient and necessary for rock glacier formation and flow.

Ice Content of Mantling Materials in Deuteronilus Mensae, Mars

AbstractLayers of ice and dust (i.e., “mantle”) are found throughout the mid‐latitudes of Mars, where they drape and mute topography and record the recent climate history. Many of the youngest Late Amazonian “latitude‐dependent mantles” have been characterized, but the full range of mantle characteristics, including variations in ice content and age, has not been well documented. To advance our understanding of these ice‐dust units, we use SHAllow RADar (SHARAD) radar sounding data to constrain the subsurface physical properties of a widespread Middle Amazonian mantle unit within Deuteronilus Mensae. This region hosts a high density of glacial landforms, which allows assessment of the stratigraphic relationship between mantle deposits and glacial ice. SHARAD reflectors at the base of the mantle units are used to determine moderate dielectric constants (∼5.75) and high loss tangents (∼0.038). These radar properties imply relatively limited ice content and high‐loss mantling materials with possible enhanced roughness or subsurface scattering. Further, we interpret deeper reflectors that are continuous with glacier basal reflectors and that suggest materials with lower dielectric constants (∼3.2–4.5) and loss tangents (0.0071–0.0120) to be due to extensions of stagnant, possibly lithic‐rich, glacier ice that are buried by a mantle layer. This study demonstrates the high variability in the origins and evolution of mid‐latitude mantling units on Mars. While the most recent latitude‐dependent mantling deposits may still be ice‐rich, older deposits, such as those preserved in Deuteronilus Mensae and elsewhere in the mid‐latitudes of Mars, may be largely desiccated or contain a more heterogeneous and limited distribution of ice.

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. We demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

94 GHz Radar Backscatter Characteristics of Alpine Glacier Ice

94 GHz Radar Backscatter Characteristics of Alpine Glacier Ice

Experimental and numerical investigation of accretion inception and heat transfer physics in ice crystal icing

Earlier research highlighted severe detrimental effects related to aircraft performance due to ice crystal icing. As part of current research impact of glaciated icing cloud on heated substrates is investigated. To that end, a series of experiments were carried out at the icing wind tunnel of Technische Universitt Braunschweig to investigate the physics of ice crystal accretion on heated substrates for various operating conditions. Dedicated transparent and metallic heatable substrates were designed for macro and microscopic investigation of ice accretion from a qualitative visualization and heat transfer physics perspective. In addition to that numerical simulations were performed at ONERA. Qualitative observations from dedicated experiments on ice accretion initiation phase resulted in advancement of numerical model to accurately capture the different stages leading up to ice accretion as well as the necessary conditions for growth of slushy and glaciated ice layers. Furthermore, a dedicated test matrix was defined to explicitly study the influence of dominant parameters such as flow velocity, heat flux, wet bulb temperature and ice water content on ice accretion process. The experiments showcased a strong influence of increasing flow velocity and ice water content yielding shorter duration required to accrete an ice layer on a heatable substrate. Further experimental investigation reflected that upon increasing the heating power of the test article the icing cloud had to overcome a larger temperature gradient resulting in longer duration required for accretion. It was found that on one hand, the test run with negative wet bulb temperature required a heating source from the substrate as the necessary condition for ice accretion resulting in glaciated ice layers. On the other hand, for positive wet bulb temperature cases natural melting of ice layer was sufficient to induce ice accretion generating slushy ice layers and the heating source from the substrate had little to no influence on the overall ice accretion growth. Numerical simulations were also performed for same conditions and were able to correctly capture the trends and orders of magnitude in comparison with experimental results. The experimental findings presented in this paper lead to development of an accretion solver based on enthalpy approach which considers accretion process as a homogeneous mixture of crystals and liquid water and not by a mere superposition of ice layer and liquid water film. The findings helped calibrate, validate and advance the numerical model for ice crystal icing and found to be more representative of accretion process than unsteady triple layer approach. The findings ensure better predictive capability resulting in improved flight safety and performance criterion.

Modelling rock glacier ice content based on InSAR-derived velocity, Khumbu and Lhotse valleys, Nepal

Abstract. Active rock glaciers are viscous flow features embodying ice-rich permafrost and other ice masses. They contain significant amounts of ground ice and serve as potential freshwater reservoirs as mountain glaciers melt in response to climate warming. However, current knowledge about ice content in rock glaciers has been acquired mainly from in situ investigations in limited study areas, which hinders a comprehensive understanding of ice storage in rock glaciers situated in remote mountains over local to regional scales. This study proposes a novel approach for assessing the hydrological value of rock glaciers in a more quantitative way and presents exploratory results focusing on a small region. We develop an empirical rheological model to infer ice content of rock glaciers using readily available input data, including rock glacier planar shape, surface slope angle, active layer thickness, and surface velocity. The model is calibrated and validated using observational data from the Chilean Andes and the Swiss Alps. We apply the model to five rock glaciers in the Khumbu and Lhotse valleys, northeastern Nepal. The velocity constraints applied to the model are derived from interferometric synthetic aperture radar (InSAR) measurements. The volume of rock glacier is estimated based on an existing scaling approach. The inferred volumetric ice fraction in the Khumbu and Lhotse valleys ranges from 70 ± 8 % to 74 ± 8 %, and the water volume equivalents lie between 1.4 ± 0.2 and 5.9±0.6×106 m3 for the coherently moving parts of individual rock glaciers. Due to the accessibility of the model inputs, our approach is applicable to permafrost regions where observational data are lacking, which is valuable for estimating the water storage potential of rock glaciers in remote areas.

Volcanogenic fluxes of iron from the seafloor in the Amundsen Sea, West Antarctica

The Amundsen Sea in the Pacific sector of West Antarctica receives meltwater from the fastest retreating Antarctic glaciers, and its coastal polynyas host the highest primary productivity per unit area observed on the Antarctic continental shelf. Polynya productivity provides the base for a robust, diverse ecosystem and is controlled primarily by light and the availability of the micronutrient iron (Fe). While the sources of Fe in the region are not yet certain, Fe could be transported within modified Circumpolar Deep Water (mCDW) that intrudes onto the retrograde shelf and into ice shelf cavities, where it gains buoyancy through the addition of glacial meltwater and is injected into the upper water column when it exits the cavity. Thus, fluxes of dissolved Fe from the seafloor into in-flowing mCDW may ultimately be a source of Fe to the euphotic zone in the Amundsen Sea. To investigate the surface sediment biogeochemistry and the potential for a significant benthic flux of Fe to the waters on the Amundsen Sea shelf, sediment cores were collected at two sites close to the calving fronts of the Pine Island and Thwaites Glacier ice shelves. Pore water was analyzed for trace element content, and sediment was analyzed for physical and chemical properties including organic carbon and trace elements. Using a novel approach based on hypothesized Fe speciation and colloidal particle radius, theoretical Fe fluxes were calculated from pore water gradients and porosity. The fluxes reveal a spatially variable Fe input to the lower water column that could ultimately fertilize primary productivity. Supported by geochemical and physical evidence, we conclude that submarine weathering of volcanic glass grains observed and quantified in seabed sediments at the Pine Island site drives nonreductive Fe fluxes that are 100-fold higher than at the Thwaites site. This study highlights the need for further investigations of benthic-pelagic coupling in the Amundsen Sea region, which will likely be impacted in coming decades by accelerating glacial melting.

The CryoGrid community model (version 1.0) – a multi-physics toolbox for climate-driven simulations in the terrestrial cryosphere

Abstract. The CryoGrid community model is a flexible toolbox for simulating the ground thermal regime and the ice–water balance for permafrost and glaciers, extending a well-established suite of permafrost models (CryoGrid 1, 2, and 3). The CryoGrid community model can accommodate a wide variety of application scenarios, which is achieved by fully modular structures through object-oriented programming. Different model components, characterized by their process representations and parameterizations, are realized as classes (i.e., objects) in CryoGrid. Standardized communication protocols between these classes ensure that they can be stacked vertically. For example, the CryoGrid community model features several classes with different complexity for the seasonal snow cover, which can be flexibly combined with a range of classes representing subsurface materials, each with their own set of process representations (e.g., soil with and without water balance, glacier ice). We present the CryoGrid architecture as well as the model physics and defining equations for the different model classes, focusing on one-dimensional model configurations which can also interact with external heat and water reservoirs. We illustrate the wide variety of simulation capabilities for a site on Svalbard, with point-scale permafrost simulations using, e.g., different soil freezing characteristics, drainage regimes, and snow representations, as well as simulations for glacier mass balance and a shallow water body. The CryoGrid community model is not intended as a static model framework but aims to provide developers with a flexible platform for efficient model development. In this study, we document both basic and advanced model functionalities to provide a baseline for the future development of novel cryosphere models.

Melt rates in the kilometer-size grounding zone of Petermann Glacier, Greenland, before and during a retreat

Warming of the ocean waters surrounding Greenland plays a major role in driving glacier retreat and the contribution of glaciers to sea level rise. The melt rate at the junction of the ocean with grounded ice-or grounding line-is, however, not well known. Here, we employ a time series of satellite radar interferometry data from the German TanDEM-X mission, the Italian COSMO-SkyMed constellation, and the Finnish ICEYE constellation to document the grounding line migration and basal melt rates of Petermann Glacier, a major marine-based glacier of Northwest Greenland. We find that the grounding line migrates at tidal frequencies over a kilometer-wide (2 to 6 km) grounding zone, which is one order of magnitude larger than expected for grounding lines on a rigid bed. The highest ice shelf melt rates are recorded within the grounding zone with values from 60 ± 13 to 80 ± 15 m/y along laterally confined channels. As the grounding line retreated by 3.8 km in 2016 to 2022, it carved a cavity about 204 m in height where melt rates increased from 40 ± 11 m/y in 2016 to 2019 to 60 ± 15 m/y in 2020 to 2021. In 2022, the cavity remained open during the entire tidal cycle. Such high melt rates concentrated in kilometer-wide grounding zones contrast with the traditional plume model of grounding line melt which predicts zero melt. High rates of simulated basal melting in grounded glacier ice in numerical models will increase the glacier sensitivity to ocean warming and potentially double projections of sea level rise.

The microbiome of the ice-capped Cayambe Volcanic Complex in Ecuador.

A major challenge in microbial ecology is to understand the principles and processes by which microbes associate and interact in community assemblages. Microbial communities in mountain glaciers are unique as first colonizers and nutrient enrichment drivers for downstream ecosystems. However, mountain glaciers have been distinctively sensitive to climate perturbations and have suffered a severe retreat over the past 40 years, compelling us to understand glacier ecosystems before their disappearance. This is the first study in an Andean glacier in Ecuador offering insights into the relationship of physicochemical variables and altitude on the diversity and structure of bacterial communities. Our study covered extreme Andean altitudes at the Cayambe Volcanic Complex, from 4,783 to 5,583 masl. Glacier soil and ice samples were used as the source for 16S rRNA gene amplicon libraries. We found (1) effects of altitude on diversity and community structure, (2) the presence of few significantly correlated nutrients to community structure, (3) sharp differences between glacier soil and glacier ice in diversity and community structure, where, as quantified by the Shannon γ-diversity distribution, the meta-community in glacier soil showed more diversity than in glacier ice; this pattern was related to the higher variability of the physicochemical distribution of variables in the former substrate, and (4) significantly abundant genera associated with either high or low altitudes that could serve as biomarkers for studies on climate change. Our results provide the first assessment of these unexplored communities, before their potential disappearance due to glacier retreat and climate change.

Historical reconstruction and future projection of mass balance and ice volume of Qiyi Glacier, Northeast Tibetan Plateau

Study regionQiyi Glacier in the Qilian Mountains, northeast Tibetan Plateau. Study focusFuture glacier changes are critical to the sustainable use of regional water resources and the survival and development of the densely populated Tibetan Plateau and surrounding regions. A distributed energy-mass balance model is applied to simulate the historical and future mass balance and ice volume of Qiyi Glacier from 1980 to 2100. The aims of this study are to illuminate the spatial and temporal patterns of mass balance and energy balance, to forecast future changes of this glacier under different scenarios of CMIP6, and to identify the timing of peak ice loss. New hydrological insights for the regionCompared to the historical reconstruction, future projections show that ablation will increase, and accumulation will decrease due to less snowfall and increased net shortwave radiation that is induced by lower albedos. The average mass balance increases by a factor of 1.08, 1.67 and 2.68 under the SSP1-2.6, SSP2-4.5 and SSP5-8.5 scenarios during 2021–2100, respectively. By the end of the 21st century, Qiyi Glacier ice volume will reduce by 107.61 × 106 m3 (95.3 %) for the SSP1-2.6 scenario, and 112.12 × 106 m3 (99.3 %) for the SSP2-4.5 scenario, and will disappear entirely by 2082 for the SSP5-8.5 scenario. In all scenarios, the glacier mass loss rate will peak around 2030–2045, then decrease in parallel with the decreasing ice volume.

Meltwater generation in ice stream shear margins: case study in Antarctic ice streams

Liquid water within glacier ice and at the glacier beds exerts a significant control on ice flow and glacier stability through a number of processes, including altering the rheology of the ice and lubricating the bed. Some of this water is generated as melt from regions of rapid deformation, including shear margins, due to heating by viscous dissipation. However, how much meltwater is generated and drained from shear margins remains unclear. Here, we apply a model that describes the evolution of ice temperature, melting, and water transport within deforming ice to estimate the flux of meltwater from shear margins in glaciers. We estimate the flux of meltwater from temperate ice zones in three Antarctic regions: Bindschadler and MacAyeal Ice Streams, Pine Island Glacier, and Byrd Glacier. We show that the flux of meltwater from shear margins in these regions may be as significant as the meltwater produced by frictional heating at the bed, with average fluxes of ∼ 0.005 – 0.1 m yr − 1 . This contribution of shear heating to meltwater flux at the bed may thus affect both the rheology of the ice as well as sliding at the bed, both key controls on fast ice flow.

A REPLICABLE OPEN-SOURCE MULTI-CAMERA SYSTEM FOR LOW-COST 4D GLACIER MONITORING

Abstract. Image-based monitoring has emerged as a prevalent technique for sensing mountain environments. Monoscopic time-lapse cameras, which rely on digital image correlation to quantify glacier motion, have limitations due to the need for a Digital Elevation Model for deriving 3D flow velocity fields. Multi-camera systems overcome this limitation, as they allow for a 3D reconstruction of the scene. This paper presents a replicable low-cost stereoscopic system designed for 4D glacier monitoring. The system consists of independent and autonomous units, built from off-the-shelves components, such as a DSLR camera, an Arduino microcontroller, and a Raspberry Pi Zero, reducing costs compared to pre-built time-lapse cameras. The units are energetically self-sufficient and resistant to harsh alpine conditions. The system was successfully tested for more than a year to monitor the northwest terminus of the Belvedere Glacier (Italian Alps). Daily stereo-pairs acquired were processed with Structure-from-Motion to derive 3D point clouds of the glacier terminus and estimate glacier retreat and ice volume loss. By combining the information about ice volume loss with ablation estimates and ice flow velocity information, e.g., derived from monoscopic-camera time series, a multi-camera system enables a comprehensive understanding of sub-seasonal glacier dynamics.