Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (104)

Search Parameters:
Keywords = rockfall hazard

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
15 pages, 3509 KiB  
Article
Damage Characteristics of Blasting Surrounding Rock in Mountain Tunnel in Fault Fracture Zones Based on the Johnson–Holmquist-2 Model
by Lizhi Cheng, Zhiquan Yang, Ping Zhao and Fengting Li
Buildings 2024, 14(11), 3682; https://doi.org/10.3390/buildings14113682 - 19 Nov 2024
Viewed by 372
Abstract
Blasting is a widely employed technique for tunnel construction in mountainous regions; however, it often causes damage to the surrounding rock mass, particularly in fault fracture zones, which can lead to hazards such as rockfalls and collapses. This study examines the characteristics of [...] Read more.
Blasting is a widely employed technique for tunnel construction in mountainous regions; however, it often causes damage to the surrounding rock mass, particularly in fault fracture zones, which can lead to hazards such as rockfalls and collapses. This study examines the characteristics of damage to surrounding rock due to tunnel blasting through fault fracture zones. Based on an actual tunnel blasting construction project, we conducted a finite element analysis using the JH-2 material model, taking into account the width of the fault fracture zone. Results indicate that as the width of the fault fracture zone increases, the disturbance effect of tunnel blasting on the surrounding rock becomes more pronounced. Compared to the arch bottom and arch waist of the tunnel, the tunnel vault primarily absorbs the slip deformation and compressive forces resulting from blasting disturbances in the fault fracture zone. The findings of this paper contribute a valuable methodology for analyzing the mechanical mechanisms in mountain tunnel blasting and provide essential theoretical parameters to inform the design and construction of tunnel blasting projects. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the distribution of fault fracture zones of Lianfeng Tunnel.</p>
Full article ">Figure 2
<p>JH-2 material intrinsic model’s strength model.</p>
Full article ">Figure 3
<p>Numerical calculation model.</p>
Full article ">Figure 4
<p>Stress distribution of surrounding rock induced by tunnel blasting in fracture area.</p>
Full article ">Figure 4 Cont.
<p>Stress distribution of surrounding rock induced by tunnel blasting in fracture area.</p>
Full article ">Figure 5
<p>Strain distribution of surrounding rock induced by tunnel blasting in fracture area.</p>
Full article ">Figure 5 Cont.
<p>Strain distribution of surrounding rock induced by tunnel blasting in fracture area.</p>
Full article ">Figure 6
<p>Shear stress cloud of X-Y section.</p>
Full article ">Figure 7
<p>Stress diagram of lining structure.</p>
Full article ">Figure 8
<p>Circumferential strain curves of the tunnel blast lining structure with the distance of 1 m, 10 m, and 20 m from the palm face under the condition of 20 m fault width.</p>
Full article ">Figure 9
<p>Characteristics of damage and destruction of the vault caused by tunnel blasting under the condition of fault width of 20 m.</p>
Full article ">
16 pages, 32403 KiB  
Article
Integrated Analysis of Rockfalls and Floods in the Jiului Gorge, Romania: Impacts on Road and Rail Traffic
by Marian Puie and Bogdan-Andrei Mihai
Appl. Sci. 2024, 14(22), 10270; https://doi.org/10.3390/app142210270 - 8 Nov 2024
Viewed by 675
Abstract
This study examines the impact of rockfalls and floods on road and rail traffic in the Jiului Gorge, Romania, a critical transportation corridor. Using Sentinel-1 radar imagery processed through ESA SNAP and ArcGIS Pro, alongside traffic detection facilitated by YOLO models, we assessed [...] Read more.
This study examines the impact of rockfalls and floods on road and rail traffic in the Jiului Gorge, Romania, a critical transportation corridor. Using Sentinel-1 radar imagery processed through ESA SNAP and ArcGIS Pro, alongside traffic detection facilitated by YOLO models, we assessed susceptibility to both rockfalls and floods. The primary aim was to enhance public safety for traffic participants by providing accurate hazard mapping. Our study focuses on the area from Bumbești-Jiu to Petroșani, traversing the Southern Carpathians. The results demonstrate the utility of integrating remote sensing with machine learning to improve hazard management and inform more effective traffic planning. These findings contribute to safer, more resilient infrastructure in areas vulnerable to natural hazards. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area.</p>
Full article ">Figure 2
<p>Sentinel 1 GRD images of study area from descending orbit (left, 20 January 2023), from ascending orbit (middle, 20 January 2023), RGB interferogram and processing software workflow.</p>
Full article ">Figure 3
<p>ESA SNAP software workflow image samples for rockfall detection from Sentinel-1 SLC product.</p>
Full article ">Figure 4
<p>Flood map of the Jiului Gorge region, illustrating the extent and severity of flood events based on Sentinel-1 GRD images.</p>
Full article ">Figure 5
<p>Rockfall map displaying incidents along National Road 66 and surrounding slopes for specified dates.</p>
Full article ">Figure 6
<p>Rockfall susceptibility map showing areas highly susceptible to rockfall, with a notable prevalence in the upper section of the gorge.</p>
Full article ">Figure 7
<p>Rockfall susceptibility map combined with affected areas from radar images, highlighting the upper part of the gorge with high susceptibility.</p>
Full article ">Figure 8
<p>Floods susceptibility map combining DEM-derived slope classifications, land cover types, rainfall, and proximity to water bodies, showing higher susceptibility in wider parts of the gorge.</p>
Full article ">Figure 9
<p>Floods susceptibility map combined with radar-detected flood areas, illustrating increased susceptibility in the central and southern parts of the gorge.</p>
Full article ">Figure 10
<p>Train detection and recognition using YOLO models, illustrating detection from a significant distance with reduced visibility.</p>
Full article ">Figure 11
<p>Road traffic element detection with greater precision due to closer camera proximity.</p>
Full article ">Figure 12
<p>Training results from YOLOv9 model, showcasing classes obtained after training.</p>
Full article ">Figure 13
<p>Detection results including several classes, highlighting various rockfall types.</p>
Full article ">Figure 14
<p>Detection results focusing on a single class, illustrating detailed rockfall identification.</p>
Full article ">
31 pages, 114861 KiB  
Article
Multitemporal Monitoring of Rocky Walls Using Robotic Total Station Surveying and Persistent Scatterer Interferometry
by Luisa Beltramone, Andrea Rindinella, Claudio Vanneschi and Riccardo Salvini
Remote Sens. 2024, 16(20), 3848; https://doi.org/10.3390/rs16203848 - 16 Oct 2024
Viewed by 708
Abstract
Rockfall phenomena are considered highly dangerous due to their rapid evolution and difficult prediction without applying preventive monitoring and mitigation actions. This research investigates a hazardous site in the Municipality of Vecchiano (Province of Pisa, Italy), characterized by vertical rock walls prone to [...] Read more.
Rockfall phenomena are considered highly dangerous due to their rapid evolution and difficult prediction without applying preventive monitoring and mitigation actions. This research investigates a hazardous site in the Municipality of Vecchiano (Province of Pisa, Italy), characterized by vertical rock walls prone to instability due to heavy fracturing and karst phenomena. The presence of anthropical structures and a public road at the bottom of the slopes increases the vulnerability of the site and the site’s risk. To create a comprehensive geological model of the area, Unmanned Aircraft System (UAS) photogrammetric surveys were conducted to create a 3D model useful in photointerpretation. In accessible and safe areas for personnel, engineering–geological surveys were carried out to characterize the rock mass and to define the portion of rock walls to be monitored. Results from nine multitemporal Robotic Total Station (RTS) measurement campaigns show that no monitoring prisms recorded significant displacement trends, both on the horizontal and vertical plane and in differential slope distance. Additionally, satellite Persistent Scatterer Interferometry (PSI) analysis indicates that the slopes were stable over the two years of study. The integration of these analysis techniques has proven to be an efficient solution for assessing slope stability in this specific rockfall-prone area. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Geological framework of the study site. (<b>A</b>) The site location in Italy; (<b>B</b>) the regional geological framework (Sheet n.273 “Pisa”) modified from [<a href="#B54-remotesensing-16-03848" class="html-bibr">54</a>]; (<b>C</b>) a subset of the geological map n.273 “Pisa” [<a href="#B54-remotesensing-16-03848" class="html-bibr">54</a>]; the red star in (<b>C</b>) indicates the precise location of the study area.</p>
Full article ">Figure 2
<p>Overall methodology flowchart.</p>
Full article ">Figure 3
<p>Stages of the GNSS (<b>A</b>) and RTS surveys (<b>B</b>).</p>
Full article ">Figure 4
<p>Photos in (<b>A</b>,<b>B</b>) show some examples of outcrops selected for the in situ engineering–geological survey; the yellow dashed lines in (<b>C</b>) indicate the location of scanlines along the slope; the red boxes in (<b>C</b>) show the areas where discontinuity sampling was performed through CloudCompare software (version 2).</p>
Full article ">Figure 5
<p>Robotic Total Station (RTS) on the iron plate anchored to a reinforced concrete curb (<b>A</b>); view of the rock walls to be monitored from the RTS position (<b>B</b>); macroprism and microprism utilized for the RTS multitemporal surveys (<b>C</b>).</p>
Full article ">Figure 6
<p>Location of RTS base (in yellow) and prisms; green colors indicate the reference prisms utilized to orient the monitoring system; red colors indicate the monitoring prisms periodically measured during the nine surveys.</p>
Full article ">Figure 7
<p>Example of error ellipse for the B4 monitoring prism.</p>
Full article ">Figure 8
<p>Metallic corner reflectors installed on the top edge above the rocky walls: location of corner reflectors in the study area (<b>A</b>). Photos in <b>CR1</b>, <b>CR2,</b> and <b>CR3</b> show detailed images of the corner reflectors.</p>
Full article ">Figure 9
<p>Satellite imagery covering the area of interest (<b>A</b>); satellite imagery spatially corresponding to the same area of interest (<b>B</b>). The yellow square shows the area of interest.</p>
Full article ">Figure 10
<p>Perspective view of the georeferenced and scaled 3D point cloud of the rocky walls; the scale bar only applies to (<b>A</b>). Georeferenced orthophotomosaic of the rocky walls and the alluvial plain (<b>B</b>).</p>
Full article ">Figure 11
<p>Stereographic projection (Schmidt equal-area method—lower hemisphere) of data collected during the in situ engineering–geological survey.</p>
Full article ">Figure 12
<p>Stereographic projection (Schmidt equal-area method—lower hemisphere) of data interpreted on the 3D point cloud.</p>
Full article ">Figure 13
<p>Slopes considered for the application of the SMR method and the following statistical kinematic stability analysis (<a href="#sec4dot4-remotesensing-16-03848" class="html-sec">Section 4.4</a>). The stereographic projections (Wulff equal-angle method—lower hemisphere) show an example of the executed kinematic analysis (ex., wedge sliding).</p>
Full article ">Figure 14
<p>Planimetric representation of multitemporal monitoring results (monitoring prisms from B1 to B15). The error ellipses for the monitoring prisms are indicated in orange. The points, differentiated by color, indicate the 9 multitemporal surveys.</p>
Full article ">Figure 15
<p>Planimetric representation of multitemporal monitoring results (monitoring prisms from B16 to B30). The error ellipses for the monitoring prisms are indicated in orange. The points, differentiated by color, indicate the 9 multitemporal surveys.</p>
Full article ">Figure 16
<p>Differential slope distance (<b>A</b>) and elevation displacement (<b>B</b>) of each prism measured in all the RTS surveys. The uncertainty thresholds for each prism are indicated by the red vertical bars.</p>
Full article ">Figure 17
<p>Differential slope distance (<b>A</b>) and elevation displacement (<b>B</b>) of each prism as computed with respect to R4. The uncertainty thresholds for each prism are indicated by the red vertical bars.</p>
Full article ">Figure 18
<p>Results of the PSI analysis in terms of LOS velocities (mm/yr) for the Sentinel-1A (<b>A</b>) and Sentinel-1B data (<b>B</b>). The yellow squares represent the position of the artificial corner reflectors installed at the top of the rocky slopes in this work. The blue triangle indicates the building where the RTS is installed.</p>
Full article ">Figure 19
<p>PSs identified by the regional LaMMa interferometric service for the study area. The red circle identifies the point FV6XKKY located near the RTS. The diagram at the bottom of the map shows the trend of this PS LOS velocity (mm/yr) considering the same interval of the RTS monitoring time span. The light blue ellipses indicate the acquisition dates of the 7th and 8th RTS surveys.</p>
Full article ">Figure 20
<p>Differential slope distance (<b>A</b>) and elevation displacement (<b>B</b>) of each prism as measured during the survey carried out on 4 March 2024.</p>
Full article ">
30 pages, 6012 KiB  
Article
A Remote-Sensing-Based Method Using Rockfall Inventories for Hazard Mapping at the Community Scale in the Arequipa Region of Peru
by Cassidy L. Grady, Paul M. Santi, Gabriel Walton, Carlos Luza, Guido Salas, Pablo Meza and Segundo Percy Colque Riega
Remote Sens. 2024, 16(19), 3732; https://doi.org/10.3390/rs16193732 - 8 Oct 2024
Viewed by 823
Abstract
Small communities in the Arequipa region of Peru are susceptible to rockfall hazards, which impact their lives and livelihoods. To mitigate rockfall hazards, it is first necessary to understand their locations and characteristics, which can be compiled into an inventory used in the [...] Read more.
Small communities in the Arequipa region of Peru are susceptible to rockfall hazards, which impact their lives and livelihoods. To mitigate rockfall hazards, it is first necessary to understand their locations and characteristics, which can be compiled into an inventory used in the creation of rockfall hazard rating maps. However, the only rockfall inventory available for Arequipa contains limited data of large, discrete events, which is insufficient for characterizing rockfall hazards at the community scale. A more comprehensive inventory would result in a more accurate rockfall hazard rating map—a significant resource for hazard mitigation and development planning. This study addresses this need through a remote method for rockfall hazard characterization at a community scale. Three communities located in geographically diverse areas of Arequipa were chosen for hazard inventory and characterization, with a fourth being used for validation of the method. Rockfall inventories of source zones and rockfall locations were developed using high-resolution aerial imagery, followed by field confirmation, and then predictions of runout distances using empirical models. These models closely matched the actual runout distance distribution, with all three sites having an R2 value of 0.98 or above. A semi-automated method using a GIS-based model was developed that characterizes the generation and transport of rockfall. The generation component criteria consisted of source zone height, slope angle, and rockmass structural condition. Transport was characterized by rockfall runout distance, estimated rockfall trajectory paths, and hazard ratings of corresponding source zones. The representative runout distance inventory model of the validation site matched that of a nearby site with an R2 of 0.98, despite inventorying less than a third of the number of rockfalls. This methodology improves upon current approaches and could be tested in other regions with similar climatic and geomorphic settings. These maps and methodology could be used by local and regional government agencies to warn residents of rockfall hazards, inform zoning regulations, and prioritize mitigation efforts. Full article
Show Figures

