The presence of expansive soils underneath pavement structures is considered one of the most common sources of pavement distresses, due to differential settlements caused by differential moisture distribution attributed to soil heterogeneity and seasonal climatic fluctuations. The cost of the repairs to the infrastructure caused by expansive soils is estimated to exceed 10 billion dollars annually in the US, as reported by Puppala and Cerato (2009). Although many studies have been developed to better understand the volume change of unsaturated soils and incorporate the effect of swelling/shrinkage behavior into pavement design procedures, current methodologies are still based on simple correlations with index properties or other empirical methods. Such solutions lead to poor or uneconomical design practices. The objective of this study was to calibrate and implement a new mechanistic, stochastic model that predicts pavement distresses caused by the presence of expansive soils. Three major tasks were completed to fulfill the objective of this study: 1) a laboratory research program performed to estimate the volume change of compacted specimens, with different expansion potential, due to the simultaneous application of suction and net normal stresses, 2) the calibration of a new mechanistic free-swell model for expansive soils tailored to pavement structures, based on elevation information collected from the Long Term Pavement Performance (LTPP) program, and 3) the incorporation and calibration of the free-swell stochastic model results into the current Pavement Mechanistic-Empirical (ME) Design procedure using the International Roughness Index (IRI) models. The results presented includes: 1) an empirical model to estimate volume change due to the coupled effect of suction, and net normal stresses, for soils with different soil index properties, 2) a calibrated model to adjust the free-swell results of the mechanistic-stochastic model developed by Olaiz et al. (2021), and 3) an updated IRI equation for asphalt concrete pavements to account for volume change fluctuations due to changes in suction stress conditions. The models presented can be easily implemented into currently available pavement design procedures and greatly improves over the existing empirical models that have been used for more than four decades.

Microbially induced desaturation (MID) via denitrification is an emerging ground improvement technique to mitigate liquefaction by stimulating the metabolic processes of native bacteria to produce biogas, biominerals and biomass. The production of biogenic gas gradually lowers the degree of saturation of treated soils, thereby dampening the pore pressure response to cyclic loading. However, the production of these metabolic products also alters the hydraulic and mechanical properties of the soil. A series of four tank tests simulating two-dimensional plane strain conditions were performed to evaluate the effectiveness of MID and the resulting changes to the hydraulic properties of the soils. Previous studies have demonstrated the mechanical response for treated homogenous granular soils at the bench scale via vertical injection methods. However, limited knowledge is available on the impact of partial desaturation on the hydraulic properties of the soil, particularly in stratified formations. Treating larger granular soil specimens via lateral injection methods is important for the up-scaling and future commercialization of the process as it may affect injection strategies, and the distribution of substrates and metabolic products. Tank tests were performed on a layered natural soil sediment collected from Richmond, British Columbia, Canada, as well as layered and unlayered laboratory grade Ottawa sands of different grain size distributions. The results demonstrated the effectiveness of treatment upon macro-scale soil properties, and showed how gas formation, migration and entrapment, and resulting degree of desaturation and hydraulic conductivity are affected by micro and macro-stratifications in granular soils.

This document presents the assessment of the swelling behavior of expansive clay stabilized with bio-based silica gel and subjected to wetting and drying cycles. The expansive clay used in this research was obtained from Anthem, Arizona. Rice husk is a rich silica by-product of rice production with commercial uses and applications in the industry. Rice husk ash from two different sources -California (named ASU) and India- were subjected to chemical characterization. Fourier Transform Infra-red Spectroscopy was used to verify the functional groups of the gel formed. Results showed differences between the ashes from different sources and confirmed the presence of silica structure bonds. X-Ray Diffraction (XRF) results showed that the ASU ash contained more amorphous silica than the Indian ash.One dimensional swell and consolidation tests were performed to investigate the volume change behavior of the untreated and silica gel treated remoulded samples. The free swell of the clay decreased from 12.3% (untreated sample) to 7.2% (ASU sample) and 11.4% (Indian sample). The effect of the wet and dry cycles on the swelling and consolidation characteristics of the untreated clay demonstrated that the treatment is irreversible after three cycles. Swelling of clay treated with ASU ash was reduced after the first cycle, while that of the clay treated with Indian ash was reduced after three cycles. This was due to the gelation time difference between treatments. Scanning Electron Microscopy images showed that the structure of the untreated clay was discontinuous, flaky and without aggregations whereas particles in the treated samples were aggregated and new bonds were created, decreasing the surface area. The X-Ray Diffraction (XRF) results showed that the main mineral responsible for expansive behavior of the clay studied was illite. The d-spacing of the illite decreased from 4.47Å for the untreated clay to 3.33Å for the treated clay. This study demonstrates a promising technique for clay swelling reduction and a more sustainable solution than that available to current practicing engineering.

