Rockfalls into forests: Analysis and simulation of rockfall trajectories — considerations with respect to mountainous forests in Switzerland

Rockfalls are a major threat to settlements and transportation routes in many places. Although the general protective effect of forests against rockfalls is currently not questioned, little is known about the ideal properties of a forest stand that provides good protection. Therefore, in this study the question was assessed of how mountainous forests may influence rockfalls of single boulders. An actual rockfall trajectory was measured, recorded, analysed and simulated with a rockfall model. Rockfalls into different forest scenarios were also modelled for the site. Results showed that the actual rockfall event can be well simulated. Furthermore, a completely forested slope reduces velocity and energy of the falling blocks much better than a sparsely forested slope. For the profile discussed in this paper, the largest effect upon falling 3 m³ blocks was obtained with a high forest containing 350 trees per ha. The results confirmed common assumptions on ideal properties of a protective forest stand against rockfalls.     

Rockfall risk management using a pre-failure deformation database

Terrestrial laser scanning (TLS) monitoring has been used to estimate the location, volume, and kinematics of a variety of small magnitude rockfalls before failure (1–1000 m³ range), and in some cases, potential failure time has been assessed through the application of inverse velocity methods. However, our current understanding of rock slope pre-failure behavior for this magnitude range and prediction ability is based on observations of a small number of failure case histories. In this study, a pre-failure deformation database was constructed for rockfall volumes exceeding 0.1 m³, observed over a 1252-day study interval at the Goldpan rock slope, British Columbia, Canada, in order to better understand the pre-failure behavior of rock slopes and provide an empirical means of estimating temporal failure ranges. Repeated TLS datasets were acquired at an average scanning interval of 2–3 months. A total of 90 rockfall events were recorded at this site, during this time period, of which 64 (71%) exhibited measurable deformation prior to failure. Classification of rockfalls by volume suggests that a scale dependency may exist, as deformation was detected for a greater proportion of rockfalls >?5 m³ (92%) than for smaller rockfalls in the range of 0.1–0.5 m³ (61%). A lower rate of pre-failure deformation detection was also reported for planar sliding failures as compared with wedge or toppling failures, suggesting that deformation was less easily detected for these failure types. This study proposes and implements a framework for rockfall assessment and forecasting that does not require continuous monitoring of deformation.     

Spatio-temporal analysis of rockfall pre-failure deformation using Terrestrial LiDAR

We present a long-term spatio-temporal analysis of rock slope evolution using a Terrestrial LiDAR aiming to improve our understanding of the link between pre-failure deformation and the spatial prediction of rockfalls. We monitored the pilot study area located at the Puigcercós cliff (Catalonia, Spain) over a period of 1,705 days and detected the deformation of nine different cliff regions together with a high rockfall activity. An exact match was observed between the progressively deformed areas and the regions recently affected by three of the highest magnitude rockfall events, demonstrating a causal relationship between pre-failure deformation and rockfall occurrence. These findings allowed us to make a forward spatial prediction of future failures, hypothesizing a high probability of failure in the six remaining regions. We observed an exponential acceleration of the deformation close to failure, in accordance with tertiary creep theory. However, the temporal analysis of the deformed areas showed a complex and variable behavior, so no exact prediction of the date of failure can yet be made. Our findings have broadened our understanding of the pre-failure behavior of rockfalls and have clear implications for the future implementation of early warning systems.     

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

There 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 m² and a maximum exposed surface of 200 m². 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.     

Rockfall Hazard Analysis for Hong Kong Based on Rockfall Inventory