Figure 1

Figure 1
<p>A map of the Arequipa region of Peru with the locations of the three study sites, Chaparra, Aplao, and Chivay, indicated by red circles. The validation site of Achoma is indicated by a blue square.</p>
Full article ">Figure 2
<p>Photographs of the three sites and the validation site, Achoma. These photos were taken by the first author, except for that of Chaparra, which was taken by a member of the research team.</p>
Full article ">Figure 3
<p>Part of the rockfall inventory for Aplao showing marked rocks, runout zones, and source zones.</p>
Full article ">Figure 4
<p>Methodology flowchart for creating a rockfall hazard map split into “generation” and “transport” processes to determine source zone and runout zone hazard ratings that are input into the final hazard map.</p>
Full article ">Figure 5
<p>Schematic of how the program identifies the bottom of a source zone with multiple source zone polygons. The red and green arrows are bearing lines extending in the direction of slope aspect from each source zone boundary point. The red arrows indicate bearing lines that have more than one intersection, and their origin points are therefore not considered the bottom of the source zone points. The green arrows indicate bearing lines that only have one intersection point and their origin points are considered the bottom of the source zone points, which are represented by the yellow circles. The cyan line that connects the yellow circles is the bottom of the source zone line that would result from using the “Points to Line” tool in ArcGIS Pro. This schematic does not show every bearing line that would be created to maintain the simplicity of the diagram.</p>
Full article ">Figure 6
<p>Drawings of joint orientation and continuity types included in <a href="#remotesensing-16-03732-t003" class="html-table">Table 3</a>.</p>
Full article ">Figure 7
<p>Schematic of the runout zone hazard based on Equation (2) where for each cell in the runout zone, the source score is multiplied by the corresponding rockfall passing percentage and the sum of rockfall trajectory paths, which results in the final scores zoned into high, medium, and low categories.</p>
Full article ">Figure 8
<p>Percentage of rockfall passing a distance from the bottom of the source zone for the sites of Aplao (blue), Chaparra (orange), and Chivay (gray), with an inset (green box) showing the sparseness of data points in the distances that are farther from the source zone.</p>
Full article ">Figure 9
<p>Rockfall hazard map for Chivay showing the area of interest that encompasses the source zone (purple) and runout zone hazard ratings of low (yellow), medium (orange), and high (red).</p>
Full article ">Figure 10
<p>Rockfall hazard map for Aplao showing the area of interest that encompasses the source zone (purple) and runout zone hazard ratings of low (yellow), medium (orange), and high (red).</p>
Full article ">Figure 11
<p>Map of Aplao showing rockfall inventory density compared to the runout zone hazard, which is depicted as yellow (low), orange (medium), and red (high) lines that indicate the extent of each runout zone.</p>
Full article ">Figure 12
<p>Rockfall hazard map for Achoma that incorporates important community landmarks (the Colca River and the Majes Canal), as well as a line with teeth for the top of the source zone to indicate which direction is downslope.</p>
Full article ">
22 pages, 19530 KiB  
Article
Cascading Landslide: Kinematic and Finite Element Method Analysis through Remote Sensing Techniques
by Claudia Zito, Massimo Mangifesta, Mirko Francioni, Luigi Guerriero, Diego Di Martire, Domenico Calcaterra and Nicola Sciarra
Remote Sens. 2024, 16(18), 3423; https://doi.org/10.3390/rs16183423 - 14 Sep 2024
Viewed by 774
Abstract
Cascading landslides are specific multi-hazard events in which a primary movement triggers successive landslide processes. Areas with dynamic and quickly changing environments are more prone to this type of phenomena. Both the kind and the evolution velocity of a landslide depends on the [...] Read more.
Cascading landslides are specific multi-hazard events in which a primary movement triggers successive landslide processes. Areas with dynamic and quickly changing environments are more prone to this type of phenomena. Both the kind and the evolution velocity of a landslide depends on the materials involved. Indeed, rockfalls are generated when rocks fall from a very steep slope, while debris flow and/or mudslides are generated by fine materials like silt and clay after strong water imbibition. These events can amplify the damage caused by the initial trigger and propagate instability along a slope, often resulting in significant environmental and societal impacts. The Morino-Rendinara cascading landslide, situated in the Ernici Mountains along the border of the Abruzzo and Lazio regions (Italy), serves as a notable example of the complexities and devastating consequences associated with such events. In March 2021, a substantial debris flow event obstructed the Liri River, marking the latest step in a series of landslide events. Conventional techniques such as geomorphological observations and geological surveys may not provide exhaustive information to explain the landslide phenomena in progress. For this reason, UAV image acquisition, InSAR interferometry, and pixel offset analysis can be used to improve the knowledge of the mechanism and kinematics of landslide events. In this work, the interferometric data ranged from 3 January 2020 to 24 March 2023, while the pixel offset data covered the period from 2016 to 2022. The choice of such an extensive data window provided comprehensive insight into the investigated events, including the possibility of identifying other unrecorded events and aiding in the development of more effective mitigation strategies. Furthermore, to supplement the analysis, a specific finite element method for slope stability analysis was used to reconstruct the deep geometry of the system, emphasizing the effect of groundwater-level flow on slope stability. All of the findings indicate that major landslide activities were concentrated during the heavy rainfall season, with movements ranging from several centimeters per year. These results were consistent with numerical analyses, which showed that the potential slip surface became significantly more unstable when the water table was elevated. Full article
Show Figures

