Unstable Rock Slopes Research Articles

Mapping rockfall hazard and detecting precursory damage in rock slopes with passive seismic: Lessons from the La Praz case-study

We conducted repeated successive passive seismic surveys on the La Praz unstable rock slope (∼13,000–15,000 m3, Savoie, France), spanning over more than ten years. The last survey was conducted two months before its complete collapse in August 2023. We used the amplitude of Horizontal to Vertical Spectral Ratio at peak (aHVSR(fp)) to discriminate stable from unstable areas and delineate the fractures with major mechanical decoupling effect. We also classified the compartments in terms of rockfall hazard: the most prone-to-fall areas showed the highest aHVSR(fp) (in the range [26–157]), compared to [7.5–19] on intermediate areas and 2.7 on adjacent stable rock mass.The frequency of HVSR main peak (fp) on three unstable rock compartments showed unambiguous −15%, −14% and −11% drops prior to the collapse. These drops significantly exceed the range of environmentally driven resonance frequency wandering observed on site and generally reported as a limitation in the literature. These drops were interpreted as progressive in-depth slope damage. Interestingly, fp losses were detected solely on the three most prone-to-fall rock compartments showing the highest aHVSR(fp). The unveiling of in-depth damage processes on such a complex, heterogenous and highly fractured rock slope using passive seismic is a first.

Engineering Geological and Geophysical Investigations to Characterise the Unstable Rock Slope of the Sopu Promontory (Gozo, Malta)

Different engineering geological and geophysical investigations were performed at the Sopu promontory in the island of Gozo (Malta), involved in an impressive lateral spreading process due to the superimposition of a stiff limestone (ULC) on a ductile clay (BC). The applied techniques include: traditional geological and engineering geological surveys, unmanned aerial vehicles (UAV) survey, electrical resistivity tomography (ERT) survey, ground-penetrating radar (GPR) investigations, single-station seismic ambient noise measurements, and array seismic ambient noise measurements. The integration of the obtained results allowed us to reconstruct a subsoil model of the promontory that includes features related to the local geology of the slope and to the landslide process, as well as to define a conceptual model that describes the main evolution phases of the expansion process. The presence of back-tilted rock blocks with no features of polarization of Rayleigh waves evidenced the different failure mechanism of the rigid UCL plateau at the Sopu promontory with respect to the Selmun promontory, located in the close island of Malta, where the lateral spreading due to the same geological setting tends to produce unstable rock blocks with a toppling mechanism. This result encourages further future observations and analyses of this topic.

Multimodal Asynchronous Kalman Filter for monitoring unstable rock slopes

Unstable rock slopes pose a hazard to inhabitants and infrastructure in their vicinity, necessitating advanced monitoring methodologies for timely risk assessment and mitigation. Recent geotechnical monitoring techniques often rely on sensor data fusion to enhance forecasting for imminent failures. Our investigation extends beyond a single sensor type to data fusion for heterogeneous sensor networks using a Multimodal Asynchronous Kalman Filter. We illustrate the application of the proposed method on a case study data set consisting of data from an on-site sensor network enriched by remote sensing data. Employing a Multimodal Asynchronous Kalman Filter, we capitalise on the distinct resolutions inherent in each sensor input. The outcome was a combined dataset with a high spatiotemporal resolution. Our approach facilitates the estimation of essential physical attributes for monitored objects, encompassing translation, rotation, velocities and accelerations. The case study site was an unstable rock section of ca. 50.000 m 3 in Aurland, Norway, which collapsed as a multi-stage failure in July 2023. Our method can be transposed to various sites with distinct sensor networks, enhancing state estimations for objects on unstable rock slopes. These estimations can significantly improve applications such as risk assessment and robust early-warning systems, enhancing predictions of critical failure points.

3D Discrete Fracture Network Modelling from UAV Imagery Coupled with Tracer Tests to Assess Fracture Conductivity in an Unstable Rock Slope: Implications for Rockfall Phenomena

