Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.

Simulation of pesticide behavior in a paddy block by a pesticide fate and transport model

A simulation model (PCPF-B) was developed based on the PCPF-1 model to predict the runoff of pesticides from paddy plots to a drainage canal in a paddy block. The block-scale model now comprises three modules: (1) a module for pesticide application, (2) a module for pesticide behavior in paddy fields, and (3) a module for pesticide concentration in the drainage canal. The PCPF-B model was first evaluated by published data in a single plot and then was applied to predict the concentration of bensulfuron-methyl in one paddy block in the Sakura river basin, Ibaraki, Japan, where a detailed field survey was conducted. The PCPF-B model simulated well the behavior of bensulfuron-methyl in individual paddy plots. It also reflected the runoff pattern of bensulfuron-methyl at the block outlet, although overestimation of bensulfuronmethyl concentrations occurred due to uncertainty in water balance estimation. Application of water management practice such as water-holding period and seepage control also affected the performance of the model. A probabilistic approach may be necessary for a comprehensive risk assessment in large-scale paddy areas.

Quantification of the inevitable: the influence of soil macrofauna on soil water movement in rehabilitated open-cut mined lands

Recolonisation of soil by macrofauna (especially ants, termites and earthworms) in rehabilitated open-cut mine sites is inevitable and, in terms of habitat restoration and function, typically of great value. In these highly disturbed landscapes, soil invertebrates play a major role in soil development (macropore configuration, nutrient cycling, bioturbation, etc.) and can influence hydrological processes such as infiltration, seepage, runoff generation and soil erosion. Understanding and quantifying these ecosystem processes is important in rehabilitation design, establishment and subsequent management to ensure progress to the desired end goal, especially in waste cover systems designed to prevent water reaching and transporting underlying hazardous waste materials. However, the soil macrofauna is typically overlooked during hydrological modelling, possibly due to uncertainties on the extent of their influence, which can lead to failure of waste cover systems or rehabilitation activities. We propose that scientific experiments under controlled conditions and field trials on post-mining lands are required to quantify (i) macrofauna–soil structure interactions, (ii) functional dynamics of macrofauna taxa,and (iii) their effects on macrofauna and soil development over time. Such knowledge would provide crucial information for soil water models, which would increase confidence in mine waste cover design recommendations and eventually lead to higher likelihood of rehabilitation success of open-cut mining land.

Depressed weir with two unequal sheetpiles in anisotropic soil

An analytical solution is presented, making use of the Schwartz-Christoffel transformation, for determining the seepage characteristics for the problem of flow under a weir having two unequal sheetpiles at the ends and embedded in an anisotropic porous medium of finite thickness. Results for several particular cases of simple hydraulic structures can be obtained from the general solution presented. Numerical results in nondimensional form have been given for quantity of seepage and exit gradient distribution for various conditions in the equivalent transformed isotropic section and, by making use of the physical parameters in the actual anisotropic plane and the set of transformation relations given, these quantities (seepage loss, exit gradient) can be interpreted in the actual anisotropic physical plane.

Interface between 2 dissimilar fluids beneath a sheetpile cofferdam

A closed form solution is presented for determining the shape and location of the interface between two dissimilar fluids (having different densities) when steady flow takes place through a homogeneous and isotropic porous medium, into a sheetpile cofferdam; the interface is assumed to be sharp and the lower fluid stationary. The solution is obtained using the inverse hodograph. Numerical results are presented in nondimensional form for various parametric conditions in the physical plane; the interface pattern, as also the seepage discharge and exit gradient distribution are shown. The critical conditions of the interface are studied.

Using a structural approach to identify relationships between soil and erosion in a semi-humid forested area, South India

Biogeochemical and hydrological cycles are currently studied on a small experimental forested watershed (4.5 km(2)) in the semi-humid South India. This paper presents one of the first data referring to the distribution and dynamics of a widespread red soil (Ferralsols and Chromic Luvisols) and black soil (Vertisols and Vertic intergrades) cover, and its possible relationship with the recent development of the erosion process. The soil map was established from the observation of isolated soil profiles and toposequences, and surveys of soil electromagnetic conductivity (EM31, Geonics Ltd), lithology and vegetation. The distribution of the different parts of the soil cover in relation to each other was used to establish the dynamics and chronological order of formation. Results indicate that both topography and lithology (gneiss and amphibolite) have influenced the distribution of the soils. At the downslope, the following parts of the soil covers were distinguished: i) red soil system, ii) black soil system, iii) bleached horizon at the top of the black soil and iv) bleached sandy saprolite at the base of the black soil. The red soil is currently transforming into black soil and the transformation front is moving upslope. In the bottom part of the slope, the chronology appears to be the following: black soil > bleached horizon at the top of the black soil > streambed > bleached horizon below the black soil. It appears that the development of the drainage network is a recent process, which was guided by the presence of thin black soil with a vertic horizon less than 2 in deep. Three distinctive types of erosional landforms have been identified: 1. rotational slips (Type 1); 2. a seepage erosion (Type 2) at the top of the black soil profile; 3. A combination of earthflow and sliding in the non-cohesive saprolite of the gneiss occurs at midslope (Type 3). Types 1 and 2 erosion are mainly occurring downslope and are always located at the intersection between the streambed and the red soil-black soil contact. Neutron probe monitoring, along an area vulnerable to erosion types 1 and 2, indicates that rotational slips are caused by a temporary watertable at the base of the black soil and within the sandy bleached saprolite, which behaves as a plane of weakness. The watertable is induced by the ephemeral watercourse. Erosion type 2 is caused by seepage of a perched watertable, which occurs after swelling and closing of the cracks of the vertic clay horizon and within a light textured and bleached horizon at the top of black soil. Type 3 erosion is not related to the red soil-black soil system but is caused by the seasonal seepage of saturated throughflow in the sandy saprolite of the gneiss occurring at midslope. (c) 2006 Elsevier B.V. All rights reserved.

