Next Article in Journal
Modelling of Metaheuristics with Machine Learning-Enabled Cybersecurity in Unmanned Aerial Vehicles
Next Article in Special Issue
Slope Stability Analysis for a Large Hydropower Station in China
Previous Article in Journal
Effects of Hydrostatic Dissolution and Seepage on the Transport and Mechanical Properties of Glauberite
Previous Article in Special Issue
Study on the Microstructure Evolution and Strength Damage Mechanism of Dolomite under Dissolution Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogeotechnical Predictive Approach for Rockfall Mountain Hazard Using Elastic Modulus and Peak Shear Stress at Soil–Rock Interface in Dry and Wet Phases at KKH Pakistan

1
Department of Civil Engineering, University of Engineering and Technology, Lahore 54890, Pakistan
2
Department of Geological Engineering, University of Engineering and Technology, Lahore 54890, Pakistan
3
Department of Structural Engineering, Future University in Egypt, New Cairo 11865, Egypt
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(24), 16740; https://doi.org/10.3390/su142416740
Submission received: 2 October 2022 / Revised: 27 October 2022 / Accepted: 9 November 2022 / Published: 14 December 2022
(This article belongs to the Special Issue Slope Stability Analysis and Landslide Disaster Prevention)

Abstract

:
Predicting the susceptibility of rockfall mountain hazards for block-in-matrix soils is challenging for critical steep cuts. This research illustrates a hydrogeotechnical approach for the prediction of rockfall triggering by performing laboratory tests on low-cohesive-matrix soil collected from steep slopes with 85° to 88° angles at the Tatta Pani site, Karakorum Highway (KKH), and then real-scale moisture-induced rockfall was conducted on site for the validation of laboratory data. Laboratory data of forty quick direct shear tests on samples collected from the field depicted a 3-fold drop in peak shear stress (PS) at the soil–soil interface and a 9.3-fold drop at the soil–rock interface by varying the moisture content from 1% (taken as dry phase) to a critical laboratory moisture content (MC)LC of 21% (taken as wet phase). Similarly, a drop in the elastic modulus (ES) was observed to be 5.7-fold at the soil–soil interface and 10-fold at the soil–rock interface for a variation of moisture content from 1 % to 21% for the matrix with a permeability (k) range of 3 × 10−4 to 5.6 × 10−4 m/s, which depicts the criticality of moisture content for the rockfall phenomenon. The critical moisture content evaluated in laboratory is validated by an innovative field-inundation method for thirty-two moisture-induced real-scale forced rockfall cases, which showed the rock-block triggering at field dry density (γd)f and the critical field moisture content (MC)FC of the matrix ranging from 1.78 g/cm3 to 1.92 g/cm3, and 1.3% to 25.4%, respectively. Hydrogeotechnical relations, i.e., MC versus PS and ES, at the soil–rock interface are developed for the prediction of rockfall triggering. The proposed correlations may be helpful in the prediction of rockfall hazards by using expected rainfall in the field for disaster warning and landslide disaster prevention at ecological geotechnical engineering projects. The results revealed that the critical (MC)FC and (MC)LC are within 20%, depicting a good confidence level of the outcomes of this research.

1. Introduction