Recent research efforts have been directed to improve the quality of pavement design procedures by considering the transient nature of soil properties due to environmental and aging effects on pavement performance. The main purpose of this research study was to investigate the existence of subgrade soil moisture changes that may have arisen due to thermal and hydraulic gradients at the Atlantic City NAPTF and to evaluate their effect on the material stiffness and the California Bearing Ratio (CBR) strength parameter of the clay subgrade materials. Laboratory data showed that at the same water content, matric suction decreases with increasing temperature; and at the same suction, hydraulic conductivity increases with increasing temperature. Models developed, together with moisture/temperature data collected from 30 sensors installed in the test facility, yielded a maximum variation of suction in field of 155 psi and changes in hydraulic conductivity from 2.9E-9 m/s at 100% saturation to 8.1E-12 at 93% saturation. The maximum variation in temperature was found to be 20.8oC at the shallower depth and decreased with depth; while a maximum variation in moisture content was found to be 3.7% for Dupont clay and 4.4% for County clay. Models developed that predicts CBR as a function of dry density and moisture content yielded a maximum variation of CBR of 2.4 for Dupont clay and 2.9 for County clay. Additionally, models were developed relating the temperature with the bulk stress and octahedral stress applied on the subgrade for dual gear, dual tandem and triple tandem gear types for different tire loads. It was found that as the temperature increases the stresses increase. A Modified Cary and Zapata model was used for predicting the resilient modulus(Mr) of the subgrade. Using the models developed and the temperature/moisture changes observed in the field, the variation of suction, bulk and octahedral stresses were estimated, along with the resilient modulus for three different gear types. Results indicated that changes in Mr as large as 9 ksi occur in the soils studied due to the combined effect of external loads and environmental condition changes.

The geotechnical community typically relies on recommendations made from numerical simulations. Commercial software exhibits (local) numerical instabilities in layered soils across soil interfaces. This research work investigates unsaturated moisture flow in layered soils and identifies a possible source of numerical instabilities across soil interfaces and potential improvement in numerical schemes for solving the Richards' equation. The numerical issue at soil interfaces is addressed by a (nonlinear) interface problem. A full analysis of the simplest soil hydraulic model, the Gardner model, identifies the conditions of ill-posedness of the interface problem. Numerical experiments on various (more advanced and practical) soil hydraulic models show that the interface problem can also be ill-posed under certain circumstances. Spurious numerical ponding and/or oscillations around soil interfaces are observed consequently. This work also investigates the impact of different averaging schemes for cell-centered conductivities on the propensity of ill-posedness of the interface problem and concludes that smaller averaging conductivities are more likely to trigger numerical instabilities. In addition, an agent-based stochastic soil model, with hydraulic properties defined at the finite difference cell level, results in a large number of interface problems. This research compares sequences of stochastic realizations in heterogeneous unsaturated soils with the numerical solution using homogenized soil parameters. The mean of stochastic realizations is not identical to the solution obtained from homogenized soil parameters.