Figure 1

Figure 1
<p>Aerial and field images of the Morino-Rendinara landslide that are representative of the impact of the landslide on the environment. (<b>a</b>) Overview of phenomenon taken from Google Earth [<a href="#B16-remotesensing-16-03423" class="html-bibr">16</a>] satellite images of 13 June 2022, from the upper sector near Morino Hamlet to the lower sector, Liri River, and deep-seated rotational slide; (<b>b</b>) Details of rockfall/avalanches sector; (<b>c</b>) Debris flow source area; (<b>d</b>) Debris flow transit zone; (<b>e</b>) Lowest debris flow transit zone; (<b>f</b>) Liri River dam; and (<b>g</b>) Effect on Liri River dam.</p>
Full article ">Figure 2
<p>Geographical location of Morino-Rendinara. Green lines indicate the regional boundaries; red lines indicate the municipality of Morino, Castronovo, and San Vincenzo Valle Roveto composing the involved municipality; the light blue square indicates the landslide and the study area.</p>
Full article ">Figure 3
<p>Geological map extract from the CARG (Geological CARtography Map n.220 Sora) Project [<a href="#B22-remotesensing-16-03423" class="html-bibr">22</a>], with indications of the geological formations and tectonic processes present in the area.</p>
Full article ">Figure 4
<p>Maps of the survey and debritic cover layer reconstruction using cross-sections to empathize the heterogeneity of deposits covering the substrate. (a) The section develops on maximum slope line. (b) The section develops on perpendicular direction.</p>
Full article ">Figure 5
<p>Conceptual flow chart of the work phases.</p>
Full article ">Figure 6
<p>Spatiotemporal baseline map of SBAS-InSAR interferometric data of ascending track (<b>a</b>) and descending track (<b>b</b>).</p>
Full article ">Figure 7
<p>The inventory map was drawn using the results of the study. The image identifies three main mechanisms: a rockfall in the upper part, a deep-seated rotational slide in the central part, and a debris flow in the lower part.</p>
Full article ">Figure 8
<p>Details of the landslide inventory map. The various fillings show the different landslide types identified in the study area: the rockfall in the upper part, the deep-seated rotational slide in the central part, and the debris flow in the lower part of the slope.</p>
Full article ">Figure 9
<p>PS velocity along the ascending (<b>a</b>) and descending (<b>b</b>) geometries from 2020 to 2023. Red dots indicate major velocity trends and instability, green and blue dots indicate minor velocity and stable sectors.</p>
Full article ">Figure 10
<p>Selected time series in ascending geometry (<b>a</b>) and descending geometry (<b>b</b>). The analyzed time series illustrates a very unstable sector represented by reflectors P106_60-61 and 102_57 in ascending geometry and P_70_141-141-136 in descending geometry. Additionally, some stable sectors are represented, such as P_83_135-143, P_85_68, and P_86_69.</p>
Full article ">Figure 11
<p>Time series vs. rainfall analyses. (<b>a</b>) Monthly cumulative rainfall for analysis period vs. one ascending and descending representative time series; (<b>b</b>) Daily cumulative rainfall for analysis period vs. one ascending and descending representative time series; (<b>c</b>) Cumulative rainfall for analysis period vs. one ascending and descending representative time series.</p>
Full article ">Figure 12
<p>Shear strain index results of the 2D numerical modelling with different pore pressure conditions. (<b>a</b>) Analysis without pore pressure; (<b>b</b>) Analysis with water table 0.5 m from ground level; (<b>c</b>) Analysis with water table 2.0 m from ground level; (<b>d</b>) Analysis with water table 6.0 m from ground level.</p>
Full article ">
18 pages, 7423 KiB  
Article
Leveraging Internet News-Based Data for Rockfall Hazard Susceptibility Assessment on Highways
by Kieu Anh Nguyen, Yi-Jia Jiang, Chiao-Shin Huang, Meng-Hsun Kuo and Walter Chen
Future Internet 2024, 16(8), 299; https://doi.org/10.3390/fi16080299 - 21 Aug 2024
Viewed by 754
Abstract
Over three-quarters of Taiwan’s landmass consists of mountainous slopes with steep gradients, leading to frequent rockfall hazards that obstruct traffic and cause injuries and fatalities. This study used Google Alerts to compile internet news on rockfall incidents along Taiwan’s highway system from April [...] Read more.
Over three-quarters of Taiwan’s landmass consists of mountainous slopes with steep gradients, leading to frequent rockfall hazards that obstruct traffic and cause injuries and fatalities. This study used Google Alerts to compile internet news on rockfall incidents along Taiwan’s highway system from April 2019 to February 2024. The locations of these rockfalls were geolocated using Google Earth and integrated with geographical, topographical, environmental, geological, and socioeconomic variables. Employing machine learning algorithms, particularly the Random Forest algorithm, we analyzed the potential for rockfall hazards along roadside slopes. The model achieved an overall accuracy of 0.8514 on the test dataset, with a sensitivity of 0.8378, correctly identifying 83.8% of rockfall locations. Shapley Additive Explanations (SHAP) analysis highlighted that factors such as slope angle and distance to geologically sensitive areas are pivotal in determining rockfall locations. The study underscores the utility of internet-based data collection in providing comprehensive coverage of Taiwan’s highway system, and enabled the first broad analysis of rockfall hazard susceptibility for the entire highway network. The consistent importance of topographical and geographical features suggests that integrating detailed spatial data could further enhance predictive performance. The combined use of Random Forest and SHAP analyses offers a robust framework for understanding and improving predictive models, aiding in the development of effective strategies for risk management and mitigation in rockfall-prone areas, ultimately contributing to safer and more reliable transportation networks in mountainous regions. Full article
Show Figures

Figure 1