The rock blocks embedded in weak soil matrix form rock blocks and soil mixtures called block-in-matrix soils [1,2]. Rockfall presents the natural event of fast down-slope freefall or the rolling movement of blocks detached from soil–rock mixtures, causing hazards to transportation corridors [3,4]. The rockfall phenomenon is attracting researchers due to its resultant sudden and frequent damage to infrastructures. Rockfall is an unstable process and severe hazard that is caused mostly by rainwater in mountainous regions, presenting an economic and social issue [5]. Soil–rock inhomogeneous mixture is frequently encountered in mountainous ranges showing weak strata prone to slope failures due to the moisture variation in strata [6,7]. The rock blocks and soil particles show different physical and strength properties in response to the loading and moisture variation in strata where the rock blocks have very high strength as compared to the soil matrix [8].
The flow of matrix material can be analyzed by topographical, hydrological, and geotechnical parameters, and it can be observed in extreme topography transition zones and regions with significant elevation differences [9]. The eco-safety of pavement infrastructures and human life in critical rockfall areas having typical extreme topography transition zones with enormous elevation differences necessitates the evaluation of the moisture strength criterion in rockfall areas.
The hydrological effect is a significant cause of the fall of boulders embedded in the matrix and the ultimate phenomenon of slope failure [10]. The flow of rock blocks and fragments, along with soil debris due to high moisture infiltration, is a complex mechanism that is still poorly understood [11]. Rockfalls are slope failures where rock and/or soil rapidly fall down, causing severe hazards to people and transport after the detachment of embedded rock blocks from the steep rock face [12]. Rockfall is the result of changes in the environmental factors along the steep cut slope portions of transport corridors, and it is considered as a significant hazard that frequently occurs on steep slopes, mainly due to rainfall [13]. Rockfall is dangerous for infrastructures and ecology along corridors and at the path of the trajectory of rock-block movement, which necessitates landslide management [10]. The initiation of rapid rockfall and rock fragments causes severe damage to structures and lifelines, where geological and geotechnical data become essential for the significant prediction of the failure of the matrix and the initiation of rock-block movement [3].
Water flow channels affect the saturation phase of matrix material, which ultimately initiates the flow of the matrix, along with the detachment of rock blocks and rock fragments embedded in the matrix [14].
The selection of predictor factors is the most significant part of the matrix failure and rockfall [15]. Rockfall is initiated by earthquakes, animals, excavations, rock slope cutting, and most frequently by the interaction of fine materials of matrix soil with high-intensity rainwater. Rainfall is the key parameter that affects the strength of the soil–rock mixture and is considered as a critical cause for disaster in landslides [16]. Rainfall is strongly related to landslides in steep-slope areas due to critical level of moisture content absorbed by soil strata, i.e., the soil matrix [17,18,19]. The threshold of critical moisture content is the most considerable aspect for the triggering of rockfall and rock-block movement from the saturated condition of the matrix [20,21]. The most critical factor that has been frequently observed worldwide in the occurrence of rockfall and the triggering of rock-block movement is rainfall, as it increases the moisture content and decreases the shear strength of the soil matrix in soil–rock mixtures to a critical level, causing the initiation of rock-block movement [22,23,24]. The pore water pressure of the soil of the slope increases due to intense rainfall, which decreases the shear strength of the soil, resulting in the detachment of rock blocks from the matrix [25,26,27].
The penetration of water in the matrix depends upon its permeability, which plays a key role for the softening of matrix material during and after rainfall in the area [28,29]. The permeability of soil–rock mixtures is around 3 × 10−6 m/s with 20% small rock blocks [16], which appreciably varies with 50% small rock blocks [30] and may reach to very low levels, i.e., 1.5 × 10−9 m/s [31]. The coefficient of permeability (k) varies with grain sizes of the soil matrix, as the permeability (k) value increases with the increase in grain size and porosity [32]; hence, the grain size and shape of the particles significantly affect the permeability of strata [16].
Soil porosity, strength, and compactness are significant factors for rockfall in block-in-matrix soil. The higher porosity favors the rapid rockfall phenomenon during rainfall. The higher level of strength and compactness of the soil strata depict the stable rock block-in-matrix soil. Hence, the evaluation of soil porosity, strength, and compactness of block-in-matrix soil are significant factors for the contact force between soil and rock and the triggering of rockfall.
The seepage of water at the soil–block interface is more critical as compared with seepage in the soil matrix during intensive rainfall [5]. The rainfall causes the softening of geomaterials of the slope at the contact surface of the blocks and initiates the gap between block and geomaterials [27]. The adhesive forces at the soil–rock interface are weakened by the strong seepage forces of water, resulting in a decrease in the friction angle between the block surface and soil [16,27,33]; hence, the block is detached from the matrix and falls freely towards the toe of the cut slope.
Rainfall intensity can be evaluated by statistical data analysis for the considered area on the basis of rainfall history for the determination of rainfall-induced slope failures and rock-block fall [5,10,34].
Water flow patterns in the area are dependent on topography of the understudied hilly area. Topography is the most important feature to be considered in the rockfall and rock-block movement on the critical slopes with debris material [9]. Steep cuts in the rock slope of a transportation corridor are considered most critical for the fall of rock blocks embedded in the soil matrix. If the steep cuts are located in the area of high-intensity rainfall, the rockfall hazard will be at a critical level, as the moisture penetration will trigger the rock-block movement downward to create a hazard [35]. The distance traveled by rock blocks detached from embedded positions can be measured by in situ surveys by identifying marks of the impacts on the toe of the slopes [36].
The permeability of the matrix material in the laboratory is evaluated for the different dry densities and moisture content. The assessment of critical moisture for the specific range of particle shapes can help prevent the failure of the steep cut face of the rock slope. Different moisture content levels result in a wide range of stress–strain curves of the matrix material collected from the field. The evaluation of critical (MC)C is the most desired parameter for the phenomenon of the triggering of rock-block movement. The shape and gradation of the matrix particles affect the water flow in the soil grains and the shear strength of soil strata [37,38]. The flowing material comprises a mixture of water and sediments of different sizes and shapes [39]. The shape of the particles of debris are mostly extremely irregular, varying from clay to boulder size [40]. The grain size of the matrix material is the crucial aspect for the movement of the rock block embedded in the soil matrix.
The failure of the matrix material in its undrained condition triggers the movement of rock blocks [36]. Matrix flow susceptibility can be assessed effectively by the use of geotechnical approaches, whereas the higher water-retention capacity of soil causes a reduction in the shear strength of the matrix material [36]. Landslides mobilize due to the increase in pore water pressure caused by rainfall in the debris flow susceptibility areas [41].
The elastic modulus (ES) of soil at the soil–rock interface, i.e., the interface of the rock block and soil matrix, plays a vital role in the triggering of rock-block movement. The ES of the soil controls the initiation of the movement of the rock block. The ES parameter depicts the elastic or recoverable strain, whereas the peak shear stress (PS) presents the failure strength. The triggering of the rock block is strongly dependent on the ES, which shows the stress mobilized before the triggering phenomena.
A peak shear strength test on fine material can be effectively performed in laboratory by a direct shear test. Stress–strain behaviors of matrix soil for PS are necessary for its characterization to assess the resistive forces to the movement of rock blocks and the soil matrix [42,43]. The PS of the matrix is the critical factor for the initiation of boulder movement, and causes the ultimate hazard for vehicles at the transport corridor. Rock-block falling is a common phenomenon occurring in cut slopes of mountainous areas where rainfall plays a vital role in the triggering of rockfall, depending upon the amount of rainfall [5]. The inundation of matrix material in lab is performed to simulate the field flooding phenomena for the evaluation of the strength of matrix soils [43]. The soil was compacted at the level of in situ/field dry density and flooded with free water. The moisture governs the residual strength of soil having low cohesion [29]. An important study for interface friction and the loss of shear strength upon wetting of the interface plane has been conducted [44] on the soil–cement mortar interface, whereas sudden failure can occur at the interface due to the sudden loss of friction [45].
The flow of matrix material with water comprises three phases, i.e., triggering, movement, and accumulation of the soil matrix with rock blocks and fragments, whereas all these phases are difficult to observe in the field [11]. Real rockfall is unlikely to be witnessed; hence, the field studies elaborating the moisture-induced forced rockfall can be carried out to simulate real rockfall events. This simulation phenomenon is very costly and risky in the context of expected hazards in the simulation process. Real-scale rockfall experimentation has been conducted by many researchers in past [46,47,48].
Past studies have addressed rainfall-induced rockfall; however, the hydrogeotechnical, i.e., moisture-strength factor for rockfall hazards has rarely been studied [49,50]. Previous studies on KKH were conducted [51,52,53] for factors such as aspect, faults, elevation, curvature, lithology, slope, topography, rainfall intensity, drainage, and geology. There is an absolute need to formulate preventive measures in most critical rockfall hazard areas, i.e., the Tatta Pani section of KKH, to correlate the moisture content to the elastic modulus (ES) and peak shear stress (PS) of matrix soil to predict the critical moisture content (MC)c causing the triggering of rockfall. In this study, an innovative, economical, and significant predictive approach involving the hydrogeotechnical factor is presented for real-scale moisture-induced forced rockfalls to manage the resources for preventive measures for the reduction in fatalities. Few studies have been seen for the novel hydrogeotechnical initiator factor of matrix softening, i.e., a substantial reduction in Es and Ps of low-cohesive-matrix soil.
The study intends to meet the following objectives:
  • Collection of rock-block samples and soil matrix from a rockfall hazard area of the Tatta Pani site, KKH, and the determination of field dry density (γd)F by a large-diameter core cutter inserted in the matrix of a steep cut slope face.
  • Laboratory evaluation of the elastic modulus (ES) and peak shear stress (PS) of the soil matrix at the in situ density condition and variable moisture contents ranging from dry to wet phases, and correlation of the moisture content of the soil matrix to the Es and Ps in laboratory for a specific range of grain sizes and constant roughness of rock, i.e., 0.5 mm for the performance of quick direct shear tests.
  • Validation of critical values of laboratory-based moisture contents, i.e., (MC)LC by performing moisture-induced real-scale forced rockfall triggering in the field by the inundation of matrix soil and determining the field critical moisture content (MC)FC.
  • Statistical analysis of the dataset.