The influence of temperature on soil engineering properties is a major concern in the design of engineering systems such as radioactive waste disposal barriers, ground source heat pump systems and pavement structures. In particular, moisture redistribution under pavement systems might lead to changes in unbound material stiffness that will affect pavement performance. Accurate measurement of thermal effects on unsaturated soil hydraulic properties may lead to reduction in design and construction costs. This thesis presents preliminary results of an experimental study aimed at determining the effect of temperature on the soil water characteristic curve (SWCC) and the unsaturated hydraulic conductivity function (kunsat). Pressure plate devices with volume change control were used to determine the SWCC and the instantaneous profile method was used to obtain the kunsat function. These properties were measured on two fine-grained materials subjected to controlled temperatures of 5oC, 25oC and 40oC. The results were used to perform a sensitivity analysis of the effect of temperature changes on the prediction of moisture movement under a covered area. In addition, two more simulations were performed where changes in hydraulic properties were done in a stepwise fashion. The findings were compared to field measured water content data obtained on the subgrade material of the FAA William Hughes test facility located in Atlantic City. Results indicated that temperature affects the unsaturated hydraulic properties of the two soils used in the study. For the DuPont soil, a soil with high plasticity, it was found that the water retention was higher at low temperatures for suction levels lower than about 10,000 kPa; while the kunsat functions at the three temperatures were not significantly different. For the County soil, a material with medium plasticity, it was found that it holds around 10% more degree of saturation at 5°C than that at 40°C for suction levels higher than about 1,000 kPa; while the hydraulic conductivity at 40°C was at least one order of magnitude higher than that at 5°C, for suction levels higher than 1,000 kPa. These properties were used to perform two types of numerical analyses: a sensitivity analysis and stepwise analysis. Absolute differences between predicted and field measured data were considered to be acceptable, ranging from 4.5% to 9% for all simulations. Overall results show an improvement in predictions when non-isothermal conditions were used over the predictions obtained with isothermal conditions.

In a laboratory setting, the soil volume change behavior is best represented by using various testing standards on undisturbed or remolded samples. Whenever possible, it is most precise to use undisturbed samples to assess the volume change behavior but in the absence of undisturbed specimens, remodeled samples can be used. If that is the case, the soil is compacted to in-situ density and water content (or matric suction), which should best represent the expansive profile in question. It is standard practice to subject the specimen to a wetting process at a particular net normal stress. Even though currently accepted laboratory testing standard procedures provide insight on how the profile conditions changes with time, these procedures do not assess the long term effects on the soil due to climatic changes. In this experimental study, an assessment and quantification of the effect of multiple wetting/drying cycles on the volume change behavior of two different naturally occurring soils was performed. The changes in wetting and drying cycles were extreme when comparing the swings in matric suction. During the drying cycle, the expansive soil was subjected to extreme conditions, which decreased the moisture content less than the shrinkage limit. Nevertheless, both soils were remolded at five different compacted conditions and loaded to five different net normal stresses. Each sample was subjected to six wetting and drying cycles. During the assessment, it was evident from the results that the swell/collapse strain is highly non-linear at low stress levels. The strain-net normal stress relationship cannot be defined by one single function without transforming the data. Therefore, the dataset needs to be fitted to a bi-modal logarithmic function or to a logarithmic transformation of net normal stress in order to use a third order polynomial fit. It was also determined that the moisture content changes with time are best fit by non-linear functions. For the drying cycle, the radial strain was determined to have a constant rate of change with respect to the axial strain. However, for the wetting cycle, there was not enough radial strain data to develop correlations and therefore, an assumption was made based on 55 different test measurements/observations, for the wetting cycles. In general, it was observed that after each subsequent cycle, higher swelling was exhibited for lower net normal stress values; while higher collapse potential was observed for higher net normal stress values, once the net normal stress was less than/greater than a threshold net normal stress value. Furthermore, the swelling pressure underwent a reduction in all cases. Particularly, the Anthem soil exhibited a reduction in swelling pressure by at least 20 percent after the first wetting/drying cycle; while Colorado soil exhibited a reduction of 50 percent. After about the fourth cycle, the swelling pressure seemed to stabilized to an equilibrium value at which a reduction of 46 percent was observed for the Anthem soil and 68 percent reduction for the Colorado soil. The impact of the initial compacted conditions on heave characteristics was studied. Results indicated that materials compacted at higher densities exhibited greater swell potential. When comparing specimens compacted at the same density but at different moisture content (matric suction), it was observed that specimens compacted at higher suction would exhibit higher swelling potential, when subjected to the same net normal stress. The least amount of swelling strain was observed on specimens compacted at the lowest dry density and the lowest matric suction (higher water content). The results from the laboratory testing were used to develop ultimate heave profiles for both soils. This analysis showed that even though the swell pressure for each soil decreased with cycles, the amount of heave would increase or decrease depending upon the initial compaction condition. When the specimen was compacted at 110% of optimum moisture content and 90% of maximum dry density, it resulted in an ultimate heave reduction of 92 percent for Anthem and 685 percent for Colorado soil. On the other hand, when the soils were compacted at 90% optimum moisture content and 100% of the maximum dry density, Anthem specimens heave 78% more and Colorado specimens heave was reduced by 69%. Based on the results obtained, it is evident that the current methods to estimate heave and swelling pressure do not consider the effect of wetting/drying cycles; and seem to fail capturing the free swell potential of the soil. Recommendations for improvement current methods of practice are provided.

