Suitability of ground-based SfM-MVS for monitoring glacial and periglacial processes

Photo-based surface reconstruction is rapidly emerging as an alternative survey technique to lidar (light detection and ranging) in many fields of geoscience fostered by the recent development of computer vision algorithms such as structure from motion (SfM) and dense image matching such as multi-view stereo (MVS). The objectives of this work are to test the suitability of the ground-based SfM-MVS approach for calculating the geodetic mass balance of a 2.1km2 glacier and for detecting the surface displacement of a neighbouring active rock glacier located in the eastern Italian Alps. The photos were acquired in 2013 and 2014 using a digital consumer-grade camera during single-day field surveys. Airborne laser scanning (ALS, otherwise known as airborne lidar) data were used as benchmarks to estimate the accuracy of the photogrammetric digital elevation models (DEMs) and the reliability of the method. The SfM-MVS approach enabled the reconstruction of high-quality DEMs, which provided estimates of glacial and periglacial processes similar to those achievable using ALS. In stable bedrock areas outside the glacier, the mean and the standard deviation of the elevation difference between the SfM-MVS DEM and the ALS DEM was-0.42 \ub1 1.72 and 0.03 \ub1 0.74 m in 2013 and 2014, respectively. The overall pattern of elevation loss and gain on the glacier were similar with both methods, ranging between-5.53 and + 3.48 m. In the rock glacier area, the elevation difference between the SfM-MVS DEM and the ALS DEM was 0.02 \ub1 0.17 m. The SfM-MVS was able to reproduce the patterns and the magnitudes of displacement of the rock glacier observed by the ALS, ranging between 0.00 and 0.48 m per year. The use of natural targets as ground control points, the occurrence of shadowed and low-contrast areas, and in particular the suboptimal camera network geometry imposed by the morphology of the study area were the main factors affecting the accuracy of photogrammetric DEMs negatively. Technical improvements such as using an aerial platform and/or placing artificial targets could significantly improve the results but run the risk of being more demanding in terms of costs and logistics

Estimating the characteristics of runoff inflow into Lake Gojal in ungauged, highly glacierized upper Hunza River Basin, Pakistan

Motivated by the potential flood outburst of Lake Gojal in the ungauged highly glacierized (27%) upper Hunza River Basin (HRB) in Pakistan that was dammed by a massive landslide on 4 January 2010, we attempt to analyze the characteristics of water inflow to the lake employing remote sensing data, two hydrological models, and sparsely observed data. One of the models (Model I) is a monthly degree-day model, while another (Model II) is the variable infiltration capacity (VIC) model. The mixture of glacier runoff output from Model I and runoff over unglacierized areas calculated by Model II has a similar seasonal variation pattern as that estimated from data recorded at a downstream station. This suggests that glacier runoff is the main source (87%) of runoff inflow into the lake. A sensitivity analysis suggests that the water inflow to the lake is highly sensitive to an increase in air temperature. Runoff in May is predicted to sharply increase by 15% to more than two-fold if the air temperature increases by 1 to 7, but it is predicted to increase only from 9% to 34% if the precipitation increases by 10% to 40%. The results suggested that the water inflow into Lake Gojal will not sharply rise even if there is heavy rain, and it needs to be in caution if the air temperature sharply increases. Analysis on long-term air temperature record indicates that the water inflow into the lake in May 2010 was probably less than average owing to the relatively low air temperature. Consequently, the flood outburst did not occur before the completion of the spillway on 29 May 2010. © 2013 China University of Geosciences and Springer-Verlag Berlin Heidelberg

Toward an imminent extinction of Colombian glaciers?