This study was focused on the Karakoram Highway (KKH), which is a part of the China–Pakistan Economic Corridor (CPEC). The main studied site was Tatta Pani, which is among the sections highly prone to rockfall hazards. Figure 1 shows embedded rock blocks in matrix soil at the studied site of Tatta Pani. The rockfall hazard value for the study section is around 757 when rated using Pierson’s Rockfall Hazard Rating System (RHRS). The average annual rainfall for Tatta Pani is 100 mm and the mean annual temperature is 16.2 °C. The average height of the cut slope face ranges from 200 ft to 225 ft from the road level. The slope angle ranges from 85° to 88°. The study section receives the most rainfall in March and July. For the last two years, severe rockfall at the study site has been observed, which has caused huge traffic blockage, even for up to one to two weeks.

2. Materials and Methods

Collection of bulk samples of matrix materials and rock blocks was performed at Tatta Pani. The matrix material collected from the steep cut slope along the road comprised mixture of clay, silt, sand, and rock blocks. The size of the matrix material ranging from 6.3 mm to 0.002 mm with embedded rock blocks was studied at the site. Disturbed matrix soil was packed in plastic bags as bulk samples and different rock blocks were collected for the cutting of rock pieces to be used in quick direct shear tests. Samples of matrix soil for the determination of field dry density and field moisture were collected separately by use of large-diameter specially fabricated core cutter.
Field and laboratory studies were performed for evaluation of matrix material. Methodology of this research comprised a field study, including (i) collection of matrix material from study area and (ii) real-scale moisture-induced rockfall at study area; and a laboratory study, comprising (i) classification of matrix material by grain size and (ii) quick direct shear tests at dry and wet conditions for (a) soil–soil and (b) soil–rock interface.
Sampling of matrix and rock blocks from the steep cut slopes is extremely difficult in mountainous areas due to danger of triggering of rock blocks by disturbance created by the sampling method adopted in field. Undisturbed samples of matrix were obtained from study area, and density was determined for use in laboratory testing. Field density tests were performed on the face of the cut slopes by specially manufactured sharp-edged cutting steel molds of 3.5″ and 5″ diameter with 4″ height each, as shown in Figure 2, with area ratio of 8.3%. Two holes of 3 mm diameter were made in the bottom of each mold wall showing the air exit path when air is compressed in the bottom of mold during hammering and insertion of mold into matrix soil. This approach is in line with the core-cutter method outlined in ASTM D2937.
The steep face of the cut slope was prepared to insert the mold in horizontal direction. Insertion of molds into the matrix required hammering to fully fill the mold with matrix soil. This technique was observed to be an effective approach for field density determination on critical steep cuts. The field density was simulated in the laboratory for the assessment of peak shear stress (Ps) and elastic modulus (ES) of the matrix material.
The study of moisture-induced rock-block triggering was performed on the blocks with volumes ranging from 0.26 ft3 to 0.31 ft3 and weights ranging from 15.7 kg to 20.3 kg, with embedment volumes of 8% to 12.3%. The embedment of rock blocks was measured by marking the boundary of matrix on the rock blocks at the embedded position of rock blocks in field.
Methods of geotechnical evaluation of matrix material included grain size analysis, Atterberg limit test, moisture content test, the direct shear (DS) test, and unconsolidated undrained (UU) triaxial compression tests for evaluation of peak shear strength at different moisture contents in unsaturated and saturated phases.
Three grain size analysis tests of the representative matrix material show grain size range from 6.3 mm to 0.002 mm. Three Atterberg limit tests on the matrix material were conducted as per ASTM D4318. Consistency parameters such as liquid limit (LL) and plasticity index (PI) were observed in the range of 27% to 28% and 7% to 9%, respectively, for the material passing #40 sieve. The classification of the material as per ASTM D2487 revealed that most of the material falls in the category of low plastic clayey sand.
Moisture content of the disturbed samples was determined by the method outlined in ASTM D2216. The material collected in the plastic jars from the field was tested for moisture content to evaluate the field moisture content of the samples recovered from the steep cut slope before and after rock-block movement. Moreover, moisture content tests were also performed for the remolded samples of matrix material for direct shear tests at different moisture content levels in saturated and unsaturated phases, i.e., wetting–drying conditions.
The quantification of resistive forces is performed by conducting direct shear tests on soil–soil and soil–rock interface. The resistive forces between the soil and rock surfaces are substantially decreased by the quick wetting process. The wetting of cohesive-matrix soil causes the detachment of rock blocks from the embedded position in matrix.
Quick direct shear tests in direct shear apparatus were performed on soil–soil interface and soil–rock interface. The samples were tested on different moisture content levels, including unsaturated and saturated phases of the matrix. The strain rate of 0.2 mm/min was selected for the shearing phase of the direct shear (DS) tests to simulate the sudden undrained shearing of the matrix in rainy season. The peak shear stress (PS) and elastic modulus (ES) of the matrix material is evaluated from stress–strain curve of shearing phase of DS tests. The granite pieces were cut in the size of (6 cm × 6 cm) with 0.5 mm roughness from the rock blocks collected from site to study the soil–rock interface resistive strength, as shown in Figure 3.
The seepage at the soil–rock contact surface is compared with the seepage through the soil particles of the matrix. For this purpose, direct shear tests were conducted for the predetermined slippage plane, i.e., the contact surface of soil–rock. The undulations were taken as rough surface of the block in field. The upper portion of the direct shear samples was taken as soil matrix. The soil matrix was taken as the representative of the thirty-two samples collected from the cut slope face selected in field for mixing of different moisture contents, i.e., 1% (taken as dry phase, simulating the field dry samples), 9%, 15%, and 21% (taken as wet phase). Four different moisture contents, i.e., 1%, 9%, 15%, and 21%, were selected in study for evaluation of PS and Es for a wide range of moistures, starting from dry phase to wet phase of field material. Initially, the remolded samples of field materials, i.e., clayey sand, were prepared, and degrees of saturation (S) were measured for each moisture content level. The moisture content level of 21% was achieved for 100% saturation level for this soil, and this moisture content was considered as the “wet phase moisture” for the experimental soil in this study. It was observed that moisture content larger than 21% revealed difficulty of sample preparation due to start of soil flow. Forty direct shear tests (twenty for soil–soil interface and twenty for soil–rock interface) were conducted, i.e., five tests each for 1%, 9%, 15%, and 21% moisture levels.
Substantial reduction in strength of soil at the soil–rock interface causes the detachment of rock blocks from the imbedded position. The soil–rock interface friction was studied in laboratory, as shown in Figure 4, by quick direct shear tests to simulate the quick undrained failure phenomena at soil–rock interface. The moisture content data show the critical moisture level for the substantial collapse of PS and Es to be used for assessment of initiation of rockfall triggering in field.
The permeability tests on the matrix material were performed to assess the infiltration of water in soil matrix collected from the steep cut slopes of Tatta Pani. Three permeability tests were conducted on the representative matrix soil as per ASTM D2434. The coefficient of permeability (k) of the matrix material was compared with literature to assess the type of soil regarding the range of permeability.