The stability of a rock slope is strongly influenced by the pattern of groundwater flow through the fracture system, which may lead to an increase in the water pressure in partly open joints and the consequent decrease in the rock wall strength. The comprehension of the fracture pattern is a challenging but vital aspect in engineering geology since the fractures’ spatial distribution, connectivity, and aperture guide both the water movement and flow quantity within the rock volume. In the literature, the most accepted methods to hydraulically characterise fractured rocks in situ are the single borehole packer test, the high-resolution flow meters for fractures, and the artificial tracer tests performed in boreholes. However, due to the high cost a borehole requires and the general absence of wells along coastal cliffs, these methods may not be appropriate in rockfall-prone areas. In this study, an unsaturated rocky cliff, strongly affected by rockfalls, was investigated by combining kinematic analysis, Discrete Fracture Network (DFN) modelling, and artificial tracer tests. The DFN model and potential rock block failure mechanisms were derived from high-resolution 3D virtual outcrop models via the Structure from Motion (SfM) photogrammetry technique. An artificial tracer was injected using a double ring infiltrometer atop the recharge zone of the slope to determine the infiltration rate and validate the DFN results. The DFN and tracer test methods are frequently used at different spatial scales and for different disciplines. However, the integration of digital photogrammetry, DFN, and tracer tests may represent a new step in rockfall and landslide studies. This approach made possible the identification of groundwater flow patterns within the fracture system and revealed about a 10-day tracer transit time from the injection area and the monitored slope, with similar conductivity values gathered from both the DFN and tracer test. Planar and wedge failures with volumes ranging from 0.1 and 1 m3 are the most probable failure mechanisms in the areas. The results were consistent with the delay between the intense rainfall and the slope failures previously documented in the study area and with their mechanisms.

Structural Characterization of a Toppling Rock Slab From Array‐Based Ambient Vibration Measurements and Numerical Modal Analysis

AbstractAccurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array‐based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross‐correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500‐m‐long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

Permafrost in monitored unstable rock slopes in Norway – new insights from temperature and surface velocity measurements, geophysical surveying, and ground temperature modelling

Abstract. The warming and subsequent degradation of mountain permafrost within alpine areas represent an important process influencing the stability of steep slopes and rock faces. The unstable and monitored slopes of Mannen (Møre and Romsdal county, southern Norway) and Gámanjunni-3 (Troms and Finnmark county, northern Norway) were classified as high-risk sites by the Norwegian Geological Survey (NGU). Failure initiation has been suggested to be linked to permafrost degradation, but the detailed permafrost distribution at the sites is unknown. Rock wall (RW) temperature loggers at both sites have measured the thermal regime since 2015, showing mean rock surface temperatures between 2.5 and −1.6 ∘C depending on site and topographic aspect. Between 2016 and 2019 we conducted 2D and 3D electrical resistivity tomography (ERT) surveys on the plateau and directly within the rock wall back scarp of the unstable slopes at both sites. In combination with geophysical laboratory analysis of rock wall samples from both sites, the ERT soundings indicate widespread permafrost areas, especially at Gámanjunni-3. Finally, we conducted 2D thermal modelling to evaluate the potential thermal regime, along with an analysis of available displacement rate measurements based on Global Navigation Satellite System (GNSS) and ground- and satellite-based interferometric synthetic aperture radar (InSAR) methods. Surface air and ground temperatures have increased significantly since ca. 1900 by 1 and 1.5 ∘C, and the highest temperatures have been measured and modelled since 2000 at both study sites. We observed a seasonality of displacement, with increasing velocities during late winter and early spring and the highest velocities in June, probably related to water pressure variations during snowmelt. The displacement rates of Gámanjunni-3 rockslide co-vary with subsurface resistivity and modelled ground temperature. Increased displacement rates seem to be associated with sub-zero ground temperatures and higher ground resistivity. This might be related to the presence of ground ice in fractures and pores close to the melting point, facilitating increased deformation. The study demonstrates and discusses the possible influence of permafrost, at least locally, on the dynamics of large rock slope instabilities.

Coseismic Stability Assessment of a Damaged Underground Ammunition Storage Chamber Through Ambient Vibration Recordings and Numerical Modelling

