This dissertation describes development of a procedure for obtaining high quality, optical grade sand coupons from frozen sand specimens of Ottawa 20/30 sand for image processing and analysis to quantify soil structure along with a methodology for quantifying the microstructure from the images. A technique for thawing and stabilizing frozen core samples was developed using optical grade Buehler® Epo-Tek® epoxy resin, a modified triaxial cell, a vacuum/reservoir chamber, a desiccator, and a moisture gauge. The uniform epoxy resin impregnation required proper drying of the soil specimen, application of appropriate confining pressure and vacuum levels, and epoxy mixing, de-airing and curing. The resulting stabilized sand specimen was sectioned into 10 mm thick coupons that were planed, ground, and polished with progressively finer diamond abrasive grit levels using the modified Allied HTP Inc. polishing method so that the soil structure could be accurately quantified using images obtained with the use of an optical microscopy technique. Illumination via Bright Field Microscopy was used to capture the images for subsequent image processing and sand microstructure analysis. The quality of resulting images and the validity of the subsequent image morphology analysis hinged largely on employment of a polishing and grinding technique that resulted in a flat, scratch free, reflective coupon surface characterized by minimal microstructure relief and good contrast between the sand particles and the surrounding epoxy resin. Subsequent image processing involved conversion of the color images first to gray scale images and then to binary images with the use of contrast and image adjustments, removal of noise and image artifacts, image filtering, and image segmentation. Mathematical morphology algorithms were used on the resulting binary images to further enhance image quality. The binary images were then used to calculate soil structure parameters that included particle roundness and sphericity, particle orientation variability represented by rose diagrams, statistics on the local void ratio variability as a function of the sample size, and the local void ratio distribution histograms using Oda's method and Voronoi tessellation method, including the skewness, kurtosis, and entropy of a gamma cumulative probability distribution fit to the local void ratio distribution.

The importance of unsaturated soil behavior stems from the fact that a vast majority of infrastructures are founded on unsaturated soils. Research has recently been concentrated on unsaturated soil properties. In the evaluation of unsaturated soils, researchers agree that soil water retention characterized by the soil water characteristic curve (SWCC) is among the most important factors when assessing fluid flow, volume change and shear strength for these soils. The temperature influence on soil moisture flow is a major concern in the design of important engineering systems such as barriers in underground repositories for radioactive waste disposal, ground-source heat pump (GSHP) systems, evapotranspirative (ET) covers and pavement systems.. Accurate modeling of the temperature effect on the SWCC may lead to reduction in design costs, simpler constructability, and hence, more sustainable structures. . The study made use of two possible approaches to assess the temperature effect on the SWCC. In the first approach, soils were sorted from a large soil database into families of similar properties but located on sites with different MAAT. The SWCCs were plotted for each family of soils. Most families of soils showed a clear trend indicating the influence of temperature on the soil water retention curve at low degrees of saturation.. The second approach made use of statistical analysis. It was demonstrated that the suction increases as the MAAT decreases. The statistical analysis showed that even though the plasticity index proved to have the greatest influence on suction, the mean annual air temperature effect proved not to be negligible. In both approaches, a strong relationship between temperature, suction and soil properties was observed. Finally, a comparison of the model based on the mean annual air temperature environmental factor was compared to another model that makes use of the Thornthwaite Moisture Index (TMI) to estimate the environmental effects on the suction of unsaturated soils. Results showed that the MAAT can be a better indicator when compared to the TMI found but the results were inconclusive due to the lack of TMI data available.