Analysis of gradually and spatially varied flow in sand-bed channels

Seepage through a sand bed affects the channel hydrodynamics, which in turn alters channel stability. Thus, the effect of seepage on its hydrodynamic parameters needs to be ascertained. The present work analyses spatially varied flow of a sand-bed channel subjected to seepage in the downward direction through a sand bed. Numerically calculated flow profiles affected by seepage have been verified using experimental observations. The present work also analyses the friction slope, velocity and bed shear stress variations along the channel for both seepage and no-seepage conditions. It was found that the downward seepage-induced channel flow has larger friction slope and bed shear stress than that of no-seepage.

Incipient motion design of sand bed channels affected by bed suction

Seepage through sand bed channels in a downward direction (suction) reduces the stability of particles and initiates the sand movement. Incipient motion of sand bed channel with seepage cannot be designed by using the conventional approach. Metamodeling techniques, which employ a non-linear pattern analysis between input and output parameters and solely based on the experimental observations, can be used to model such phenomena. Traditional approach to find non-dimensional parameters has not been used in the present work. Parameters, which can influence the incipient motion with seepage, have been identified and non-dimensionalized in the present work. Non-dimensional stream power concept has been used to describe the process. By using these non-dimensional parameters; present work describes a radial basis function (RBF) metamodel for prediction of incipient motion condition affected by seepage. The coefficient of determination, R-2 of the model is 0.99. Thus, it can be said that model predicts the phenomena very well. With the help of the metamodel, design curves have been presented for designing the sand bed channel when it is affected by seepage. (C) 2010 Elsevier B.V. All rights reserved.

Analysis of flow near a dug well in an unconfined aquifer

A numerical analysis of flow to a dug well in an unconfined aquifer is made, taking into account well storage, elastic storage release, gravity drainage, anisotropy, partial penetration, vertical flow and seepage surface at the well face, and treating the water table in the aquifer and water level in the well as unknown boundaries. The pumped discharge is maintained constant. The solution is obtained by a two-level iterative scheme. The effects of governing parameters on the drawdown, development of seepage surface and contribution from aquifer flow to the total discharge are discussed. The degree of anisotropy and partial penetration are found to be the parameters which affect the flow characteristics most significantly. The effect of anisotropy on the development of seepage surface is very pronounced.

Permeability and Compressibility Behavior of Bentonite-Sand/Soil Mixes

As a seepage barrier slurry trench material should have a relatively low coefficient of permeability, in the range of 10(-7) cm/s, and at the same time should be compatible with surrounding material with regard to compressibility. Although bentonite-sand/soil mixes are used widely, there is no specific engineering approach to proportion these mixes that satisfies the above practical requirements. In this paper, a generalized approach is presented for predicting the permeability and compressibility characteristics of mixes with minimum input parameters. This approach will be helpful in proportioning mixes and predicting corresponding changes in engineering behavior. It is possible to proportion a mix to arrive at the required compressibility without affecting the permeability. This is explained using an illustrative example.

Permeability and consolidation behavior of fly ashes

A knowledge of permeability and consolidation is essential in a number of engineering problems such as settlement, seepage, and stability of the structures. Since fly ash is used very widely for several geotechnical applications, there is a need to understand its permeability and consolidation behavior. This paper presents a detailed study conducted on two Indian fly ashes. It brings out the role of chemical composition (free lime) on the permeability and consolidation behavior of fly ashes. It is found that the permeability values computed based on grain-size distribution agree well with those obtained based on test data.

Optimal crop planning and conjunctive use of water resources in a coastal river basin

