The Status of Earth Observation Techniques in Monitoring High Mountain Environments at the Example of Pasterze Glacier, Austria: Data, Methods, Accuracies, Processes, and Scales

: Earth observation offers a variety of techniques for monitoring and characterizing geomorphic processes in high mountain environments. Terrestrial laserscanning and unmanned aerial vehicles provide very high resolution data with high accuracy. Automatic cameras have become a valuable source of information—mostly in a qualitative manner—in recent years. The availability of satellite data with very high revisiting time has gained momentum through the European Space Agency’s Sentinel missions, offering new application potential for Earth observation. This paper reviews the status of recent techniques such as terrestrial laserscanning, remote sensed imagery, and synthetic aperture radar in monitoring high mountain environments with a particular focus on the impact of new platforms such as Sentinel-1 and -2 as well as unmanned aerial vehicles. The study area comprises the high mountain glacial environment at the Pasterze Glacier, Austria. The area is characterized by a highly dynamic geomorphological evolution and by being subject to intensive scientiﬁc research as well as long-term monitoring. We primarily evaluate landform classiﬁcation and process characterization for: (i) the proglacial lake; (ii) icebergs; (iii) the glacier river; (iv) valley-bottom processes; (v) slope processes; and (vi) rock wall processes. We focus on assessing the potential of every single method both in spatial and temporal resolution in characterizing different geomorphic processes. Examples of the individual techniques are evaluated qualitatively and quantitatively in the context of: (i) morphometric analysis; (ii) applicability in high alpine regions; and (iii) comparability of the methods among themselves. The ﬁnal frame of this article includes considerations on scale dependent process detectability and characterization potentials of these


Introduction
Glaciers and their changes are well recognized as crucial indicators for climate change [1][2][3][4]. Glacier retreat has many, potentially severe, impacts on human life, as exemplified by its influence on the availability of freshwater [5,6] or its role in an increase of hazardous events [7]. Glacier fluctuations cause massive impacts on glacio-hydrological or geomorphological process systems across various scales. Changes within a cryospheric environment in the form of glacier retreat result, e.g., in local hazard events [8,9], changes in regional water cycle systems [10][11][12], and sea level rise on a global scale [13,14]. The impacts of glacier retreat in a socio-economic context are far-reaching [15,16], affecting different areas from energy supply [15] to tourism [17,18]. As a consequence, glacier retreat induces hazards due to changing conditions within and different resilience of process systems [8]. Concluding insights from glacier monitoring are able to raise people's awareness to the importance of glaciers for the society [19].
High mountain environments are undergoing major changes due to the impact of the ongoing climate change [2,20]. A large variety of processes-often showing accelerating magnitudes and rates in the last two decades [21]-have been reshaping high mountain environments in recent years [22,23]. Especially areas extensively covered by glaciers-such as the European Alps-show a fast transition from glacially dominated to pro-and paraglacial landscapes since the 1970s [24,25]. Information about, e.g., changes in glacier length, area, and volume are high-confidence indicators of climate change [26,27]. In the European Alps, monitoring the cryosphere, and in particular glaciers, has been performed since the end of the 19th century [28,29]. Since then, there has been an increase in both the number of glaciers observed and the number of measurements per glacier. As a result, countries within the Alps with a long tradition of monitoring systems (such as Austria, France, and Switzerland) have valuable information of glacier change [19,30].
Glaciers are a crucial part of a geomorphological process system and should not be analyzed separately. Different processes occurring with different magnitudes on different spatial and temporal scales are a challenge for establishing a comprehensive monitoring system. Therefore, careful evaluation of every single remote sensing method is essential. This paper aims to present a quantitative and qualitative evaluation of remote sensing methods for monitoring geomorphic processes in a cryospheric environment. Geomorphic processes under investigation are characterized by different occurrence frequency and different superordinate process types (glacial vs. proglacial processes and valley bottom to rock wall processes). Consequently, particular focus lies on the evaluation of the observation of the dynamics of: (i) the proglacial lake; (ii) icebergs; (iii) the glacier river; (iv) valley-bottom processes; (v) slope processes; and (vi) rock wall processes.
Therefore, this work seeks to provide answers to the following aspects: • What are the specifications, uncertainties, and constraints of Earth observation techniques in monitoring geomorphic processes in an cryospheric environment? • How do the Earth observation techniques presented prove to be appropriate for monitoring glacial, and paraglacial processes/landforms of different magnitudes and scales?
This paper provides a review of a variety of Earth observation techniques and demonstrates their applicability in monitoring different geomorphic processes. This extensive overview results in the following breakdown: Section 2 gives a broad overview of the use of Earth observation techniques in the characterization of geomorphic processes in high alpine regions, with special focus on the Alpine region. Section 3 provides an outline of the geomorphic processes and landforms investigated in the study area. Section 4 describes Earth observation methods specifically used for monitoring these processes at the Pasterze Glacier area (in terms of technical specifications, monitoring configurations, and quality assessment). Section 5 presents the quantitative results of Earth observation techniques in monitoring certain geomorphic processes and landforms in the study area. Section 6 discusses the applicability of every single method for characterizing the analysed geomorphic processes. Section 7 provides a classification of Earth observation methods comparing data acquisition specification and processes characteristics. Finally, Section 8 evaluates the applied Earth observation techniques with regard to their suitability for monitoring the processes and landforms, as well as their practical application in the study area.

Earth Observation Techniques for Characterizing Processes in Cryospheric Environments
Earth observation techniques provide valuable databases by area wide data acquisition in order to characterize geomorphic processes in a cryospheric environment. Typically, these include high resolution optical data (aerial images and satellite-borne multi-spectral images) as well as data which are subsequently processed to high resolution digital terrain models (DTMs). In Austria, aerial images have been used to, e.g., delineate the extent of glaciers since the 1950s leading to the first so-called Glacier cadaster in 1969 [29] in high mountain applications. Beginning around 2000, airborne (ALS) [31] as well as terrestrial (TLS) [32] laserscanning became a crucial technique to acquire very precise and high resolution surface data. This section provides a short technical description of every Earth observation technique used in this work, and addresses main applications in monitoring geomorphic processes.

Terrestrial Laserscanning
TLS is a time-of-flight system which allows measurements with a range accuracies of a few centimeters [32]. Laserscanning combines the specifications and hence advantages of laser (directional nature of the rays) and radar (location) [33]. Principles of TLS are summarized in [34,35]. The resulting point cloud is registered with the aid of (reflective) objects such as spheres or cylinders representing known coordinates. Carefully registered point clouds enable precise comparison of multi-temporal measurements. TLS uses the near-infrared section of the spectrum at different wavelength such as approximately 1000 nm for snow and ice applications (e.g., Riegl LPM-i800HA [36] and Riegl VZ6000 [37]), in which the wavelength coincides with the measurement range [38] and approximately 1500 nm for other applications (e.g., Riegl LMS Z620 [24]).
In the last decade, massive glacier retreat, especially of large valley glaciers, has resulted in an increased formation of glacial lakes [57]. Due to natural hazards related to glacial lakes such as glacier lake outburst floods (GLOFs), monitoring these lakes has become a crucial task in cryosphere research [58,59]. Since GLOFs are not among the most frequent hazards in the Alps, monitoring the evolution of glacier lakes using TLS was limited to work on Brenva Glacier [54]. For TLS, methodological consideration discussing monitoring concepts, accuracies and uncertainties were made by Ingensand et al. [60], and in a geomorphological context by Fey and Wichmann [61].

Radar Satellite Application: Backscatter
SAR instruments generally provide backscatter measurements which are influenced by the terrain structure (surface roughness). High backscatter values are caused by surfaces with higher roughness as incoming radar pulses are scattered in all directions (diffuse reflection) [62]. Contrarily, calm water surfaces show a very smooth surface, reflecting the radar pulse away from the sensor (specular reflection). Therefore, water surfaces typically show lower backscatter values than adjacent surface types and can be measured applying a simple threshold approach. A widely used approach for threshold selection was provided by Otsu [63]. This method, however, relies on a bimodal histogram with a clear minimum, dividing areas of water and land surface. Manual classification was shown applicable to delineate glacial lakes [64] as well as extracting buffered polygons of the lake area in order to obtain a bimodal histogram [65], or recently by level-set segmentation [66].