A method for evaluating the integrity of geosynthetic elements of a waste containment system subject to seismic loading is developed using a large strain finite difference numerical computer program. The method accounts for the effect of interaction between the geosynthetic elements and the overlying waste on seismic response and allows for explicit calculation of forces and strains in the geosynthetic elements. Based upon comparison of numerical results to experimental data, an elastic-perfectly plastic interface model is demonstrated to adequately reproduce the cyclic behavior of typical geomembrane-geotextile and geomembrane-geomembrane interfaces provided the appropriate interface properties are used. New constitutive models are developed for the in-plane cyclic shear behavior of textured geomembrane/geosynthetic clay liner (GMX/GCL) interfaces and GCLs. The GMX/GCL model is an empirical model and the GCL model is a kinematic hardening, isotropic softening multi yield surface plasticity model. Both new models allows for degradation in the cyclic shear resistance from a peak to a large displacement shear strength. The ability of the finite difference model to predict forces and strains in a geosynthetic element modeled as a beam element with zero moment of inertia sandwiched between two interface elements is demonstrated using hypothetical models of a heap leach pad and two typical landfill configurations. The numerical model is then used to conduct back analyses of the performance of two lined municipal solid waste (MSW) landfills subjected to strong ground motions in the Northridge earthquake. The modulus reduction "backbone curve" employed with the Masing criterion and 2% Rayleigh damping to model the cyclic behavior of MSW was established by back-analysis of the response of the Operating Industries Inc. landfill to five different earthquakes, three small magnitude nearby events and two larger magnitude distant events. The numerical back analysis was able to predict the tears observed in the Chiquita Canyon Landfill liner system after the earthquake if strain concentrations due to seams and scratches in the geomembrane are taken into account. The apparent good performance of the Lopez Canyon landfill geomembrane and the observed tension in the overlying geotextile after the Northridge event was also successfully predicted using the numerical model.

In this work, the vapor transport and aerobic bio-attenuation of compounds from a multi-component petroleum vapor mixture were studied for six idealized lithologies in 1.8-m tall laboratory soil columns. Columns representing different geological settings were prepared using 20-40 mesh sand (medium-grained) and 16-minus mesh crushed granite (fine-grained). The contaminant vapor source was a liquid composed of twelve petroleum hydrocarbons common in weathered gasoline. It was placed in a chamber at the bottom of each column and the vapors diffused upward through the soil to the top where they were swept away with humidified gas. The experiment was conducted in three phases: i) nitrogen sweep gas; ii) air sweep gas; iii) vapor source concentrations decreased by ten times from the original concentrations and under air sweep gas. Oxygen, carbon dioxide and hydrocarbon concentrations were monitored over time. The data allowed determination of times to reach steady conditions, effluent mass emissions and concentration profiles. Times to reach near-steady conditions were consistent with theory and chemical-specific properties. First-order degradation rates were highest for straight-chain alkanes and aromatic hydrocarbons. Normalized effluent mass emissions were lower for lower source concentration and aerobic conditions. At the end of the study, soil core samples were taken every 6 in. Soil moisture content analyses showed that water had redistributed in the soil during the experiment. The soil at the bottom of the columns generally had higher moisture contents than initial values, and soil at the top had lower moisture contents. Profiles of the number of colony forming units of hydrocarbon-utilizing bacteria/g-soil indicated that the highest concentrations of degraders were located at the vertical intervals where maximum degradation activity was suggested by CO2 profiles. Finally, the near-steady conditions of each phase of the study were simulated using a three-dimensional transient numerical model. The model was fit to the Phase I data by adjusting soil properties, and then fit to Phase III data to obtain compound-specific first-order biodegradation rate constants ranging from 0.0 to 5.7x103 d-1.

