scholarly journals Investigation of the Influence of Rainfall Runoff on Shallow Landslides in Unsaturated Soil Using a Mathematical Model

Water ◽  
2019 ◽  
Vol 11 (6) ◽  
pp. 1178 ◽  
Author(s):  
Yen-Yu Chiu ◽  
Hung-En Chen ◽  
Keh-Chia Yeh

Infiltration and groundwater have been widely considered as the main factors that cause shallow landslides; however, the effect of runoff has received less attention. In this study, an in-house physical-process-based shallow landslide model is developed to demonstrate the influence of runoff. The model is controlled by coupling the shallow water equation (dynamic) and Richards’ equation. An infinite slope stability analysis is applied to evaluate the possibility of regional landslides. A real, small catchment topography is adopted as a demonstration example. The simulation illustrates the variations of runoff and the factor of safety (FS) during a storm. The results indicate that, after the surface becomes saturated, the FS may keep varying due to the increasing pressure head, which is caused by increasing surface water depth. This phenomenon most likely occurs downstream where the slopes easily accumulate water. The depth of the surface water may also be a factor of slope failure. Therefore, it is essential to increase the accuracy of calculating the runoff depth when assessing regional shallow landslides.

2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Jun Kameda ◽  
Atsushi Okamoto

AbstractDestructive landslides were triggered by the 6.7 Mw Eastern Iburi earthquake that struck southern Hokkaido, Japan, on 6 September 2018. In this study, we carried out 1-D inversion analysis of one of the shallow landslides near the epicenter using a Bing debris-flow model. At this site, the slope failure comprised cover soil with an initial down-slope length of ~ 80 m and a thickness of ~ 7 m on a slope with < 20° dip. The landslide moved southeastward with a run-out distance of ~ 100 m. Inversion analysis of the post-failure deposit geometry was conducted with the Markov Chain Monte Carlo method (MCMC) to optimize the Bingham rheological parameters of the debris. The analysis reproduced several features of the deposit geometry with a yield stress of ~ 1500 Pa and dynamic viscosity of 800–3000 Pa s. The results suggest that the shallow landslide can be approximated by the flow of a viscoplastic fluid with high-mobility debris and a maximum frontal velocity of 6–9 m/s, with a flow duration of 2–4 min.


2021 ◽  
Author(s):  
Enrico D'Addario ◽  
Leonardo Disperati ◽  
Josè Luis Zezerè ◽  
Raquel Melo ◽  
Sergio Cruz Oliveira

&lt;p&gt;Landsliding is a complex phenomenon and its modelling aimed at predicting where the processes are most likely to occur is a tricky issue to be performed. Apart the chosen modelling approach, for both data-driven and physically-based models, paying adequate attention to the predisposing and triggering factors, as well as the input parameters is no less important. Generally, shallow landslides mobilize relatively small volumes of material sliding along a nearly planar rupture surface which is assumed to be roughly parallel to the ground surface. In the literature it is also widely accepted that shallow landslides involve only unconsolidated slope deposits (i.e., the colluvium), then the rupture surface corresponds to the discontinuity between the bedrock and the overlying loose soil. In this work, based on systematic field observations, we highlight that shallow landslides often involve also portions of the sub-surface bedrock showing different levels of weathering and fracturing. Then, we show that the engineering geological properties of slope deposits, as well as those related to the underlying bedrock, must be considered to obtain more reliable shallow landslides susceptibility assessment. As a first task, a multi-temporal shallow landslide inventory was built by photointerpretation of aerial orthoimages. Then, a new fieldwork-based method is proposed and implemented to acquire, process and spatialize the engineering geological properties of both slope deposits and bedrock. To support the regional scale approach, field observations were collected within, in the neighbour and far from the shallow landslide areas. Finally, both physically-based and data-driven methods were implemented to assess and compare shallow landslide susceptibility at regional scale, as well as to analyse the role of spatial distribution of rock mass quality for shallow slope failure development. The results highlight that, according to geology, structural setting and morphometric conditions, bedrock properties spatially change, defining clusters influencing both the distribution and characters of shallow landslides. As a consequence, the physically-based modelling provides better prediction accuracy when two possible rupture surfaces are analysed, the shallower one located at the slope deposit / bedrock discontinuity, and the deeper one located at the bottom of the fractured and weathered bedrock horizon. Even though the physically-based and data-driven models provide similar results in terms of ROC curves, the resulting susceptibility maps highlight quite substantial differences.&lt;/p&gt;


