Abstract. Knowledge about the occurrence and characteristics of surge-type glaciers is crucial due to the impact of surging on glacier melt and glacier-related hazards. One of the super-clusters of surge-type glaciers is High Mountain Asia (HMA). However, no consistent region-wide inventory of surge-type glaciers in HMA exists. We present a regionally resolved inventory of surge-type glaciers based on their behaviour across High Mountain Asia between 2000 and 2018. We identify surge-type behaviour from surface velocity, elevation and feature change patterns using a multi-factor remote sensing approach that combines yearly ITS_LIVE velocity data, DEM differences and very-high-resolution imagery (Bing Maps, Google Earth). Out of the ≈95 000 glaciers in HMA, we identified 666 that show diagnostic surge-type glacier behaviour between 2000 and 2018, which are mainly found in the Karakoram (223) and the Pamir regions (223). The total area covered by the 666 surge-type glaciers represents 19.5 % of the glacierized area in Randolph Glacier Inventory (RGI) V6.0 polygons in HMA. Only 68 glaciers were already identified as “surge type” in the RGI V6.0. We further validate 107 glaciers previously labelled as “probably surge type” and newly identify 491 glaciers, not previously reported in other inventories covering HMA. We finally discuss the possibility of self-organized criticality in glacier surges. Across all regions of HMA, the surge-affected area within glacier complexes displays a significant power law dependency with glacier length.

This paper presents a survey of georeferenced point clouds. Concentration is, on the one hand, put on features, which originate in the measurement process themselves, and features derived by processing the point cloud. On the other hand, approaches for the processing of georeferenced point clouds are reviewed. This includes the data structures, but also spatial processing concepts. We suggest a categorization of features into levels that reflect the amount of processing. Point clouds are found across many disciplines, which is reflected in the versatility of the literature suggesting specific features.

Highlights d Glaciers around Mt. Everest have thinned by more than 100 m since the 1960s d The rate of ice mass loss has consistently accelerated over the past six decades d Glacier thinning has occurred at above 6,000 masl d Surge-type glacier behavior has been identified for the first time in the region

An active landslide in Doren, Austria, has been studied by multitemporal airborne and terrestrial laser scanning from 2003 to 2012. To evaluate the changes, we have determined the 3D motion using the range flow algorithm, an established method in computer vision, but not yet used for studying landslides. The generated digital terrain models are the input for motion estimation; the range flow algorithm has been combined with the coarse-to-fine resolution concept and robust adjustment to be able to determine the various motions over the landslide. The algorithm yields fully automatic dense 3D motion vectors for the whole time series of the available data. We present reliability measures for determining the accuracy of the estimated motion vectors, based on the standard deviation of components. The differential motion pattern is mapped by the algorithm: parts of the landslide show displacements up to 10 m, whereas some parts do not change for several years. The results have also been compared to pointwise reference data acquired by independent geodetic measurements; reference data are in good agreement in most of the cases with the results of range flow algorithm; only some special points (e.g., reflectors fixed on trees) show considerably differing motions.

Planet Labs have recently launched a large constellation of small satellites (3U cubesats) capable of imaging the whole Earth landmass everyday. These small satellites capture multiple images of an area on consecutive days or sometimes on the same day with a spatial resolution of 3–4 m. Planet Labs endeavors to operate the constellation in a nadir pointing mode, however, the view angle of these satellites currently varies within a few degrees from the nadir leading to varying B/H ratio for overlapping image pairs. Due to relatively small scene footprint and small off-nadir angle, the baseline to height ratio (B/H) of the overlapping PlanetScope images is often less than 1:10, which is not ideal for 3D reconstruction. Therefore, this paper explores the potential of Digital Elevation Model generation from this multi-date, multi-satellite PlanetScope imagery. The DEM generation from multiple PlanetScope images is achieved using a volumetric stereo reconstruction technique, which applies semi global matching in georeferenced object space. The results are evaluated using a LiDAR based DEM (5 m) over Mount Teide (3718 m) in Canary Islands and the ALOS (30 m) DEM on rugged terrain of the Nanga Parbat massif (8126 m) in the western Himalaya range. The proposed methodology is then applied on images from two PlanetScope satellites overpasses within a couple of minutes difference to compute the DEM of the Khurdopin glacier in the Karakoram range, known for its recent surge. The quantitative assessment of the generated elevation models is done by comparing statistics of the elevation differences between the reference LiDAR and ALOS DEM and the PlanetScope DEM. The Normalized Median of Absolute Deviation (NMAD) of the elevation differences between the computed PlanetScope DEM and LiDAR DEM is 4.1 m and the elevation differences for the ALOS DEM over stable terrain is 3.9 m. The results show that PlanetScope imagery can lead to sufficient quality DEM even with a small baseline to height ratio. Therefore, the daily PlanetScope imagery is a valuable data source and the DEM generated from this imagery can potentially be employed in numerous applications requiring multi temporal DEMs.