3. Results

Presented in Figure 5 is the average grain size analysis curve for the matrix material for the study area. The curve shows that the grain size of the Tatta Pani site ranges from 6 mm to 0.075 mm and 87 % sand, 5% silt, and 8% clay. The gradation curve in Figure 5 shows a large portion of sand (coarser portion of filed soil), which depicts high porosity and permeability. The sand fraction of the field soil governs the flow of water at the soil–rock interface [32]. The porosity of the soil was measured as per the method described by [32] to access the absorption of rainfall during the experimentation of real-scale rockfall. The plasticity index (PI) ranges from 5 to 7, depicting the classification of matrix material as per ASTM D2487 to be SC (clayey sand).
Typical stress–strain curves are presented in Figure 6 for the soil matrix. The peak shear stress (PS) and elastic modulus (ES) were evaluated by performing 40 direct shear tests at different moisture content levels of matrix soil. The PS and ES values were observed to be in the range of 88.83 kPa to 29.5 kPa and 11.83 MPa to 2.1 MPa, respectively, from the soil–soil interface resistance tests for the variation of moisture content from 1% to 21%, whereas the PS and ES were observed to be in the range of 67.2 kPa to 7.2 kPa and 8.7 MPa to 1.0 MPa, respectively, in the case of soil–rock interface resistance tests, respectively, at the moisture content range of 1% to 21%. The collapse of PS and ES was observed at 21% moisture content, i.e., laboratory-based critical moisture content (MC)LC. The correlations of the PS and ES are high in the dry phase due to the bonding of the clay fraction (cohesive material) with the sand fraction (non-cohesive material). Clay shows high bonding strength in its dry condition and depicts collapsing behavior in its wet condition [54].
The following correlations of moisture content (MC) with the peak shear stress (PS) (Equations (1) and (2)) and elastic modulus (ES) (Equations (3) and (4)) are developed (Figure 7 and Figure 8) at the matrix density ranging from 1.78 g/cm3 to 1.92 g/cm3 (equivalent to field density) for the soil–soil (Equations (1) and (3)) and soil–rock interfaces (Equations (2) and (4)):
P S = 0.1841 ( MC ) 2 + 1.4277 ( MC ) + 84.009
P S = 0.1838 ( MC ) 2 + 1.4443 ( MC ) + 60.851
Es = 0.0221 ( MC ) 2 + 0.0652 ( MC ) + 11.061
Es = 0.0165 ( MC ) 2 + 0.002 ( MC ) + 8.268
The developed equations (correlations) in the present study have been compared with the literature, and a resemblance of the drop in the shear stress and elastic modulus was observed when compared with the study of [54], which presents a substantial drop in the shear stress and elastic modulus from 15% to 23% moisture content, validating the data of present study; however, the difference in values is attributed to the difference in material properties.
The substantial collapse of the PS and ES causes the triggering of rockfall. The critical values of the PS and ES can predict the critical rainfall intensity in the area expected to cause rockfall. Laboratory data on samples collected from the field from forty quick direct shear tests depicted a 3-fold drop in the PS at the soil–soil interface and a 9.3-fold drop at the soil–rock interface by varying the moisture content from 1% (taken as dry phase) to a critical laboratory moisture content (MC)LC of 21% (taken as wet phase). Similarly, a drop in the ES was observed to be 5.7-fold at the soil–soil interface and 10-fold at the soil–rock interface, causing a triggering in rockfall in the clayey sand with blocks, i.e., block-in-matrix soil (BIM soil). Hence, the PS and ES values for the triggering of rockfall have been mentioned in the text for the soil–rock interface, in addition to the single factor, i.e., critical moisture content.
The quick direct shear tests and unconsolidated undrained triaxial tests were performed in laboratory, due to the fact that the rockfall is a rapid process. The collapse occurrence time was observed in a field simulation test, which helped in the setting of the total time of failure in the direct shear tests, i.e., the strain rate assessment for the direct shear tests.
Physical studies of rockfall are significant for real data analysis and the precise prediction of hazards [55]. A novel approach was applied to saturate the surrounding material of the large blocks for intended fall of boulders. The present study is intended to simulate rockfall at Tatta Pani and the laboratory studies of the moisture strength relation of the matrix material of Tatta Pani. To devise the field experimentation of moisture-induced rockfall, the goals of the field study were set. As the porosity of the soil is a very important factor in study of moisture-induced rockfall, the porosity of the field soil material was determined in the laboratory having an average value of 0.31. The water was showered on the periphery of the embedded rock blocks to make the matrix soil in a fully saturated condition. Rock blocks of about 12–16 inches were selected for forced fall by showering water from upside of the rock blocks. It was observed that the water penetrated the soil–rock interface for a duration of 2 to 3 h and forced the triggering of rock-block movement.
Figure 9 and Figure 10 show the phenomena of rockfall from the embedded position in the soil matrix, i.e., the case of block-in-matrix soils. In the field, it can be observed that the rock blocks embedded in the soil matrix are retained in place due to the peak shear strength and non-exceeding critical moisture (MC)C of the matrix material. The rockfall is mainly caused due to the reduction in the shear strength of the matrix material by the increase in moisture content; in particular, when the shear strength of the soil mass drops to a level lower than the sliding force during slope failure, it ultimately results in a disaster with fast soil flow, i.e., debris flow in hilly areas [3,56,57]. Due to rainfall in the mountainous area, the water infiltrates into the matrix material and pore water pressure increases, resulting in a 90% chance of loss of shear strength in the soil–rock mixture (SRM), causing landslide occurrence [13,58].
The failure pattern of the matrix in the field depicts the movement of different grain size particles of the matrix in the field. Water flow patterns in the area depict the availability of the critical quantity of water required for the initiation of boulder movement. The water flow channels and the rainfall intensity of the area affect the erosion of the matrix material of the steep cut slopes in the study areas. The rock-block movement pattern observed in the field shows that the fine material at the surface of the rock-block loses its shear strength and collapses, resulting in the initiation of rock-block movement, and a rockfall hazard takes place in the field.
Thirty-two real-scale moisture-induced rockfalls were performed in the field (after stopping traffic on both sides). The simulated watering in the field was applied near the upper line of the soil–rock contact to assure the total water absorption at the soil–rock contact surface. The rate of water falling was related to the rainfall intensity for rockfall expected in the area. The critical moisture contents that caused rockfall in Figure 9 and Figure 10 were observed to be 21.5% and 22.7%, respectively, whereas the corresponding critical rainfall intensities were found to be 99 mm and 102 mm, respectively. Soil moisture sensors were used to monitor the field moisture content during the moisture-induced forced rockfall at the site. Furthermore, the samples of soil in the dry and wet phase (at the rockfall stage) were checked for moisture by use of a speedy moisture tester for validation in the field. In addition, the soil samples taken from the soil contact surface after rockfall were packed in airtight plastic containers for moisture content tests in laboratory. The critical moisture content range of the soil matrix collected from the field, i.e., (MC)FC at the rockfall phase, were observed to be 20% to 23%, which is very near to the critical moisture content (MC)LC determined in laboratory, i.e., (MC)FC is within 20 % of (MC)LC. Hence, it is inferred that direct shear tests conducted in laboratory for the soil–rock interface resistance can be effectively used for the prediction of the collapse of shear strength and the elastic modulus in the case of rockfall for block in-matrix soils. It was observed in the field that the wetting of the contact surface caused the triggering of rockfall instead of the additional wetting of the surrounding area of the rock block. Hence, the moisture content of the soil in field was determined for the central portion of the soil patch, i.e., the slippage surface for the rock block, as shown in Figure 9f and Figure 10f.
Analysis of variance (ANOVA) of datasets of the field and laboratory results was conducted, and p-values showed that data were found in good health, as no outlier was found in the data, including laboratory and field experimentation.
The field research in rockfall areas is always risky and costly. This study covered only the extremely critical steep cut slopes to study the fall of boulders or rock blocks. Almost the same sizes of blocks ranging from 0.26 ft3 to 0.31 ft3 and weights ranging from 15.7 kg to 20.3 kg were selected to analyze the variation of PS and ES with the increase in moisture content of the soil–rock interface in the field. The same phases of unsaturated and saturated matrix were used in laboratory. However, the large-volume blocks were not covered in this study due to expected severe hazard creation on the adjacent transportation corridor; however, these aspects should be covered in future research. It is also noteworthy that the steep slope at the Tatta Pani site with a different matrix soil mixture should also be studied in future to cover a full range of matrix materials for a comprehensive rockfall hazard protection model.
Although the slope angle, weight of rock block, size of rock block, and matrix density at the soil–rock interface play a role in the quantification of rockfall, the critical moisture content is observed to be a significant factor in the triggering of rockfall. In this study, the real-scale moisture-induced forced rockfall was triggered by inducing moisture at the soil–rock interface by selecting a very narrow range of slope angles, the weight of the rock block, the size of the rock block, and the matrix density at the soil–rock interface to critically observe the effect of single parameter, i.e., the critical moisture content (MC)C of the matrix at the soil–rock interface. Although the developed correlations are parameter-specific and site-specific, the proposed innovative methodology of field inundation in comparison with laboratory data can be used to develop predictive correlations to prevent/reduce fatalities near vertical rock cut slopes in KKH.
A thirty-liter-capacity container was arranged for the showering test. Showering was conducted by a small pot of known volume. The water showering was applied on a 2 ft × 2 ft (4 sft) area, including the face of the rock block. Showering was continued until the fall of the rock block. The net volume of water was noted as the critical moisture causing rockfall. The total volume of the water used in showering for each of the 32 rockfall-triggering cases was converted into millimeters of rainfall, having a range from 97 mm–107 mm of water (simulated rainfall amount is shown in Figure 9 and Figure 10), which falls in the annual precipitation of the Tatta Pani site, KKH, i.e., 112 mm. Keeping in view the historical event-based precipitation in the study area, the simulated rainfall can be considered too high. This hydrogeotechnical study indicates that the critical moisture (MC)C (converted in mm of rainfall) can well predict the anticipated triggering of rockfall, which can be assessed by the rainfall prediction of the site so that the time/day for stoppage of traffic may be decided to prevent fatalities.

