Evaluation of maximum rockfall dimensions based on probabilistic assessment of the penetration of the sliding planes into the slope

There is intrinsic difficulty in the investigation of the largest volume of rockfalls that is expected in an area, which lies in the small number of large events, in registrable times. The maximum credible rockfall size has been associated with the properties of the rock mass discontinuities, as they delimit detachable rock blocks, and in particular with the penetration of those discontinuities that comprise rockfall sliding planes. In highly fractured rock masses, the evaluation of the penetration remains an issue. A probabilistic methodology is proposed, to measure the penetration of potential sliding planes into the interior of a rocky slope. The main hypothesis of the method is that the sliding plane persistence is interrupted along its two directions, at the intersection with two lateral discontinuity sets, as the latter displaces the former. Due to the displacement, the sliding planes are formed by quasi-planes that contain a maximum number of spacings of the intersecting joints, hence their size is restricted. The methodology requires as an input the spacing of the intersecting joint sets. Its application to a granodiorite slope confirms that for the study site, there is a maximum volume of rockfalls, excluding the possibility of large stepped failures and rupture of rock bridges. The maximum calculated persistence for the two existing sliding planes in the study site is, respectively, 28.0 m and 48.5 m. The maximum calculated sliding plane surfaces are, accordingly, 282.5 m2 and 289.3 m2. These results are compared against the observed scar dimensions at the study site, which have been retrieved alternatively, by processing a LiDAR point cloud. The results from the two alternative sources are similar, indicating that the methodology can be efficiently used to assess the sliding plane persistence and the expected maximum rockfall magnitude at the study site.Peer ReviewedPostprint (published version

Comparing rockfall scar volumes and kinematically detachable rock masses

Scenario-based risk assessment for rockfalls, requires assumptions for different scenarios of magnitude (volume). The magnitude of such instabilities is related to the properties of the jointed rock mass, with the characteristics of the existing unfavourably dipping joint sets playing a major role. The critical factors for the determination of the maximum credible rockfall volume in a study site, the Forat Negre in Andorra, are investigated. The results from two previous analyses for the rockfall size distribution at this site are discussed. The first analysis provides the observed size distribution of the rockfall scars, and it is an empirical evidence of past rockfalls. The second one, calculates the kinematically detachable rock masses, indicating hypothetical rockfalls that might occur in the future. The later gives a maximum rockfall volume, which is one order of magnitude higher, because the persistence of the basal planes is overestimated. The tension cracks and lateral planes interrupt systematically the basal planes, exerting a control over their persistence, and restricting the formation of extensive planes and large rockfall failures. Nonetheless, the formation of basal planes across more than one spacings of tension cracks is possible and small step-path failures have been observed too. Concluding, the key factor for the determination of the maximum credible volume at the study-site is the maximum realistic length of the basal planes, penetrating into the rock mass, their spacing, and, if applied, the contribution of the rock bridges to the overall rock mass resistance.Peer ReviewedPostprint (author's final draft

A fractal fragmentation model for rockfalls

The final publication is available at Springer via http://dx.doi.org/10.1007/s10346-016-0773-8The impact-induced rock mass fragmentation in a rockfall is analyzed by comparing the in situ block size distribution (IBSD) of the rock mass detached from the cliff face and the resultant rockfall block size distribution (RBSD) of the rockfall fragments on the slope. The analysis of several inventoried rockfall events suggests that the volumes of the rockfall fragments can be characterized by a power law distribution. We propose the application of a three-parameter rockfall fractal fragmentation model (RFFM) for the transformation of the IBSD into the RBSD. A discrete fracture network model is used to simulate the discontinuity pattern of the detached rock mass and to generate the IBSD. Each block of the IBSD of the detached rock mass is an initiator. A survival rate is included to express the proportion of the unbroken blocks after the impact on the ground surface. The model was calibrated using the volume distribution of a rockfall event in Vilanova de Banat in the Cadí Sierra, Eastern Pyrenees, Spain. The RBSD was obtained directly in the field, by measuring the rock block fragments deposited on the slope. The IBSD and the RBSD were fitted by exponential and power law functions, respectively. The results show that the proposed fractal model can successfully generate the RBSD from the IBSD and indicate the model parameter values for the case study.Peer ReviewedPostprint (author's final draft

Magnitude and frequency relations: are there geological constraints to the rockfall size?

The final publication is available at Springer via http://dx.doi.org/10.1007%2Fs10346-017-0910-zThere exists a transition between rockfalls, large rock mass failures, and rock avalanches. The magnitude and frequency relations (M/F) of the slope failure are increasingly used to assess the hazard level. The management of the rockfall risk requires the knowledge of the frequency of the events but also defining the worst case scenario, which is the one associated to the maximum expected (credible) rockfall event. The analysis of the volume distribution of the historical rockfall events in the slopes of the Solà d’Andorra during the last 50 years shows that they can be fitted to a power law. We argue that the extrapolation of the F-M relations far beyond the historical data is not appropriate in this case. Neither geomorphological evidences of past events nor the size of the potentially unstable rock masses identified in the slope support the occurrence of the large rockfall/rock avalanche volumes predicted by the power law. We have observed that the stability of the slope at the Solà is controlled by the presence of two sets of unfavorably dipping joints (F3, F5) that act as basal sliding planes of the detachable rock masses. The area of the basal sliding planes outcropping at the rockfall scars was measured with a terrestrial laser scanner. The distribution of the areas of the basal planes may be also fitted to a power law that shows a truncation for values bigger than 50 m2 and a maximum exposed surface of 200 m2. The analysis of the geological structure of the rock mass at the Solà d’Andorra makes us conclude that the size of the failures is controlled by the fracture pattern and that the maximum size of the failure is constrained. Two sets of steeply dipping faults (F1 and F7) interrupt the other joint sets and prevent the formation of continuous failure surfaces (F3 and F5). We conclude that due to the structural control, large slope failures in Andorra are not randomly distributed thus confirming the findings in other mountain ranges.Peer ReviewedPostprint (author's final draft

Rockfall occurrence and fragmentation

Rockfalls are very rapid and damaging slope instability processes that affect mountainous regions, coastal cliffs and slope cuts. This contribution focuses on fragmental rockfalls in which the moving particles, particularly the largest ones, propagate following independent paths with little interaction among them. The prediction of the occurrence and frequency of the rockfalls has benefited by the rapid development of the techniques for the detection and the remote acquisition of the rock mass surface features such as the 3D laser scanner and the digital photogrammetry. These techniques are also used to monitor the deformation experienced by the rock mass before failure. The quantitative analysis of the fragmental rockfalls is a useful approach to assess risk and for the design of both stabilization and protection measures. The analysis of rockfalls must consider not only the frequency and magnitude of the potential events but also the fragmentation of the detached rock mass. The latter is a crucial issue as it affects the number, size and the velocity of the individual rock blocks. Several case studies of the application of the remote acquisition techniques for determining the size and frequency of rockfall events and their fragmentation are presented. The extrapolation of the magnitude-frequency relationships is discussed as well as the role of the geological factors for constraining the size of the largest detachable mass from a cliff. Finally, the performance of a fractal fragmentation model for rockfalls is also discussed.Peer ReviewedPostprint (published version

Metodología integrada para la evaluación de riesgos en la Red de Carreteras de Gipuzkoa

Se presenta una metodología de carácter integradorpara evaluar el nivel de riesgo de fallo o rotura en estructuras y taludes anclados,terraplenes, desprendimientos en laderas y desmontes, corrientes de derrubios y deslizamientos que ha sido aplicada a la Red de Carreteras de Gipuzkoa. En la presente comunicación se proporcionan los detalles del cálculo del riesgo para desprendimientos rocosos y estructuras de contención ancladas.Postprint (published version

Methodology to evaluate rock slope stability under seismic conditions at Solà de Santa Coloma, Andorra

An analytical methodology is presented to evaluate rock slope stability under seismic conditions by considering the geomechanical and topographic properties of a slope. The objective is to locate potential rockfall source areas and evaluate their susceptibility in terms of probability of failure. For this purpose, the slope face of a study area is discretized into cells having homogenous aspect, slope angle, rock properties and joint set orientations. A pseudostatic limit equilibrium analysis is performed for each cell, whereby the destabilizing effect of an earthquake is represented by a horizontal force. The value of this force is calculated by linear interpolation between the peak horizontal ground acceleration PGA at the base and the top of the slope. The ground acceleration at the top of the slope is increased by 50% to account for topographic amplification. The uncertainty associated with the joint dip is taken into account using the Monte Carlo method. The proposed methodology was applied to a study site with moderate seismicity in Sol`a de Santa Coloma, located in the Principality of Andorra. The results of the analysis are consistent with the spatial distribution of historical rockfalls that have occurred since 1997. Moreover, the results indicate that for the studied area, 1) the most important factor controlling the rockfall susceptibility of the slope is water pressure in joints and 2) earthquake shaking with PGA of 0.16 g will cause a significant increase in rockfall activity only if water levels in joints are greater than 50% of the joint height.Postprint (published version

Calculation of the rockwall recession rate of a limestone cliff, affected by rockfalls, using cosmogenic chlorine-36 : case study of the Montsec Range (Eastern Pyrenees, Spain)

Cliff erosion may be a major problem in settled areas affecting populations and producing economic and ecological losses. In this paper we present a procedure to calculate the long-term retreat rate of a cliff affected by rockfalls in the Montsec Range, Eastern Pyrenees (Spain). It is composed of low, densely fractured limestones; and the rockwall is affected by rockfalls of different sizes. The rockfall scars are clearly distinguishable by their regular boundaries and by their orange colour, which contrast with the greyish old reference surface (S0) of the cliff face. We have dated different stepped surfaces of the rockwall, including S0, using cosmogenic 36Cl. The total amount of material released by rockfall activity was calculated using a high definition point cloud of the slope face obtained with a terrestrial laser scanner (TLS). The present rockwall surface has been subtracted from the reconstructed old cliff surface. This has allowed the calculation of the total volume released by rockfalls and of the retreat rate. The latter ranges from 0.31 to 0.37 mm·a−1 . This value is of the same order of magnitude as that obtained by other researchers in neighbouring regions in Spain, having similar geology and affected by rockfalls

Size distribution for potentially unstable rock masses and in situ rock blocks using LIDAR-generated digital elevation models

In this paper, two analytical procedures which are independent from the existence of empirical data are presented for the calculation of (1) the size distribution of potentially unstable rock masses that expresses the potential rockfall size distribution, including big volumes corresponding to potential rare events with low susceptibility of failure and (2) the in situ block distribution on the slope face. Two approaches are, respectively, used. The first one involves the detection of kinematically unstable surfaces on a digital elevation model (DEM) and on orthophotos and the calculation of the volumes resting on them. For the second one the in situ block volumes formed by the intersection of the existing discontinuity sets are calculated using a high-resolution DEM. The procedures are presented through an application example at the country of Andorra and in particular at the chute of Forat Negre. The results from the first procedure indicate that it is kinematically possible to have mobilized volumes of some thousands of cubic meters; however, these are considered rare events with low susceptibility of failure. The size distribution of potentially unstable rock masses for big volume events was well fitted by a power law with an exponent of -0.5. The in situ block distribution on the slope face from the second procedure, assuming three types of intersection between the joints of the existing discontinuity sets and two extreme cases of discontinuity persistence, was also found to follow a power law, but with an exponent of -1.3. The comparison with the observed in the field block volume distribution on the slope face indicates that in reality discontinuities have a very high persistence and that considering only their visible trace length overestimates volumes, which is conservative.Peer ReviewedPostprint (author’s final draft

A methodology to obtain the block size distribution of fragmental rockfall deposits

Peer ReviewedPostprint (author's final draft