Figure 1
<p>The highway system in Taiwan, highlighting the locations of rockfall and non-rockfall events.</p>
Full article ">Figure 2
<p>Predictive variables used for rockfall analysis: (<b>a</b>) elevation, (<b>b</b>) slope angle in degrees, (<b>c</b>) NDVI, (<b>d</b>) river system, (<b>e</b>) urban areas, (<b>f</b>) fault lines, (<b>g</b>) geologically sensitive areas, (<b>h</b>) potential debris flow areas, (<b>i</b>) average yearly rainfall, (<b>j</b>) average monthly temperature, (<b>k</b>) average monthly wind speed, (<b>l</b>) average monthly humidity, (<b>m</b>) average number of rainy days, (<b>n</b>) strata, (<b>o</b>) epoch, (<b>p</b>) population density.</p>
Full article ">Figure 2 Cont.
<p>Predictive variables used for rockfall analysis: (<b>a</b>) elevation, (<b>b</b>) slope angle in degrees, (<b>c</b>) NDVI, (<b>d</b>) river system, (<b>e</b>) urban areas, (<b>f</b>) fault lines, (<b>g</b>) geologically sensitive areas, (<b>h</b>) potential debris flow areas, (<b>i</b>) average yearly rainfall, (<b>j</b>) average monthly temperature, (<b>k</b>) average monthly wind speed, (<b>l</b>) average monthly humidity, (<b>m</b>) average number of rainy days, (<b>n</b>) strata, (<b>o</b>) epoch, (<b>p</b>) population density.</p>
Full article ">Figure 3
<p>The ROC curve of the Random Forest model, highlighting the AUC value of 0.9317, which indicates excellent model performance.</p>
Full article ">Figure 4
<p>Rockfall susceptibility map of a section of Provincial Highway 7, showing the classification into five susceptibility levels. Black dots indicate the locations of past rockfalls.</p>
Full article ">Figure 5
<p>Variable importance derived from the Random Forest model, highlighting the significant predictors across geographical, topographical, environmental, geological, and socioeconomic categories.</p>
Full article ">Figure 6
<p>Mean SHAP values indicating the average impact of each feature on model predictions. Features are ranked by importance, with <tt>slope_deg</tt> and <tt>dist_geo_sens_area</tt> showing the highest impact.</p>
Full article ">Figure 7
<p>SHAP values illustrating the impact of each feature on the model’s output. Red dots represent high feature values, while blue dots represent low feature values.</p>
Full article ">
21 pages, 34311 KiB  
Case Report
Drone-Borne LiDAR and Photogrammetry Together with Historical Data for Studying a Paleo-Landslide Reactivated by Road-Cutting and Barrier Construction outside Jerusalem
by Yaniv Darvasi, Ben Laugomer, Ido Shicht, John K. Hall, Eli Ram and Amotz Agnon
Geotechnics 2024, 4(3), 786-806; https://doi.org/10.3390/geotechnics4030041 - 9 Aug 2024
Viewed by 713
Abstract
Assessment of landslide hazards often depends on the ability to track possible changes in natural slopes. To that end, historical air photos can be useful, particularly when slope stability is compromised by visible cracking. Undocumented landsliding rejuvenates a paleo-landslide on a busy motorway [...] Read more.
Assessment of landslide hazards often depends on the ability to track possible changes in natural slopes. To that end, historical air photos can be useful, particularly when slope stability is compromised by visible cracking. Undocumented landsliding rejuvenates a paleo-landslide on a busy motorway connecting Jerusalem to a small Jewish settlement. Recently, a plan for broadening the motorway was approved, and we were asked to study the hazards of the road by Israeli NGOs and Palestinian residents of the area. We captured high-resolution topography around the unstable slope using drone-borne photogrammetry and LiDAR surveys. The modern data allow us to analyze historic air photos and topo maps to assess the level of sliding prior to and during modern landscaping. Our results indicate horizontal offsets of ~0.9–1.8 m and vertical offsets of 1.54–2.95 m at selected sites. We next assess the possible role of anthropogenic versus natural factors in compromising slope stability. We analyze monthly rain records together with seismic catalogs spanning several decades. Shortly after the motorway construction in 1995, a January 1996 rainstorm triggered a massive rockfall. The rockfall blocked traffic with up to 4 m-diameter boulders. We found that while a certain level of rain is a necessary condition for mobilizing the rock mass, it is the anthropogenic intervention that caused the rockfall in this site. We conclude that the recent plan for broadening the motorway jeopardizes the lives of vehicle passengers and the lives of future residents should the development materialize. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location map (the source of the two upper maps is Open Topo Map). The aerial photo (bottom) was taken by a drone, and the map was compiled with QGIS software 3.22).</p>
Full article ">Figure 2
<p>The rockfall of 24 January 1996, shortly after Motorway 385 was completed. (<b>A</b>) From a news release on the daily Yedi’ot A’hronot; (<b>B</b>) slightly later, after the 3 m-wide boulder in the foreground (rolled through the protective concrete wall) was shattered for removal (photo by Yossi Leshem). The newspaper (from the day following the fall) noted that a real disaster was prevented when a four-metric-ton rock rolled onto the road, as traffic had been stopped earlier when smaller rocks fell; residents and Society for Nature Protection people had warned that, due to the haste to evacuate Bethlehem, the road was poorly constructed and was dangerous for traffic.</p>
Full article ">Figure 3
<p>1:50,000 Geological map [<a href="#B4-geotechnics-04-00041" class="html-bibr">4</a>]. Stratigraphic and lithological symbols are given in <a href="#geotechnics-04-00041-f004" class="html-fig">Figure 4</a>. The blue quadrangle demarcates the slide site. Motorway 385 is red. White circles mark rock fall blocks mapped by Katz et al., 2011 [<a href="#B1-geotechnics-04-00041" class="html-bibr">1</a>]. the inscriptions in Hebrew symbolize ancient sites, springs and streams.</p>
Full article ">Figure 4
<p>Stratigraphy [<a href="#B4-geotechnics-04-00041" class="html-bibr">4</a>].</p>
Full article ">Figure 5
<p>(<b>A</b>). A LiDAR point cloud of an open crevasse with depth reaching up to ~6.5 m. (<b>B</b>,<b>C</b>) An open crack in a dirt road, cut for installing a barrier around 2012; (<b>B</b>) A person halfway within the crevasse. (<b>C</b>) Zoom-in with a 0.4 m geological pick for scale. Notice the paleo-boulder within the breccia rock mass, split open by the new crevasse (photos taken 15 December 2020).</p>
Full article ">Figure 6
<p>LiDAR sensors: (<b>Left</b>), GeoSlam Zeb Horizon LiDAR system mounted on a Matrice 600 Pro UAV. (<b>Right</b>), DJI L1 LiDAR system mounted on a Matrice 300 RTK UAV.</p>
Full article ">Figure 7
<p>Orthophotos since 1970. The 1970 photo shows mild terracing, a traditional method practiced by farmers for soil accumulation and preservation (e.g., [<a href="#B6-geotechnics-04-00041" class="html-bibr">6</a>]).</p>
Full article ">Figure 8
<p>Contour maps: The yellow ellipse marks the area where the slide has modified the shapes of the contour map. (<b>A</b>) Contours traced from British Mandate (pre-1948) maps by Hall (2008) (25 m DTM) [<a href="#B5-geotechnics-04-00041" class="html-bibr">5</a>], superimposed on an orthophoto acquired by the drone. (<b>B</b>,<b>C</b>) are from <a href="https://amudanan.co.il/" target="_blank">https://amudanan.co.il/</a> (1:50,000). Contours in the more recent map (<b>C</b>) show inverted concavity where the slide intercepted the surface. This concavity controlled by the landslide head scarp is more pronounced on the 5 m contour map of the Israel Mapping Center.</p>
Full article ">Figure 9
<p>Orthophoto with the five main areas where movements have been recognized.</p>
Full article ">Figure 10
<p>LiDAR 2021/2022 comparison of areas 1 and 2 indicate the difference in topography (location in <a href="#geotechnics-04-00041-f009" class="html-fig">Figure 9</a>).</p>
Full article ">Figure 11
<p>Photogrammetry 2021/2022 comparison of main differences (location in <a href="#geotechnics-04-00041-f009" class="html-fig">Figure 9</a>).</p>
Full article ">Figure 12
<p>LiDAR 2021/2022 comparison of area 3 indicates the difference in topography (location in <a href="#geotechnics-04-00041-f009" class="html-fig">Figure 9</a>).</p>
Full article ">Figure 13
<p>LiDAR 2021/2022 comparison of area 5 indicates the difference in topography (location in <a href="#geotechnics-04-00041-f009" class="html-fig">Figure 9</a>).</p>
Full article ">Figure A1
<p>The amount of monthly rain (only for rainy months) measured in Jerusalem between 1967 and 2023. Median, average, and average + standard deviation are marked green, orange, and violet, respectively. The storm preceding the rockfall is outlined (Jan 1996, 157.5 mm). While this storm was significant, it was far from unusual (as compared to the Pinatubo winter marked—1991/2). We suggest that the combination of the storm together with the landscaping and road cutting for Motorway 385, with possible contribution of a regional earthquake, triggered the rock fall. Data from Israel Meteorological Service (IMS).</p>
Full article ">Figure A2
<p>Annual average of rain measured in Jerusalem. The winter of 1995/6 was slightly below average. Also shown (top) are the years for which we present air photos (<a href="#geotechnics-04-00041-f007" class="html-fig">Figure 7</a>). Data from Israel Meteorological Service (IMS).</p>
Full article ">Figure A3
<p>Epicenters of earthquakes with a magnitude larger than 5.0 around the study site (yellow star) between 1900 and 2023. The largest earthquake (M7.2) hit the Gulf of Aqaba on 22 November 1995. Data from the Geological Survey of Israel <a href="https://eq.gsi.gov.il/earthquake" target="_blank">https://eq.gsi.gov.il/earthquake</a> (accessed on 2 March 2024).</p>
Full article ">Figure A4
<p>(<b>Left</b>): A portion of a map (1:50,000) marking intermediate potential for ground amplification during earthquakes [<a href="#B10-geotechnics-04-00041" class="html-bibr">10</a>], with the landslide site marked by a blue star. The potential was calculated based on slope and lithology. The semi-circular band about the site is formed by the marly-clayey Motza Formation (<a href="#geotechnics-04-00041-f004" class="html-fig">Figure 4</a>). The blue square and star, respectively, mark the sloping hill and the slide. (<b>Right</b>): A portion of the geological map (<a href="#geotechnics-04-00041-f003" class="html-fig">Figure 3</a>) where the blue square is marked for ease of orientation between panels. The inscriptions in Hebrew symbolize ancient sites, springs and streams.</p>
Full article ">Figure A5
<p>A traditional agricultural installation carved in the carbonate cemented breccia (likely an olive or grape press), Coord. 31.731274/35.148377. A crack is cutting through the stone-carved installation (yellow ellipse). We suggest that this is a part of the ongoing cracking triggered by the landscaping related to Motorway 385 and the pedestrian barrier.</p>
Full article ">Figure A6
<p>Field photos.</p>
Full article ">
24 pages, 24217 KiB  
Article
Evaluating the Impact of DEM Spatial Resolution on 3D Rockfall Simulation in GIS Environment
by Maria P. Kakavas, Paolo Frattini, Alberto Previati and Konstantinos G. Nikolakopoulos
Geosciences 2024, 14(8), 200; https://doi.org/10.3390/geosciences14080200 - 29 Jul 2024
Cited by 1 | Viewed by 883
Abstract
Rockfalls are natural geological phenomena characterized by the abrupt detachment and freefall descent of rock fragments from steep slopes. These events exhibit considerable variability in scale, velocity, and trajectory, influenced by the geological composition of the slope, the topography, and other environmental conditions. [...] Read more.
Rockfalls are natural geological phenomena characterized by the abrupt detachment and freefall descent of rock fragments from steep slopes. These events exhibit considerable variability in scale, velocity, and trajectory, influenced by the geological composition of the slope, the topography, and other environmental conditions. By employing advanced modeling techniques and terrain analysis, researchers aim to predict and control rockfall hazards to prevent casualties and protect properties in areas at risk. In this study, two rockfall events in the villages of Myloi and Platiana of Ilia prefecture were examined. The research was conducted by means of HY-STONE software, which performs 3D numerical modeling of the motion of non-interacting blocks. To perform this modeling, input files require the processing of base maps and datasets in a GIS environment. Stochastic modeling and 3D descriptions of slope topography, based on Digital Elevation Models (DEMs) without spatial resolution limitations, ensure multiscale analysis capabilities. Considering this capability, seven freely available DEMs, derived from various sources, were applied in HY-STONE with the scope of performing a large number of multiparametric analyses and selecting the most appropriate and efficient DEM for the software requirements. All the necessary data for the multiparametric analyses were generated within a GIS environment, utilizing either the same restitution coefficients and rolling friction coefficient or varying ones. The results indicate that finer-resolution DEMs capture detailed terrain features, enabling the precise identification of rockfall source areas and an accurate depiction of the kinetic energy distribution. Further, the results show that a correct application of the model to different DEMs requires a specific parametrization to account for the different roughness of the models. Full article
(This article belongs to the Special Issue Earth Observation by GNSS and GIS Techniques)
Show Figures