No abstract

Glacial lakes mapping using satellite remote sensing data are important for studying the effects of climate change as well as for the mitigation and risk assessment of a Glacial Lake Outburst Flood (GLOF). The 3U cubesat constellation of Planet Labs offers the capability of imaging the whole Earth landmass everyday at 3–4 m spatial resolution. The higher spatial, as well as temporal resolution of PlanetScope imagery in comparison with Landsat-8 and Sentinel-2, makes it a valuable data source for monitoring the glacial lakes. Therefore, this paper explores the potential of the PlanetScope imagery for glacial lakes mapping with a focus on the Hindu Kush, Karakoram and Himalaya (HKKH) region. Though the revisit time of the PlanetScope imagery is short, courtesy of 130+ small satellites, this imagery contains only four bands and the imaging sensors in these small satellites exhibit varying spectral responses as well as lower dynamic range. Furthermore, the presence of cast shadows in the mountainous regions and varying spectral signature of the water pixels due to differences in composition, turbidity and depth makes it challenging to automatically and reliably extract surface water in PlanetScope imagery. Keeping in view these challenges, this work uses state of the art deep learning models for pixel-wise classification of PlanetScope imagery into the water and background pixels and compares the results with Random Forest and Support Vector Machine classifiers. The deep learning model is based on the popular U-Net architecture. We evaluate U-Net architecture similar to the original U-Net as well as a U-Net with a pre-trained EfficientNet backbone. In order to train the deep neural network, ground truth data are generated by manual digitization of the surface water in PlanetScope imagery with the aid of Very High Resolution Satellite (VHRS) imagery. The created dataset consists of more than 5000 water bodies having an area of approx. 71km2 in eight different sites in the HKKH region. The evaluation of the test data show that the U-Net with EfficientNet backbone achieved the highest F1 Score of 0.936. A visual comparison with the existing glacial lake inventories is then performed over the Baltoro glacier in the Karakoram range. The results show that the deep learning model detected significantly more lakes than the existing inventories, which have been derived from Landsat OLI imagery. The trained model is further evaluated on the time series PlanetScope imagery of two glacial lakes, which have resulted in an outburst flood. The output of the U-Net is also compared with the GLakeMap data. The results show that the higher spatial and temporal resolution of PlanetScope imagery is a significant advantage in the context of glacial lakes mapping and monitoring.

The Corona KH-4 reconnaissance satellite missions acquired panoramic stereo imagery with high spatial resolution of 1.8-7.5 m from 1962-1972. The potential of 800,000+ declassified Corona images has not been leveraged due to the complexities arising from handling of panoramic imaging geometry, film distortions and limited availability of the metadata required for georeferencing of the Corona imagery. This paper presents the Corona Stereo Pipeline (CoSP): A pipeline for processing of Corona KH-4 stereo panoramic imagery. CoSP utilizes deep learning based feature matcher SuperGlue to automatically match features point between Corona KH-4 images and recent satellite imagery to generate Ground Control Points (GCPs). To model the imaging geometry and the scanning motion of the panoramic KH-4 cameras, a rigorous camera model consisting of modified collinearity equations with time-dependent exterior orientation parameters is employed. Using the entire frame of the Corona image, bundle adjustment with well-distributed GCPs results in an average standard deviation or σ0 of less than two pixels. We evaluate fiducial marks on the Corona films and show that preprocessing the Corona images to compensate for film bending improves the 3D reconstruction accuracy. The distortion pattern of image residuals of GCPs and y-parallax in epipolar resampled images suggest that film distortions due to long-term storage likely cause systematic deviations of up to six pixels. Compared to the SRTM DEM, the Corona DEM computed using CoSP achieved a Normalized Median Absolute Deviation of elevation differences of ≈ 4 m over an area of approx. 4000 km 2 after a tilebased fine coregistration of the DEMs. We further assess CoSP on complex scenes involving high relief and glacierized terrain and show that the resulting DEMs can be used to compute long-term glacier elevation changes over large areas.