Water ◽  
2020 ◽  
Vol 12 (5) ◽  
pp. 1229
Author(s):  
Hung-En Chen ◽  
Yen-Yu Chiu ◽  
Tung-Lin Tsai ◽  
Jinn-Chuang Yang

To analyze the effect of runoff on shallow landslides, a model coupling one-dimensional rainfall–runoff and two-dimensional infiltration was established to simulate rainfall, infiltration, and runoff processes. Based on Bishop’s limit equilibrium method, the slope failure of a hypothetical footslope was studied. First, conditions with and without inflow were compared. The results reveal a remarkable difference in factors of safety (FS) between the two conditions, suggesting that considering the effect of runoff is crucial for landslide modeling. In terms of a series of tests of the various magnitudes, durations, lag-time, and peak position of the hydrograph, analyses show that larger inflow leads to more accumulated infiltration and triggers landslides earlier. A long-term duration inflow decreases the stability more than short intensive inflow does. With subsequent surface inflow, slope failure may occur after rainfalls stop, owing to the inflow, and the shape of inflow hydrographs could slightly affect the variance in FS. Results also indicate the necessity of considering the surface runoff when using a numerical model to analyze landslide, particularly on a footslope.


2018 ◽  
Vol 18 (7) ◽  
pp. 1919-1935 ◽  
Author(s):  
Teresa Salvatici ◽  
Veronica Tofani ◽  
Guglielmo Rossi ◽  
Michele D'Ambrosio ◽  
Carlo Tacconi Stefanelli ◽  
...  

Abstract. In this work, we apply a physically based model, namely the HIRESSS (HIgh REsolution Slope Stability Simulator) model, to forecast the occurrence of shallow landslides at the regional scale. HIRESSS is a physically based distributed slope stability simulator for analyzing shallow landslide triggering conditions during a rainfall event. The modeling software is made up of two parts: hydrological and geotechnical. The hydrological model is based on an analytical solution from an approximated form of the Richards equation, while the geotechnical stability model is based on an infinite slope model that takes the unsaturated soil condition into account. The test area is a portion of the Aosta Valley region, located in the northwest of the Alpine mountain chain. The geomorphology of the region is characterized by steep slopes with elevations ranging from 400 m a.s.l. on the Dora Baltea River's floodplain to 4810 m a.s.l. at Mont Blanc. In the study area, the mean annual precipitation is about 800–900 mm. These features make the territory very prone to landslides, mainly shallow rapid landslides and rockfalls. In order to apply the model and to increase its reliability, an in-depth study of the geotechnical and hydrological properties of hillslopes controlling shallow landslide formation was conducted. In particular, two campaigns of on site measurements and laboratory experiments were performed using 12 survey points. The data collected contributed to the generation of an input map of parameters for the HIRESSS model. In order to consider the effect of vegetation on slope stability, the soil reinforcement due to the presence of roots was also taken into account; this was done based on vegetation maps and literature values of root cohesion. The model was applied using back analysis for two past events that affected the Aosta Valley region between 2008 and 2009, triggering several fast shallow landslides. The validation of the results, carried out using a database of past landslides, provided good results and a good prediction accuracy for the HIRESSS model from both a temporal and spatial point of view.


2021 ◽  
Vol 13 (12) ◽  
pp. 2385
Author(s):  
Iuliana Armaș ◽  
Mihaela Gheorghe ◽  
George Cătălin Silvaș

A multi-temporal satellite radar interferometry technique is used for deriving the actual surface displacement patterns in a slope environment in Romania, in order to validate and improve a landslide susceptibility map. The probability the occurrence of future events is established using a deterministic approach based on a classical one-dimension infinite slope stability model. The most important geotechnical parameters for slope failure in the proposed study area are cohesion, unit weight and friction angle, and the triggering factor is a rapid rise in groundwater table under wetting conditions. Employing a susceptibility analysis using the physically based model under completely saturated conditions proved to be the most suitable scenario for identifying unstable areas. The kinematic characteristics are assessed by the Small BAseline Subsets (SBAS) interferometry technique applied to C-band synthetic aperture radar (SAR) Sentinel-1 imagery. The analysis was carried out mainly for inhabited areas which present a better backscatter return. The validation revealed that more than 22% of the active landslides identified by InSAR were predicted as unstable areas by the infinite slope model. We propose a refinement of the susceptibility map using the InSAR results for unravelling the danger of the worst-case scenario.