4. Conclusions

A hydrogeotechnical study for the prediction of rockfall triggering was conducted at a critical cut slope face of the rockfall hazard area of the Tatta Pani site, KKH, with cut slope angles ranging from 85° to 88°, rock-block volumes ranging from 0.26 ft3 to 0.31 ft3, and weights ranging from 15.7 kg to 20.3 kg for moisture-induced forced rockfall. The following conclusions can be drawn from this research:
  • Quick direct shear laboratory data showed a 3-fold collapse in peak shear stress (PS) at the soil–soil interface and a 9.3-fold drop at the soil–rock interface by varying the moisture content from 1% to a critical laboratory moisture content (MC)LC of 21%. Similarly, there was a 5.7-fold drop in the elastic modulus (ES) at the soil–soil interface and a 10-fold decrease at the soil–rock interface for the variation of moisture content from 1 % to 21% for the low-cohesive-soil matrix with a permeability range of 3 × 10−4 to 5.6 × 10−4 m/s. Hence, the (MC)LC of 21% depicted the substantial loss of PS and Es at the soil–rock interface.
  • Hydrogeotechnical relations, i.e., moisture content (MC) versus the peak shear stress (PS) and elastic modulus (ES), at the soil–rock interface are developed for the prediction of rock by using MC as the critical factor, which causes the collapse of PS and ES in the saturated phase of soil at the soil–rock interface.
  • A proposed laboratory approach for the collapse of PS and ES is validated by an innovative field inundation of the soil–rock interface to cause moisture-induced real-scale rockfall triggering at field dry density (γd)f and a critical field moisture content (MC)FC of the matrix ranging from 1.78 g/cm3 to 1.92 g/cm3 and 1.3% to 25.4%.
  • Variation in (MC)FC was within 20% of (MC)LC, which depicts a good confidence level for using proposed correlations in matching the lithology of block-in-matrix soil. It was also observed that the ES and PS values of the soil–rock interface are substantially lower as compared with the soil–soil interface.
  • The outcome of hydrogeotechnical studies for 32 rockfall cases also revealed that the triggering of rockfall may start at rainfall ranges from 97–107 mm, falling in the intense rainfalls of the monsoon and torrential period of the Tatta Pani site, KKH.

Author Contributions