In the past decade, ambient vibration measurements found numerous applications on unstable rock slopes and developed into a powerful tool for site characterization of slope instabilities. In this study, for the first time ambient vibration measurements were applied to a rock mass strongly disturbed and damaged by subsurface explosions. The site above the ammunition storage chamber at Mitholz (Switzerland) is especially interesting because the subsurface geology below the seismic array is well known, including the location of the caverns, and the degree of degradation caused by the subsurface explosions in 1947 of around 40 t TNT of ammunition. Measurement data were analyzed using current state-of-the-art seismic single-station and array methods, focusing on surface-wave dispersion analysis, wave field polarization, wave amplification using site-to-reference spectral ratios and analysis of normal mode behavior. The results allow for calibrating the elastic properties of a 2D numerical rock mechanical model which was used to simulate the stability of the disturbed rock mass during seismic loading. Therefore, ambient vibration measurements can contribute not only to a better understanding of the subsurface, but also for the assessment of earthquake risk.

Monitoring the changing seismic site response of a fast-moving rockslide (Brienz/Brinzauls, Switzerland)

SUMMARYSeismic measurements on unstable rock slopes are a complementary tool to surface displacement surveys to characterize and monitor landslides. A key parameter is seismic amplification, which tends to scale with the degree of rock mass degradation. Amplification also provides a direct measure of how the wavefield is intensified during seismic loading, eventually leading to coseismic failure. Here we present the dynamic response of the fast-moving Brienz/Brinzauls rock slope instability in Switzerland (10 $ \times $ 106 to 25 $ \times $ 106 m3), which threatens settlements and infrastructure in the area. The rockslide shows strong seismic amplification at two resonant frequencies with factors of up to 11 and wavefield polarization influenced by the local fracture network orientation. We monitored the dynamic response over a period of 30 months using ambient vibrations and regional earthquake recordings. We observed a change in wavefield polarization of up to 50°, coinciding with a rotation of the relative surface displacement vector field measured by geodetic systems, highlighting the linkage between wavefield polarization and stress field (i.e. rock mass kinematics). For the analysis of secondary, relative surface displacements, we propose a singular value filtering of the displacement field to remove the principal component of landslide motion. In addition, we found increased seismic amplification values after periods of strong precipitation, providing empirical field evidence that the local precipitation history is a key parameter for assessing the hazard of earthquake-induced slope failure.

Monitoring the Preonzo Rock Slope Instability Using Resonance Mode Analysis

AbstractReliable monitoring of unstable rock slopes is a prerequisite for successful mitigation of landslide hazards. However, most state‐of‐the art techniques rely on measuring the local surface displacement in the potential release area. In contrast, recording ambient vibration data allows for analyzing structural dynamic parameters of the unstable slope, such as resonance frequency, polarization of vibration, and energy dissipation. These parameters can be linked to properties of the instability, for example, to rock stiffness and fracture network orientation. We developed a processing method for continuous seismic data based on enhanced frequency domain decomposition modal analysis and applied it to the unstable rock slope Preonzo in Switzerland (∼140,000 m3). Four years of ambient vibration data recorded at two permanent seismometers on the instability were analyzed, providing the resonance frequency, damping ratio, and normal mode shapes of the fundamental (∼3.5 Hz) and the first higher (∼4.2 Hz) vibrational mode. We found that modal analysis can be reliably used to monitor the dynamic response of an unstable rock slope. We observed annual changes of all parameters with a damping ratio varying between 6.0% and 9.7% for the fundamental mode. The dynamic parameters appear to be primarily driven by temperature and only secondarily by opening and closing of fractures. No large slope failure was registered during the observation period. However, the data provide a baseline model for ongoing slope monitoring to recognize structural changes before a future collapse. The setup proposed builds a complementary monitoring system to displacement‐based surveying.

Movements, failure and climatic control of the Veslemannen rockslide, Western Norway

