Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.

An elasto-plastic critical state framework for unsaturated soil.

Data from a series of controlled suction triaxial tests on samples of compacted speswhite kaolin were used in the development of an elasto–plastic critical state framework for unsaturated soil. The framework is defined in terms of four state variables: mean net stress, deviator stress, suction and specific volume. Included within the proposed framework are an isotropic normal compression hyperline, a critical state hyperline and a state boundary hypersurface. For states that lie inside the state boundary hypersurface the soil behaviour is assumed to be elastic, with movement over the state boundary hypersurface corresponding to expansion of a yield surface in stress space. The pattern of swelling and collapse observed during wetting, the elastic–plastic compression behaviour during isotropic loading and the increase of shear strength with suction were all related to the shape of the yield surface and the hardening law defined by the form of the state boundary. By assuming that constant–suction cross–sections of the yield surface were elliptical it was possible to predict test paths for different types of triaxial shear test that showed good agreement with observed behaviour. The development of shear strain was also predicted with reasonable success, by assuming an associated flow rule.

Aspects of Mechanical Behavior and Modeling of a Tropical Unsaturated Soil

The research studies the applicability of two elastoplastic models for the collapse prediction of the lateritic soil profile from Southeastern Brazil. These tropical soils have peculiar geotechnical behavior, due to their mineralogical composition and porous structure coming from intense process of formation. Two elastoplastic models were analyzed: the Barcelona Basic Model (BBM) and another one based on BBM, however developed for tropical soils. Oedometric tests with suction control were performed at three distinct depths of the soil profile. The BBM was not suitable for the upper layer of the soil profile, because BBM considers the compressible behavior of the soil in function of the reduction of the elastoplastic compressibility index with the increase of the matric suction. The model developed for tropical soils showed better suited to the compressible behavior of the soil profile, resulting in good prediction of the collapse potential, mainly by accepting increasing values of the elastoplastic compressibility index of the soil profile with the matric suction rise. © 2013 Springer Science+Business Media Dordrecht.

Modeling field load tests in lateritic unsaturated soil

Many models for unsaturated soil have been developed in the last years, accompanying the development of experimental techniques to deal with such soils. The benchmark of the models for unsaturated soil can be assigned to the Barcelona Basic Model (BBM) now incorporated in some codes such as the CODE_BRIGHT. Most of those models were validated considering limited laboratory test results and not much validation is available considering real field problems. This paper presents modeling results of field plate load tests performed under known suction on a lateritic unsaturated soil. The required input data were taken from laboratory tests performed under suction control. The modeling nicely reproduces field tests allowing appreciating the influence of soil suction on the stress-settlement curve. In addition, wetting induced or collapse settlements were calculated from field tests and were nicely duplicated by the numerical analysis performed.

Small strain elastic behaviour of unsaturated soil investigated by bender/extender element testing

The aim of this project was to investigate very small strain elastic behaviour of soils under unsaturated conditions, using bender/extender element (BEE) testing. The behaviour of soils at very small strains has been widely studied under saturated conditions, whereas much less work has been performed on very small strain behaviour under unsaturated conditions. A suction-controlled double wall triaxial apparatus for unsaturated soil testing was modified to incorporate three pairs of BEEs transmitting both shear and compression waves with vertical and horizontal directions of wave transmission and wave polarisation. Various different techniques for measuring wave travel time were investigated in both the time domain and the frequency domain and it was concluded that, at least for the current experimental testing programme, peak-to-first-peak in the time domain was the most reliable technique for determining wave travel time. An experimental test programme was performed on samples of compacted speswhite kaolin clay. Two different forms of compaction were employed (i.e. isotropic and anisotropic). Compacted kaolin soil samples were subjected to constant suction loading and unloading stages at three different values of suction, covering both unsaturated conditions (s= 50kPa and s= 300kPa) and saturated conditions (s=0). Loading and unloading stages were performed at three different values of stress ratio (η=0, η=1 and η=-1 ). In some tests a wetting-drying cycle was performed before or within the loading stage, with the wetting-drying cycles including both wetting-induced swelling and wetting-induced collapse compression. BEE tests were performed at regular intervals throughout all test stages, to measure shear wave velocity Vs and compression wave velocity Vp and hence to determine values of shear modulus G and constrained modulus M. The experimental test programme was designed to investigate how very small strain shear modulus G and constrained modulus M varied with unsaturated state variables, including how anisotropy of these parameters developed either with stress state (stress-induced anisotropy) or with previous straining (strain-induced anisotropy). A new expression has been proposed for the very small strain shear modulus G of an isotropic soil under saturated and unsaturated conditions. This expression relates the variation of G to only mean Bishop’s stress p* and specific volume v, and it converges to a well-established expression for saturated soils as degree of saturation approaches 1. The proposed expression for G is able to predict the variation of G under saturated and unsaturated conditions at least as well as existing expressions from the literature and it is considerably simpler (employing fewer state variables and fewer soil constants). In addition, unlike existing expressions from the literature, the values of soil constants in the proposed new expression can be determined from a saturated test. It appeared that, in the current project at least, any strain-induced anisotropy of very small strain elastic behaviour was relatively modest, with the possible exception of loading in triaxial extension. It was therefore difficult to draw any firm conclusion about evolution of strain-induced anisotropy and whether it depended upon the same aspects of soil fabric as evolution of anisotropy of large strain plastic behaviour. Stress-induced anisotropy of very small strain elastic behaviour was apparent in the experimental test programme. An attempt was made to extend the proposed expression for G to include the effect of stress-induced anisotropy. Interpretation of the experimental results indicated that the value of shear modulus was affected by the values of all three principal Bishop’s stresses (in the direction of wave transmission, the direction of wave polarisation and the third mutually perpendicular direction). However, prediction of stress-induced anisotropy was only partially successful, and it was concluded that the effect of Lode angle was also significant.