Figure 1

Figure 1
<p>Pictures of the rockfall events in the Myloi and Platiana region in relation to the Hellenic region.</p>
Full article ">Figure 2
<p>Pictures of the rockfall in the Myloi region illustrate the final position of rock blocks. Red arrows indicate the final positions of the rock masses at the conclusion of the rockfall event.</p>
Full article ">Figure 3
<p>Picture taken from the Platiana region after the rockfall event. The slope picture is taken from the Greek Cadastral (2008).</p>
Full article ">Figure 4
<p>Reforestation procedures in the slope next to the Platiana village. The slope pictures are taken from the Greek Cadastral.</p>
Full article ">Figure 5
<p>Myloi area before (<b>left</b> image) and after (<b>right</b> image) removing the vegetation through Cloud Compare Software.</p>
Full article ">Figure 6
<p>HY-STONE results for the site of Myloi using parameters (E<sub>N</sub>, E<sub>T</sub>, and A<sub>T</sub>) calibrated on the UAV DEM. The transit frequencies of blocks are shown for (<b>a</b>) UAV DEM, (<b>b</b>) Greek Cadastral DEM, (<b>c</b>) ALOS AW3D30 DEM, (<b>d</b>) ASTER GDEM, and (<b>e</b>) SRTM30 DEM.</p>
Full article ">Figure 7
<p>HY-STONE results for the site of Myloi using parameters (E<sub>N</sub>, E<sub>T</sub>, and A<sub>T</sub>) calibrated on the UAV DEM. The maximum translation kinetic energies are shown for (<b>a</b>) UAV DEM, (<b>b</b>) Greek Cadastral DEM, (<b>c</b>) ALOS AW3D30 DEM, (<b>d</b>) ASTER GDEM, and (<b>e</b>) SRTM30 DEM.</p>
Full article ">Figure 8
<p>Slope profiles from (<b>a</b>) UAV DSE, (<b>b</b>) Greek Cadastral DEM, (<b>c</b>) ALOS AW3D30 DEM, (<b>d</b>) ASTER GDEM, and (<b>e</b>) SRTM30 DEM.</p>
Full article ">Figure 9
<p>HY-STONE results, in terms of transit frequency for the site of Platiana. (<b>a</b>) Greek Cadastral DEM, (<b>b</b>) ALOS AW3D30 DEM, (<b>c</b>) ASTER GDEM, (<b>d</b>) SRTM30 DEM, (<b>e</b>) SRTM90 DEM, and (<b>f</b>) TanDEM_X.</p>
Full article ">Figure 10
<p>HY-STONE results in terms of the maximum translation kinetic energy for the site of Platiana. (<b>a</b>) Greek Cadastral DEM, (<b>b</b>) ALOS AW3D30 DEM, (<b>c</b>) ASTER GDEM, (<b>d</b>) SRTM30 DEM, (<b>e</b>) SRTM90 DEM, and (<b>f</b>) TanDEM_X.</p>
Full article ">Figure 11
<p>Slope profiles from (<b>a</b>) Greek Cadastral DEM, (<b>b</b>) ALOS AW3D30 DEM, (<b>c</b>) ASTER GDEM, (<b>d</b>) SRTM30 DEM, (<b>e</b>) SRTM90 DEM, and (<b>f</b>) TanDEM_X.</p>
Full article ">Figure 12
<p>HY-STONE results by using the Greek Cadastral DEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t003" class="html-table">Table 3</a> for the site of Myloi: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 13
<p>HY-STONE results by using the ALOS AW3D30 DEM with optimal coefficients shown in <a href="#geosciences-14-00200-t004" class="html-table">Table 4</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 14
<p>HY-STONE results by using the ASTER GDEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t005" class="html-table">Table 5</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 15
<p>HY-STONE results by using the SRTM30 DEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t006" class="html-table">Table 6</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 16
<p>HY-STONE results by using the SRTM90 DEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t007" class="html-table">Table 7</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 17
<p>HY-STONE results by using the TanDEM_X with the optimal coefficients shown in <a href="#geosciences-14-00200-t008" class="html-table">Table 8</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 18
<p>The effect of spatial resolution in the normal (E<sub>N</sub>) and tangential (E<sub>T</sub>) restitutions and the rolling friction (A<sub>T</sub>) coefficients. For 30 m and 90 m, three and two DEMS are available, respectively.</p>
Full article ">
18 pages, 9240 KiB  
Article
Identification and Analysis of the Geohazards Located in an Alpine Valley Based on Multi-Source Remote Sensing Data
by Yonglin Yang, Zhifang Zhao, Dingyi Zhou, Zhibin Lai, Kangtai Chang, Tao Fu and Lei Niu
Sensors 2024, 24(13), 4057; https://doi.org/10.3390/s24134057 - 21 Jun 2024
Viewed by 1005
Abstract
Geohazards that have developed in densely vegetated alpine gorges exhibit characteristics such as remote occurrence, high concealment, and cascading effects. Utilizing a single remote sensing datum for their identification has limitations, while utilizing multiple remote sensing data obtained based on different sensors can [...] Read more.
Geohazards that have developed in densely vegetated alpine gorges exhibit characteristics such as remote occurrence, high concealment, and cascading effects. Utilizing a single remote sensing datum for their identification has limitations, while utilizing multiple remote sensing data obtained based on different sensors can allow comprehensive and accurate identification of geohazards in such areas. This study takes the Latudi River valley, a tributary of the Nujiang River in the Hengduan Mountains, as the research area, and comprehensively uses three techniques of remote sensing: unmanned aerial vehicle (UAV) Light Detection and Ranging (LiDAR), Small Baseline Subset interferometric synthetic aperture radar (SBAS-InSAR), and UAV optical remote sensing. These techniques are applied to comprehensively identify and analyze landslides, rockfalls, and debris flows in the valley. The results show that a total of 32 geohazards were identified, including 18 landslides, 8 rockfalls, and 6 debris flows. These hazards are distributed along the banks of the Latudi River, significantly influenced by rainfall and distribution of water systems, with deformation variables fluctuating with rainfall. The three types of geohazards cause cascading disasters, and exhibit different characteristics in the 0.5 m resolution hillshade map extracted from LiDAR data. UAV LiDAR has advantages in densely vegetated alpine gorges: after the selection of suitable filtering algorithms and parameters of the point cloud, it can obtain detailed terrain and geomorphological information on geohazards. The different remote sensing technologies used in this study can mutually confirm and complement each other, enhancing the capability to identify geohazards and their associated hazard cascades in densely vegetated alpine gorges, thereby providing valuable references for government departments in disaster prevention and reduction work. Full article
(This article belongs to the Topic Advanced Risk Assessment in Geotechnical Engineering)
Show Figures

Figure 1