Radar Satellite Application: Differential Interferometric Synthetic Aperture Radar (DInSAR)
DInSAR offers a range of approaches to detect small surface deformations with sub-centimeter accuracy. Repeat acquisition of the same constellation can be used to identify small changes in range direction through measured differences in phase. Consequently, it is possible to measure ground deformations in the magnitude of a fraction of the used wavelength [67]. To determine this particular component of the phase (caused by surface displacement), other interfering aspects have to be separated, such as error components due to atmosphere, orbital errors, or phase noise [68]. This can be achieved by using pixels of small phase-noise, which is realized by two reflector types: (i) Permanent or Persistent Scatterers (PS), which consist of a dominating scatterer persisting over time; and (ii) Distributed Scatterers (DS), which feature constant signals caused by different small scatterers [68]. In applying Persistent Scatterer Interferometry (PSI), PS within a time series are used by calculating interferograms in relation to one single master scene [69,70]. Small Baseline Subset (SBAS) is the second multi-temporal InSAR method, which particularly incorporates DS. In contrast to PSI, the retrieval of DS in SBAS becomes increasingly unlikely for interferograms with larger temporal baselines. Therefore, SBAS is not based on one single master scene. Rather, multiple master-slave combinations with small baselines are used in the calculation of the interferograms [71].
In high mountain environments, DInSAR using ERS data was applied to detect slope movements in the Swiss Alps, allowing the implementation of an inventory of mass movement types [72]. However, limitations were identified for steep rock walls and northern and southern facing slopes due to partial illumination by the sensor [72,73]. Slope movements showing varying deformation patterns were further investigated using TerraSAR-X data applying PSI and SBAS [74]. To assess rock glacier movement rates, DInSAR was successfully applied to determine displacement rates [75][76][77]. Landslide, rock glacier movement rates, and, e.g., surface displacement mapping using DInSAR was determined in order to assess hazards related to GLOFs [78].
DInSAR was also used to measure glacier movement at several study areas [79,80]. However due to snowfall, snowdrift, and melting, which leads to temporal decorrelation, mainly data with short repeat-pass were used [81]. To estimate surface velocities of glaciers, DInSAR using Sentinel-1 data was applied, exhibiting low deformation rates and therefore low decorrelation [82].

Multi-Spectral Satellite Data
The widespread use of multi-spectral satellite data to monitor and characterize high mountain processes started about 30 years ago. Landsat 5 [83], Landsat ETM+ [84], or ASTER [85,86] were the basis of space-borne glacier mapping mainly using the visible and near-infrared section of the spectrum. Creation of glacier inventories is mainly based on multi-spectral satellite data using automatic procedures [87]. For the Alps, the use of multi-spectral data for monitoring high mountain areas was presented, e.g., for the example of Switzerland in a synoptic view by Huggel et al. [88] as well as the perspectives for a worldwide assessment e.g., using satellite data by Gärtner-Roer et al. [19]. A comprehensive review of global glacier characterization using space-borne sensors was presented by Kääb et al. [89]. For characterizing glacier lake dynamics, Landsat 8 images were used by Li and Sheng [90]. To create a nation-wide glacier lake inventory, Sentinel-2 data were the basis for mapping more than 400 glacier lakes in Norway by Nagy and Andreassen [91] or a classification of glacier lakes by Verma and Ghosh [92].
Sentinel-2 is a two-satellite mission: Sentinel-2A was launched on 23 June 2015 and Sentinel-2B on 07 March 2017 [93]. This constellation yields a revisiting time of five days at the equator showing a substantial improvement in temporal resolution to other missions such as Landsat 8 [94]. Sentinel-2 satellites are equipped with a multi-spectral instrument providing 13 bands from the visible light (VIS), near (NIR), and short wave infrared (SWIR) with up to 10 m spatial resolution [95]. Compared to Landsat 8, the spatial resolution of Sentinel-2 (10 m VIS) is nine times better than Landsat-8 (30 m VIS).