On September 5, 2019, the Veslemannen unstable rock slope (54,000 m3) in Romsdalen, Western Norway, failed catastrophically after 5 years of continuous monitoring. During this period, the rock slope weakened while the precursor movements increased progressively, in particular from 2017. Measured displacement prior to the failure was around 19 m in the upper parts of the instability and 4–5 m in the toe area. The pre-failure movements were usually associated with precipitation events, where peak velocities occurred 2–12 h after maximum precipitation. This indicates that the pore-water pressure in the sliding zones had a large influence on the slope stability. The sensitivity to rainfall increased greatly from spring to autumn suggesting a thermal control on the pore-water pressure. Transient modelling of temperatures suggests near permafrost conditions, and deep seasonal frost was certainly present. We propose that a frozen surface layer prevented water percolation to the sliding zone during spring snowmelt and early summer rainfalls. A transition from possible permafrost to a seasonal frost setting of the landslide body after 2000 was modelled, which may have affected the slope stability. Repeated rapid accelerations during late summers and autumns caused a total of 16 events of the red (high) hazard level and evacuation of the hazard zone. Threshold values for velocity were used in the risk management when increasing or decreasing hazard levels. The inverse velocity method was initially of little value. However, in the final phase before the failure, the inverse velocity method was useful for forecasting the time of failure. Risk communication was important for maintaining public trust in early-warning systems, and especially critical is the communication of the difference between issuing the red hazard level and predicting a landslide.

Reclassification of Microseismic Events through Hypocenter Location: Case Study on an Unstable Rock Face in Northern Italy

Passive seismic methods are increasingly used in monitoring unstable rock slopes that are likely to cause rockfalls. Event classification is a basic step in microseismic monitoring. However, the classification of events generated by the propagation of fractures and rockfalls is still uncertain due to their similar features in the time and frequency domains. Hypocenter localization might be a powerful tool to distinguish events generated by fracture propagation from those caused by rockfalls. In this study, a classification procedure based on hypocenter location was validated using a selected subset of high-quality data recorded by a five-geophone network installed on a steep rock slope in Northern Italy. Considering the complexity and heterogeneity of the rock mass, a 3D velocity model that was derived from a tomographic experiment was used. We performed the localization using the equal differential time method. The location results fairly fit our expectations on suspected rockfall events because most signals were located near the rock face. However, only 4 out of 20 suspected fracture events were unquestionably confirmed as fractures being located inside the rock mass and far enough from the rock face. Further improvements in location accuracy are still necessary to distinguish suspected fracture events located close to the rock face from rockfalls. This study demonstrates that hypocenter location is a promising method to improve the final classification of microseismic events.

Evolution of events before and after the 17 June 2017 rock avalanche at Karrat Fjord, West Greenland – a multidisciplinary approach to detecting and locating unstable rock slopes in a remote Arctic area

Abstract. The 17 June 2017 rock avalanche in the Karrat Fjord, West Greenland, caused a tsunami that flooded the nearby village of Nuugaatsiaq and killed four people. The disaster was entirely unexpected since no previous records of large rock slope failures were known in the region, and it highlighted the need for better knowledge of potentially hazardous rock slopes in remote Arctic regions. The aim of the paper is to explore our ability to detect and locate unstable rock slopes in remote Arctic regions with difficult access. We test this by examining the case of the 17 June 2017 Karrat rock avalanche. The workflow we apply is based on a multidisciplinary analysis of freely available data comprising seismological records, Sentinel-1 spaceborne synthetic-aperture radar (SAR) data, and Landsat and Sentinel-2 optical satellite imagery, ground-truthed with limited fieldwork. Using this workflow enables us to reconstruct a timeline of rock slope failures on the coastal slope here collectively termed the Karrat Landslide Complex. Our analyses show that at least three recent rock avalanches occurred in the Karrat Landslide Complex: Karrat 2009, Karrat 2016, and Karrat 2017. The latter is the source of the abovementioned tsunami, whereas the first two are described here in detail for the first time. All three are interpreted as having initiated as dip-slope failures. In addition to the recent rock avalanches, older rock avalanche deposits are observed, demonstrating older (Holocene) periods of activity. Furthermore, three larger unstable rock slopes that may pose a future hazard are described. A number of non-tectonic seismic events confined to the area are interpreted as recording rock slope failures. The structural setting of the Karrat Landslide Complex, namely dip slope, is probably the main conditioning factor for the past and present activity, and, based on the temporal distribution of events in the area, we speculate that the possible trigger for rock slope failures is permafrost degradation caused by climate warming. The results of the present work highlight the benefits of a multidisciplinary approach, based on freely available data, to studying unstable rock slopes in remote Arctic areas under difficult logistical field conditions and demonstrate the importance of identifying minor precursor events to identify areas of future hazard.