Figure 1
<p>Overview of the geographical location of the Latudi River.</p>
Full article ">Figure 2
<p>Geological overview of the Latudi River valley. (<b>a</b>) The bedrock along the valley bottom is exposed; (<b>b</b>) development of fissures in the rock formations.</p>
Full article ">Figure 3
<p>Aerial photo of the lower reaches of the Latudi River in May 2014. Photo from Yunnan Bureau of Surveying, Mapping and Geoinformation.</p>
Full article ">Figure 4
<p>SBAS-InSAR main processing flow. (<b>a</b>) Coverage area of each acquired Sentinel-1A image; (<b>b</b>) SBAS-InSAR spatial baseline map; (<b>c</b>) phase interference map; (<b>d</b>) phase unwrapping map.</p>
Full article ">Figure 5
<p>Processing process of point cloud data.</p>
Full article ">Figure 6
<p>Diagram of IPTD filtering principle.</p>
Full article ">Figure 7
<p>Deformation velocity of Latudi River valley based on SBAS-InSAR.</p>
Full article ">Figure 8
<p>Identified geohazard areas A, B, and C. (<b>a</b>) Deformation rate map of areas A, B, and D; (<b>b</b>) UAV image of areas A, B, and D; (<b>c</b>) deformation rate map of area C; (<b>d</b>) UAV image of area C.</p>
Full article ">Figure 9
<p>Identification results and examples of geohazards in the Latudi River valley.</p>
Full article ">Figure 10
<p>Comparison of identification results from different remote sensing methods. (<b>a</b>) Identification results from SBAS-InSAR. (<b>b</b>) Identification results from LiDAR. (<b>c</b>) Identification results from UAV optical image.</p>
Full article ">Figure 11
<p>Landslide features. (<b>a</b>) Hillshade map; (<b>b</b>) UAV image; (<b>c</b>) field survey photo—landslide rear wall; (<b>d</b>) field survey photo—landslide front edge.</p>
Full article ">Figure 12
<p>Rockfall features. (<b>a</b>) Three-dimensional hillshade map; (<b>b</b>) UAV image; (<b>c</b>) three-dimensional image; (<b>d</b>) field survey photo—landslide site; (<b>e</b>) field survey photo—collapsed debris; (<b>f</b>) field survey photo—rockfall surface.</p>
Full article ">Figure 13
<p>Debris flow features. (<b>a</b>) Hillshade map; (<b>b</b>) UAV image; (<b>c</b>) 3D scene.</p>
Full article ">Figure 14
<p>The relationship between deformation points P1, P2, and P3 and monthly average rainfall.</p>
Full article ">Figure 15
<p>Based on the IPTD algorithm, the filtering effects of steep terrain obtained with different angle and distance thresholds are displayed. (<b>a</b>) Angle/distance: 55°/1.5 m; (<b>b</b>) angle/distance: 30°/0.8 m; (<b>c</b>) angle/distance: 30°/1.6 m; (<b>d</b>) angle/distance: 30°/1.6 m and 45°/1.5 m.</p>
Full article ">Figure 16
<p>Cascading relationship of geohazards in the Latudi River valley.</p>
Full article ">
24 pages, 14076 KiB  
Article
The Interconnection between Climate Cycles and Geohazards in Urban Areas of the Tourist Island of Mallorca, Spain
by Juan A. Luque-Espinar, Rosa M. Mateos, Roberto Sarro, Cristina Reyes-Carmona and Mónica Martínez-Corbella
Sustainability 2024, 16(12), 4917; https://doi.org/10.3390/su16124917 - 7 Jun 2024
Viewed by 998
Abstract
The island of Mallorca has experienced major interventions and transformations of the territory, with unprecedented urban development related to growing tourism activity. In this paper, we present a spatio-temporal analysis—by using spectral analysis techniques—of climate cycles on the island of Mallorca (Spain) and [...] Read more.
The island of Mallorca has experienced major interventions and transformations of the territory, with unprecedented urban development related to growing tourism activity. In this paper, we present a spatio-temporal analysis—by using spectral analysis techniques—of climate cycles on the island of Mallorca (Spain) and their correlation with the occurrences of landslides and flash floods. Both geohazards are closely related to wet periods, which are controlled by different, well-known natural cycles: ENSO, the NAO, sunspot, etc. Geostatistical methods are used to map the distribution of rainfall, as well as a spatial representation of the spectral confidence of the different natural cycles, to define the hazardous areas on the island. The cycles with the greatest influence on rainfall in Mallorca are El Niño–Southern Oscillation (ENSO) (5.6 y and 3.5 y), the North Atlantic Oscillation (NAO) (7.5 y), and Quasi-Biennial Oscillation (QBO). Recorded events of both rockfalls and flash floods exhibit a strong correlation with the climate indices of QBO, ENSO, the NAO, and sunspot activity. This correlation is particularly pronounced with QBO, as this cycle has a higher frequency than the others, and QBO is observed as part of the other cycles in the form of increases and decreases during periods of higher ENSO, NAO, and sunspot values. However, the impact of flash floods is also significant in the southeast part of the island, despite its lower levels of rainfall. The most dangerous episodes are related to ENSO (6.4 y) and the NAO. The validation of the methodology employed is strengthened by incorporating information from the flash flood data, as it offers comprehensive coverage of the entire island, compared to the landslide database, which is confined to the Serra de Tramuntana region. The study reveals that the city of Palma and the municipality of Calvià, as well as the central and eastern urban areas of the island, are the most vulnerable regions to intense rainfall and its consequences. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) Location of Mallorca in the western Mediterranean. (<b>B</b>) Rainfall stations used in the present work. (<b>C</b>) The main geomorphological domains on the island, highlighting the Tramuntana range in the northwestern part. The territory designated for urban use in Mallorca is indicated.</p>
Full article ">Figure 2
<p>Workflow of the methodology.</p>
Full article ">Figure 3
<p>Distribution of the 423 landslide events registered during the past 30 years (1991–2020) in the Tramuntana range of Mallorca and the temporal location of the five largest rockfalls (&gt;28,000 m<sup>3</sup>). There are no records for the first 2 years. Photos of two of the largest rockfalls that took place during the rainy cold period 2008–2010. Right: The Son Cocó rock avalanche in December 2008, the largest rockfall (300,000 m<sup>3</sup>). Left: The Gorg Blau rockfall (30,000 m<sup>3</sup>), also in December 2008.</p>
Full article ">Figure 4
<p>(<b>a</b>) Correlation matrix of the variables in <a href="#sustainability-16-04917-t003" class="html-table">Table 3</a>. (<b>b</b>) Logistic regression model of max_rainfall and QBO.</p>
Full article ">Figure 5
<p>Rainfall records and results of the spectral analysis. Blue line &gt; 99% statistical confidence. Green line &gt; 95% spectral confidence. Orange line &gt; 90% statistical confidence. Red line &lt; 90% statistical confidence.</p>
Full article ">Figure 6
<p>Main estimated climate cycles and spectral confidence values.</p>
Full article ">Figure 7
<p>Spatial estimation of the spectral confidence of the climatic cycles with the greatest influence on the island of Mallorca using geostatistical methods. (<b>a</b>–<b>d</b>) have been estimated by Ordinary Kriging and (<b>e</b>–<b>g</b>) have been estimated by Indicator Kriging. (<b>a’</b>,<b>b’</b>,<b>e’</b>,<b>g’</b>) show the estimation error. The rockfalls are located and classified by volume. Rainfall: data from weather stations and associated with flash flood events recorded between 1932 and 2010 are presented in <a href="#sustainability-16-04917-t003" class="html-table">Table 3</a>.</p>
Full article ">Figure 8
<p>Theoretical variograms fitted for each of the dates selected for the spatial estimation of precipitation. Red cross: experimental variogram of wet-type year. Blue cross: experimental variogram of dry-type year.</p>
Full article ">Figure 9
<p>Estimated rainfall in mm in the year 2008 (<b>a</b>) and wet-type year (<b>b</b>), with similar values to those of the wet-type year. The landslides are located and classified by volume.</p>
Full article ">Figure 10
<p>Superposition over urban areas of the most influential climatic cycles on the island (QBO; 3.5 y, 7.5 y, and 5.6 y cycles), where boundaries are defined according to the value of the third quartile of the estimated data range. The landslides are located (classified by volume), as well as extreme rainfall events recorded during the spanning period 2015–2022.</p>
Full article ">
21 pages, 10718 KiB  
Article
A Comprehensive Approach to Quantitative Risk Assessment of Rockfalls on Buildings Using 3D Model of Rockfall Runout
by Mohammad Al-Shaar, Pierre-Charles Gerard, Ghaleb Faour, Walid Al-Shaar and Jocelyne Adjizian-Gérard
J 2024, 7(2), 183-203; https://doi.org/10.3390/j7020011 - 30 May 2024
Viewed by 1277
Abstract
Rockfalls are incidents of nature that take place when rocks or boulders break from a steep slope and fall to the ground. They can pose considerable threats to buildings placed in high-risk zones. Despite the fact that the impact of a rockfall on [...] Read more.
Rockfalls are incidents of nature that take place when rocks or boulders break from a steep slope and fall to the ground. They can pose considerable threats to buildings placed in high-risk zones. Despite the fact that the impact of a rockfall on a building can cause structural and non-structural damage, few studies have been undertaken to investigate the danger associated with this event. Most of these studies indicated that the risk resulting from rockfall hazards is hard to forecast and assess. A comprehensive quantitative risk assessment approach for rockfalls on buildings is developed and described in this paper and applied for the Mtein village in Mount Lebanon. This method employs a 3D model to simulate the rockfall trajectories using a combination of digital elevation data, field surveys, and orthorectified aerial photographs. The spatial and temporal probability of rockfalls were evaluated using the analysis of historical data in two triggering-factor scenarios: earthquake and precipitation. The findings show that, during the period of 1472 years between the years 551 (the first observed large earthquake in Lebanon) and the current year of the study (2023), the temporal probability will potentially be equal to 0.002 and 0.105 in the cases of earthquake- and rainfall-triggered rockfalls, respectively, while the maximal damage values are expected to be 232 USD and 10,511 USD per year, respectively. The end result is a final map presenting the risk values assigned to each building that could be damaged by rockfalls. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart representing the methodology.</p>
Full article ">Figure 2
<p>Map of the study area showing RF distribution, RF release areas, fallen rock sizes, maximum observed reaching point, and location of buildings (built before and after the year 2000).</p>
Full article ">Figure 3
<p>Reach probability in earthquake case.</p>
Full article ">Figure 4
<p>Reach probability in rainfall case.</p>
Full article ">Figure 5
<p>Mean kinetic energy of rockfalls—the case of earthquake.</p>
Full article ">Figure 6
<p>Mean kinetic energy of rockfalls—the case of rainfall.</p>
Full article ">Figure 7
<p>Risk of rockfalls on buildings, in USD per annum—case of earthquake.</p>
Full article ">Figure 8
<p>Risk of rockfalls on buildings, in USD per annum—case of rainfall.</p>
Full article ">
30 pages, 59008 KiB  
Article
Managing Rockfall Hazard on Strategic Linear Stakes: How Can Machine Learning Help to Better Predict Periods of Increased Rockfall Activity?
by Marie-Aurélie Chanut, Hermann Courteille, Clara Lévy, Abdourrahmane Atto, Lucas Meignan, Emmanuel Trouvé and Muriel Gasc-Barbier
Sustainability 2024, 16(9), 3802; https://doi.org/10.3390/su16093802 - 30 Apr 2024
Viewed by 1435
Abstract
When rockfalls hit and damage linear stakes such as roads or railways, the access to critical infrastructures (hospitals, schools, factories …) might be disturbed or stopped. Rockfall risk management often involves building protective structures that are traditionally based on the intensive use of [...] Read more.
When rockfalls hit and damage linear stakes such as roads or railways, the access to critical infrastructures (hospitals, schools, factories …) might be disturbed or stopped. Rockfall risk management often involves building protective structures that are traditionally based on the intensive use of resources such as steel or concrete. However, these solutions are expensive, considering their construction and maintenance, and it is very difficult to protect long linear stakes. A more sustainable and effective risk management strategy could be to account for changes on rockfall activity related to weather conditions. By integrating sustainability principles, we can implement mitigation measures that are less resource-intensive and more adaptable to environmental changes. For instance, instead of solely relying on physical barriers, solutions could include measures such as restriction of access, monitoring and mobilization of emergency kits containing eco-friendly materials. A critical step in developing such a strategy is accurately predicting periods of increased rockfall activity according to meteorological triggers. In this paper, we test four machine learning models to predict rockfalls on the National Road 1 at La Réunion, a key road for the socio-economic life of the island. Rainfall and rockfall data are used as inputs of the predictive models. We show that a set of features derived from the rainfall and rockfall data can predict rockfall with performances very close and almost slightly better than the standard expert model used for operational management. Metrics describing the performance of these models are translated in operational terms, such as road safety or the duration of road closings and openings, providing actionable insights for sustainable risk management practices. Full article
Show Figures

