Spaceborne radar interferometry is a powerful tool to characterize landslides at local and regional scales. However, its application to very slow rock slope deformations in alpine environments (displacement rates < 5 cm/year) remains challenging, mainly due to low signal to noise ratio, atmospheric disturbances, snow cover effects, and complexities resulting from heterogeneous displacement in space and time. Here we combine SqueeSAR^TM data, targeted multi-temporal baseline DInSAR, GPS data, and detailed field morpho-structural mapping, to unravel the kinematics, internal segmentation, and style of activity of the Mt. Mater deep-seated gravitational slope deformation (DSGSD) in Valle Spluga (Italy). We retrieve slope kinematics by performing 2D decomposition (2D InSAR) of SqueeSAR^TM products derived from Sentinel-1 data acquired in ascending and descending orbits. To achieve a spatially-distributed characterization of DSGSD displacement patterns and activity, we process Sentinel-1 A/B images (2016-2019) with increasing temporal baselines (ranging from 24-days to 1-year) and generate several multi-temporal interferograms. Unwrapped displacement maps are validated using ground-based GPS data. Interferograms derived with different temporal baselines reveal a strong kinematic and morpho-structural heterogeneity and outline nested rockslides and active sectors, that arise from the background displacement signal of the main DSGSD. Seasonal interferograms, supported by GPS displacement measurements, reveal non-linear displacement trends suggesting a complex response of different slope sectors to rainfall and snowmelt. Our analyses clearly outline a composite slope instability with different nested sectors possibly undergoing different evolutionary trends towards failure. The results herein outline the potential of a targeted use of DInSAR for the detailed investigation of very slow rock slope deformations in different geological and geomorphological settings.

The reliable in situ quantification of rock mass fracturing and engineering quality is critical for slope stability, surface mining and rock engineering applications, yet it remains difficult due to the heterogeneous nature of fracture networks. We propose a method to quantify and map the slope-scale geomechanical quality of fractured rock masses using infrared thermography (IRT). We use the Mt. Gorsa quarry (Trentino, Italy) as a field laboratory to upscale a physics-based approach, which was developed in the laboratory, to in situ conditions, including the effects of fracture heterogeneity, environmental conditions and IRT limitations. We reconstructed the slope in 3D using UAV photogrammetry, characterized the rock mass quality in the field at selected outcrops in terms of the Geological Strength Index (GSI) and measured their cooling behavior through 18h time-lapse IRT surveys. With ad hoc field experiments, we developed a novel procedure to correct IRT data in outdoor environments with complex topography. This allowed for a spatially distributed quantification of the rock mass surface cooling behavior in terms of a Curve Shape Parameter (C_SP). Using non-linear regression, we established a quantitative C_SP-GSI relationship, which allowed for the C_SP to be translated into GSI maps. Our results demonstrate the possibility of applying infrared thermography to the slope-scale mapping of rock mass fracturing based on a physics-based experimental methodology.

The fracturing state of rocks is a fundamental control on their hydro-mechanical properties. It can be quantified in the laboratory by non-destructive geophysical techniques that are hardly applicable in situ, where biased mapping and statistical sampling strategies are usually exploited. We explore the suitabilty of infrared thermography (IRT) to develop a quantitative, physics-based approach to predict rock fracturing starting from laboratory scales and conditions. To this aim, we performed an experimental study on the cooling behaviour of pre-fractured gneiss and mica schist samples, whose 3D fracture networks were reconstructed using Micro-CT and quantified by unbiased fracture abundance measures. We carried out cooling experiments in both controlled (laboratory) and natural (outdoor) environmental conditions and monitored temperature with a thermal camera. We extracted multi-temporal thermograms to reconstruct the spatial patterns and time histories of temperature during cooling. Their synthetic description show statistically significant correlations with fracture abundance measures. More intensely fractured rocks cool at faster rates and outdoor experiments show that differences in thermal response can be detected even in natural environmental conditions. 3D FEM models reproducing laboratory experiments outline the fundamental control of fracture pattern and convective boundary conditions on cooling dynamics. Based on a lumped capacitance approach, we provided a synthetic description of cooling curves in terms of a Curve Shape Parameter, independent on absolute thermal boundary conditions and lithology. This provides a starting point toward the development of a quantitative methodology for the contactless in situ assessment of rock mass fracturing.