This study documents the current state of glacier coverage in the Colombian Andes, the glacier shrinkage over the twentieth century and discusses indication of their disappearance in the coming decades. Satellite images have been used to update the glacier inventory of Colombia reflecting an overall glacier extent of about 42.4 ± 0.71 km2 in 2016 distributed in four glacierized mountain ranges. Combining these data with older inventories, we show that the current extent is 36% less than in the mid-1990s, 62% less than in the mid-twentieth century and almost 90% less than the Little Ice Age maximum extent. Focusing on Nevado Santa Isabel (Los Nevados National Park), aerial photographs from 1987 and 2005 combined with a terrestrial LiDAR survey show that the mass loss of the former ice cap, which is nowadays parceled into several small glaciers, was about −2.5 m w.e. yr−1 during the last three decades. Radar measurements performed on one of the remnant glaciers, La Conejeras glacier, show that the ice thickness is limited (about 22 m in average in 2014) and that with such a mass loss rate, the glacier should disappear in the coming years. Considering their imbalance with the current climate conditions, their limited altitudinal extent and reduced accumulation areas, and in view of temperature increase expected in future climate scenarios, most of the Colombian glaciers will likely disappear in the coming decades. Only the largest ones located on the highest summits will probably persist until the second half of the twenty-first century although very reducedThis study was conducted in the context of the project Capacity Building and Twinning for Climate Observing Systems (CATCOS) supported by the Federal Office of Meteorology and Climatology MeteoSwiss [contract no. 7F-08114.1], between the Swiss Agency for Development and Cooperation (SDC) and MeteoSwiss, by the Swiss State Secretariat for Economic Affairs (SECO). This work was also supported by SNO GLACIOCLIM; LMI GREAT ICE (IRD); Labex OSUG@2020, Investissements d’avenir: [Grant Number ANR10 LABX56]. M. Ménégoz is supported by the project VOLCADEC funded by the Spanish programme Retos (MINECO/FEDER, ref. CGL2015–70177-R).Peer ReviewedPostprint (author's final draft

Including debris cover effects in a distributed model of glacier ablation

Distributed glacier melt models generally assume that the glacier surface consists of bare exposed ice and snow. In reality, many glaciers are wholly or partially covered in layers of debris that tend to suppress ablation rates. In this paper, an existing physically based point model for the ablation of debris-covered ice is incorporated in a distributed melt model and applied to Haut Glacier d’Arolla, Switzerland, which has three large patches of debris cover on its surface. The model is based on a 10 m resolution digital elevation model (DEM) of the area; each glacier pixel in the DEM is defined as either bare or debris-covered ice, and may be covered in snow that must be melted off before ice ablation is assumed to occur. Each debris-covered pixel is assigned a debris thickness value using probability distributions based on over 1000 manual thickness measurements. Locally observed meteorological data are used to run energy balance calculations in every pixel, using an approach suitable for snow, bare ice or debris-covered ice as appropriate. The use of the debris model significantly reduces the total ablation in the debris-covered areas, however the precise reduction is sensitive to the temperature extrapolation used in the model distribution because air near the debris surface tends to be slightly warmer than over bare ice. Overall results suggest that the debris patches, which cover 10% of the glacierized area, reduce total runoff from the glacierized part of the basin by up to 7%

Terminus-driven retreat of a major southwest Greenland tidewater glacier during the early 19th century : insights from glacier reconstructions and numerical modelling

Peer reviewedPublisher PD

Suitability of a constant air temperature lapse rate over an Alpine glacier: testing the Greuell and Böhm model as an alternative

Near-surface air temperature, typically measured at a height of 2 m, is the most important control on the energy exchange and the melt rate at a snow or ice surface. It is distributed in a simplistic manner in most glacier melt models by using constant linear lapse rates, which poorly represent the actual spatial and temporal variability of air temperature. In this paper, we test a simple thermodynamic model proposed by Greuell and Böhm in 1998 as an alternative, using a new dataset of air temperature measurements from along the flowline of Haut Glacier d'Arolla, Switzerland. The unmodified model performs little better than assuming a constant linear lapse rate. When modified to allow the ratio of the boundary layer height to the bulk heat transfer coefficient to vary along the flowline, the model matches measured air temperatures better, and a further reduction of the root-mean-square error is obtained, although there is still considerable scope for improvement. The modified model is shown to perform best under conditions favourable to the development of katabatic winds – few clouds, positive ambient air temperature, limited influence of synoptic or valley winds and a long fetch – but its performance is poor under cloudy conditions

Simulations of snow distribution and hydrology in a mountain basin

We applied a version of the Regional Hydro‐Ecologic Simulation System (RHESSys) that implements snow redistribution, elevation partitioning, and wind‐driven sublimation to Loch Vale Watershed (LVWS), an alpine‐subalpine Rocky Mountain catchment where snow accumulation and ablation dominate the hydrologic cycle. We compared simulated discharge to measured discharge and the simulated snow distribution to photogrammetrically rectified aerial (remotely sensed) images. Snow redistribution was governed by a topographic similarity index. We subdivided each hillslope into elevation bands that had homogeneous climate extrapolated from observed climate. We created a distributed wind speed field that was used in conjunction with daily measured wind speeds to estimate sublimation. Modeling snow redistribution was critical to estimating the timing and magnitude of discharge. Incorporating elevation partitioning improved estimated timing of discharge but did not improve patterns of snow cover since wind was the dominant controller of areal snow patterns. Simulating wind‐driven sublimation was necessary to predict moisture losses

Accurate simulation of ice and snow runoff for the mountainous terrain of the Kunlun Mountains, China

While mountain runoff provides great potential for the development and life quality of downstream populations, it also frequently causes seasonal disasters. The accurate modeling of hydrological processes in mountainous areas, as well as the amount of meltwater from ice and snow, is of great significance for the local sustainable development, hydropower regulations, and disaster prevention. In this study, an improved model, the Soil Water Assessment Tool with added ice-melt module (SWATAI) was developed based on the Soil Water Assessment Tool (SWAT), a semi-distributed hydrological model, to simulate ice and snow runoff. A temperature condition used to determine precipitation types has been added in the SWATAI model, along with an elevation threshold and an accumulative daily temperature threshold for ice melt, making it more consistent with the runoff process of ice and snow. As a supplementary reference, the comparison between the normalized difference vegetation index (NDVI) and the quantity of meltwater were conducted to verify the simulation results and assess the impact of meltwater on the ecology. Through these modifications, the accuracy of the daily flow simulation results has been considerably improved, and the contribution rate of ice and snow melt to the river discharge calculated by the model increased by 18.73%. The simulation comparison of the flooding process revealed that the accuracy of the simulated peak flood value by the SWATAI was 77.65% higher than that of the SWAT, and the temporal accuracy was 82.93% higher. The correlation between the meltwater calculated by the SWATAI and the NDVI has also improved significantly. This improved model could simulate the flooding processes with high temporal resolution in alpine regions. The simulation results could provide technical support for economic benefits and reasonable reference for flood prevention

Glacier motion estimation using SAR offset-tracking procedures

Two image-to-image patch offset techniques for estimating feature motion between satellite synthetic aperture radar (SAR) images are discussed. Intensity tracking, based on patch intensity cross-correlation optimization, and coherence tracking, based on patch coherence optimization, are used to estimate the movement of glacier surfaces between two SAR images in both slant-range and azimuth direction. The accuracy and application range of the two methods are examined in the case of the surge of Monacobreen in Northern Svalbard between 1992 and 1996. Offset-tracking procedures of SAR images are an alternative to differential SAR interferometry for the estimation of glacier motion when differential SAR interferometry is limited by loss of coherence, i.e., in the case of rapid and incoherent flow and of large acquisition time intervals between the two SAR images. In addition, an offset-tracking procedure in the azimuth direction may be combined with differential SAR interferometry in the slant-range direction in order to retrieve a two-dimensional displacement map when SAR data of only one orbit configuration are available