Figure 1

Figure 1
<p>Two approaches for rockfall risk management: the main chosen approach, using protective works (<b>A</b>) and the proposed approach, using temporary mitigation measures (<b>B</b>).</p>
Full article ">Figure 2
<p>(<b>a</b>) Localization of the RN1 road (red line) along the coast of La Reunion Island, linking the cities of Saint Denis and La Possession. (<b>b</b>) Photograph of the RN1 road after a rockfall that reached the road on the 09/03/2015, breaking steel nets and gabions in its path (image from report BRGM/RP-64556-FR). (<b>c</b>) Zoom in on the RN1 settings, with indications of several landmarks (KP) along the road. Two important ravines reach the road: Jacques, at KP 7.5 km, and Grande Chaloupe, at KP 8.5 km. The yellow rectangles correspond to the figures presenting the detailed geology of the site in <a href="#app1-sustainability-16-03802" class="html-app">Appendix A</a>. (<b>d</b>) Photograph of the cliff around KP 7 km, with the indication of the limits between several geological units and other hydro-geomorphological features. Pictures (<b>e</b>,<b>f</b>) represent outcrops of the superior unit, with variable basalt flow thickness and roughness. Plot (<b>f</b>) also shows an instable rock volume.</p>
Full article ">Figure 3
<p>Rainfall (blue) and rockfall (red star) time series (respectively, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mi>t</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> </mrow> <mrow> <mi>t</mi> </mrow> </msub> </mrow> </semantics></math>) from 2000 to 2018.</p>
Full article ">Figure 4
<p>Cross-correlation between daily rockfalls and rainfall delayed for several delays: (<b>a</b>) for the 2000–2006 period and (<b>b</b>) for the 2009–2018 period.</p>
Full article ">Figure 5
<p>Rainfall distribution on days with rockfalls: (<b>a</b>) from 2000 to 2006, (<b>b</b>) from 2009 to 2018.</p>
Full article ">Figure 6
<p>Rainfall conditions on days with rockfalls depending on the number of events per days, as observed in (<b>a</b>) the 2000–2006 period and (<b>c</b>) the 2009–2018 period. Additionally, they are dependent on the rockfall mass in (<b>b</b>) the 2000–2006 period and (<b>d</b>) the 2009–2018 period.</p>
Full article ">Figure 7
<p>Selected features for ML models.</p>
Full article ">Figure 8
<p>Prediction rate of the DNN over 2000–2006 and bagged tree over 2009–2018 compared with the number of rockfalls per days (<b>a</b>) and the mass of rockfalls (<b>b</b>).</p>
Full article ">Figure 9
<p>Comparison of the distribution of the number of rockfalls per day over the period 2000–2006 (<b>a</b>) and 2009–2018 (<b>c</b>) and comparison of the distribution of the mass of rockfalls over the period 2000–2006 (<b>b</b>) and 2009–2018 (<b>d</b>).</p>
Full article ">Figure 10
<p>Distribution of the duration of periods with consecutive days with rockfalls (<b>left</b>) and the range between these periods (<b>right</b>) during 2000–2006.</p>
Full article ">Figure A1
<p>Geological map of the area around the RN1 road. For the location of the different sub-figures (<b>a</b>–<b>c</b>), see <a href="#sustainability-16-03802-f002" class="html-fig">Figure 2</a>. After [<a href="#B37-sustainability-16-03802" class="html-bibr">37</a>,<a href="#B38-sustainability-16-03802" class="html-bibr">38</a>,<a href="#B39-sustainability-16-03802" class="html-bibr">39</a>,<a href="#B40-sustainability-16-03802" class="html-bibr">40</a>].</p>
Full article ">Figure A2
<p>(<b>a</b>) Drapery systems set up in 1990 and (<b>b</b>) deflector net set up in 2016 along the national road 1 in La Réunion.</p>
Full article ">Figure A3
<p>Extract from the raw rockfall database (source: RN1 Management Services).</p>
Full article ">Figure A4
<p>Principle of the four tested ML models: (<b>a</b>) nearest neighbor algorithm, (<b>b</b>) bagged decision trees, (<b>c</b>) logistic regression, (<b>d</b>) dense neural network with 2 hidden layers.</p>
Full article ">Figure A4 Cont.
<p>Principle of the four tested ML models: (<b>a</b>) nearest neighbor algorithm, (<b>b</b>) bagged decision trees, (<b>c</b>) logistic regression, (<b>d</b>) dense neural network with 2 hidden layers.</p>
Full article ">Figure A5
<p>Binary confusion matrix: 1 stands for at least one rockfall, 0 for no rockfall. TN = true negative, FP = false positive.</p>
Full article ">
17 pages, 3609 KiB  
Article
Real-Time Dynamic Intelligent Image Recognition and Tracking System for Rockfall Disasters
by Yu-Wei Lin, Chu-Fu Chiu, Li-Hsien Chen and Chao-Ching Ho
J. Imaging 2024, 10(4), 78; https://doi.org/10.3390/jimaging10040078 - 26 Mar 2024
Viewed by 1831
Abstract
Taiwan, frequently affected by extreme weather causing phenomena such as earthquakes and typhoons, faces a high incidence of rockfall disasters due to its largely mountainous terrain. These disasters have led to numerous casualties, government compensation cases, and significant transportation safety impacts. According to [...] Read more.
Taiwan, frequently affected by extreme weather causing phenomena such as earthquakes and typhoons, faces a high incidence of rockfall disasters due to its largely mountainous terrain. These disasters have led to numerous casualties, government compensation cases, and significant transportation safety impacts. According to the National Science and Technology Center for Disaster Reduction records from 2010 to 2022, 421 out of 866 soil and rock disasters occurred in eastern Taiwan, causing traffic disruptions due to rockfalls. Since traditional sensors of disaster detectors only record changes after a rockfall, there is no system in place to detect rockfalls as they occur. To combat this, a rockfall detection and tracking system using deep learning and image processing technology was developed. This system includes a real-time image tracking and recognition system that integrates YOLO and image processing technology. It was trained on a self-collected dataset of 2490 high-resolution RGB images. The system’s performance was evaluated on 30 videos featuring various rockfall scenarios. It achieved a mean Average Precision (mAP50) of 0.845 and mAP50-95 of 0.41, with a processing time of 125 ms. Tested on advanced hardware, the system proves effective in quickly tracking and identifying hazardous rockfalls, offering a significant advancement in disaster management and prevention. Full article
(This article belongs to the Section Image and Video Processing)
Show Figures