2008 ◽  
Vol 10 (3) ◽  
pp. 227-244 ◽  
Author(s):  
Olaf Kolditz ◽  
Jens-Olaf Delfs ◽  
Claudius Bürger ◽  
Martin Beinhorn ◽  
Chan-Hee Park

In this paper we present an object-oriented concept for numerical simulation of multi-field problems for coupled hydrosystem analysis. Individual (flow) processes modelled by a particular partial differential equation, i.e. overland flow by the shallow water equation, variably saturated flow by the Richards equation and saturated flow by the groundwater flow equation, are identified with their corresponding hydrologic compartments such as land surface, vadose zone and aquifers, respectively. The object-oriented framework of the compartment approach allows an uncomplicated coupling of these existing flow models. After a brief outline of the underlying mathematical models we focus on the numerical modelling and coupling of overland flow, variably saturated and groundwater flows via exchange flux terms. As each process object is associated with its own spatial discretisation mesh, temporal time-stepping scheme and appropriate numerical solution procedure. Flow processes in hydrosystems are coupled via their compartment (or process domain) boundaries without giving up the computational necessities and optimisations for the numerical solution of each individual process. However, the coupling requires a bridging of different temporal and spatial scales, which is solved here by the integration of fluxes (spatially and temporally). In closing we present three application examples: a benchmark test for overland flow on an infiltrating surface and two case studies – at the Borden site in Canada and the Beerze–Reusel drainage basin in the Netherlands.


2018 ◽  
Author(s):  
Anna Botto ◽  
Enrica Belluco ◽  
Matteo Camporese

Abstract. Data assimilation has been recently the focus of much attention for integrated surface-subsurface hydrological models, whereby joint assimilation of water table, soil moisture, and river discharge measurements with the ensemble Kalman filter (EnKF) have been extensively applied. Although the EnKF has been specifically developed to deal with nonlinear models, integrated hydrological models based on the Richards equation still represent a challenge, due to strong nonlinearities that may significantly affect the filter performance. Thus, more studies are needed to investigate the capabilities of the EnKF to correct the system state and identify parameters in cases where the unsaturated zone dynamics are dominant, as well as to quantify possible tradeoffs associated with assimilation of multi-source data. Here, the model CATHY (CATchment HYdrology) is applied to reproduce the hydrological dynamics observed in an experimental two-layered hillslope, equipped with tensiometers, water content reflectometer probes, and tipping bucket flow gages to monitor the hillslope response to a series of artificial rainfall events. Pressure head, soil moisture, and subsurface outflow are assimilated with the EnKF in a number of scenarios and the challenges and issues arising from the assimilation of multi-source data in this real-world test case are discussed. Our results demonstrate that the EnKF is able to effectively correct states and parameters even in a real application characterized by strong nonlinearities. However, multi-source data assimilation may lead to significant trade-offs: the assimilation of additional variables can lead to degradation of model predictions for other variables that were otherwise well reproduced. Furthermore, we show that integrated observations such as outflow discharge cannot compensate for the lack of well-distributed data in heterogeneous hillslopes.


2015 ◽  
Vol 3 (1) ◽  
pp. 791-836 ◽  
Author(s):  
B.-G. Chae ◽  
J.-H. Lee ◽  
H.-J. Park ◽  
J. Choi