Semi-empirical prediction of dam height and stability of dams formed by rock slope failures in Norway

Abstract. Based on an inventory of 69 dams formed by rock slope failures in southwestern Norway and published inventories from other parts of the world, we developed semi-empirical relationships linking the maximum dam height (HD.max in metres) to dam volume (VD in 106 m3) and other relevant parameters such as valley width (WV in metres) or dam area (AD in square kilometres). Power laws are obtained for HD.max=f(VD) and HD.max=f(VD, WV), while a linear relationship links HD.max to the ratio VD∕AD. For dams in southwestern Norway, the linear relationship HD.max=1.75×VD/AD has the least uncertainties and provides the best results when comparing predicted dam heights with a validation dataset composed of existing dams in northern Norway and numerically modelled dams for possible rock slope failures. To assess the stability of future dams, we use the predicted dam heights in the dimensionless blockage index (DBI) and relating this index to the probability of dam failure derived from our dataset and other published databases on landslide dams. This study underlines the potential of semi-empirical relationships for assessing dam height and stability that needs to be included in preliminary hazard and risk assessment for unstable rock slopes, because damming of a river is an important secondary effect of landslides due to upstream flooding and possible outburst floods in the case of dam failure.

Measurements and modeling of the post-failure micro-deformations and tilts of the Preonzo unstable slope, Alpe di Roscioro, Switzerland

The Alpe di Roscioro, is a site of an unstable rock slope located above the village of Preonzo in southern Switzerland. There have been numerous rock-slope failures with a major event that happened in May 2012. After that event, several seismic stations were set up at the site of the remaining highly fractured and unstable rock mass. The analysis of continuous seismic recordings has shown a high number of permanent micro-deformations and tilts, recorded by a seismometer located on top of the unstable part of the slope. In terms of the ground motion, the directions of these disturbances are parallel to the mean strike of the mapped fracture network, implying a connection to the ongoing deformations of the rock mass. The static deformation field of the fractured rock mass is modeled by finite-difference method (FD). The FD method enables us to apply a reciprocal approach by tilting the affected body of block 1, closest to the stable massif, by measured values and observing the deformations below the seismometer and in the neighboring blocks. The numerical model explains the following features observed in situ: feature 1. The combination of the horizontal displacement and tilt observed approximately at the same measurement point; feature 2. The block 1 center of rotation is shifted towards SSE by ~10 m off its center; feature 3. The sensors placed on the stable massif did not record any micro-deformation and tilt signals above the background noise.In the light of the numerical modeling, the weak coupling between the micro-tilt active block 1 and neighboring media explains the absence of detected events in the reference stations placed on the stable massif. The simplified “shear and bending” force system FD model, indicates that the aftermath of the May 2012 main failure consisted of numerous episodic elementary relaxations of the SSE part of block 1 after the collapse of its NNW compartment, which was previously connected to it.

Monitoring and analysis of active rockslide-glacier interactions (Moosfluh, Switzerland)

Valley glaciers have traditionally been expected to significantly influence the stability and movement rates of adjacent paraglacial landslides. However, detailed studies related to the mechanical and displacement interactions between glacier ice and unstable rock slopes are very rare. Here we present a detailed in-situ investigation of the spatial variations of the displacement field of the Great Aletsch Glacier in the surroundings of a large active instability, i.e., the Moosfluh Landslide, a Deep-Seated Gravitational Slope Deformation, with superimposed large (1–5 million m3) secondary rockslides formed during fall 2016. We performed repeat UAV-based photogrammetric surveys during ~3 days (74 h) and applied Digital Image Correlation techniques to record high-resolution surface displacement vector fields of the landslide, stable slopes and adjacent glacier. Our results show that the secondary rockslide adjacent to the glacier is composed of two parts of 1.5 and 2.8 million m3 volume respectively, showing significant differences in mean displacement velocities (0.4 and 0.9 m in 74 h respectively, excluding rapid movements from detached blocks). Both rockslide compartments induce clear deflections of the glacier flow field, moving with a maximum velocity of about 0.3 to 0.4 m in 74 h. This influence is highest near the ice-contact boundary and decreases within a distance of about 100–200 m from the rock slope instability. We investigate the viscous forces at the landslide/glacier contact using a straightforward analytical model for an incompressible rockslide block sliding along a planar, cohesionless surface into ice. These forces are then applied to a slope stability model based on the limit equilibrium concept, representing the real geometry at the interface boundary to quantitatively explore the true buttressing effects of valley glaciers on an already moving slope instability. We show that the viscous ice deformation plays an important role in mediating the displacement velocities of landslides in unstable conditions, while, on the other hand, the slope support from the valley glacier has very little influence on the stability of the investigated rockslides.As most valley glaciers are currently strongly retreating due to global warming, uncovering significant numbers of pre-LIA slope instabilities, this detailed investigation provides important hints on their potential displacement behavior. Understanding the factors controlling landslide velocity is of great importance for hazard analyses and early warning. Hence, this study has implications beyond academic interest, e.g., for the planning and operation of alpine infrastructure, such as cable cars or hydropower systems.