Conceptualization, E.M.; methodology, I.R. and A.T.; formal analysis, E.M.; investigation, E.M.; resources, A.M.E.; data curation, K.F.; writing—original draft, E.M.; writing—review and editing, A.M.E.; visualization, F.A.; supervision, I.R. and K.F.; project administration, I.R.; funding acquisition, A.M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kalender, A.; Sonmez, H.; Medley, E.; Tunusluoglu, C.; Kasapoglu, K.E. An approach to predicting the overall strengths of unwelded bimrocks and bimsoils. Eng. Geol. 2014, 183, 65–79. [Google Scholar] [CrossRef]
  2. Wang, Y.; Li, C.H.; Hu, Y.Z. Use of X-ray computed tomography to investigate the effect of rock blocks on meso-structural changes in soil-rock mixture under triaxial deformation. Constr. Build Mater. 2018, 164, 386–399. [Google Scholar] [CrossRef]
  3. Castelli, F.; Lentini, V.; Venti, A.D. Evaluation of Unsaturated Soil Properties for a Debris- Flow Simulation. Geosciences 2021, 11, 64. [Google Scholar] [CrossRef]
  4. Yang, Y.; Chen, T.; Wu, W.; Zheng, H. Modelling the stability of a soil-rock- mixture slope based on the digital image technology and strength reduction numerical manifold method. Eng. Anal. Bound. Elem. 2021, 126, 45–54. [Google Scholar] [CrossRef]
  5. Melillo, M.; Gariano, S.L.; Peruccacci, S.; Sarro, R.; Mateos, R.M.; Brunetti, M.T. Rainfall and rockfalls in the Canary Islands: Assessing a seasonal link. Nat. Hazards Earth Syst. Sci. 2020, 20, 2307–2317. [Google Scholar] [CrossRef]
  6. Medley, E.; Lindquist, E.S. The engineering significance of the scale-independence of some Franciscan Melanges in California, USA. In Proceedings of the 35th US Rock Mechanics Symposium (USRMS), Reno, NV, USA, 5–7 June 1995. [Google Scholar]
  7. Wang, Y.; Li, X.; Wu, Y.F.; Lin, C.; Zhang, B. Experimental study on meso-damage cracking characteristics of RSA by CT test. Environ. Earth Sci. 2015, 73, 5545–5558. [Google Scholar] [CrossRef]
  8. Wang, Y.; Li, X.; Hu, R.L.; Li, S.D.; Wang, J.Y. Experimental study of the ultrasonic and mechanical properties of SRM under compressive loading. Environ. Earth Sci. 2015, 74, 5023–5037. [Google Scholar] [CrossRef]
  9. Jiang, H.; Zou, Q.; Zhou, B.; Hu, Z.; Li, C.; Yao, S.; Yao, H. Susceptibility Assessment of Debris Flows Coupled with Ecohydrological Activation in the Eastern Qinghai-Tibet Plateau. Remote Sens. 2022, 14, 1444. [Google Scholar] [CrossRef]
  10. Emeka, O.J.; Nahazanan, H.; Kalantar, B.; Khuz e aimah, Z.; Sani, O.S. Evaluation of the effect of hydroseeded vegetation for slope reinforcement. Land 2021, 10, 995. [Google Scholar] [CrossRef]
  11. Xia, C.; Haoyong, T. A Quasi-Single-Phase Model for Debris Flows Incorporating Non-Newtonian Fluid Behavior. Water 2022, 14, 1369. [Google Scholar] [CrossRef]
  12. Highland, L.M.; Bobrowsky, P. The Landslide Handbook: A Guide to Understanding Landslides; US Geological Survey: Reston, VA, USA, 2008; p. 29. [Google Scholar]
  13. Zhao, F.; Meng, X.; Zhang, Y.; Chen, G.; Su, X.; Yue, D. Landslide susceptibility mapping of karakorum highway combined with the application of SBAS- InSAR technology. Sensors 2019, 19, 2685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Haque, U.; Blum, P.; Da Silva, P.F.; Andersen, P.; Pilz, J.; Chalov, S.R.; Malet, J.-P.; Auflič, M.J.; Andres, N.; Poyiadji, E. Fatal landslides in Europe. Landslides 2016, 13, 1545–1554. [Google Scholar] [CrossRef]
  15. Dias, H.; Gramani, M.; Grohmann, C.; Bateira, C.; Vieira, B.C. Statistical-based Shallow Landslide Susceptibility Assessment for a Tropical Environment: A Case Study in the Southeastern Brazilian Coast. Nat. Hazards 2021, 108, 205–223. [Google Scholar] [CrossRef]
  16. Wang, Y.; Li, X.; Zheng, B.; Ma, C.F. An experimental investigation of the flow–stress coupling characteristics of soil–rock mixture under compression. Transp. Porous Media 2016, 112, 429–450. [Google Scholar] [CrossRef]
  17. Zhao, B.; Dai, Q.; Zhuo, L.; Mao, J.; Zhu, S.; Han, D. Accounting for satellite rainfall uncertainty in rainfall-triggered landslide forecasting. Geomorphology 2022, 398, 108051. [Google Scholar] [CrossRef]
  18. Emberson, R.; Kirschbaum, D.B.; Amatya, P.; Tanyas, H.; Marc, O. Insights from the topographic characteristics of a large global catalog of rainfall-induced landslide event inventories. Nat. Hazards Earth Syst. Sci. 2022, 22, 1129–1149. [Google Scholar] [CrossRef]
  19. Sun, T.; Sun, D.Y.; Wang, X.K.; Ma, Q.; Gourbesville, P.; Nohara, D. Numerical analysis of landslide-generated debris flow on July 3, 2021 in Izu Mountain area, Shizuoka County, Japan. J. Mt. Sci. 2022, 19, 1738–1747. [Google Scholar] [CrossRef]
  20. Segoni, S.; Piciullo, L.; Gariano, S.L. A review of the literature on rainfall thresholds for landslide occurrence. Landslides 2018, 15, 1483–1501. [Google Scholar] [CrossRef]
  21. Emadi-Tafti, M.; Ataie-Ashtiani, B. A modeling platform for landslide stability: A hydrological approach. Water 2019, 11, 2146. [Google Scholar] [CrossRef] [Green Version]
  22. Melillo, M.; Brunetti, M.T.; Peruccacci, S.; Gariano, S.L.; Guzzetti, F. Rainfall thresholds for the possible landslide occurrence in Sicily (southern Italy) based on the automatic reconstruction of rainfall events. Landslides 2016, 13, 165–172. [Google Scholar] [CrossRef]
  23. Collins, B.D.; Stock, G.M. Rockfall triggering by cyclic thermal stressing of exfoliation fractures. Nat. Geosci. 2016, 9, 395–400. [Google Scholar] [CrossRef]
  24. Contino, A.; Bova, P.; Esposito, G.; Giuffré, I.; Monteleone, S. Historical analysis of rainfall-triggered rockfalls: The case study of the disaster of the ancient hydrothermal Sclafani Spa (Madonie Mts, northern-central Sicily, Italy) in 1851. Nat. Hazards Earth Syst. Sci. 2017, 17, 2229–2243. [Google Scholar] [CrossRef] [Green Version]
  25. Wu, W.A.; Yang, Y.T.; Zheng, H. Hydro-mechanical simulation of the semi-saturated porous soil-rock mixtures using the numerical manifold method. Comput. Methods Appl. Mech. 2020, 370, 113238. [Google Scholar] [CrossRef]
  26. Guo, J.; Yi, S.; Yin, Y.; Cui, Y.; Qin, M.; Li, T.; Wang, C. The effect of topography on landslide kinematics: A case study of the Jichang town landslide in Guizhou, China. Landslides 2020, 17, 959–973. [Google Scholar] [CrossRef]
  27. He, K.; Ma, G.; Hu, X.; Liu, B.; Han, M. The July 2, 2017, Lantian landslide in Leibo, China: Mechanisms and mitigation measures. Geomech. Eng. 2022, 28, 283–298. [Google Scholar]
  28. Li, J.; Liu, Z.; Wang, R.; Zhang, X.; Liu, X.; Yao, Z. Analysis of Debris Flow Triggering Conditions for Different Rainfall Patterns Based on Satellite Rainfall Products in Hengduan Mountain Region, China. Remote Sens. 2022, 14, 2731. [Google Scholar] [CrossRef]
  29. Liu, Y.; Qiu, H.; Yang, D.; Liu, Z.; Ma, S.; Pei, Y.; Zhang, J.; Tang, B. Deformation responses of landslides to seasonal rainfall based on InSAR and wavelet analysis. Landslides 2022, 19, 199–210. [Google Scholar] [CrossRef]
  30. Shakoor, A.; Cook, B.D. The effect of stone content, size, and shape on the engineering properties of a compacted silty clay. Bull. Assoc. Eng. Geol. 1990, 27, 245–253. [Google Scholar] [CrossRef]
  31. Proctor, R. Fundamental principles of soil compaction. Eng. News-Rec. 1933, 111. [Google Scholar]
  32. Sun, J.; Huang, Y. Modeling the Simultaneous Effects of Particle Size and Porosity in Simulating Geo-Materials. Materials 2022, 15, 1576. [Google Scholar] [CrossRef]
  33. Xu, Q.; Fan, X.; Dong, X. Characteristics and formation mechanism of a catastrophic rainfall–induced rock avalanche–mud flow in Sichuan, China, 2010. Landslides 2011, 9, 143–154. [Google Scholar] [CrossRef]
  34. Greenwood, J.R. SLIP4EX–A program for routine slope stability analysis to include the effects of vegetation, reinforcement and hydrological changes. Geotech. Geol. Eng. 2006, 24, 449. [Google Scholar] [CrossRef]
  35. Leyva, S.; Cruz-Pérez, N.; Rodríguez-Martín, J.; Miklin, L.; Santamarta, J.C. Rockfall and Rainfall Correlation in the Anaga Nature Reserve in Tenerife (Canary Islands, Spain). Rock Mech. Rock Eng. 2022, 55, 2173–2181. [Google Scholar] [CrossRef]
  36. Budetta, P. Rockfall-induced impact force causing a debris flow on a volcanoclastic soil slope: A case study in southern Italy. Nat. Hazards Earth Syst. Sci. 2010, 10, 1995–2006. [Google Scholar] [CrossRef]
  37. Cagnoli, B.; Romano, G.P. Effect of grain size on mobility of dry granular flows of angular rock fragments: An experimental determination. J. Volcanol. Geotherm. Res. 2010, 193, 18–24. [Google Scholar] [CrossRef]
  38. Ren, Z.; Wang, K.; Yang, K.; Zhou, Z.H.; Tang, Y.J.; Tian, L.; Xu, Z.M. The grain size distribution and composition of the Touzhai rock avalanche deposit in Yunnan, China. Eng. Geol. 2018, 234, 97–111. [Google Scholar] [CrossRef]
  39. Takahashi, T. Debris Flow Mechanics, Prediction and Countermeasures; Taylor & Francis: London, UK, 2007. [Google Scholar]
  40. Hu, K.; Cui, P.; Li, P. Debris flow dynamic models and numerical computation. Chin. J. Nat. 2014, 36, 313–318. [Google Scholar]
  41. Iverson, R.M.; Reid, M.E.; LaHusen, R.G. Debris-flow mobilization from landslides. Annu. Rev. Earth Planet. Sci. 1997, 25, 85–138. [Google Scholar] [CrossRef]
  42. Blatz, J.; Graham, J. A system for controlled suction in triaxial tests. Géotechnique 2000, 50, 465–469. [Google Scholar] [CrossRef]
  43. Armanini, A.; Capart, H.; Fraccarollo, L.; Larcher, M. Rheological stratification in experimental free-surface flows of granular liquid mixtures. J. Fluid Mech. 2005, 532, 269–319. [Google Scholar] [CrossRef] [Green Version]
  44. Li, M.; Li, Y.; Islam, M.R. Effects of water content and interface roughness on the shear strength of silt–cement mortar interface. Soils Found. 2021, 61, 1615–1629. [Google Scholar] [CrossRef]
  45. Li, M.; Li, Y.; Islam, M.R.; Zhang, D. Experimental Study of the Shear Behaviors of the Cement Mortar–Concrete Interface. J. Mater. Civ. Eng. 2022, 34, 4022287. [Google Scholar] [CrossRef]
  46. Ushiro, T.; Shinohara, S.; Tanida, K.; Yagi, N. A Study on the Motion of Rockfalls on Slopes. In Proceedings of the 5th Symposium on Impact Problems in Civil Engineering, Tokyo, Japan, June 2000; pp. 91–96. [Google Scholar]
  47. Huang, R.Q.; Liu, W.; Zhou, J.P.; Pei, X. Rolling tests on movement characteristics of rock blocks. Chin. J. Geotech. Eng. 2007, 29, 1296–1302. [Google Scholar]
  48. Hu, J.; Li, S.; Li, L.; Shi, S.; Zhou, Z.; Liu, H.; He, P. Field, experimental, and numerical investigation of a rockfall above a tunnel portal in southwestern China. Bull. Eng. Geol. Environ. 2017, 77, 1365–1382. [Google Scholar] [CrossRef]
  49. Herrera, C.; Custodio, E. Groundwater low in a relatively old oceanic volcanic island: The Betancuria area, Fuerteventura Island, Canary Islands, Spain. Sci. Total Environ. 2014, 496, 531–550. [Google Scholar] [CrossRef]
  50. Rodríguez-Losada, J.A.; Hernández-Gutiérrez, L.E.; Olalla, C.; Perucho, A.; Serrano, A.; Ef-Darwich, A. Geomechanical parameters of intact rocks and rock masses from the Canary Islands: Implications on their lank stability. J. Volcanol. Geotherm. Res. 2009, 182, 67–75. [Google Scholar] [CrossRef]
  51. Kanwal, S.; Atif, S.; Shafiq, M. GIS based landslide susceptibility mapping of northern areas of Pakistan, a case study of Shigar and Shyok Basins. Geomat. Nat. Hazards Risk 2016, 5705, 1–19. [Google Scholar] [CrossRef] [Green Version]
  52. Khan, H.; Shafique, M.; Khan, M.A.; Bacha, M.A.; Shah, S.U.; Calligaris, C. Landslide susceptibility assessment using Frequency Ratio, a case study of northern Pakistan, Egypt. J. Remote Sens. Space Sci. 2018, 22, 11–24. [Google Scholar]
  53. Ali, S.; Biermanns, P.; Haider, R.; Reicherter, K. Landslide susceptibility mapping by using a geographic information system (GIS) along the China–Pakistan Economic Corridor (Karakoram Highway), Pakistan. Nat. Hazards Earth Syst. Sci. 2019, 19, 999–1022. [Google Scholar] [CrossRef] [Green Version]
  54. Tankiewicz, M.; Strózyk, J.; Zieba, Z. Undrained Shear Strength cu and Undrained Elastic Modulus Eu of Anthropogenic Soils from Laboratory Tests. IOP Conf. Ser. Earth Environ. Sci. 2021, 906, 12118. [Google Scholar] [CrossRef]
  55. Crosta, G.B.; Agliardi, F. A methodology for physically based rockfall hazard assessment. Nat. Hazards Earth Syst. Sci. 2003, 3, 407–422. [Google Scholar] [CrossRef] [Green Version]
  56. Asch, T.; Buma, J.; Beek, L. A View on Some Hydrological Triggering Systems in Landslides. Geomorphology 1999, 30, 25–32. [Google Scholar] [CrossRef]
  57. Liu, B.; Li, S.; Zhang, L. Experimental and Discrete Element Numerical Analysis of Side Slope Instability Induced by Fissure Water Underlying Impervious Bed. Sci. China Technol. Sci. 2005, 48, 65–80. [Google Scholar]
  58. Wienhofer, J.; Lindenmaier, F.; Zehe, E. Challenges in Understanding the Hydrologic Controls on the Mobility of Slow-Moving Landslides. Vadose Zone J. 2011, 10, 496–511. [Google Scholar] [CrossRef]