Due to increasing trend of intensive rice cultivation in a coastal river basin, crop planning and groundwater management are imperative for the sustainable agriculture. For effective management, two models have been developed viz. groundwater balance model and optimum cropping and groundwater management model to determine optimum cropping pattern and groundwater allocation from private and government tubewells according to different soil types (saline and non-saline), type of agriculture (rainfed and irrigated) and seasons (monsoon and winter). A groundwater balance model has been developed considering mass balance approach. The components of the groundwater balance considered are recharge from rainfall, irrigated rice and non-rice fields, base flow from rivers and seepage flow from surface drains. In the second phase, a linear programming optimization model is developed for optimal cropping and groundwater management for maximizing the economic returns. The models developed were applied to a portion of coastal river basin in Orissa State, India and optimal cropping pattern for various scenarios of river flow and groundwater availability was obtained.

Analytical solutions for protective filters based on soil-retention and permeability criteria with respect to the phenomenon of soil boiling

For the successful performance of a granular filter medium, existing design guidelines, which are based on the particle size distribution (PSD) characteristics of the base soil and filter medium, require two contradictory conditions to be satisfied, viz., soil retention and permeability. In spite of the wider applicability of these guidelines, it is well recognized that (i) they are applicable to a particular range of soils tested in the laboratory, (ii) the design procedures do not include performance-based selection criteria, and (iii) there are no means to establish the sensitivity of the important variables influencing performance. In the present work, analytical solutions are developed to obtain a factor of safety with respect to soil-retention and permeability criteria for a base soil - filter medium system subjected to a soil boiling condition. The proposed analytical solutions take into consideration relevant geotechnical properties such as void ratio, permeability, dry unit weight, effective friction angle, shape and size of soil particles, seepage discharge, and existing hydraulic gradient. The solution is validated through example applications and experimental results, and it is established that it can be used successfully in the selection as well as design of granular filters and can be applied to all types of base soils.

Inverse models to analyze the spatiotemporal variations of chemical weathering fluxes in a granito-gneissic watershed: Mule Hole, South India

Water-rock reactions are driven by the influx of water, which are out of equilibrium with the mineral assemblage in the rock. Here a mass balance approach is adopted to quantify these reactions. Based on field experiments carried out in a granito-gneissic small experimental watershed (SEW), Mule Hole SEW (similar to 4.5 km(2)), quartz, oligoclase, sericite, epidote and chlorite are identified as the basic primary minerals while kaolinite, goethite and smectite are identified as the secondary minerals. Observed groundwater chemistry is used to determine the weathering rates, in terms of `Mass Transfer Coefficients' (MTCs), of both primary and secondary minerals. Weathering rates for primary and secondary minerals are quantified in two steps. In the first step, top red soil is analyzed considering precipitation chemistry as initial phase and water chemistry of seepage flow as final phase. In the second step, minerals present in the saprolite layer are analyzed considering groundwater chemistry as the output phase. Weathering rates thus obtained are converted into weathering fluxes (Q(weathering)) using the recharge quantity. Spatial variability in the mineralogy observed among the thirteen wells of Mule Hole SEW is observed to be reflected in the MTC results and thus in the weathering fluxes. Weathering rates of the minerals in this silicate system varied from few 10 mu mol/L (in case of biotite) to 1000 s of micromoles per liter (calcite). Similarly, fluxes of biotite are observed to be least (7 +/- 5 mol/ha/yr) while those of calcite are highest (1265 791 mol/ha/yr). Further, the fluxes determined annually for all the minerals are observed to be within the bandwidth of the standard deviation of these fluxes. Variations in these annual fluxes are indicating the variations in the precipitation. Hence, the standard deviation indicated the temporal variations in the fluxes, which might be due to the variations in the annual rainfall. Thus, the methodology adopted defines an inverse way of determining weathering fluxes, which mainly contribute to the groundwater concentration. (C) 2011 Elsevier B.V. All rights reserved.

Pseudo-Static and Pseudo-Dynamic Stability Analysis of Tailings Dam Under Seismic Conditions

In this paper the seismic slope stability analyses are performed for a typical section of 44 m high water retention type tailings earthen dam located in the eastern part of India, using both the conventional pseudo-static and recent pseudo-dynamic methods. The tailings earthen dam is analyzed for different upstream conditions of reservoir like filled up with compacted and non-compacted dumped waste materials with different water levels of the pond tailings portion. Phreatic surface is generated using seepage analysis in geotechnical software SEEP/W and that same is used in the pseudo-static and pseudo-dynamic analyses to make the approach more realistic. The minimum values of factor of safety using pseudo-static and pseudo-dynamic method are obtained as 1.18 and 1.09 respectively for the chosen seismic zone in India. These values of factor of safety show clearly the demerits of conventional pseudo-static analysis compared to recent pseudo-dynamic analysis, where in addition to the seismic accelerations, duration, frequency of earthquake, body waves traveling during earthquake and amplification effects are considered.