Tomographic Experiments for Defining the 3D Velocity Model of an Unstable Rock Slope to Support Microseismic Event Interpretation

To monitor the stability of a mountain slope in northern Italy, microseismic monitoring technique has been used since 2013. Locating microseismic events is a basic step of this technique. We performed a seismic tomographic survey on the mountain surface above the rock face to obtain a reliable velocity distribution in the rock mass for the localization procedure. Seismic travel-time inversion showed high heterogeneity of the rock mass with strong contrast in velocity distribution. Low velocities were found at shallow depth on the top of the rock cliff and intermediate velocities were observed in the most critical area of the rock face corresponding to a partially detached pillar. Using the 3D velocity model obtained from inversion, localization tests were performed based on the Equal Differential Time (EDT) localization method. The results showed hypocenter misfits to be around 15 m for the five geophones of the microseismic network and the error was significantly decreased compared to the results produced by a constant velocity model. Although the localization errors are relatively large, the accuracy is sufficient to distinguish microseismic events occurring in the most critical zone of the monitored rock mass from microseismic events generated far away. Thus, the 3D velocity model will be used in future studies to improve the classification of the recorded events.

Probabilistic hazard analysis of impulse waves generated by multiple subaerial landslides and its application to Wu Gorge in Three Gorges Reservoir, China

Subaerial landslides can generate impulse waves in reservoirs, fjords, and other water bodies and cause severe damage to the shipping industry, infrastructure, and human communities along the shorelines. Compared with deterministic hazard models, probabilistic hazard analysis is an important tool for assessing the intensity and exceedance probability of impulse waves generated by multiple landslide sources, which are essential for evidence-based risk mitigation and contingency planning. This paper proposes a methodological procedure for systematic probabilistic hazard analysis of impulse waves generated by subaerial landslides considering all relevant landslide sources and multiple uncertainties (i.e., landslide mechanisms, triggering scenarios, and material parameters). The main steps include: (1) definition of the source parameters by parameterizing the epistemic and aleatory uncertainties. In this respect, the logic tree, event tree, and Monte Carlo models were used to organize the uncertainties of landslide mechanism models, parametric scenarios, and material parameters, respectively; (2) estimation of generated waves based on a landslide dynamic analytical model and previously validated empirical wave models derived from prototype scaled experiments; and (3) aggregation of the hazard analysis results in the form of hazard maps and probability maps under parametric scenarios and a specific scenario. Considering the Wu Gorge in the Three Gorges Reservoir as a case study, soil slides and unstable rock slopes were analyzed to quantify the wave hazard along the reservoir under parametric and specific scenarios while considering rainfall and water level fluctuations. The results show considerable wave hazard in the Wu Gorge. For a 20-year return period, the maximum expected height of propagation waves is 14.41 m. The probability of experiencing propagation waves higher than 2 m in 20 years is over 50% in the river segment with high landslide density. The results of the probabilistic wave hazard analysis could be aggregated and disaggregated to accommodate specific requirements of wave risk management.

Geological Survey and Unstable Rock Block Movement Monitoring of a Post-Earthquake High Rock Slope Using Terrestrial Laser Scanning