Summary ¶This paper compiles and analyzes the rockfall data in Hong Kong in the last fifty years. A simple rockfall hazard analysis is presented based on this rockfall inventory. A frequency-magnitude relation, which is analogous to the Gutenberg-Richter relation for earthquake occurrence, is proposed for rockfall, and direct correlation between rockfall frequency and the daily rainfall is observed. Data analysis shows that a threshold daily rainfall of about 150–200mm is expected in order to trigger rockfall events in Hong Kong. Among the 368 rockfall events in the 13 year period from 1984 to 1996 in Hong Kong, 35% of the incidents lead to blockage of or damage to roads, 22% lead to damages or evacuation of squatter huts, 21% lead to blockage of pedestrian pavement and footpath, and 15% affect buildings, such as housing apartments and schools. Only 15% of these rockfalls fell onto open space and caused negligible effects on human activities. Most of these rockfall events occurred during heavy rain and when landslide warning should has been issued by the Hong Kong Observatory, thus only 6% of these events led to injury or casualty, car damages, and damages to public utilities. Rockfall hazard zonation maps in terms of the spatial distribution of previous rockfalls are proposed for both Kowloon Peninsula and Hong Kong Island.     

Calibration and validation of rockfall modelling at regional scale: application along a roadway in Mallorca (Spain) and organization of its management

The Tramuntana range, in the northwestern sector of the island of Mallorca (Balearic Islands, Spain), is frequently affected by rockfalls which have caused significant damage, mainly along the road network. In this work, we present the procedure we have applied to calibrate and validate rockfall modelling in this region, using 103 cases of the available detailed rockfall inventory (630 rockfall events collected since the eighteenth century). We have exploited STONE (Guzzetti et al. 2002), a GIS-based rockfall simulation software which computes 2D and 3D rockfall trajectories starting from a DTM and maps of the dynamic rolling friction coefficient and of the normal and tangential energy restitution coefficients. The appropriate identification of these parameters determines the accuracy of the simulation. To calibrate them, we have selected 40 rockfalls along the range which include a wide variety of outcropping lithologies. Coefficients values have been changed in numerous attempts in order to select those where the extent and shape of the simulation matched the field mapping. Best results were summarized with the average statistical values for each parameter and for each geotechnical unit, determining that mode values represent more precisely the data. Initially, for the validation stage, 10 well-known rockfalls exploited in the calibration phase have been selected. Confidence tests have been applied to their modelling results taking into account not only the success but also the mistakes. The best accuracy is obtained when the rockfall has a preferential trajectory and worse results when the rockfall follows two or more trajectories. Additionally, the greater the rockfall runout length, the less precise the simulation is. We have further validated the calibrated parameters along the Ma-road (111 km), the main transportation corridor in the range, using 63 rockfall events that occurred during the past 18 years along the road. Of the rockfalls where source areas were properly identified, 81.5 % are well represented by STONE modelling, as the travel paths and the depositional areas are successfully ascertained. Results of the analysis have been used by the Road Maintenance Service of Mallorca to assess hazard and risk posed by rockfall at regional scale to design the road management plan.     

融合区域潜在落石源区识别的峡谷区落石危险性评价

山西壶关太行山大峡谷景区为中国最美十大峡谷之一,但景区落石灾害频发,严重威胁景区安全运营。本文基于高精度地形信息与岩土体强度特性,采用坡度角分布方法开展区域尺度潜在落石源区识别,并引入岩体破坏敏感性指标定量描述潜在落石源区失稳概率。然后,利用经验模型Flow-R模拟落石运动扩散过程,获取落石的传播概率与能量分布情况。最后,提出落石危险性双因子评价模型实现落石危险性定量评估。获得主要结论如下:(1)研究区内潜在落石源区面积为25.7 km2(35.7%),主要以条带状分布于峡谷两侧陡壁。其中岩体破坏高敏感性区为3.3 km2。(2)研究区落石高危险区面积达3.22 km2,主要威胁景区内游客集散地与交通线路,尤其在S327荫林线红豆峡入口处落石危险性最高。(3)野外调查验证结果表明了应用坡度角分布方法识别潜在落石源区的高效性与准确性,提出的双因子评价模型可为峡谷区落石危险性评估提供快速解决方案。本文提出的“区域落石源区识别-源区失稳概率分析-落石危险性评估”的一整套技术方案能够为类似的高山峡谷区落石灾害早期识别及风险防控提供技术参考。     