Abstract. Most landslides in Korea are classified as shallow landslides with an average depth of less than 2 m. These shallow landslides are associated with the advance of a wetting front in the unsaturated soil due to rainfall infiltration, which results in an increase in water content and a reduction in the matric suction in the soil. Therefore, this study presents a modified equation of infinite slope stability analysis based on the concept of the saturation depth ratio to analyze the slope stability change associated with the rainfall on a slope. A rainfall infiltration test in unsaturated soil was performed using a column to develop an understanding of the effect of the saturation depth ratio following rainfall infiltration. The results indicated that the rainfall infiltration velocity due to the increase in rainfall in the soil layer was faster when the rainfall intensity increased. In addition, the rainfall infiltration velocity tends to decrease with increases in the unit weight of soil. The proposed model was applied to assess its feasibility and to develop a regional landslide susceptibility map using a Geographic Information System (GIS). For that purpose, the spatial databases for input parameters were constructed and landslide locations were obtained. In order to validate the proposed approach, the results of the proposed approach were compared with the landslide inventory using ROC (Receiver Operating Characteristics) graph. In addition, the results of the proposed approach were compared with the previous approach used steady state hydrological model. Consequently, the approach proposed in this study displayed satisfactory performance in classifying landslide susceptibility and showed better performance than the steady state approach.


2013 ◽  
Vol 10 (8) ◽  
pp. 10997-11033 ◽  
Author(s):  
C. Massari ◽  
L. Brocca ◽  
S. Barbetta ◽  
C. Papathanasiou ◽  
M. Mimikou ◽  
...  

Abstract. Floods are one of the most dangerous natural hazards in Mediterranean regions. Flood forecasting tools and early warning systems can be very beneficial to reduce flood risk. Event-based rainfall runoff models are frequently employed for operational flood forecasting purposes because of their simplicity and the reduced number of parameters involved with respect to continuous models. However, the advantages that are related with the reduced parameterization face against the need for a correct initialization of the model, especially in areas affected by strong climate seasonality. On the other hand, the use of continuous models may be very problematic in poorly gauged areas. This paper introduces a simplified continuous rainfall-runoff model, which uses globally available soil moisture retrievals to identify the initial wetness condition of the catchment, and, only event rainfall data to simulate discharge hydrographs. The model calibration involves only 3 parameters. For soil moisture, beside in situ and modelled data, satellite products from the Advanced SCATterometer (ASCAT) and the Advanced Microwave Scanning Radiometer for Earth observation (AMSR-E) sensors are employed. Additionally, the ERA-LAND reanalysis soil moisture product of the European Centre for Medium Range Weather Forecasting (ECMWF) is used. The model was tested in the small catchment of Rafina, 109 km2 located in the Eastern Attica region, Greece. Specifically, fifteen rainfall-runoff events were modelled by considering different configurations for the initial soil moisture conditions. Comparing the performance of the different soil moisture products, it was found that all global indicators allow reproducing fairly well the selected flood events providing much better results than the situation where a constant initial condition is provided. ERA-LAND slightly outperforms the satellite soil moisture products and in general, all the indicators give the same performance obtained by ground and continuously simulated soil moisture data. Due to the wide diffusion of globally available soil moisture retrievals and the small amount of parameters used, the proposed modelling approach is very suitable for runoff prediction in poorly gauged areas.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-25
Author(s):  
Thanh Son Nguyen ◽  
Kuo-Hsin Yang ◽  
Chia-Chun Ho ◽  
Feng-Chi Huang

Although the mechanisms of slope failure caused by rising groundwater have been widely investigated, the kinematic behavior of landslides in the postfailure stage, which contains essential information for hazard mitigation and risk assessment, has not yet been fully studied. Thus, in this study, a series of numerical simulations using the material point method (MPM) were conducted to analyze the kinematic behavior and soil movement of shallow landslides (infinite slope problems). First, the proposed MPM formulation was validated in a full-scale landslide flume test. The simulated results of final slope profile, runout distance, deposit height, shear band development, slope displacement, and velocity accorded with the experimental results, suggesting that the MPM can quantitatively simulate large deformations. A parametric study of shallow slopes with various hydrological conditions and soil hydraulic and soil mechanical parameters was then performed to assess the influence of the aforementioned factors on landslide kinematics. The simulation results indicated that mechanical behavior at the slope toe is complex; the multiple plastic shear bands generated at the slope toe were due to a combination of shearing and compression. The deposition profile of the slopes was significantly influenced by all input parameters. Among the aforementioned parameters, soil cohesion, location of the groundwater table, and saturated soil permeability most greatly affected runout distance in the sensitivity assessment. Soil friction angle had a minor influence on the kinematic behavior of the slope.


Sign in / Sign up

Export Citation Format

Share Document