Urbanization and infrastructure development often brings dramatic changes in the surface and groundwater regimes. These changes in moisture content may be particularly problematic when subsurface soils are moisture sensitive such as expansive soils. Residential foundations such as slab-on ground may be built on unsaturated expansive soils and therefore have to resist the deformations associated with change in moisture content (matric suction) in the soil. The problem is more pronounced in arid and semi arid regions with drying periods followed by wet season resulting in large changes in soil suction. Moisture content change causes volume change in expansive soil which causes serious damage to the structures. In order to mitigate these ill effects various mitigation are adopted. The most commonly adopted method in the US is the removal and replacement of upper soils in the profile. The remove and replace method, although heavily used, is not well understood with regard to its impact on the depth of soil wetting or near-surface differential soil movements. In this study the effectiveness of the remove and replace method is studied. A parametric study is done with various removal and replacement materials used and analyzed to obtain the optimal replacement depths and best material. The depth of wetting and heave caused in expansive soil profile under climatic conditions and common irrigation scenarios are studied for arid regions. Soil suction changes and associated soil deformations are analyzed using finite element codes for unsaturated flow and stress/deformation, SVFlux and SVSolid, respectively. The effectiveness and fundamental mechanisms at play in mitigation of expansive soils for remove and replace methods are studied, and include (1) its role in reducing the depth and degree of wetting, and (2) its effect in reducing the overall heave potential, and (3) the effectiveness of this method in pushing the seat of movement deeper within the soil profile to reduce differential soil surface movements. Various non-expansive replacement layers and different surface flux boundary conditions are analyzed, and the concept of optimal depth and soil is introduced. General observations are made concerning the efficacy of remove and replace as a mitigation method.

Trenchless technology is a group of techniques whose utilization allows for the installation, rehabilitation, and repair of underground infrastructure with minimal excavation from the ground surface. As the built environment becomes more congested, projects are trending towards using trenchless technologies for their ability to quickly produce a quality product with minimal environmental and social costs. Pilot tube microtunneling (PTMT) is a trenchless technology where new pipelines may be installed at accurate and precise line and grade over manhole to manhole distances. The PTMT process can vary to a certain degree, but typically involves the following three phases: jacking of the pilot tube string to achieve line and grade, jacking of casing along the pilot bore and rotation of augers to excavate the borehole to a diameter slightly larger than the product pipe, and jacking of product pipe directly behind the last casing. Knowledge of the expected productivity rates and jacking forces during a PTMT installation are valuable tools that can be used for properly weighing its usefulness versus competing technologies and minimizing risks associated with PTMT. This thesis outlines the instrumentation and monitoring process used to record jacking frame hydraulic pressures from seven PTMT installations. Cyclic patterns in the data can be detected, indicating the installation of a single pipe segment, and enabling productivity rates for each PTMT phase to be determined. Furthermore, specific operations within a cycle, such as pushing a pipe or retracting the machine, can be observed, allowing for identification of the critical tasks associated with each phase. By identifying the critical tasks and developing more efficient means for their completion, PTMT productivity can be increased and costs can be reduced. Additionally, variations in depth of cover, drive length, pipe diameter, and localized ground conditions allowed for trends in jacking forces to be identified. To date, jacking force predictive models for PTMT are non-existent. Thus, jacking force data was compared to existing predictive models developed for the closely related pipe jacking and microtunneling methodologies, and the applicability of their adoption for PTMT jacking force prediction was explored.