郑万高铁宜万段边坡危岩崩落破坏特征

危岩是山区常见的地质灾害之一。以往研究缺少对危岩整体破坏导致危岩解体方面的关注,而危岩在失稳崩落过程中的解体行为却是预测危岩影响范围和防治成效的关键所在。为此,文章以郑万（郑州—万州）高铁宜万段沿线隧道洞口边坡危岩为研究对象,从结构面角度出发,对危岩崩落破坏特征进行研究。通过对15个隧道洞口边坡的调查,首先从边坡坡度、岩性组合、相对高差三个方面总结了研究区危岩发育分布规律;然后根据边坡岩体结构特征,分析了危岩失稳模式,并基于边坡上部危岩和下部落石的体积和形状对应关系,进一步探讨了边坡危岩崩落破坏演化过程;在此基础上,利用Rockfall模拟软件对落石运动特征进行预测分析。结果表明:（1）研究区边坡呈上陡下缓地形,上部基岩裸露,坡度基本上≥70°;危岩主要发育于弱风化的灰岩和白云岩中;边坡高差在150~300 m之间。（2）边坡上部危岩将呈阶梯状方式逐渐沿基底结构面滑移或者沿后缘结构面拉裂坠落。（3）研究区危岩崩落破坏模式主要为边坡上部岩体沿结构面解体破坏。（4）大部分隧道洞口边坡落石危险性较大,严重威胁隧道洞口的安全,需要采取相应的防治措施。研究成果可为在建的郑万高铁宜万段隧道边坡危岩的有效防治提供参考。     

Rockfall susceptibility zoning at a large scale: From geomorphological inventory to preliminary land use planning

Rockfall, up to several hundreds of cubic meters, is a frequent and rapid landslide which menaces extensive areas in mountainous territories. Rockfall susceptibility zoning map at a large scale (1:5000–1:25 000) can be the first tool for land use planning in order to manage rockfall risk. A methodology allowing to analyze susceptibility in extensive areas with optimum cost/benefit relationship is needed. This work analyzes rockfall susceptibility in an extensive rocky mountain of the Principality of Andorra (Pyrenees Mountains), first on the rock slope and then on the exposed area located below. The rockfall record, obtained by means of geomorphological analysis, supplies the main data to analyze the susceptibility on the rock slope. An additional historical inventory verifies the accuracy of rockfall sizes recorded by means of the geomorphological analysis. According to the classification recommended by the Guidelines of Joint Technical Committee, the density of rockfall features on the rock slope assesses susceptibility in four levels. Subsequently, susceptibility on exposed areas has been analyzed by means of reach probability of rock blocks analysis using empirical models. Data acquired from thirteen recent events, from 1999 to 2004, have been used to verify the accuracy of the two empirical models mainly used (reach angle and shadow angle). Five reach probability limits (1, 0.5, 0.25, 0.01, and 0) establish boundaries between susceptibility levels. The resulting rockfall susceptibility zoning map allows: (a) to identify land areas and human elements exposed to rockfalls and, (b) to establish several exposition levels. This map can be a useful and cost-effective tool for administrations responsible to manage natural risk in order to guide urban grow in extensive areas or decide upon work programs based on in-depth analysis (hazard and risk).     

青海玉树断裂带地震落石的地震地质意义

2010年青海玉树Ms 7.1级大地震引发了一系列次生地质灾害,其中地震落石是除地震滑坡外沿断裂带及其邻侧最常见的现象。对玉树震区落石的调查发现,该区多处存在非常典型的多期地震落石分布现象,指示该区地震落石的发育与其他古地震现象类似,具有多期性和一定的原地复发性。实地调查表明,该区地震落石分布的主要特征为:多集中发育在活动断裂带附近的陡峭基岩斜坡下方,分布零散,且滚动较远,并常与古地震滑坡相伴生。初步获得的8个地震落石钙膜U系测年结果分布在距今6030±300a BP、4720±210a BP、3530±490~3560±280a BP、2010±160a BP、1090±70a BP、760±20a BP和230±20a BP年龄段,与该区古地震探槽和滑坡反映的地震事件比较吻合,进一步揭示玉树断裂带附近在全新世中晚期发生过多期可导致地表产生地震落石的事件。同时也说明,地震落石及其钙膜测年是特别值得进一步探索的潜在古地震研究方法或途径。     