Figure 1

Figure 1
<p>This figure summarizes the entire optimized motion history analysis process, including image overlay, image differencing, binary thresholding, and the acquisition of motion trajectories from short-term motion variations to long-term historical tracks.</p>
Full article ">Figure 2
<p>Blue rectangular boxes indicating the trajectory of the falling object.</p>
Full article ">Figure 3
<p>The size definition on the image in pixel values.</p>
Full article ">Figure 4
<p>This flowchart illustrates the process starting from the input of the initial image, which undergoes preprocessing using the OpenCV module. It then proceeds to target detection through the YOLO deep learning model, followed by data processing involving coordinate tracking and template matching on the images. This process aims to achieve the final tracking and detection of rockfall, ultimately outputting warnings based on the detection of rockfall.</p>
Full article ">Figure 5
<p>The experimental setup diagram. The structure is divided into three components: the camera, the rock target, and the background.</p>
Full article ">Figure 6
<p>In the actual experiment, one person will release stones in compliance with the field conditions to simulate high-speed free-falling stones dropping from outside the scene into the scene.</p>
Full article ">Figure 7
<p>Results of <span class="html-italic">F</span><sub>1</sub> score with scores of 0.84.</p>
Full article ">Figure 8
<p>Results of <span class="html-italic">mAP</span> at 0.5 with scores of 0.845.</p>
Full article ">Figure 9
<p>Results of precision with values of 0.93.</p>
Full article ">Figure 10
<p>Results of recall with values of 0.91.</p>
Full article ">Figure 11
<p>Results of the Precision–Recall figure.</p>
Full article ">Figure 12
<p>The complete and well-defined trajectory of a falling rock from (<b>a</b>–<b>s</b>) corresponds to frames 82 to 100 in the experimental sequence, which is the main focus of the experiment.</p>
Full article ">Figure 12 Cont.
<p>The complete and well-defined trajectory of a falling rock from (<b>a</b>–<b>s</b>) corresponds to frames 82 to 100 in the experimental sequence, which is the main focus of the experiment.</p>
Full article ">Figure 12 Cont.
<p>The complete and well-defined trajectory of a falling rock from (<b>a</b>–<b>s</b>) corresponds to frames 82 to 100 in the experimental sequence, which is the main focus of the experiment.</p>
Full article ">
3 pages, 458 KiB  
Abstract
An Autonomous Multi-Technological LoRa Sensor Network for Landslide Monitoring
by Mattia Ragnoli, Paolo Esposito, Gianluca Barile, Giuseppe Ferri and Vincenzo Stornelli
Proceedings 2024, 97(1), 11; https://doi.org/10.3390/proceedings2024097011 - 13 Mar 2024
Viewed by 786
Abstract
Hazards like landslides have significant economic and societal repercussions; hence, the issue of remote structure health monitoring has grown in significance for geologic applications. Wireless sensor networks (WSNs) stand out among the new sensing architectures as a particularly well-suited solution, thanks to the [...] Read more.
Hazards like landslides have significant economic and societal repercussions; hence, the issue of remote structure health monitoring has grown in significance for geologic applications. Wireless sensor networks (WSNs) stand out among the new sensing architectures as a particularly well-suited solution, thanks to the versatility they offer. This research, necessary for safety reasons, predictive maintenance and emergency evacuation, presents a WSN-based landslide monitoring system with multi-technology sensor implementation. Its goal is to track the land movements on a hillside. The network is composed of long range (LoRa) sensor nodes connected using a LoRaWAN media access control (MAC) layer. The nodes are several and of different natures and help monitor land movements, hydric parameters and rockfall events, and they also offer a camera view of the landslide in case of an emergency. The system is built on an Internet of Things (IoT) framework, enabling online access to data and reports. The final work will include a system description of the hardware and functionality of all the devices, a description of the web section for remote monitoring, a power analysis and statistics from actual scenarios. Full article
(This article belongs to the Proceedings of XXXV EUROSENSORS Conference)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) System application scheme; (<b>b</b>) an installation point for one of the sensing stations.</p>
Full article ">
21 pages, 2923 KiB  
Article
Dynamic Response of PCCP under the Rockfall Impact Based on the Continuous–Discontinuous Method: A Case Study
by Chunhui Ma, Ying Tu, Yonglin Zhou, Jie Yang and Lin Cheng
Water 2024, 16(6), 801; https://doi.org/10.3390/w16060801 - 7 Mar 2024
Cited by 1 | Viewed by 1338
Abstract
Rockfalls are major geological hazards threatening prestressed concrete cylinder pipes (PCCPs) in water diversion projects. To accurately assess the impact of large deformation movements of rockfalls on PCCPs, this study utilized the continuous–discontinuous method to investigate the dynamic response of a PCCP under [...] Read more.
Rockfalls are major geological hazards threatening prestressed concrete cylinder pipes (PCCPs) in water diversion projects. To accurately assess the impact of large deformation movements of rockfalls on PCCPs, this study utilized the continuous–discontinuous method to investigate the dynamic response of a PCCP under a rockfall. The impact mode of rockfalls, the mechanical characteristics of PCCP, and the nonlinear-contact characteristics between soil and PCCP were considered in this study. The advantages of continuous and discontinuous numerical simulation methods were utilized to establish a continuous and discontinuous coupling model of “tube-soil-rock” considering the interaction of soil and structure. The impact mechanism and process of PCCP under the rockfall were investigated by simulating the rockfall process and analyzing its spatiotemporal evolution. The influence of PCCP under rockfalls with different heights and radii was studied to clarify the effects of these two parameters on the PCCP. Combined with a practical application example of large-scale water transfer projects, there is a tendency of center flattening under static load and dynamic impact load, and the PCCP part directly below the impact point is the most dangerous. This investigation provided a comprehensive understanding of the impact mechanism of the PCCPs under rockfall. The findings of this study have significant implications for the design of the protection engineering of PCCPs and ensuring the safe operation of water diversion projects. Full article
Show Figures

Figure 1

Figure 1
<p>The principle of DEM.</p>
Full article ">Figure 2
<p>Linear contact model.</p>
Full article ">Figure 3
<p>Transmission of continuous–discontinuous coupling boundary force.</p>
Full article ">Figure 4
<p>Continuous–discontinuous coupling boundary.</p>
Full article ">Figure 5
<p>The continuous–discontinuous “Pipe-Soil-Rock” model: (<b>a</b>) the overall model of rockfall’s impact on PCCP, (<b>b</b>) semi-section of the coupling model, and (<b>c</b>) local profile of the coupling model.</p>
Full article ">Figure 6
<p>Construction model of a PCCP based on the FDM.</p>
Full article ">Figure 7
<p>Position coordinate diagram of a rockfall.</p>
Full article ">Figure 8
<p>Ball–pipe contact during an impact process.</p>
Full article ">Figure 9
<p>Displacement time–history curve of concrete and mortar coating under impact loads: (<b>a</b>) concrete–midspan, (<b>b</b>) concrete–spigot joint, (<b>c</b>) concrete–socket joint, (<b>d</b>) coating mortar–midspan, (<b>e</b>) coating mortar–spigot joint, and (<b>f</b>) coating mortar–socket joint.</p>
Full article ">Figure 10
<p>Displacement time–history curve of prestressing wire and steel cylinder under impact loads: (<b>a</b>) prestressing wire–midspan, (<b>b</b>) prestressing wire–spigot joint, (<b>c</b>) prestressing wire–socket joint, (<b>d</b>) steel cylinder–midspan, (<b>e</b>) steel cylinder–spigot joint, and (<b>f</b>) steel cylinder–socket joint.</p>
Full article ">Figure 11
<p>Cloud diagram of the displacement of each PCCP material: (<b>a</b>) the vertical displacement changes in the prestressing wire before and after impact, (<b>b</b>) the vertical displacement changes in the steel cylinder before and after impact, (<b>c</b>) the horizontal displacement changes in the prestressing wire before and after impact, and (<b>d</b>) the horizontal displacement changes in the steel cylinder before and after impact.</p>
Full article ">Figure 12
<p>Displacement time history of each component of the PCCP under rockfall impacts with different radii.</p>
Full article ">Figure 13
<p>The displacement–rockfall radius relationship of each material of the PCCP under rockfall impact.</p>
Full article ">Figure 14
<p>Displacement time history of each component of a PCCP under rockfall impacts from different heights.</p>
Full article ">Figure 15
<p>The displacement–rockfall height relationship of each material of a PCCP under rockfall impacts.</p>
Full article ">
Back to TopTop