Effects of dry density and grain size distribution on soil-water characteristic curves of sandy soils

The unsaturated soil mechanics is receiving increasing attention from researchers and as well as from practicing engineers. However, the requirement of sophisticated devices to measure unsaturated soil properties and time consumption have made the geotechnical engineers keep away from implication of the unsaturated soil mechanics for solving practical geotechnical problems. The application of the conventional laboratory devices with some modifications to measure unsaturated soil properties can promote the application of unsaturated soil mechanics into engineering practice. Therefore, in the present study, a conventional direct shear device was modified to measure unsaturated shear strength parameters at low suction. Specially, for the analysis of rain-induced slope failures, it is important to measure unsaturated shear strength parameters at low suction where slopes become unstable. The modified device was used to measure unsaturated shear strength of two silty soils at low suction values (0 ~ 50 kPa) that were achieved by following drying path and wetting path of soil-water characteristic curves (SWCCs) of soils. The results revealed that the internal friction angle of soil was not significantly affected by the suction and as well as the drying-wetting SWCCs of soils. The apparent cohesion of soil increased with a decreasing rate as the suction increased. Further, the apparent cohesion obtained from soil in wetting was greater than that obtained from soil in drying. Shear stress-shear displacement curves obtained from soil specimens subjected to the same net normal stress and different suction values showed a higher initial stiffness and a greater peak stress as the suction increased. In addition, it was observed that soil became more dilative with the increase of suction. A soil in wetting exhibited slightly higher peak shear stress and more contractive volume change behaviour than that of in drying at the same net normal stress and the suction.

Laboratory measurement of hydraulic conductivity functions of two unsaturated sandy soils during drying and wetting processes

The importance of applying unsaturated soil mechanics to geotechnical engineering design has been well understood. However, the consumption of time and the necessity for a specific laboratory testing apparatus when measuring unsaturated soil properties have limited the application of unsaturated soil mechanics theories in practice. Although methods for predicting unsaturated soil properties have been developed, the verification of these methods for a wide range of soil types is required in order to increase the confidence of practicing engineers in using these methods. In this study, a new permeameter was developed to measure the hydraulic conductivity of unsaturated soils using the steady-state method and directly measured suction (negative pore-water pressure) values. The apparatus is instrumented with two tensiometers for the direct measurement of suction during the tests. The apparatus can be used to obtain the hydraulic conductivity function of sandy soil over a low suction range (0-10 kPa). Firstly, the repeatability of the unsaturated hydraulic conductivity measurement, using the new permeameter, was verified by conducting tests on two identical sandy soil specimens and obtaining similar results. The hydraulic conductivity functions of the two sandy soils were then measured during the drying and wetting processes of the soils. A significant hysteresis was observed when the hydraulic conductivity was plotted against the suction. However, the hysteresis effects were not apparent when the conductivity was plotted against the volumetric water content. Furthermore, the measured unsaturated hydraulic conductivity functions were compared with predictions using three different predictive methods that are widely incorporated into numerical software. The results suggest that these predictive methods are capable of capturing the measured behavior with reasonable agreement.

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.

Relative performances of textural models in estimating soil moisture characteristic

The soil moisture characteristic (SMC) forms an important input to mathematical models of water and solute transport in the unsaturated-soil zone. Owing to their simplicity and ease of use, texture-based regression models are commonly used to estimate the SMC from basic soil properties. In this study, the performances of six such regression models were evaluated on three soils. Moisture characteristics generated by the regression models were statistically compared with the characteristics developed independently from laboratory and in-situ retention data of the soil profiles. Results of the statistical performance evaluation, while providing useful information on the errors involved in estimating the SMC, also highlighted the importance of the nature of the data set underlying the regression models. Among the models evaluated, the one possessing an underlying data set of in-situ measurements was found to be the best estimator of the in-situ SMC for all the soils. Considerable errors arose when a textural model based on laboratory data was used to estimate the field retention characteristics of unsaturated soils.