Application of a long-range Terrestrial Laser Scanner to a detailed rockfall study at Vall de Núria (Eastern Pyrenees, Spain)

In this study we show the application of a long-range Terrestrial Laser Scanner (TLS) to a detailed rockfall study in a test zone at Vall de Núria, located in the Eastern Pyrenees. Data acquisition was carried out using TLS-Ilris3D, the new generation of reflector-less laser scanners with a high range, accuracy and velocity of measurements. Eight scans were performed at 3 stations to acquire coordinates of almost 4 million points. The results from the acquired data are a high accuracy Digital Elevation Model (DEM) and the reconstruction of the joint geometry. The former is used for inventory of rockfalls and for more accurate rockfall simulation (trajectories and velocities). The latter allows us to model the geometry and volume of the source area in recent rockfalls. Our findings suggest that TLS technology could be a tool of reference in rockfall studies in the near future.     

An analysis of failure mechanism constraints on pre-failure rock block deformation using TLS and roto-translation methods

Rockfall hazards increase the risk of train derailment along railway corridors in western Canada. In this study, repeated terrestrial laser scanning (TLS) datasets were collected every 2–3 months at three different sites along the Thompson and Fraser River corridors in British Columbia, referred to as the Goldpan, White Canyon, and Mile 109 sites. A total of 207 rockfall events occurring across all three sites between November 11, 2014 and October 18, 2016 were recorded in a database. For each of these rockfalls, pre-failure deformation was measured using a method of three-dimensional roto-translation block tracking. Each rockfall was classified by its deformation behaviour and further categorised based on failure mechanism, volume, lithology, and the roughness condition of the failure plane. Results reveal that detectable levels of deformation were measured in 33% of the total number of rockfall events using the present methods. Rotation deformation was most commonly observed in toppling failures with relatively steep joint orientations. Conversely, planar sliding blocks generally exhibited the least measurable deformation, with the majority not showing any precursory translation or rotation. It is postulated that overhanging rockfall configurations may suppress the expression of deformation in rockfall source blocks, though additional research is required to confirm this.     

Seismic-induced rockfalls and landslide dam following the October 30, 2016 earthquake in Central Italy

On October 30, 2016, a seismic event and its aftershocks produced diffuse landslides along the SP 209 road in the Nera River Gorge (Central Italy). Due to the steep slopes and the outcropping of highly fractured and bedded limestone, several rockfalls were triggered, of which the main event occurred on the slope of Mount Sasso Pizzuto. The seismic shock acted on a rock wedge that, after an initial slide, developed into a rockfall. The debris accumulation blocked the SP 209 road and dammed the Nera River, forming a small lake. The river discharge was around 3.6 m³/s; the water overtopped the dam and flooded the road. By a preliminary topographic survey, we estimated that the debris accumulation covers an area of about 16,500 m², while the volume is around 70,000 m³. The maximum volume occupied by the pre-existing talus mobilized by the rockfall is about 20% of the total volume. Besides blocking the road, the rockfall damaged a bridge severely, while, downstream of the dam, the water flow caused erosion of a road embankment. A rockfall protection gallery, a few hundred meters downstream of the dam, was damaged during the event. Other elastic nets and rigid barriers were not sufficient to protect the road from single-block rockfalls, with volumes around 1–2 m³. Considering the geological and geomorphological conditions, as well as the high seismicity and the socioeconomic importance of the area, a review of the entire rockfall protection systems is required to ensure protection of critical infrastructure and local communities.     

Progressive failure leading to the 3 December 2013 rockfall at Puigcercós scarp (Catalonia,Spain)