Local topography and structural features of discontinuities have a great effect on the movement characteristics and failure mechanisms of the unstable rock blocks of the rock slope. By taking full advantage of the geometric information contained in the multi-temporal terrestrial laser scanning (TLS) point clouds, a technique based on a roto-translation method has been developed to track and monitor the movement behavior of blocks in unstable rock slopes. The technique can be completed in three stages: (a) a preliminary change detection is used to distinguish the unstable area by means of a shortest distance (SD) algorithm; (b) an automatic discontinuity identification algorithm is implemented to visually identify the discontinuity sets and to exactly delineate rock blocks in unstable areas; and (c) finally a roto-translation movement monitoring method is applied to track and monitor the unstable rock blocks. This technique permits the user to visually identify the discontinuity sets and to determine their orientation, and therefore provides detailed information of geometric compounding relationships between the slope and joint systems as well as the forming rock blocks which are the key impacts to stability of the rock slope. It permits the user to exactly and visually delineate the rock blocks in the unstable area, and thus the completely tracking and monitoring of lock block movement in unstable rock slopes can be effectively conducted. It overcomes the limited precision of the single TLS points and provides an actual 3D movement measuring method for the blocks rather than a shortest distance. The technique was first used in the movement behavior monitoring of the unstable rock blocks in a post-earthquake high rock slope. The results show that the rock block movement monitoring technique is well suited for providing high-quality data in the assessment of rock failure hazards.

Reliability and Uncertainties of the Analysis of an Unstable Rock Slope Performed on RPAS Digital Outcrop Models: The Case of the Gallivaggio Landslide (Western Alps, Italy)

A stability investigation based on Digital Outcrop Models (DOMs) acquired in emergency conditions by photogrammetric surveys based on Remote Piloted Aerial System (RPAS) was conducted on an unstable rock slope near Gallivaggio (Western Alps, Italy). The predicted mechanism of failure and volume of the unstable portion of the slope were successively verified on the DOMs acquired after the rockfall that effectively collapsed the May 29th, 2018. The comparison of the pre- and post-landslide 3D models shows that the estimated mode of failure was substantially correct. At the same time, the predicted volume of rock involved in the landslide was overestimated by around 10%. To verify if this error was due to the limited accuracy of the models georeferenced in emergency considering only the Global Navigation Satellite System/Inertial Measurement Unit (GNSS/IMU)-information of RPAS, several Ground Control Points (GCPs) were acquired after the failure. The analyses indicate that the instrumental error in the volume calculation due to the direct-georeferencing method is only of the 1.7%. In contrast, the significant part is due to the geological uncertainty in the reconstruction of the real irregular geometry of the invisible part of the failure surface. The results, however, confirm the satisfying relative accuracy of the direct-georeferenced DOMs, compatible with most geological and geoengineering purposes.

Innovative Directional Controlled Blasting Technique for Excavation of Unstable Slopes Along a Busy Transportation Route: a Case Study of Konkan Railway in India

Many cases of derailments, train accidents and traffic interruptions due to the slope failure and rockfall have been reported in the past on Konkan Railway trackline because of unstable rock slopes. The authorities installed stability enhancement measures such as wire netting, rock bolting and shotcreting which were found insufficient. The failed geotechnical measures prompted authorities to excavate the unstable slopes. Rock blasting is generally not preferred for the excavation of slopes along the transportation route. Stabilisation and widening of slopes by blasting need closure of the route for a considerable time. Discontinuing traffic on Konkan Railway route was undesirable as it connects two important port cities Mumbai and Mangalore. The paper presents an innovative directional controlled blasting (DCB) technique for excavating three unstable slopes (i.e. cuttings) in hilly terrain without disturbing the traffic. New empirical formulae for small geometry blasting to estimate burden and throw have been developed. A unique sequence of firing and excavation with powder factor varying between 0.04 and 0.39 kg/m3 irrespective of the actual free face was used for rock blasting. The process controlled the throw of blasted material over the track and resulted in a smooth and stable wall. The slope angle has been changed from 80 to 82° to 45–47° and probability of wedge, planar and toppling failures was eliminated. About 100 such accident-prone slopes (52 km cumulative length, 7% of total route) situated at discrete locations along Konkan Railway route were stabilised using the technique. No incident of rockfall and slope failures along these slopes has been observed after the stabilisation.