Figure 1. A view of scattering pattern of rock blocks in matrix at Tatta Pani, KKH.
Figure 1. A view of scattering pattern of rock blocks in matrix at Tatta Pani, KKH.
Sustainability 14 16740 g001
Figure 2. (a) Large-diameter samplers (b) 3.5″ dia. sampling (c) 5″ dia. sampling.
Figure 2. (a) Large-diameter samplers (b) 3.5″ dia. sampling (c) 5″ dia. sampling.
Sustainability 14 16740 g002
Figure 3. The cut pieces of rock blocks for direct shear tests of soil–rock interface resistance.
Figure 3. The cut pieces of rock blocks for direct shear tests of soil–rock interface resistance.
Sustainability 14 16740 g003
Figure 4. Direct shear test for (a) soil–soil interface resistance and (b) soil–rock interface resistance.
Figure 4. Direct shear test for (a) soil–soil interface resistance and (b) soil–rock interface resistance.
Sustainability 14 16740 g004
Figure 5. Grain size curve of matrix material collected from Tatta Pani.
Figure 5. Grain size curve of matrix material collected from Tatta Pani.
Sustainability 14 16740 g005
Figure 6. Typical stress–strain curves of direct shear test in dry and wet phases of soil–rock interface.
Figure 6. Typical stress–strain curves of direct shear test in dry and wet phases of soil–rock interface.
Sustainability 14 16740 g006
Figure 7. PS of soil−soil and soil−rock interface of Tatta Pani site.
Figure 7. PS of soil−soil and soil−rock interface of Tatta Pani site.
Sustainability 14 16740 g007
Figure 8. ES of soil−soil and soil−rock interface at Tatta Pani site.
Figure 8. ES of soil−soil and soil−rock interface at Tatta Pani site.
Sustainability 14 16740 g008
Figure 9. Stages of real-scale rockfall at Tatta Pani section KKH (Location 1) (a) Rock-block (b) Matrix collection (c) inundation started (d) initiation of detachment of rock-block (e) widening of soil/rock-block gap (f) occurrence of rockfall.
Figure 9. Stages of real-scale rockfall at Tatta Pani section KKH (Location 1) (a) Rock-block (b) Matrix collection (c) inundation started (d) initiation of detachment of rock-block (e) widening of soil/rock-block gap (f) occurrence of rockfall.
Sustainability 14 16740 g009
Figure 10. Stages of real-scale rockfall at Tatta Pani section KKH (Location 2) (a) Rock-block (b) inundation started (c) Initiation of rock-block detachment (d) Enlargement of soil/ rock-block gap (e) Widened soil/rock-block gap (f) Occurrence of rockfall.
Figure 10. Stages of real-scale rockfall at Tatta Pani section KKH (Location 2) (a) Rock-block (b) inundation started (c) Initiation of rock-block detachment (d) Enlargement of soil/ rock-block gap (e) Widened soil/rock-block gap (f) Occurrence of rockfall.
Sustainability 14 16740 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mehmood, E.; Rashid, I.; Ahmed, F.; Farooq, K.; Tufail, A.; Ebid, A.M. Hydrogeotechnical Predictive Approach for Rockfall Mountain Hazard Using Elastic Modulus and Peak Shear Stress at Soil–Rock Interface in Dry and Wet Phases at KKH Pakistan. Sustainability 2022, 14, 16740. https://doi.org/10.3390/su142416740