Due to the lack of understanding of soil thermal behavior, rules-of-thumb and generalized procedures are typically used to guide building professionals in the design of ground coupled heat pump systems. This is especially true when sizing the ground heat exchanger (GHE) loop. Unfortunately, these generalized procedures often encourage building engineers to adopt a conservative design approach resulting in the gross over-sizing of the GHE, thus drastically increasing their installation cost. This conservative design approach is particularly prevalent for buildings located in hot and arid climates, where the soils are often granular and where the water table tends to exist deep below the soil surface. These adverse soil conditions reduce the heat dissipation efficiency of the GHE and have hindered the adoption of ground coupled heat pump systems in such climates. During cooling mode operation, heat is extracted from the building and rejected into the ground via the GHE. Prolonged heat dissipation into the ground can result in a coupled flow of both heat and moisture, causing the moisture to migrate away from the GHE piping. This coupled flow phenomenon causes the soil near the GHE to dry out and results in the degradation of the GHE heat dissipation capacity. Although relatively simple techniques of backfilling the GHE have been used in practice to mitigate such coupled effects, methods of improving the thermal behavior of the backfill region around the GHE, especially in horizontal systems, have not been extensively studied. This thesis presents an experimental study of heat dissipation from a horizontal GHE, buried in two backfill materials: (1) dry sand, and (2) wax-sand composite mixture. The HYDRUS software was then used to numerically model the temperature profiles associated with the aforementioned backfill conditions, and the influence of the contact resistance at the GHE-backfill interface was studied. The modeling strategy developed in HYDRUS was proven to be adequate in predicting the thermal performance of GHE buried in dry sand. However, when predicting the GHE heat dissipation in the wax-sand backfill, significant discrepancies between model prediction and experimental results still exist even after calibrating the model by including a term for the contact resistance. Overall, the thermal properties of the backfill were determined to be a key determinant of the GHE heat dissipation capacity. In particular, the wax-sand backfill was estimated to dissipate 50-60% more heat than dry sand backfill.

Expansive soils impose challenges on the design, maintenance and long-term stability of many engineered infrastructure. These soils are composed of different clay minerals that are susceptible to changes in moisture content. Expansive clay soils wreak havoc due to their volume change property and, in many cases, exhibit extreme swelling and shrinking potentials. Understanding what type of minerals and clays react in the presence of water would allow for a more robust design and a better way to mitigate undesirable soil volume change. The relatively quick and widely used method of X-ray Diffraction (XRD) allows identifying the type of minerals present in the soil. As part of this study, three different clays from Colorado, San Antonio Texas, and Anthem Arizona were examined using XRD techniques. Oedometer-type testing was simultaneously preformed in the laboratory to benchmark the behavior of these soils. This analysis allowed performing comparative studies to determining if the XRD technique and interpretation methods currently available could serve as quantitative tools for estimating swell potential through mineral identification. The soils were analyzed using two different software protocols after being subjected to different treatment techniques. Important observations include the formation of Ettringite and Thaumasite, the effect of mixed-layer clays in the interpretation of the data, and the soils being subject to Gypsification. The swelling data obtained from the oedometer-type laboratory testing was compared with predictive swelling functions available from literature. A correlation analysis was attempted in order to find what index properties and mineralogy parameters were most significant to the swelling behavior of the soils. The analysis demonstrated that Gypsification is as important to the swelling potential of the soil as the presence of expansive clays; and it should be considered in the design and construction of structures in expansive soils. Also, the formation of Ettringite and Thaumasite observed during the treatment process validates the evidence of Delayed Ettringite Formation (DEF) reported in the literature. When comparing the measured results with a proposed method from the University of Texas at Arlington (UTA), it was found that the results were somewhat indicative of swell potential but did not explain all causes for expansivity. Finally, it was found that single index properties are not sufficient to estimate the free swell or the swell pressure of expansive soils. In order to have a significant correlation, two or more index properties should be combined when estimating the swell potential. When properties related to the soil mineralogy were correlated with swell potential parameters, the amount of Gypsum present in the soil seems to be as significant to the swell behavior of the soil as the amount of Smectite found.

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.