A detailed analysis of the pre-failure behavior of the 3 December 2013 rockfall (1,012 m³) occurred on Puigcercós pilot study area (Catalonia, Spain) is presented. The exact date of failure was obtained based on a photographic monitoring performed every 4 h. The long-term monitoring (2,217 days) of the rock slope carried out by a Terrestrial LiDAR allowed the early detection of both pre-failure deformation and precursory rockfalls preceding the final failure. By means of the analysis of the pre-failure deformation, four different deformed areas were detected and the tertiary creep phase was observed in three of them. An attempt to predict the time to failure was performed using the Fukuzono’s (1985) method. Furthermore, the temporal evolution of the precursory rockfalls occurred in those four areas during the progressive failure showed a close resemblance with the exponential pattern of the cumulated displacements at tertiary creep stage. Finally, the study of the meteorological conditions did not show any single triggering factor associated with the final failure. Reversely, the increase in the occurrence of precursory rockfalls on several areas of the slope together with the observed acceleration on the deformation pattern reinforce the role of a progressive degradation of the stability conditions, which ultimately leaded to the 3 December rockfall event.     

Evaluating seismically induced mass movement hazard in Caramanico Terme (Italy)

Mass movements and earthquakes represent two major geological hazards in the municipal territory of Caramanico Terme (south-central Apennines). Available records revealed the contemporaneous occurrence of earthquakes and slope failures on four occasions in the last four centuries (1627, 1706, 1933, and 1984). These events, with local intensities ranging from VI to IX, generated mass movements varying from a rotational slope failure to rock/block falls. All occurred in the southern periphery of the town and involved a thick carbonate megabreccia caprock and coarse colluvia which overlie a clayey substratum. Field investigation and review of historical records helped to delimit the areas susceptible to seismically triggered rockfalls. The mapping of historic and pre-historic rockfall deposits revealed their dispersal patterns and provided the basis for a determination of potential hazard zones. We approximate the temporal hazard assessment by relating the rockfall occurrence to the probability of earthquake triggering. Considering the VI degree triggering threshold indicated by local historical data, the statistical analysis of the regional seismic activity shows that events capable of inducing rockfalls have an approximately decennial recurrence in Caramanico. The approach presented could be readily applied to other potential risk areas of Italy by exploiting the rich long-term record of historical seismicity. In general, temporal hazard estimates at relatively low intensity levels will be possible even where the seismic history of the site is only well documented for a relatively limited time interval, provided that this interval is much longer than the recurrence time of the events exceeding the threshold considered.     

贺兰口岩画载体落石危害及治理研究

因暴雨、地震和岩羊踩踏等不确定因素,贺兰口沟口至“水关”处发育的34处孤石极易失稳滚落形成落石对岩画和游人安全构成严重威胁。通过对34处孤石数量、成因、分布、岩性和几何特征全面调查与分析的基础上,运用Rocfall软件分析其失稳后的运动过程和能量变化进行分析,得出如下结论:(1)薄层状变质砂岩风化内凹的临空面、岩层层面与两组节理组合切割巨厚层砂岩失稳后堆积于坡面;(2)将落石划分为基岩裸露型、孤石带型和坡面泥石流型3种不同类型;(3)数值模拟结果显示,失稳后弹跳落点有52%~95%的概率在岩画赋存区,所具有能量为301~8412 kJ;最终落点有51%~95%的概率在游客栈道,所具有能量为299~7698 kJ;(4)监测结果表明,水剂比为0.32、平行布孔、间排距20 cm和抵抗线20 cm布置方案为最优,静态破碎结合人工搬运为岩画保护区落石源头治理的最佳方案。为贺兰口岩画保护区落石灾害的源头治理提供了新的途径与参考。     

无人机影像在高陡边坡危岩体调查中的应用