AMA Style

Mehmood E, Rashid I, Ahmed F, Farooq K, Tufail A, Ebid AM. Hydrogeotechnical Predictive Approach for Rockfall Mountain Hazard Using Elastic Modulus and Peak Shear Stress at Soil–Rock Interface in Dry and Wet Phases at KKH Pakistan. Sustainability. 2022; 14(24):16740. https://doi.org/10.3390/su142416740

Chicago/Turabian Style

Mehmood, Ehtesham, Imtiaz Rashid, Farooq Ahmed, Khalid Farooq, Akbar Tufail, and Ahmed M. Ebid. 2022. "Hydrogeotechnical Predictive Approach for Rockfall Mountain Hazard Using Elastic Modulus and Peak Shear Stress at Soil–Rock Interface in Dry and Wet Phases at KKH Pakistan" Sustainability 14, no. 24: 16740. https://doi.org/10.3390/su142416740

APA Style

Mehmood, E., Rashid, I., Ahmed, F., Farooq, K., Tufail, A., & Ebid, A. M. (2022). Hydrogeotechnical Predictive Approach for Rockfall Mountain Hazard Using Elastic Modulus and Peak Shear Stress at Soil–Rock Interface in Dry and Wet Phases at KKH Pakistan. Sustainability, 14(24), 16740. https://doi.org/10.3390/su142416740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop