Rockfall inventories are essential to quantify a rockfall activity and characterize the hazard. Terrestrial laser scanning and advancements in processing algorithms have resulted in three-dimensional (3D) semi-automatic extraction of rockfall events, permitting detailed observations of evolving rock masses. Currently, multiscale model-to-model cloud comparison (M3C2) is the most widely used distance computation method used in the geosciences to evaluate 3D changing features, considering the time-sequential spatial information contained in point clouds. M3C2 operates by computing distances using points that are captured within a projected search cylinder, which is locally oriented. In this work, we evaluated the effect of M3C2 projection diameter on the extraction of 3D rockfalls and the resulting implications on rockfall volume and shape. Six rockfall inventories were developed for a highly active rock slope, each utilizing a different projection diameter which ranged from two to ten times the point spacing. The results indicate that the greatest amount of change is extracted using an M3C2 projection diameter equal to, or slightly larger than, the point spacing, depending on the variation in point spacing. When the M3C2 projection diameter becomes larger than the changing area on the rock slope, the change cannot be identified and extracted. Inventory summaries and illustrations depict the influence of spatial averaging on the semi-automated rockfall extraction, and suggestions are made for selecting the optimal projection diameter. Recommendations are made to improve the methods used to semi-automatically extract rockfall from sequential point clouds.

This study presents a point cloud de-noising and calibration approach that takes advantage of point redundancy in both space and time (4D). The purpose is to detect displacements using terrestrial laser scanner data at the sub-mm scale or smaller, similar to radar systems, for the study of very small natural changes, i.e., pre-failure deformation in rock slopes, small-scale failures or talus flux. The algorithm calculates distances using a multi-scale normal distance approach and uses a set of calibration point clouds to remove systematic errors. The median is used to filter distance values for a neighbourhood in space and time to reduce random type errors. The use of space and time neighbours does need to be optimized to the signal being studied, in order to avoid smoothing in either spatial or temporal domains. This is demonstrated in the application of the algorithm to synthetic and experimental case examples. Optimum combinations of space and time neighbours in practical applications can lead to an improvement of an order or two of magnitude in the level of detection for change, which will greatly improve our ability to detect small changes in many disciplines, such as rock slope pre-failure deformation, deformation in civil infrastructure and small-scale geomorphological change.

Laser scanning is routinely being used for the characterization and management of rockfall hazards. A key component of many studies is the ability to use the high-resolution topographic datasets for detailed volume estimates. 2.5-Dimensional (2.5D) approaches exist to estimate the volume of rockfall events; however these approaches require rasterization of the point cloud. These 2.5D volume estimates are therefore sensitive to picking an appropriate cell size to preserve resolution while minimizing interpolation, especially for lower volume rockfall events. To overcome the limitations of working with 2.5D raster datasets, surface reconstruction methods originating from the field of computational geometry can be implemented to assess the volume of rockfalls in 3D. In this technical note, the authors address the methods and implications of how the surface of 3D rockfall objects, derived from sequential terrestrial laser scans (TLS), are reconstructed for volumetric analysis. The Power Crust, Convex Hull and Alpha-shape algorithms are implemented to reconstruct a synthetic rockfall object generated in Houdini, a procedural modeling and animation software package. The reconstruction algorithms are also implemented for a selection of three rockfall cases studies which occurred in the White Canyon, British Columbia, Canada. The authors find that there is a trade-off between accurate surface topology reconstruction and ensuring the mesh is watertight manifold; which is required for accurate volumetric estimates. Power Crust is shown to be the most robust algorithm, however, the iterative Alpha-shape approach introduced in the study is also shown to find a balance between hole-filling and loss of detail.

Key to the quantification of rockfall hazard is an understanding of its magnitude-frequency behaviour. Remote sensing has allowed for the accurate observation of rockfall activity, with methods being developed for digitally assembling the monitored occurrences into a rockfall database. A prevalent challenge is the quantification of rockfall volume, whilst fully considering the 3D information stored in each of the extracted rockfall point clouds. Surface reconstruction is utilized to construct a 3D digital surface representation, allowing for an estimation of the volume of space that a point cloud occupies. Given various point cloud imperfections, it is difficult for methods to generate digital surface representations of rockfall with detailed geometry and correct topology. In this study, we tested four different computational geometry-based surface reconstruction methods on a database comprised of 3668 rockfalls. The database was derived from a 5-year LiDAR monitoring campaign of an active rock slope in interior British Columbia, Canada. Each method resulted in a different magnitude-frequency distribution of rockfall. The implications of 3D volume estimation were demonstrated utilizing surface mesh visualization, cumulative magnitude-frequency plots, power-law fitting, and projected annual frequencies of rockfall occurrence. The 3D volume estimation methods caused a notable shift in the magnitude-frequency relations, while the power-law scaling parameters remained relatively similar. We determined that the optimal 3D volume calculation approach is a hybrid methodology comprised of the Power Crust reconstruction and the Alpha Solid reconstruction. The Alpha Solid approach is to be used on small-scale point clouds, characterized with high curvatures relative to their sampling density, which challenge the Power Crust sampling assumptions.

Rockfall constitutes a major threat to the safety and sustainability of transport corridors bordered by rocky cliffs. This research introduces a new approach to rockfall susceptibility modeling for the identification of potential rockfall source zones. This is achieved by developing a data-driven model to assess the local slope morphological attributes with respect to the rock slope evolution processes. The ability to address “where” a rockfall is more likely to occur via the analysis of historical event inventories with respect to terrain attributes and to define the probability of a given area producing a rockfall is a critical advance toward effective transport corridor management. The availability of high-quality digital volumetric change detection products permits new developments in rockfall assessment and prediction. We explore the potential of simulating the conceptualization of slope-scale rockfall susceptibility modeling using computer power and artificial intelligence (AI). We employ advanced 3D computer vision algorithms for analyzing point clouds to interpret high-resolution digital observations capturing the rock slope evolution via long-term, LiDAR-based 3D differencing. The approach has been developed and tested on data from three rock slopes: two in Canada and one in the UK. The results indicate clear potential for AI advances to develop local susceptibility indicators from local geometry and learning from recent rockfall activity. The resultant models produce slope-wide rockfall susceptibility maps in high resolution, producing up to 75% agreement with validated occurrences.

Reality capture technologies, also known as close-range sensing, have been increasingly popular within the field of engineering geology and particularly rock slope management. Such technologies provide accurate and high-resolution n-dimensional spatial representations of our physical world, known as 3D point clouds, that are mainly used for visualization and monitoring purposes. To extract knowledge from point clouds and inform decision-making within rock slope management systems, semantic injection through automated processes is necessary. In this paper, we propose a model that utilizes a segmentation procedure which delivers segments ready to classify and be retained or rejected according to complementary knowledge-based filter criteria. First, we provide relevant voxel-based features based on the local dimensionality, orientation, and topology and partition them in an assembly of homogenous segments. Subsequently, we build a decision tree that utilizes geometrical, topological, and contextual information and enables the classification of a multi-hazard railway rock slope section in British Columbia, Canada into classes involved in landslide risk management. Finally, the approach is compared to machine learning integrating recent featuring strategies for rock slope classification with limited training data (which is usually the case). This alternative to machine learning semantic segmentation approaches reduces substantially the model size and complexity and provides an adaptable framework for tailored decision-making systems leveraging rock slope semantics.