Structure from Motion-Unmanned Aerial Vehicles
UAV-based images are processed using structure-from-motion (SfM) photogrammetry [96]. SfM photogrammetry uses images captured from different perspectives, automatically assembled to point clouds using image matching techniques. This matching uses the identification of interest points and is based on the Scale-Invariant Feature Transform algorithm [97]. In combination with multi-view stereo (MVS) techniques, SfM photogrammetry allows simultaneously reconstructing dense 3D models, camera positions, and orientations [98,99]. Currently, SfM-MVS photogrammetry is increasingly used for generating ortho-images and DTMs for different applications (e.g., [100][101][102][103]. All applications differ in both scaling and flight altitude. Thus, to achieve the necessary image resolution for every particular application, different heights over ground are necessary; e.g., for sediment analysis, low height over ground is required (limit of approximately H = 12 m [104]) compared to large scale applications (e.g., H = 100 m for mapping river sections [105]).
UAV applications at high mountain environments are challenging due to the following reasons: (1) limited accessibility to unstable surfaces due to hazardous conditions causes constraints in the acquisition of ground data (for registration and accuracy assessment); (2) changing meteorological conditions (temperature, wind, fog); and (3) surface conditions: e.g., proglacial lake water surface covers large parts of the investigation area, which cannot be used in photogrammetry. Analysis using UAV-based SfM-MVS often lack of information about survey design and image measurement and processing [106][107][108]. In addition, the SfM-MVS approach is somewhat limited, and results such as DTMs and ortho-images can be erroneous [109].
Recently, unmanned aerial vehicles (UAV) have become a well-established platform in acquiring high resolution data, covering various fields of application, e.g., surveying [110][111][112]. Currently, UAV-based glacio-morphological analyses are increasingly conducted, e.g., the monitoring of glacier surface changes [113][114][115][116][117][118]. The use of UAVs is also increasing for monitoring fluvial systems [119]. UAV are the basis for data collection for requested bathymetry of river sections (for numerical modeling) and the observation of morphodynamic processes, such as quantifying erosion and accumulation or stream bank failure [105,120]. On a smaller scale, UAV-based SfM analysis are used to determine: (i) the sediment composition and distribution; and (ii) hydraulic parameters, such as grain roughness, etc. [104,121]. For the description of channel characteristics and follow-up hydraulic processes, however, both are necessary: the geometry (DTM) as a basis for hydraulic modeling and the characteristic grain sizes, grain roughness, etc. as its input parameters. For this reason, it is usually required in river engineering to conduct additional flights of characteristic sections at different flight altitudes.

Automatic Cameras
Due to essential developments in cameras for terrestrial photogrammetry, automatic time-laps cameras have been frequently used in monitoring geomorphological processes in the last decade.
High spatiotemporal resolution, cost efficiency, and independent operating in remote places [122][123][124] are valuable improvements. In contrast to optical remote sensing techniques, cloud coverage is a by far smaller problem. In the Austrian Alps, a yearly average of 60% of pixels are hidden by clouds [125]. The application of automatic time-lapse cameras as a monitoring tool in mountain regions is versatile: selected examples are glacier terminus position [126], snow cover monitoring [127][128][129][130][131], monitoring mass movements [132], supraglacial lakes drainage events [133], snow melt, and vegetation phenology [134].

Study Area: Relevant Pre-Work and Geomorphological Setting
Starting in 1893, Pasterze Glacier was one of the first glaciers constantly monitored in annual measurements, leading to the longest record of measurements of a single glacier [135]. Next to linear and punctual information, the complexity of a glacier system revealed the need for measurements showing more spatial significance such as spatially well distributed measurements. Multi-spectral satellite data were used to quantify changes in glacier extent changes using images from Landsat MSS (1976), Landsat TM (5 during 1984-1992), Landsat ETM+ (2000), and Ikonos (2000) by Hall et al. [136]. Methodological considerations of this work were presented by Hall et al. [137]. DInSAR was only used to determine flow velocity patterns of the Pasterze Glacier by Kaufmann et al. [138] using five ERS-1/2 image pairs between 1995 and 2001. Aerial images were widely used to characterize several processes and impacts: both using multi-temporal ortho-images, Kellerer-Pirklbauer et al. [139] analyzed the influence of supra-glacial debris cover at the Pasterze Glacier tongue and Kaufmann et al. [102] gave a comprehensive analysis of Pasterze Glacier retreat between 2003 and 2009. At Mittlerer Burgstall mountain (MBUG), a first quantification of the large rockfall event in 2007 and possible relations to climate change was given by Kellerer-Pirklbauer et al. [140]. Kaufmann [141] provided a detailed determination using high resolution aerial images. At Pasterze Glacier terminus area, a first assessment of sub-surface ice and glacier lake evolution using TLS and UAV was given by Kellerer-Pirklbauer et al. [142] and Kellerer-Pirklbauer et al. [143].
The study area comprises the catchment area of Pasterze glacier, which includes the maximum extent of Pasterze Glacier of 1851 (Little Ice Age (LIA) glaciation ( Figure 1) [102]. Since the LIA-maximum, Pasterze Glacier has undergone a constant retreat, which accelerates since the 1990s [135]. Pasterze Glacier lost 37% of its area (a decrease from 26.5 to 16.6 km 2 ) and 63% of its volume (a decrease from 3.10 to 1.16 km 3 ) between 1852 and 2012 [135].

•
Connection between Pasterzenboden and Glacier tongue ( Figure 1A(2)): The so-called Hufeisenbruch summarizes a terrain step from approximately 2700 to 2300 m a.s.l. connecting the latter areas.
In 2019, only two small connections are left such as one approximately 50-m-wide serac-area from the so-called Teufelskampkees and one larger connection SW from the MBUG area.  (6)) is mainly composed of mica-schist and therefore large debris covered. There is minor glaciation at the Burgstall area, relief up to 3331 m a.s.l. Upper slopes are very prone to rock falls, lower slope sections are dominated by consolidation after glacier retreat and linear processes (fluvial erosion). Footslope areas are mainly characterized by dead ice degradation and fluvial erosion.
• Glacier lakes ( Figure 1A (7)): Historically, the Pasterze Glacier system includes only the so-called Sandersee Lake since 1958 at approximately 2060 m a.s.l., with a maximum extent in 1979 (121,500 km 2 ). From 2004 onwards, a braided river system evolved into a second lake (Pasterzensee), which was established in 2010. A constant decrease of dead-ice islands sections and therefore subsequently prevailing water surface developed since then. Lake Pasterzensee is a main object of interest in this work.

Application of Earth Observation Techniques at Pasterze Glacier: Data Basis, Spatial and Temporal Variability, Quality Assessment
Based on the technical description presented in Section 2, this section comprises the particular applications of the single monitoring configurations at Pasterze Glacier area indicating: (i) technical specifications of used instruments, and satellite systems; and (ii) monitoring configuration. Furthermore, we present the characteristics of the single datasets and applied approaches (including quality assessment). To ensure interpretability of measurements, monitoring of glacial and proglacial landscape should always be carried out as close as possible to the end of the hydrological year. By definition, the hydrological year ends with the ablation period before the onset of winter (at alpine glaciers mostly in September). Therefore, for the Pasterze Glacier area, at least one measurement is carried out in September designated as annual measurement campaigns.

Terrestrial Laserscanning
Annual TLS measurements have been carried out since 2001 at the scanning positions, Franz-Josef-Höhe and Hofmanns-Hütte (FJH and HH, Figure 1C). Two devices have been used since 2001: from 2001 to 2009 the Riegl LPM-2k system and from 2009 on the Riegl LMS-Z620 system ( Table 1). Technical specifications can be found in [24,49]. A comprehensive overview of previous work including scanning geometries at Pasterze Glacier tongue area is also given in [24]. In the upper part of the Pasterze Glacier area, a massive rockfall at MBUG occurred in 2007 [140]. Therefore, the TLS monitoring network 'Burgstall rock fall area' was established in 2010, comprising the eastern part of the rock fall area of MBUG and the S-face of the Hohe Burgstall Mountain (HBUG, Figure 1B). The scanning position Burgstall (BUG) is located at the eastern margin of the Wasserfallwinkelkees glacier at a bedrock ridge in 2800.34 m a.s.l. TLS acquisition at Burgstall area was conducted in a distance of 500-750 m at MBUG and a distance of 100-300 m at HBUG ( Figure 1B) using six stable, permanent reference points [144]. Thus, the respective ground sampling distances (GSD) at MBUG was 0.25 m (at 650 m distance) and 0.10 m (at 100 m distance) at HBUG. At MBUG, for special analysis such as void size and density for geological interpretation, detail scans with a GSD of 0.15 m were acquired to cover very active rock fall areas.
TLS data processing (e.g., registration) was performed in Riegl RiScan; the rectified point cloud was afterwards exported to Golden Software Surfer 15 to calculate respective DTMs to obtain the area-wide vertical elevation changes and volumetric information in reasonable calculation time. DTMs with different spatial resolutions were calculated: 1 m for the area-wide assessment of vertical surface changes (e.g., for mass balances) and 0.5 m for geomorphological interpretation and process characterization of specific areas of interest [24]. The quality of a TLS measurements is influenced by four factors: instrument calibration, atmospheric conditions, object properties, and scan geometry [145]. As we compare measurements using the same instrument and meteorological conditions are integrated in the processing chain, we focus on the influence of distance measurements and the incidence angle on the accuracy of TLS measurements. Inaccuracies of measurements are assessed by calculating euclidian distances between two point clouds (using CloudCompare) with respect to measurement distance and incidence angle ( Figure 2). Uncertainties of distance measurement show sufficiently small values for incidence angles larger than 50 • : at distances of around 1150 m, measurements on rock walls (mean incidence angle 86 • ) show rather small uncertainties between 0.091 and 0.154 m. The influence of the incidence angle on distance measurements was shown at several stable areas such as 0.042 and 0.049 m (mean incidence angle 50 • ) and 0.188 and 0.203 m (mean incidence angle 50 • ) at two adjacent stable areas in a distance of 377 and 390 m. For the single stable areas, uncertainty values were stable over the observation period ( Figure 2).

Radar Satellite Application-Backscatter
To quantify the extent of Lake Pasterzensee using backscatter information, Sentinel-1 Single Look Complex (SLC) data were processed using Sentinel Application Platform (SNAP). As precise geolocation is a crucial precondition of comparability, precise orbit files were applied. Pixel values were radiometrically calibrated to derive physical units. Speckle filtering and multi-looking was applied to reduce noise with a subsequent SAR-simulation terrain correction for geometric adjustment. To further reduce noise, the mean of VV and VH data was used to delineate lake extents. Thresholds were set for each scene individually between −14 and −17 dB, and lake extents were calculated per year for every scene taken between June and October. Quality assessment mainly covers co-registration of images on one master image to avoid any existing shift to assure spatial comparability.
Being a mountainous study area, Pasterze Glacier area is widely affected by topographical limitations such as shadowing, foreshortening and layover. which limit a threshold selection based on the image histogram.

Radar Satellite Application-DInSAR
The DInSAR analysis carried out at Pasterze Glacier is a pilot study for the applicability of Sentinel-1 DInSAR for surface deformation assessment in the entire Grossglockner area. The analysis is based on 149 SLC images of Sentinel-1A and 1B taken in Interferometric Wide swath (IW) between late spring 2017 (2017-06-04) and late fall 2019 (2019-10-28). It covers swaths of 250 km with a spatial resolution of 5 m by 20 m and incidence angles from 29 • to 46 • with a repeat cycle of six days. Three different DInSAR approaches were investigated and compared. While the underlying data basis is the same in terms of raw input data and temporal resolution, the three methods differ with respect to (aggregated) spatial resolution. Due to decorrelation effects in the winter seasons (snow cover), only summer scenes were used. The area is covered by ascending orbits 44 and 117, and descending orbit 95.
The first approach is based on a SNAP-StaMPS (Stanford Method for Persistent Scatterers)workflow to perform persistent scatterer interferometry (PSI) using the Stanford Method for Persistent Scatterers [146]. Sentinel-1 pre-processing for PSI was performed utilizing the software SNAP. At first, the optimal master scene was selected for ascending and descending orbits, respectively. The sub-swath covering the region of interest was identified and images were split accordingly. After application of an orbit correction, the images were co-registered and interferograms were computed. Potential PS-points were pre-selected by using the amplitude dispersion, and phase stability was estimated using phase analysis. PS pixels were subsequently filtered and dropped if they were too noisy or if they were influenced by neighboring elements. The wrapped phase was then corrected for spatially-uncorrelated look angle errors followed by phase unwrapping. Eventually, the spatially-correlated look angle error was calculated. Compared to regular grid of aggregated SBAS pixels, the PSI point distribution is more irregular with a point distance of roughly 3 m × 14 m.
The second approach tested within this study is the P-SBAS (Parallel SBAS) service of European Space Agency (ESA) Geohazards Thematic Exploitation Platform (GEP). This online service provides an unsupervised implementation of the P-SBAS algorithm, which is a parallel computing implementation of the SBAS approach [147,148]. The spatial resolution is specified with 90 m, yet actual point intervals were measured with approximately 60 × 90 m.
The third approach for assessing surface deformations derived via DInSAR is also based on small baselines and includes the developments in [71,[149][150][151]. For the pre-processing of the Sentinel-1 stack, the joint azimuth shift estimation as in [152] was applied, which is slightly different from the method used in [147]. The analysis was conducted using the Remote Sensing software Graz (RSG [151]) software suite. Results are aggregated to a spatial resolution of 80 × 80 m.
To assess quality and accuracy of the processed DInSAR results, 24 stable areas were manually identified by delineating areas where no deformations are expected ( Figure 1C). Selection of these areas was based on a geomorphological assessment. These areas comprise 19 stable terrain areas (bedrock) and 5 polygons related to artificial structures (buildings and parking lots around TLS scanning position FJH). Areas were chosen considering a trade-off between the size of the polygons and the uniformity with respect to aspect and slope angle. As the area of each single polygon is still comparably small-particularly with respect to the SBAS pixel spacing-deformation results were aggregated the main two categories 'bedrock' and 'infrastructure' in order to provide more robust indication on mean annual (pseudo-)deformation rates. Resulting deformation rates are computed as a linear regression of deformation rates on the date. To quantify variability, we use the standard deviation (SD) of the residuals.
(Pseudo-)deformation rates of the ascending orbit (44) show a higher variability with respect to the (linear) trend across all DInSAR methods and for both the bedrock and the infrastructure cluster (Table 2). This might be explained by the slope exposition, as almost all areas are located on slopes facing southwest. In fact, foreshortening effects are more prominent in ascending orbits than in descending ones. Compared to the bedrock cluster, the standard deviation of residuals is slightly lower for the infrastructure cluster, which is consistent with the expectation that particularly PSI should perform better on PS than on DS [153].
Notably, PSI shows a systematic trend on both orbit directions with a comparably high variability. However, results are consistent across both orbits and both clusters. Arguably, SBAS-based methods seem to be better suited in (high) alpine environments because of their capability to handle DS, too. Both SBAS methods show rather small trends with residual standard deviation of only a few millimeters. Nevertheless one has to keep in mind, that all time series are not optimum as we had to mask out all 'winter-acquisitions' leading to long time periods without data. This has severe negative impacts on, e.g., the removal of atmospheric effects.

Multi-Spectral Satellite Data
Sentinel-2 data were used to enhance inter-annual data availability due to potential high revisiting time. ESA's Sentinel-2 mission consists of two polar-orbiting satellites flying on the same orbit, phased at 180 • to each other. Its wide-swath, high-resolution, multi-spectral imaging satellites comprise 13 spectral bands, with four bands at 10 m, six bands at 20 m, and three bands at 60 m spatial resolution covering the Visible and Near Infra-Red (VNIR) and Short Wave Infra-Red (SWIR). The mission has a high revisit time of five days at the equator for both satellites with an orbital swath width of 290 km.
To map glacier lake extension accurately, respective satellite images should be chosen showing the glacier lake's maximum extent. This is ensured when using images without snow and ice cover to avoid misinterpretations. Furthermore, cloud coverage and shadowed areas due to the relief have to be considered. In detail, 27 Sentinel-2 Top of atmosphere Level 1C images (Table 3) were available in the observation period 2015-2019. Due to a misregistration of more than one pixel (>10 m) of multi-temporal Sentinel-2 acquisitions [154], all scenes had to be co-registered to the acquisition of 2016-08-07, which matched satisfactorily with the corresponding TLS dataset. For co-registration the software SNAP was used, which computes the offset between master and slave images by maximizing the cross-correlation within sub-images. Subsequently, Sentinel-2 data were classified applying the Normalized Difference Water Index (NDWI), using the green and near-infrared bands [88,155,156]. For this analysis, the respective Bands 3 and 8 of Sentinel-2 images were used: For water bodies, the reflectivity of the green light is maximized, while the near infrared reflectivity is typically low and therefore minimized. Water features exhibit positive values, while while soil and vegetation show lower values due to higher reflectivity of near infrared than green light [157]. Thresholding of the resulting NDWI-maps has widely been used [155,158] but the individual values are dependent on the specific application. In delineating glacier lakes, the threshold value must be high enough to distinguish between glacier ice and water but also low enough for the omission of the water pixels [91].

Structure from Motion
Structure from Motion (SfM) using UAV for monitoring glacier and proglacial area was carried out at the tongue of the Pasterze Glacier, the proglacial area of Pasterze Glacier, the Burgstall rockfall area, and the Oberer Pasterzenboden (Figure 1). Data from Burgstall represent a first flight campaign in order to monitor the entire Burgstall mountains for subsequent geological analysis. Data acquired at the Oberer Pasterzenboden are the basis for glacier mass balance measurements.
The glacial/proglacial transition zone was covered by two UAV surveys in September 2016 and June 2019 (Table 4, VB/GL). Using a fixed-wing UAV (a Quest UAV) in the 2016-11 survey, the consumer grade camera Sony α 6000 was used (E 16 mm F2.8 lens, sensor size is 23.5 × 15.6 mm, resolution 6000 × 4000 px). During the campaign in June 2019, topographic conditions ( Figure 1) and platform specifications furthermore necessitated that the southeastern part of the glacier lake was measured using a multi-rotary UAV (DJI Phantom 4. integrated camera (resolution 4000 × 3000 px)). The SfM-MVS photogrammetry processing was based on ground control points (GCPs), which were used for indirect georeferencing of the UAV imagery. GCPs were measured by GNSS solution in real-time kinematic (RTK) mode (EPOSA). The SfM-MVS processing was conducted using the software Agisoft Photoscan (1.2.5 build 2735; 1.3.4 build 5067) with a default key point limit of 40,000 and a maximum of 4000 tie points (sparse point cloud generation). Thereafter, a bundle adjustment and camera self-calibration were conducted followed by a dense point cloud processing.
To independently assess the vertical accuracy of the photogrammetrically processed DTMs, the same geodetic method was used to measure so-called independent check points (ICPs) for validation of geolocation. The mean vertical differences between the DTMs and ICPs was 0.08-0.13 m with a standard deviation of from ±0.12 to ±0.15 m ( Table 4). The resulting ortho-images were visually interpreted and the analysis of DTMs focused on the quantification of vertical changes using DTM differencing.
SfM for glacier river characterization comprised the documentation of the channel evolution downstream in autumn 2018 with regard to estimate the future sediment (bed load) input into the reservoir Margaritze [159]. The investigation area (proglacial river, Figure 1C), between the glacier terminus and the delta area of the Lake Pasterzensee (100 m wide and 1000 m long), was mapped with an UAV (Hexacopter model KR615), equipped with compact camera (ILCE-6000; 6000 × 4000 px; Table 4). In the case of the high image resolution for sediment analysis and the depth of the incised river channel, an additional flight for the inaccessible sections was done at surrounding terrain level (around 15 m above water level).
The calculated DTM serves as a basis for the 1D hydraulic model of the proglacial part of the river Möll. The DTM was generated using the program Agisoft [160] and registered using GCPs measured with a GNSS device with RTK provider (APOS). Some of the GCPs (20-30%) serve as ICPs to quantify uncertainties of the dataset. The resulting sparse point cloud was cleaned up by depth filtering and the removing of points with high reprojection error. Finally, the dense point cloud contains around 478 million points (computation time of 14 days, 10 h). The resulting point cloud shows a point density of around 4000 points/m 2 , a GSD of 1.59 cm/px and a RMSE (X,Y,Z) of 2.55, 4.38, and 2.44 cm (Table 4). This is a significantly higher accuracy than the one, e.g., presented by Vázquez-Tarrío et al. [121] (GSD = 2 cm/px) and more detailed mapping of the topography is shown. Table 4. UAV based surveys for the focus of proglacial river (GR) and valley-bottom/glacier lake assessment (VB/GL). Information about the height over ground (HoG), ground sampling distance (GSD), number of images, the root mean square error (RMSE) in both the horizontal plane (xz) and the vertical plane (z), as well as vertical quality of the DTM are indicated.

Automatic Cameras
In the Pasterze Glacier area, six automatic cameras are installed in order to monitor mainly glaciological processes with a very high temporal resolution and, as a side-effect, to validate other methods qualitatively (configuration: Table 5; location: Figure 1). The main focus of cameras operated by ZAMG (Table 5, 1-5) is the assessment of the snow cover extent. Due to complex terrain, different camera locations were selected for a maximum coverage of the glacier. Next to snow cover, individual cameras were also applied for flow velocity or determination of the glacier outline. The panoramic camera at Franz-Josefs-Höhe area ( Figure 1 and Table 5, 6) was installed in 2010, which is a valuable source for scientific work and has been used for time-lapse studies of the entire glacier tongue. The five cameras operated by ZAMG are Canon EOS1200D digital single lens reflex cameras. The two cameras at the site Pasterze have a EF-S10-22mm f/3.5-4.5 USM lens. The tree other cameras are equipped with a EF-S 18-55mm III lens.

Results: Monitoring of Processes and Landforms
In the following, the quantitative results of every single Earth observation techniques with respect to particular geomorphic process systems are presented. Geomorphic process systems were distinguished into glacial processes (directly at the glacier) and paraglacial processes (in the vicinity of the glacier since the LIA maximum of 1851).

Glacial Processes
At Pasterze Glacier, the glacial process group include glacier lake evolution (fluctuations), and the identification of icebergs.

Glacier Lake Evolution
The evolution of Lake Pasterzensee at the terminus area of Pasterze Glacier was analyzed using all three sensor types: TLS, Sentinel-1 for backscatter analysis, and Sentinel-2 (using different bands for calculating the NDWI). TLS delineated the largest lake extent compared to other methods, except for the year 2015, where the value was even slightly lower than for the same month of the preceding year.
Quantitatively, for Sentinel-1, the lake extent nearly doubled in the time from 2015 to 2019. For TLS and Sentinel-2, the increase of extent was larger due to lower values in 2015 and larger values in 2019 compared to Sentinel-1. The lake areas delineated from Sentinel-1 and -2 are similar for similar epochs with some differences in 2016 and in the early season of 2019 ( Figure 3). Results exhibit an inter-annual pattern, which is visible in each year: Lake extent increases to a certain level, reaches a peak around early September, and decreases again towards the end of the season. This might be an interesting aspect to be further investigated. Results of qualitative assessment revealed that TLS data provided very accurate results in delineating the glacier lake extent: rhe mean GSD of TLS raw data is between 0.12 and 0.54 m [24] with a mean error in single point geolocation of approximately 0.15 m (Figure 2). Data gaps due to shadowing are of minor importance at the dead ice zones at the hillslope area, but point density at the transition between the glacier terminus area and the lake is very small. Summing up, uncertainties affecting the delineation of water bodies using TLS data are much smaller than the increase of the extent of Lake Pasterzensee (Figure 4).
For the period 2015-2019, in total, 27 Sentinel-2 scenes (at least three scenes per year) could be used to delineate the extent of Lake Pasterzensee (Table 3). Sentinel-2 images were classified using NDWI, where thresholds of 0.12-0.4 in the different scenes were used to classify water. NDWI values between 0.09 and 0.23 typically indicate ice; water areas are classified with NDWI values < 0.60 [91]. However, lakes with strong turbidity are classified as glacier surface [93], which is proved at Lake Pasterzensee. Glacier lake extent at Lake Pasterzensee show a clear increase of area in the observation period from 136 × 10 3 (±12 × 10 3 ) to 288 × 10 3 (±8 × 10 3 ) m 2 (Figure 3). To assess uncertainties of automatic glacier lake extent qualitatively, we used: (i) high resolution ortho-images (2018); and (ii) a combination of TLS-based delineation of lake extent and visual interpretation of the automatic camera at FJH (at a 5 min basis) (see Figure 4). Assuming that 0.5 pixels are an average error in image classification, we applied buffers of ±5 m that indicated a mean accuracy of ±34 × 10 3 m 2 for water body classification.
The glacial lake delineation based on Sentinel-1 data was conducted using scenes from the same orbit with a revisiting time of 12 days for the years 2015 and 2016 (Sentinel 1A) and 6 days for the years 2017, 2018, and 2019 (Sentinel-1 A and B). For the study area, between 10 and 25 relevant scenes per year could be analyzed during the snow-free period (June to October) between 2015 and 2019. In total, four acquisitions (in June and October) had to be excluded due to partial freezing of the lake. Although mostly weather independent, wind and rain may cause roughening of the water surface leading to a reduction in contrast between water and land surfaces. Furthermore, the noise-like effect of speckle may decrease the accuracy of the delineated lake extent. To reduce this effect, the mean of the available polarizations VV (vertically transmitted and received radiation) and VH (vertically transmitted and horizontally received radiation) was used and overall threshold values of 14-17 dB were applied. To estimate classification accuracies, the change in lake extent was assessed depending on threshold value. A change in threshold of 0.5 dB led to a change of 3-7% depending on contrast conditions between water and terrain.

Proglacial River
In addition to the challenges of using UAV applications in high mountain environments, there are also problems in recording channel geometries. On the one hand, the recording of moving elements, such as the water surface level, provides blanks in recording and no points in the model, and on the other hand the river bed topography cannot be mapped as there is no detection of the bathymetry below the water surface; a general problem of UAV applications for mapping river geometries. Thus, the river bed geometry must be obtained from terrestrial surveying (e.g., GNSS devices) and will finally be combined with the UAV data to get the whole river geometry for high water stages. Thus, it is inevitable to intersect elements with different point density to a final DTM, linear cross-sections with planar information from UAVs. This 'sampling problem' applies to TLS mapping as well.
These results are calculated based on the UAV derived DTM, whose positional accuracy is determined by a point using a GNSS device (dx = 4 cm, dy = 0 cm). Due to the high mountain environment, just one point was available. Based on the chosen flight altitude, camera and modeling settings, this DTM has a resolution of 1.59 cm/px (Table 4), necessary for the photogrammetric sediment analyses.

Icebergs
Icebergs are a frequent phenomenon at Pasterze Glacier terminus. Contrary to other glaciers, icebergs at Pasterze Glacier seem to develop more from sub-water dead-ice areas than directly from the glacier part itself. Hence, at Pasterze Glacier, we combine analysis of dead-ice areas and icebergs at the glacier lake as dependent processes. The limited lifecycle of icebergs complicate the detection in remote sensing techniques ( Figure 5). TLS and multispectral satellite images only detect icebergs randomly and therefore a detection does not support a substantial analysis. The only method providing valuable information is the use of temporal high-frequent automatic cameras. However, since information is provided only qualitatively (occurrence and trajectories), a quantification of iceberg characteristics (e.g., size) cannot be conducted satisfactorily with the available monitoring system at Pasterze Glacier (using automatic cameras).
TLS provides the most accurate delineation of icebergs however in our case at the lowest temporal frequency. In addition, with Sentinel-2 icebergs are clearly distinguishable, although the resolution of 10 m often leads to mixed pixel information for these oftentimes small and highly dissected features. In many cases, clusters of icebergs may seemingly merge into one or be connected to the lake border. Smaller icebergs-clearly detectable with TLS-may also be misclassified. To some extent, icebergs are also distinguishable with Sentinel-1, however the speckled nature of SAR data often hinders a confident delineation. Furthermore, islands, similar to sandbanks, are hardly distinguishable due to their smooth surface and therefore backscatter values similar to water.   14:00. During this process, the iceberg also rotates around its own vertical axis (C,D). To automatically track the movement of icebergs, high temporal resolution of the data basis is required. Additionally, due to rotational movements of the iceberg, such automatic monitoring must also be able to process changes in the geometry of the iceberg. Images provided by GROHAG.

Valley-Bottom and Slope Processes
At the valley-bottom area of Pasterze Glacier terminus area, TLS was conducted with measurement distances between 300 (dead-ice terraces) and 800 m (debris covered glacier area) to the scanning position. This leads to distance-dependent GSD of 0.12-0.20 m at 300 m and 0.33-0.54 m at 800 m. Depending on the analysis of different process groups, calculated DTMs were resampled to the following spatial resolution: 0.25 m for visual geomorphological interpretation, 0.5 for sectoral, and 1.0 m for area wide calculations [24]. To quantify uncertainties, 13 stable (bedrock) areas beneath the scanning position FJH (Figure 1) were used at distances of 345-1186 m (Figure 2). Calculated uncertainties are 0.083-0.192 m for mean Euclidean distance error and −0.025 to 0.033 m for dz.
TLS based analysis revealed several interesting areas based on the work of [24]. For the period 2018-2019, two areas were exemplarily chosen for further analysis. Figure 6A shows the vertical elevation patterns of a sliding area towards the Lake Pasterzensee. This landform covers an area of 16,750 m 2 with a mean surface subsidence of −0.32 m of the entire area. Within this entire landform, patterns of different vertical surface elevation changes can be detected comprising distinct sub areas such as with maximum dz of −0.92 m Figure 6A(1)), maximum dz of −0.81 m (Figure 6A(2)), and maximum dz of −1.28 m (Figure 6A(3)). Figure 6B shows typical erosion process with downslope channel structure ( Figure 6B(3)) covering an area of 14,500 m 2 . In terms of mass balance, erosional ( Figure 6B(1)) and deposition processes ( Figure 6B(2)) nearly cancel each other out with a calculative dV of 290 m 3 (corresponding to a mean surface subsidence of −0.02 m).
Using Sentinel-1 for DInSAR applications, no significant deformation rates indicating valley bottom or slope processes can be detected in the main study area Pasterze Glacier. However, close to the main study area, DInSAR analysis delivered one particularly interesting process. Ongoing erosional processes are clearly detectable in PSI and SBAS analysis for the landslide Guttal (Figure 7) in the detachment area. We measured deformation values of −6.1 cm/a (away from the satellite) for P-SBAS in orbit 117 and −4.6 cm/a for PSI in orbit 44.

Rock Fall Processes
Rock wall processes are monitored at the Burgstall Mountains mainly using TLS. Quality considerations were conducted similar to the glacier terminus area using stable areas Figure 2.

Discussion
In this chapter, we analyze the applicability and the constraints of single methods used in monitoring particular geomorphological processes at the Pasterze Glacier area.

Glacial Lakes
TLS was conducted once a year at Pasterze Glacier area, only in 2017 inter-annual TLS data were acquired three times. Measurements at Pasterze Glacier area are limited to three regular measurement campaigns due to logistics and as a consequence costs. However, spatial resolution show remarkable results both in high dynamic (e.g., evolution of Lake Pasterzensee) as well as processes showing small magnitudes (e.g., slope processes) depending on the measurement distances [161]. TLS provides very dense point clouds but lacks in reflectivity especially in close vicinity to the Pasterze glacier terminus.
Water saturated sander and dead-ice bodies as well debris covered glacier do not reflect laser impulses sufficiently and therefore margins and transition zones from water to terrain are not represented accurately. Water at Lake Pasterzensee exhibits high turbidity and diurnal variations of water level which influences signal reflection and interpretability of results. Topographical situations leading to unfavorable scanning geometry (small incidence angle close to the glacier terminus) yield different extents or data gaps, which subsequently have to be interpreted carefully.
The selection of appropriate Sentinel-2 images proved to be difficult due to frequent cloud coverage. Therefore, of the total of about 30 potentially available images taken between June and October, only the images in 2015 and eight images in 2018 are suitable for further analysis. However, in terms of assessing the inter-annual variability, Sentinel-2 images do increase the temporal resolution of data availability, e.g., for delineating glacier lakes. Analysis on a one-to two-monthly basis (in the summer period) can be enabled (Figure 3), which is already sufficient for a single test region with a few lakes and would be a valuable support for area wide quantification of glacier lake extent.
The spatial resolution of 10 m of Sentinel-2 images is a valuable improvement for area-wide assessment compared with available multi-spectral data so far (e.g., Landsat 8: 30 m spatial resolution, panchromatic 15 m). However, in detail, small scale dynamics such as degradation of ice-terraces or the dynamics of sander areas are hardly detectable with Sentinel-2 images. Spatial resolution of Sentinel-2 and reflectivity characteristics of the surface cause misclassifications (Figure 4). Especially shallow water areas (Figure 4 A-D(1,5)) are often classified as terrain and small-scale features are not represented accurately ( Figure 4A-D(3)). At Pasterze Glacier, the water level shows distinct diurnal variations of max. 0.5 m especially in August. A comparison of small-scale landforms in shape and size based on TLS and Sentinel-2 data is very problematic due to the rapidly varying water levels ( Figure 4A-D(1-3)).
Sentinel-1 provides data with high temporal resolution with a revisiting time of six days (for Sentinel-1 A and B). Since results are hardly affected by cloud cover, this yields clear benefits compared to optical remote sensing for frequent mapping of glacial lakes [64,65]. Because Sentinel-1 is a side-looking sensor, terrain correction with a high resolution DTM is essential in order to achieve correct geolocation in mountainous terrain. The necessity of an appropriate DTM for glacier lake delineation, which is essential in cases of rapidly developing lakes [162,163], was also pointed out by Strozzi et al. [64] and Wangchuk et al. [65]. In this study, a DTM with a spatial resolution of 10 m was used. DTMs with lower resolution, such as SRTM-3 (3 s, i.e. approximately 90 m spatial resolution) would lead to distortions of several tens of meters.
However, compared to other sensors used for glacial lake delineation, SAR data are subject to higher noise levels. Especially the noise-like effect of speckle may lead to classification uncertainties. In addition, special care must be taken regarding acquisitions with lower contrast, when the water surface is roughened due to wind or rain. In this type of situation, the delineation of lake extent or areas which exhibit lower backscatter due to topographic effects may be difficult. Consequently, manual correction of delineation results is necessary. The challenges of speckle and the effects of wind and waves was also pointed out by Strozzi et al. [64], who therefore also preferred a manual classification of glacial lake outlines. Topographic effects may also lead to classification challenges. The slopes facing away from the satellite naturally exhibit lower backscatter values, as the density of scatter points is lower. Very steep slopes may show low values comparable to water surfaces. We therefore chose an orbit with a steep incidence angle in order to minimize these topographic influences. Furthermore, Sentinel-1 provides the same spatial resolution as Sentinel-2 (10 × 10 m), also leading to mixed pixel information. This mainly affects small island detection in the lake which potentially leads to misclassification. Another source of misclassification are areas of wet sand or wet snow, which exhibit very low backscatter values similar to water [64]. Sander islands (cf. Figure 4A,B) are therefore mostly misclassified due to their low backscatter intensity.
Summarizing glacier lake evolution (at Lake Pasterzensee) is one particular process which benefits from the variety of data availability. TLS and UAV provides information with very high spatial resolution and accuracy in a minimum annual and maximum three times a year temporal resolution. Multi-spectral data Glacier lake extent in late summer season undergoes diurnal variations, which is a crucial indication for subsequent process related data interpretation.

Icebergs
Because of a short lifecycle, icebergs on Lake Pasterzensee cannot be identified with most remote sensing techniques thus far. Icebergs can be detected using TLS due to the high spatial resolution and accuracy. However, to characterize icebergs, a drawback of using TLS is the low temporal resolution (maximum three times per year). Thus, the detection of icebergs and their general occurrence at the time of measurement and the validity of interpretation of long-term is rather random. The usage of multi-spectral is strongly influenced by turbidity and shallow water areas, which delivers mixed-pixel information. Icebergs can be detected due to different reflectivity but delineation causes inaccuracies due to spatial resolution of the multi-spectral satellite data and the small size of icebergs. However, at Pasterze Glacier, the usage of the automatic camera is a very valuable support for geomorphic interpretation. The automatic camera gives qualitative information about the occurrence and disappearing, the movement ( Figure 5), and special dynamics of icebergs as a high dynamic process (e.g., tilting of icebergs [142]). Figure 5 also shows the high dynamics in the movement of several other icebergs within this period of 3 h. Automatic tracking is limited due to sudden rotating and even tilting of the iceberg and accordingly leading to a changing geometry of the iceberg for tracking purposes [143].
To automatically track the movement of this landform, a high temporal resolution is required. In addition, due to the rotational movements of the iceberg in the water, such automatic monitoring must also be able to process changes in the geometry of the iceberg

River Processes
To determine the potential of erosional processes, the quantification of fluvial processes particularly of the proglacial channel system was based on UAVs and TLS. To describe processes accurately, grain sizes of the sediment have to be determined by using sampling information from several pixels [104]. In detail, roughness coefficients are subsequently deduced from calculated DTMs.
To achieve all requirements from a process point of view, the resulting point density of the point cloud has to be in the range of 15-20 points/m 2 . By analyzing statistical parameters (e.g., 2σ) from the vertical elevation differences (dz) of the calculated DTM, correlations between the calculated statistical values and reference samples (from field work) can be observed [121,164]. Consequently characteristic grain sizes can be defined, but no grain size distributions can be determined.
Despite low GSD, the fines fraction is strongly underestimated in both SfM and TLS-based DTMs. Thus, these two applications are limited for characterizing diamictic sediment [121,164]. In [159], the characteristic grain diameter d90 (90% of the particles in the sediment sample are finer than the respective d90 grain size), necessary for hydraulic 1D modeling, could be determined photogrammetrically with the achieved GSD (1.59 cm/px), similar to what Hauer and Pulg [165] obtained during fieldwork. The in-situ determination of any characteristic grain sizes in the proglacial river at the Pasterze, however, is not be possible due to the described characteristics of the torrent.
However, the photogrammetric grain size determination has a major uncertainty: the orientation of each grain may not be determined exactly due to the top view. Thus, the important b-axis for grain size distribution curves probably cannot be measured. Nevertheless, for the determination of large characteristic grain sizes (e.g., d90), photogrammetric evaluation was found to be applicable.
Next to the challenges of using SfM applications in high mountain environments, there are further problems in measuring channel geometries. On the one hand, the recording of, e.g., the water level provides data gaps in the model. To avoid limitations of high water level in assessing river bed topography, low water level situations support area wide determination of bathymetry.
If water level is high and only mono-temporal surveys are possible, the river bed topography must be obtained from terrestrial surveying (e.g., GNSS devices) and will then finally be combined with the SfM-based data to get the entire river geometry. Thus, it is essential to include the information of these linear transverse cross-sections in the 3D information of the SfM-based DTM.
This 'sampling problem' applies to TLS mapping as well. In addition, such high (up to 15 m) and steep slopes, such as the studied canyon-like channel in [159], lead to big 'scan shadows', which, however, can be minimized by changing the position of the TLS device several times. This in turn implies a good accessibility of the channel.

Valley Bottom and Slope Processes
In high mountain environments, slope processes in the vicinity of ongoing glacier retreat show a large variety of process magnitudes [166]. At Pasterze Glacier slope, processes can be distinguished in areas: (i) reworked by downwasting of dead ice; and (ii) areas mainly stabilized, but partly reworked by gullying and debris faces [24].
For this work, valley bottom and slope processes were analyzed exemplarily comprising changes in dead-ice areas, erosional processes, and small landslides ( Figure 6). Due to very fast glacier retreat at Pasterze Glacier these process types show an accentuated temporal succession (e.g., [24]). In the vicinity of the Lake Pasterzensee dead-ice degradation with a vertical subsidence of some decimeters per year are detectable, e.g., using TLS ( Figure 6A). High accuracy of data enable the determination of process patterns leading to the fragmentation of landform due to their activity status ( Figure 6A(1-3)). Erosional process chains showing lower magnitudes are also clearly detectable indicating areas of erosion and areas of deposition ( Figure 6B).
Using multi-spectral images, the information received from quantitative assessment of slope processes is limited. The change of size of dead-ice areas bordering the the glacial lake is the only application which is possible due to significant difference of the spectral information and the spatial resolution of data. Thus, the monitoring of slope processes using multi-spectral surface information is not possible due to the coarse spatial resolution of Sentinel-2 in regard to process magnitudes.
The applicability of DInSAR techniques in characterizing high mountain processes has to be discussed diversely. To measure surface deformations, the approach of using stable areas for testing PSI and SBAS applications provides interesting insights into using different orbits and software applications (StaMPS, GEP, and RSG). However, results do not deliver satisfying results in terms of the required accuracy needed for the characterization of processes at this stage of analysis. Results at slope areas close to the Pasterze Glacier do not show consistent patterns of deformation values, which can be validated with other techniques.
DInSAR analyses show only minor drifts over time. However, deviations of single epochs from the trend line are higher than reported in other publications. This is probably mainly attributable to insufficiently modeled atmospheric turbulences. For the correction of atmospheric phase delays, the methods of Cong et al. [167] are available. However, first experiments with ECMWF (European Centre for Medium-Range Weather Forecasts) ERA5 model parameters confirmed larger inaccuracies of such corrections, especially in the 'turbulent' summer season. Therefore, only annual mean deformation rates were estimated for all DInSAR methods applied in this study. Furthermore, results are based on an assessment in line-of-sight only. This may cause systematic effects in the vertical (up/down) or horizontal (east/west) deformation rate estimation [168].
One particular process, which is clearly detectable with DInSAR techniques, is surface deformation in the upper part of the landslide Guttal, some 4.3 km east of Franz-Josefs-Höhe. This area shows significantly higher deformation values than the surrounding surface in every analysis (cf. Figure 7). Due to the lack of reference data, the validity of the magnitude of the single deformation values cannot be stated in this early stage of using Sentinel-1 data for high mountain applications. There is neither additional information about movement patterns of the Guttal landslide using other techniques, nor recent quantifications in prework [169,170].

Rockfall Processes
For the characterization of rockfall processes, TLS provides point clouds as a valuable basis for further analysis such as geological interpretation. At both test sites (MBUG and HBUG), analysis delivered information about distinct areas of rock fall, larger block falls, and the changes in the accumulation area of the large rock fall of 2007. Rock fall processes produce small to large blocks, thus the minimum detectable object (MDO) of the methods is a crucial factor. Substantially different measurement distances cause different reasonable GSDs leading to different MDOs. At MBUG, the mean GSD varies from 15 (for special areas) to 40 cm for practical reasons (measurement time and data handling). Consequently, rocks larger than 15 cm are potentially detectable in the accumulation area, as well as, interestingly, changes in the detachment area.
Due to small measurement distances (100-250 m), joints and planes of bedding are detectable in data from HBUG. At MBUG, we adapted the scanning increment the achieve the similar GSD in order to measure special areas of rock falls effectively. In the deposition area of the rock fall in 2007, the positive vertical changes deduce rock fall activity and area-wide negative vertical changes are most likely evidence of block consolidation and/or subsidence due to underlying glacier retreat. In some situations, scanning geometry and topography lead to pseudo-deformations due to misalignments. As geolocation of particular point clouds are not equal, breaklines such as ridges are represented differently.
Qualitative assessment of rock fall activity especially at HBUG east-face was conducted using the automatic camera of Fuscherkarkopf (FKK in Figure 1). Visual inspections of the rock faces were conducted during measurement campaigns to avoid misinterpretations due to inaccurate representation of the surfaces in the models. In the last four years, rock fall activity was not limited to the south-face of HBUG. Frequent rock falls were also reported from the east-face affecting the track to the Oberwalder Hütte, which is one of the main training bases of the Austrian Alpine Club. TLS scanning sectors do not cover this areal thus, for subsequent monitoring of possible rock faces, the automatic camera was a valuable support. In August 2019, a first UAV-based campaign was conducted covering both MBUG and HBUG as reference measurements. The quality assessment of the resulting DTM is still in progress. Annual SfM derived DTMs will be the basis for a comprehensive geological analysis in upcoming years.

Synopsis: Processes/Landforms and Data Acquisition
All the consideration made in this section are condensed into a final synopsis of spatial and temporal scales for both methods and processes. To draw conclusions about the applicability of Earth observation techniques for monitoring high alpine environments, both the processes with respect to landforms and the required data are analyzed in spatial and temporal resolution. For this purpose, the processes or consequent landforms are classified with respect to their: (i) persistence of occurrence (lifecycle); and (ii) the extent of occurrence (including the determined morphometric parameter(s)). The data acquisition was quantified with respect to the minimum spatial and temporal resolution in order to monitor these processes and landforms to their changes independently of the sensor specification ( Table 6).
The occurrence of Lake Pasterzensee is persistent over a year. The landform itself and the changes are detectable over the parameter area. Changes over one year are in the order of >100 m 2 and thus data acquisition with a spatial resolution of 10 m is sufficient to monitor the glacier lake accurately. If the temporal resolution is increased, information about the water body size delivers further information about, e.g., dead-ice degradation, iceberg dynamics, and sediment transport.
Icebergs show pronounced dynamics in lifecycle of days to weeks. These fast changes occurring in size and shape are quantified with the parameter area and qualified with movement tracks.
Proglacial rivers are characterized by a high spatiotemporal variability, whereby there may be big differences in the longitudinal direction of supply and sorting of sediments. Near the glacier terminus, lateral channel changes may occur within days to weeks. The already developed and incised part of the channel, however, seems to be stable for months up to a year. Both sections were determined morphometrically over the parameters area (=position), length (=elongation), and depth (dz). For the upper section, the short persistence of occurrence leads to minimum observation periods of approximately one day to one month. In terms of sediment budgets, acquisition should be conducted at least monthly close to the glacier terminus with spatial resolution of 1-2 cm to determine grain sizes (b-axis = second longest one). A similar high temporal resolution is necessary for the measurement of large roughness elements in the developed section.
The persistence of occurrence for dead-ice landforms strongly relates to the distance to the glacier. Close to the glacier, they are characterized by a high temporal variation with lifecycles of about weeks. With increasing distance to the glacier, dead-ice landforms persist for an observation year. The determinate morphometric parameters are area and vertical elevation changes (dz). To detect the area of these landforms, a minimum spatial resolution of 0.5-10 m is necessary (depending on the distance to the glacier). Temporal resolution of one year is sufficient to calculate area-wide characterization of dead-ice areas. Detailed process characterization needs higher temporal resolution due to higher process rates, e.g., lateral melting of dead-ice peninsulas. Qualitative considerations can be done using automatic cameras but they lack quantification so far.
Slope processes were analyzed exemplarily for the processes erosion and landslides. Erosional processes have a lifecycle ranging from a few weeks to a year. They are primarily detectable using the parameter area and vertical elevation changes (dz). A minimum spatial resolution in the horizontal plane for area of 0.5-1 m is necessary to detect process related landforms and changes. A temporal resolution of at least six months is sufficient for area-wide characterization. An increased temporal resolution provides the differentiation of spontaneous and continuous slope processes (supported by automatic cameras). Landslides show a lifecycle of minimum one year and are detectable by the parameters area, dz, and dr. Rock fall is a spontaneous and extreme rapid process. Rock falls range in size from single stones to large failures, involving several 100,000 m 3 (see Rockfall Burgstall). Due to their infrequent nature, it is almost impossible to monitor the process trigger in detail-except with automatic cameras. Nevertheless, the accumulation of debris over time is persistent and thus the process and its changes is detectable over the parameters volume and dz with one year.
Overall, the DInSAR results demonstrate the feasibility to monitor deformation even in the challenging alpine environment. Further analysis is required to better understand the effects of DS and PS characteristics, the influence of slope aspect, and exposition and especially best suited atmospheric corrections that can deal with the very local meteorological conditions and large missing data in the analyzed time series. Table 6. Schematic classification of processes and data acquisition: quantification of lifecycle (persistence), morphometry, spatial resolution, and temporal resolution. Codes: area, area of landform; dz, vertical difference; dr, differential range → euclidean distance.

Conclusions
The characterization of the dynamics of geomorphological processes is always a crucial task in assessing geohazards. The evaluation of the specification, uncertainties, and limitations of monitoring techniques provide valuable information about their applicability. Amongst other aspects, we therefore focused on the comparison of available spatial resolution and feasible, effective temporal resolution.
This article comprises a variety of approaches monitoring geomorphological processes in different magnitudes and temporal as well as spatial scales in the test site Pasterze Glacier area in the observation period 2015-2019. The synoptic usage of data acquired by different remote sensing sensors proves to be a reasonable approach in monitoring geomorphic processes. As different sensors becomes more easily accessible and the usage more frequent in recent years, a closer look at accuracies and effective availability is crucial in terms of interpretability of the geomorphic processes (Table 7). Table 7. Comparison of applied methods for monitoring geomorphological processes in terms of technical specification and quality characteristics: Codes: coverage, spatial extent of dataset; completion time. time period from the survey to the usable dataset; accuracy. both geolocation and pixel value (elevation for TLS and UAV; reflectivity for backscatter; phase for DInSAR reflectivity for multi-spectral; no additional pixel information for automatic camera); TRL, technology readiness level for the application monitoring geomorphic processes. Nearly all applied techniques were able to quantify glacier lake extent sufficiently. TLS and UAV delivered high resolution datasets with very high accuracy. One drawback is the comparable low temporal resolution due to logistics and costs. Automatic classification of multi-spectral data and radar backscatter analysis provided quantifications with a very high temporal resolution but show some problems in accuracy due to spatial resolution. Summarizing glacier lake extent assessment using this pool of techniques delivered very interesting details in lake evolution.
Conducting sediment analysis such as in the example of river bed characterization, the most efficient method seems to be UAV applications due to: (i) low water depth and high geometric (canyon) and roughness features (boulders); and (ii) limiting errors due to lacking penetration of the water surface. Moreover, terrestrial surveying (e.g., TLS) is limited due to: (i) problems with the scanning geometry (shadowing, GSD); (ii) the dangerous characteristics of torrents (e.g., pronounced riverbed structures and high flow velocities); and (iii) the very steep slopes of torrents.
To quantify valley bottom, slope, and rock wall processes, TLS and UAV also provided precise data with very high spatial resolution. Even slow processes with small magnitudes can be detected over the observation period of one year. UAV could be utilized for large areas consistently, whereas TLS shows drawbacks of small incidence angle or shadowing. To cope with these challenges several scanning positions are necessary, which is not possible at Pasterze Glacier area due to accessibility and time effectiveness. DInSAR using Sentinel-1 data is a rather new application in monitoring geomorphic processes. Spatial patterns of deformation values do not deliver explicit information. We assume that accuracy is significantly dependent on atmospheric influences, which are difficult to correct for, particularly in high alpine environments, although larger deformation values of, e.g., landslide processes are clearly discriminable from adjacent terrain.
Today, automatic cameras become a valuable source of information mostly in a qualitative manner. Especially the availability of public, high quality cameras increasingly improve scientific work. Temporal resolution is a major upgrade in order to better understand processes or to validate processes dynamics. Further development of camera system (e.g., resolution) and available software (registration) will enhance the applicability in upcoming years.

Acknowledgments:
We gratefully appreciate the support of the students of University of Graz, Graz University of Technology, and University of Natural Resources and Life Sciences, Vienna. Furthermore, we thank Grossglockner Hochalpenstraßen AG (GROHAG) as a long-term partner for valuable support in terms of logistics, transport and energy supply.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: