Analysis of Deep Moisture Barriers in Expansive Soils. I: Constitutive Model Formulation

Expansive soils cause important economical losses in many arid or semiarid countries in the world. Considering the large economic impact, relatively few efforts have been devoted to develop analytical methods that may help practitioner engineers to adequately design civil infrastructure on this type of soil. A rational design method should be able to quantify the heave or subsidence of the soil associated with the suction changes during water diffusion, as well as the contact pressures on soil-structure interfaces. Accordingly, in this and in a companion paper, the problem of volume changes due to nonpermanent water flow in expansive soils is studied and applied to the case of vertical moisture barriers. In this paper, a constitutive model for expansive soils is proposed. This model is an extension of that developed by Alonso et al. in 1990, in the sense that it can take into account the behavior of expansive soils. The advanced model is evaluated by comparing the numerical results with experimental data.     

Modeling Expansive Soil Movements Beneath Structures

This paper describes the use of an existing numerical model which is capable of generating continuous records over time of 3D soil suction profiles beneath a structure. The model uses recorded climatic data and representative soil properties. The model was used to obtain data required for the simulation of ground movements and the resulting structural response in a separate soil/structure interaction model (the results of the subsequent soil/structure interaction simulations are not reported). The factors influencing the soil moisture distribution beneath a structure were identified and careful consideration was given to quantifying the variability in these factors. The model explicitly captures the long-term moisture redistribution occurring beneath a structure as a result of introducing a ground cover. It also captures the effect of different construction dates on soil moisture conditions. If a range of construction dates are selected at random and 3D simulations over appropriate time periods are conducted, the variability in ground movements due to seasonal and long-term climatic effects and due to the choice of construction date can be quantified.     

Method for Evaluation of Depth of Wetting in Residential Areas

Field sampling was performed at a number of residential structures in the Denver metropolitan area for the purpose of assessing the extent of wetting below residential structures after construction and commencement of landscape irrigation. Total suction measurements using filter paper methods were undertaken on undisturbed samples. A similar data set from sites without previous development or irrigation was used to estimate the suction profile before the imposition of residential construction and landscape irrigation. Comparison of the postdevelopment profile measured at the residential structures to the predevelopment profile estimated using a site-specific procedure was used to assess the depth of wetting at each structure. Cumulative probability curves are presented for the depth of wetting arising from common residential landscape schemes and site drainage for the Denver area. The Denver-area data set was used to develop a method for assessment of depth of wetting. This method can be used to assess depth of wetting from residential development for other regions having different climatic conditions and landscape practices.     

Critical Review of the Methodologies Employed for Soil Suction Measurement

Modeling the behavior of unsaturated soils necessitates the measurement of soil suction and the establishment of its variation with the water content, which is commonly known as the soil-water characteristic curve (SWCC). Several methodologies have been developed for measuring either total suction ψ (sum of matric suction ψm and osmotic suction ψo) or ψm. While employing different methodologies for suction measurement, there is a possibility that various factors (viz., type of the soil, measurement methodology, range of the suction measurement, equilibration time, and presence of salts or contaminants in the soil) may influence the results and hence the SWCC. Therefore, it is essential to investigate the uniqueness of SWCC, determined by using some commonly adopted suction measurement methodologies. This study indicates that the SWCC established by adopting different methodologies may not be unique and is primarily influenced by the range of suction measurement. As such, it is essential to highlight the range of suction values involved for establishing the SWCC, to facilitate unambiguous modeling and to precisely understand the behavior of unsaturated soil.     

Soil Moisture Flow Modeling with Water Uptake by Plants (Wheat) under Varying Soil and Moisture Conditions

A variably saturated soil moisture flow model is developed for planted soils with depth varying properties by incorporating a nonuniform macroscopic root water uptake function. The model includes spatial and temporal variation of the root density with dynamic root growth for simulating water uptake by plants along with the impact of soil moisture availability. The governing partial differential moisture flow equation integrated over the depth with a plant water uptake term is solved numerically by the implicit finite difference method using an iterative scheme. The model is first tested for barren soils for two profiles considering constant and depth varying soil characteristics under constant inflow condition. The results obtained are later tested with experimental data available in the literature. A nonuniform plant water uptake term is subsequently incorporated in the model and water uptake by wheat plants under different soil moisture availability conditions is studied. Finally, the moisture flow model is validated with field data of rain fed wheat (Triticum aestivum) using a dynamic root growth model for a layered root zone soil profile. The simulated soil moisture regime of the layered root zone shows a reasonably good agreement with the observed data.     

Soil-Water Characteristic Curve Equation with Independent Properties

The soil-water characteristic curve (SWCC) has traditionally been represented using equations whose fitting parameters do not individually correspond to clearly defined soil properties or to features of the curve. As a result, unique sets of parameters are often nonexistent, and sensitivity analyses and statistical assessments of SWCC parameters become difficult. In order to overcome these difficulties, a new class of equations to represent unimodal and bimodal SWCCs is proposed. The chosen fitting parameters are the air-entry value, the residual suction, the residual degree of saturation, and a parameter that controls the sharpness of the curvatures. The physical meaning for the soil parameters is discussed for different soil types. A unique relation between each of the equation parameters and the individual features of SWCCs is assured. The proposed equations are fitted to data corresponding to a variety of soil types and a good fit is observed.     

Study of Expansive Soils and Residential Foundations on Expansive Soils in Arizona

Construction on expansive soils is challenging and thus prone to some problems and litigation. The engineering community makes extensive use of local experience and empirical procedures to address these problems. Although there has been extensive study of expansive soils and foundations on expansive soils, data related to performance of residential structures are limited in general and limited in the Phoenix area, in particular. In this study, an overview of the Phoenix Valley, Arizona, geotechnical practice and foundation performance related to residential structures on expansive clays, was developed through surveys and interviews with geotechnical engineers, structural engineers, and homebuilders. Using data obtained from files of Phoenix area geotechnical firms and government agencies, the existing Natural Resource Conservation Service map showing expansive soil locations throughout the Phoenix region was updated through the use of correlation developed in this study relating expansion index to common soil index properties such as Atterberg limits and percent passing the No. 200 sieve. Files of forensic investigations linked to expansive soil regions were made available for this study by several geotechnical engineering firms, and Phoenix Valley areas where forensic investigations have been identified, were mapped for comparison to regions identified in the updated map as having expansive soils. Comparison of the forensic investigation map to the updated map of expansive clay locations revealed that most of the forensic investigations were in regions identified with clays labeled as high to moderately high expansion potential, with a few forensic investigations in regions of medium expansion potential. Finally, unsaturated flow analyses were conducted for an Arizona expansive clay profile for two very different landscaped conditions of well-irrigated turf and desert landscape. The results of the numerical analyses were consistent with the reported observations and modes of failure identified through the surveys and interviews conducted with engineering and homebuilder professionals, including the finding that site drainage was found to be extremely important to good foundation performance, regardless of the type of landscape selected.     

Experimental Study on Moisture Migration in Unsaturated Loess under Effect of Temperature

In this paper, moisture migration in loess considering temperature effect is studied by tests on unsaturated loess samples with different densities and initial moisture contents. Test results reveal that obvious changes in moisture content distribution in a loess sample can be observed after temperature difference is exerted on the two ends of the sample. Moisture content at the cold end increases and that at the hot end decreases. Under the effect of temperature difference, moisture content difference at the two ends of a soil sample is related to the initial moisture content, soil density, and magnitude of the temperature difference. Generally speaking, larger temperature differences and smaller soil densities result in more obvious moisture migration and larger moisture content differences at the two ends of the soil sample. When the initial moisture content is large, the moisture content difference caused by a temperature difference is small; when the initial moisture content is small, the moisture content difference caused by a temperature difference is also small; when the initial moisture content is moderate, the moisture content difference caused by a temperature difference is large. After the analysis of test results, taking the soil density and moisture content into account, a formula is obtained to determine the moisture content gradient resulting from the temperature gradient. Reliability of the formula is verified by comparing the measured and calculated data. Because of the reverse migration of liquid water and water vapor at the end of the experiment, it is difficult to determine the thermal potential and matrix potential. Based on the experimental data, this paper probes into the water potential equation that can be used for stability analysis. The equation considers the comprehensive impact of soil density, temperature gradient, moisture content, and moisture content gradient on water potential. It only applies to analyze stable distributions of temperature and does not apply to unstable temperature distributions.     

Profiles of Steady-State Suction Stress in Unsaturated Soils

Application of the effective stress principle in unsaturated geotechnical engineering problems often requires explicit knowledge of the stress acting on the soil skeleton due to suction pore water pressure. This stress is defined herein as the suction stress. A theoretical formulation of suction stress profiles, based on the soil water characteristics curve, the soil permeability characteristic curve, and previous shear strength experimental verification, is developed. The theory provides a general quantitative way to calculate vertical suction stress profiles in various unsaturated soils under steady flow rate in the form of infiltration or evaporation.     

Soil-Water Characteristic Curves of Stabilized Expansive Soils

The engineering properties of expansive soils are conventionally improved through the use of additives such as fly ash, lime, and chemical additives. Such soils are often referred to as stabilized or modified or treated expansive soils. The soil-water characteristic curves (SWCC) of two expansive soils from Texas were measured both in natural and stabilized conditions using the pressure plate apparatus in the suction range of 0-1,000 kPa. The SWCC results are used to interpret the expansive soil behavior due to stabilizer treatment. In addition, relationships were developed between the basic soil and stabilizer properties such as water content, dry density, liquid limit, plastic limit, and stabilizer dosages and the model constants of the SWCC formulation of Fredlund and Xing via multiple linear regression analysis. The analysis showed that higher coefficients of correlations can be achieved by using six independent soil properties. The comparisons between the predicted and measured volumetric water contents are within ±20% for ash-treated expansive soils, and within ±15% for combined ash- and fiber-treated expansive soils. The research data and interpretation analysis presented here can be extended to understand volume change behaviors of other stabilized expansive soils using the SWCC test data.     

Impacts of Unsaturated Zone Soil Moisture and Groundwater Table on Slope Instability

The combined effect of soil moisture in unsaturated soil layers and pore-water pressure in saturated soil layers is critical to predict landslides. An improved infinite slope stability model, that directly includes unsaturated zone soil moisture and groundwater, is derived and used to analyze the factor of safety’s sensitivity to unsaturated zone soil moisture. This sensitivity, the change in the factor of safety with respect to variable unsaturated zone soil moisture, was studied at local and regional scales using an active landslide region as a case study. Factors of safety have the greatest sensitivity to unsaturated zone soil moisture dynamics for shallow soil layers (<2?m) and comparatively deep groundwater tables (1 m). For an identical groundwater table, the factor of safety for a 1 m thick soil mantle was four times more sensitive to soil moisture changes than a 3-m thick soil. At a regional scale, the number of unstable areas increases nonlinearly with increasing unsaturated zone soil moisture and with moderately wet slopes exhibiting the greatest sensitivity.     

Analysis of Deep Moisture Barriers in Expansive Soils. II: Water Flow Formulation and Implementation

Several alternatives have been proposed to prevent damage to civil infrastructure founded on expansive soils. For example, deep moisture barriers have been used in highways and buildings. However, in some cases, the protected lanes or structures degrade to similar levels as the unprotected ones, although at a smaller rate. In spite of these poor results, relatively few efforts have been devoted to the development of analytical methods for rational designs on expansive soils. This paper couples the constitutive model for expansive soils developed in the companion paper with flow equations in a deforming medium and a finite element code is developed. The resulting numerical tool has the capacity of computing soil suction and volumetric strain changes of expansive soils under a defined wetting–drying regime. To verify the capabilities of this computer code, a laboratory barrier model was built. The model was instrumented to measure soil suction changes and the corresponding surface displacements. The experimental and theoretical results were compared. Finally, the numerical model was applied to a design example of a deep moisture barrier.     

Three-Dimensional Field-Scale Coupled Thermo-Hydro-Mechanical Modeling: Parallel Computing Implementation

An approach for the simulation of three-dimensional field-scale coupled thermo-hydro-mechanical problems is presented, including the implementation of parallel computation algorithms. The approach is designed to allow three-dimensional large-scale coupled simulations to be undertaken in reduced time. Owing to progress in computer technology, existing parallel implementations have been found to be ineffective, with the time taken for communication dominating any reduction in time gained by splitting computation across processors. After analysis of the behavior of the solver and the architecture of multicore, nodal, parallel computers, modification of the parallel algorithm using a novel hybrid message passing interface/open multiprocessing (MPI/OpenMP) method was implemented and found to yield significant improvements by reducing the amount of communication required. This finding reflects recent enhancements of current high-performance computing architectures. An increase in performance of 500% over existing parallel implementations on current processors was achieved for the solver. An example problem involving the Prototype Repository experiment undertaken by the Swedish Nuclear Fuel and Waste Management Co. [Svensk K?rnbr?nslehantering AB (SKB)] in ?sp?, Sweden, has been presented to demonstrate situations in which parallel computation is invaluable because of the complex, highly coupled nature of the problem.     

Inclusion of Some Aspects of Chemical Behavior of Unsaturated Soil in Thermo/Hydro/Chemical/Mechanical Models. II: Application and Transport of Soluble Salts in Compacted Bentonite

This paper presents an application of a coupled thermo/hydro/chemical/mechanical model via simulation of a laboratory experiment in order to investigate the transport behavior of ions in bentonite pore water. Chemical reactions considered include ion exchange reactions involving major cations (Na+, K+, Mg2+, and Ca2+) and precipitation-dissolution of trace minerals (calcite, dolomite, anhydrite, and halite). The following conclusions are drawn based on the numerical results. The development of both the temperature and moisture fields was captured by simulation, and a good correlation with the experimental water uptake results was observed. For all ions, the model showed a good qualitative and reasonable quantitative agreement with the experimental results.     

Centrifuge Permeameter for Unsaturated Soils. II: Measurement of the Hydraulic Characteristics of an Unsaturated Clay

This paper presents the hydraulic characteristics of an unsaturated, compacted clay, including its soil-water retention curve (SWRC) and hydraulic conductivity function (K function), determined using a new centrifuge permeameter developed at the University of Texas at Austin. A companion paper describes the apparatus, its instrumentation layout, and data reduction procedures. Three approaches are evaluated in this study to define the SWRC and K function of the compacted clay under both drying and wetting paths, by varying the inflow rate, the g level, or both. For imposed inflow rates ranging from 20 to 0.1 mL/h and g levels ranging from 10 to 100 g, the measured matric suction ranged from 5 to 70 kPa, the average volumetric water content ranged from 23 to 33%, and the hydraulic conductivity ranged from 2×10?7 to 8×10?11?m/s. The SWRCs and K functions obtained using the three different testing approaches were very consistent, and yielded suitable information for direct determination of the hydraulic characteristics. The approaches differed in the time required to complete a testing stage and in the range of measured hydraulic conductivity values. The g level had a negligible effect on the measured hydraulic characteristics of the compacted clay. The SWRCs and K functions defined using the centrifuge permeameter are consistent with those obtained using pressure chamber and column infiltration tests. The K functions defined using the centrifuge permeameter follow the same shape as those obtained from predictive relationships, although the measured and predicted K functions differ by two orders of magnitude at the lower end of the volumetric water content range.     

Elastic Model for Partially Saturated Granular Materials

This paper presents the development of an elastic model for partially saturated granular materials based on micromechanical factor consideration. A granular material is considered as an assembly of particles. The stress-strain relationship for an assembly can be determined by integrating the behavior at all interparticle contacts and by using a static hypothesis, which relates the average stress of the granular assembly to a mean field of particle contact forces. As for the nonsaturated state, capillary forces at grain contacts are added to the contact forces created by an external load. These are then calculated as a function of the degree of saturation, depending on the grain size distribution and on the void ratio of the granular assembly. Hypothesizing a Hertz-Mindlin law for the grain contacts leads to an elastic nonlinear behavior of the particulate material. The prediction of the stress-strain model is compared to experimental results obtained from several different granular materials in dry, partially saturated and fully saturated states. The numerical predictions demonstrate that the model is capable of taking into account the influence of key parameters, such as degree of saturation, void ratio, and mean stress.     

Unsaturated Infinite Slope Stability Considering Surface Flux Conditions

A slope stability model is derived for an infinite slope subjected to unsaturated infiltration flow above a phreatic surface. Closed form steady state solutions are derived for the matric suction and degree of saturation profiles. Soil unit weight, consistent with the degree of saturation profile, is also directly calculated and introduced into the analyzes, resulting in closed-form solutions for typical soil parameters and an infinite series solution for arbitrary soil parameters. The solutions are coupled with the infinite slope stability equations to establish a fully realized safety factor function. In general, consideration of soil suction results in higher factor of safety. The increase in shear strength due to the inclusion of soil suction is analogous to making an addition to the cohesion, which, of course, increases the factor of safety against sliding. However, for cohesive soils, the results show lower safety factors for slip surfaces approaching the phreatic surface compared to those produced by common safety factor calculations. The lower factor of safety is due to the increased soil unit weight considered in the matric suction model but not usually accounted for in practice wherein the soil is treated as dry above the phreatic surface. The developed model is verified with a published case study, correctly predicting stability under dry conditions and correctly predicting failure for a particular storm.     

Inclusion of Some Aspects of Chemical Behavior of Unsaturated Soil in Thermo/Hydro/Chemical/Mechanical Models. I: Model Development

This paper presents an investigation of the inclusion of some aspects of chemical behavior within a model of coupled thermo/hydro/chemical/mechanical behavior of unsaturated soils. In particular, multicomponent reactive chemical transport behavior is addressed. The chemical transport model is based on the advection/dispersion/reaction equation, while geochemical reactions are considered via coupling with an established geochemical speciation model. A numerical solution of the governing differential equations is achieved by the use of the Galerkin-weighted residual method for spatial discretization and an implicit backward Eulerian finite-difference method for temporal discretization. The solution of the geochemical reactions is achieved externally to the main solution procedure. Coupling between the chemical transport and geochemical models is achieved via the implementation of both sequential iterative and sequential noniterative techniques. Three application problems are then presented to demonstrate the capability of the coupled model.     

Evaporation Theory for Deformable Soils

Desiccation of a deformable soil is basically the removal of water by evaporation, which is controlled by evaporativity and evaporability. Surface evaporation improves the trafficability, which is essential for the access of construction equipment in areas reclaimed with soft clay. The existing traditional methods for evaluating evaporation cannot account for the deformation of soils during evaporation. Therefore, a theoretical model for predicting the rate of evaporation from the surface of a deformable material is proposed. The model is based on a system of equations for coupled heat and mass transfer in unsaturated soils. The modified pressure plate extractor test and glass desiccator test were carried out to obtain the soil-water characteristic curve for a deformable soil. A column-drying test was conducted to investigate one-dimensional water flow, heat flow, and evaporation in the surface. A finite difference program was developed to solve the coupled nonlinear partial differential equations, which permits the study of liquid, diffusive vapor, and heat flows in the deformable soil. Comparison between measured and simulated values shows good agreement.     

Finite Difference Method for Computation of 1D Pollutant Migration through Saturated Homogeneous Soil Media

The physical processes such as advection, dispersion, and diffusion and interaction between the solution and the soil solids such as sorption, biodegradation, and retention processes have been considered in the governing equation used in the present study. Finite difference method has been adopted herein to solve the one-dimensional contaminant transport model to predict the pollutant migration through soil in waste landfill. In the finite difference technique, the velocity field is first determined within a hydrologic system, and these velocities are then used to calculate the rate of contaminant migration by solving the governing equation. A total of seven contaminants have been chosen for analysis to represent a wide variety of wastes both organic and inorganic. A computer software CONTAMINATE has been developed for solution of the contaminant transport model. Results of this study have been compared with existing analytical solution for validation of the proposed solution technique. Design charts for liners have also been developed to facilitate the designers. The liner thickness has been optimized by considering the effect of velocity of advection, dispersion coefficient, and geochemical reactions for all the contaminants of this study.