The dissimilatory reduction of nitrate, or denitrification, offers the potential of a sustainable, cost effective method for the non-disruptive mitigation of earthquake-induced soil liquefaction. Worldwide, trillions of dollars of infrastructure are at risk for liquefaction damage in earthquake prone regions. However, most techniques for remediating liquefiable soils are either not applicable to sites near existing infrastructure, or are prohibitively expensive. Recently, laboratory studies have shown the potential for biogeotechnical soil improvement techniques such as microbially induced carbonate precipitation (MICP) to mitigate liquefaction potential in a non-disruptive manner. Multiple microbial processes have been identified for MICP, but only two have been extensively studied. Ureolysis, the most commonly studied process for MICP, has been shown to quickly and efficiently induce carbonate precipitation on particle surfaces and at particle contacts to improve the stiffness, strength, and dilatant behavior of liquefiable soils. However, ureolysis also produces copious amounts of ammonium, a potentially toxic byproduct. The second process studied for MICP, denitrification, has been shown to precipitate carbonate, and hence improve soil properties, much more slowly than ureolysis. However, the byproducts of denitrification, nitrogen and carbon dioxide gas, are non-toxic, and present the added benefit of rapidly desaturating the treated soil. Small amounts of desaturation have been shown to increase the cyclic resistance, and hence the liquefaction resistance, of liquefiable soils. So, denitrification offers the potential to mitigate liquefaction as a two-stage process, with desaturation providing short term mitigation, and MICP providing long term liquefaction resistance. This study presents the results of soil testing, stoichiometric modeling, and microbial ecology characterization to better characterize the potential use of denitrification as a two-stage process for liquefaction mitigation.

Up to 25 percent of the operating budget for contaminated site restoration projects is spent on site characterization, including long-term monitoring of contaminant concentrations. The sensitivity, selectivity, and reproducibility of analytical methods have improved to the point where sampling techniques bear the primary responsibility for the accuracy and precision of the data. Most samples represent discrete concentrations in time and space; with sampling points frequently limited in both dimensions, sparse data sets are heavily extrapolated and the quality of data further limited.

Methods are presented for characterizing contaminants in water (groundwater and surface waters) and indoor air. These techniques are integrative, providing information averaged over time and/or space, as opposed to instantaneous point measurements. Contaminants are concentrated from the environment, making these methods applicable to trace contaminants. These methods have the potential to complement existing techniques, providing the practitioner with opportunities to reduce costs and improve the quality of the data used in decision making.

A conceptual model for integrative sampling of environmental waters is developed and a literature review establishes an advantage in precision for active samplers. A programmable sampler was employed to measure the concentration of chromate in a shallow aquifer exhibiting time-dependent contaminant concentrations, providing a unique data set and sustainability benefits. The analysis of heat exchanger condensate, a waste stream generated by air conditioning, is demonstrated in a non-intrusive method for indoor air quality assessment. In sum, these studies present new opportunities for effective, sustainable environmental characterization.

Nanotechnology has been applied to many areas such as medicine, manufacturing, catalysis, food, cosmetics, and energy since the beginning 21st century. However, the application of nanotechnology to geotechnical engineering has not received much attention. This research explored the technical benefits and the feasibility of applying nanoparticles in geotechnical engineering. Specific studies were conducted by utilizing high-pressure devices, axisymmetric drop shape analysis (ADSA), microfluidics, time-lapse technology, Atomic Force Microscopy (AFM) to develop experiments. The effects of nanoparticle on modifying interfacial tension, wettability, viscosity, sweep efficiency and surface attraction forces were investigated. The results show that nanoparticles mixed in water can significantly reduce the interfacial tension of water in CO2 in the applications of nanofluid-CO2 flow in sediments; nanoparticle stabilized foam can be applied to isolate contaminants from clean soils in groundwater/soil remediation, as well as in CO2 geological sequestration or enhanced oil/gas recovery to dramatically improve the sweep efficiency; nanoparticle coatings are capable to increase the surface adhesion force so as to capture migrating fine particles to help prevent clogging near wellbore or in granular filter in the applications of oil and gas recovery, geological CO2 sequestration, geothermal recovery, contaminant transport, groundwater flow, and stormwater management system.