在高陡边坡危岩体的调查中,复杂的地形条件经常限制工作的正常开展,如何快速准确地获取地质灾害信息一直是地质灾害调查研究中的难点之一。以往的研究中对无人机遥感技术在黄土、高原等地区应用有所报道,但对西南地区高陡边坡危岩体灾害调查的研究尚无报道。文章以锦屏二级水电站出线场边坡落石灾害所在区域为例,将无人机摄影测量技术应用于高陡边坡危岩体调查中,通过无人机倾斜摄影获取高分辨率遥感影像,开展遥感影像三维建模,进行地质灾害遥感解译,总结了无人机遥感系统在高陡边坡危岩体调查的技术流程。通过三维实景模型,精确地分析了落石灾害的空间分布、失稳模式及演化过程,查明了区域内危岩隐患点的分布特征;基于三维点云模型,提取出地质灾害体的属性信息,测得落石方量为11.7 m3,采用最小二乘法进行平面拟合,得到落石两组主控结构面产状为275.4°∠31.2°、103.5°∠63.3°。实践表明,无人机遥感技术在高陡边坡地区落石灾害调查中具有明显的可行性和优越性,可以较好地应用于高陡边坡危岩体调查中。     

Evaluation of rock falls in an urban area: the case of Boğaziçi (Erzincan/Turkey)

Disasters caused by events such as earthquake, flooding, rock falls, landslides are often encountered. However, generally, the reasons for the destructive and devastating effects of these nature events are that settlement locations were chosen without site investigation studies, or that available studies were inadequate. Such inadequacies in the field are related to inappropriate settlement location and the resulting damage caused by rock falls. This study evaluated rockfall risk in a settlement that developed in a similar manner. The study was carried out in Bo?aziçi village of Kemah (Erzincan/Turkey), which is located in a very important tectonic zone. The study site is located on the lower sections of an area with very steep cliffs and 50–75° slopes. This cliff, which is the source of rockfalls, has a slope dip of approximately 90°. The cliff comprises 25–30 m high, fractured and cracked basaltic volcanic mass. To determine block size in the study area, scanline survey measurements and block size measurements were performed on blocks that loosened and fell from the cliff face. It was found that block sizes reached 6 m³. Rockfall analyses were performed along the selected profiles using the Rockfall V.4.0 software. Kinetic energy, bounce height, horizontal location of rock end-points, and velocity of the rocks along each section were evaluated separately for each profile. This data were used to produce distribution maps for each profile and the settlement was evaluated in terms of rockfall risk. The results indicate that the study area was at risk of future rockfalls and that it would be appropriate to relocate one part of the settlement.     

Rockfall detachment susceptibility map in El Hierro Island, Canary Islands, Spain

A semi-quantitative heuristic methodology is developed to map a rockfall detachment susceptibility zonation of El Hierro Island (Canary Archipelago). The rationalized procedure, which we called non-weighted bounded indicators, is based on overlapping thematic maps of conditioning factors to mass movement, which are appropriately and individually rescaled and then composed by addition to obtain a susceptibility numerical index through a GIS. As the consistency of the geomorphological analysis depends on the expert subjective criteria and the appropriate interpretation of the landscape, the use of this methodology reduces subjectivity and quantifies the degree of susceptibility. The main factors affecting the mass movement phenomena (rockfalls events), also recognized in the field and, therefore, considered in the presented GIS arrangement, are slope, profile curvature, lithology, vegetation cover and dykes density. To calculate the slope threshold or minimum angle characteristic of rockfall source areas, mixed Gaussian slope frequency decomposition is used. The curvature index reveals stepwise areas. Qualitative geomechanical characteristics are linked to a quantitative index according to a volcanic lithological-complexes classification. Both destabilization (root-wedging) and stabilization effects are considered into the vegetation cover index. The dyke density index incorporates the bearing rock capacity decrease produced in the halo around a dyke network intrusion. Slope, curvature and vegetation indexes thresholds have been fitted following field observations. A rockfall detachment susceptibility map is obtained and classified based on the histogram maxima. The rockfall inventory, based on rockfall events reported within the island, was used for the model validation. A 12?% of the whole island shows medium to very high susceptibility.