Increased land use change near fragile ecosystems can affect the ecosystem energy balance leading to increased global warming. One component of surface energy balance is soil storage heat flux. In past work, a complex thermal behavior was noticed in the shrink-swell sediment of the upper Nueces Delta (upper Rincon) during summer months as it dried. Soil storage heat flux was found to first increase, then decrease, as the soil dried. It was suggested that the complex behavior was due to the relationship between thermal diffusivity and soil moisture, where thermal diffusivity increases to a local maximum before decreasing with respect to decreasing soil moisture. This study explores the observed phenomenon in a controlled laboratory environment by relating the sediment shrinkage curve to changing heat transfer properties. Due to the complicated nature of the drying-shrinking sediment, it was necessary to measure the sediment shrinkage curve and heat transfer properties in separate experiments. The shrinkage curve was found by correlating measured sample volume with gravimetric moisture content. Heat transfer properties were found using a single needle heat pulse probe. A normalized gravimetric moisture content was used as a common variable to relate the shrinkage curve and heat transfer data. Data suggests that the shrink-swell Rincon sediment portrays different behavior in drying than that which occurs for a non-shrink-swell soil. For the shrink-swell Rincon sediment, thermal conductivity is seen to increase with decreasing moisture, the suggested mechanism being increased surface area contact between particles as the shrinking sediment dries.

Problems associated with shrink-swell soils are well known geotechnical problems that have been studied and researched by many geotechnical researchers for many decades. Potentially shrink-swell soils can be found almost anywhere in the world especially in the semi-arid regions of the tropical and temperate climate. Foundation slabs on grade on shrink-swell soils are one of the most efficient and inexpensive solutions for this kind of problematic soil. It is commonly used in residential foundations or any light weight structure on shrink-swell soils. Many design methods have been established for this specific problem such as Building Research Advisory Board (BRAB), Wire Reinforcement Institute (WRI), Post- Tensioning Institute (PTI), and Australian Standards (AS 2870) design methods. This research investigates most of these methods, and then, proposes a moisture diffusion soil volume change model, a soil-weather interaction model, and a soil-structure interaction model. The proposed moisture diffusion soil volume change model starts with proposing a new laboratory test to determine the coefficient of unsaturated diffusivity for intact soils. Then, it introduces the development of a cracked soil diffusion factor, provides a chart for it, and explains a large scale laboratory test that verifies the proposed moisture diffusion soil volume change model. The proposed soil-weather interaction model uses the FAO 56-PM method to simulate a weightless cover performance for six cities in the US that suffer significantly from shallow foundation problems on shrink-swell soils due to seasonal weather variations. These simulations provide more accurate weather site-specific parameters such as the range of surface suction variations. The proposed weather-site specific parameters will be input parameters to the soil structure models. The proposed soil-structure interaction model uses Mitchell (1979) equations for moisture diffusion under covered soil to develop a new closed form solution for the soil mound shape under the foundation slab. Then, it presents a parametric study by carrying out several 2D finite elements plane strain simulations for plates resting on a semiinfinite elastic continuum and resting on different soil mounds. The parametric study outcomes are then presented in design charts that end with a new design procedure for foundation slabs on shrink-swell soils. Finally, based on the developed weather-soil-structure interaction models, this research details two procedures of a proposed new design method for foundation slabs on grade on shrink-swell soils: a suction based design procedure and a water content based design procedure.

A study of 12 expansive soils in four major physiographic provinces in Virginia was initiated to examine and quantify the relationship between shrink-swell potential, shrink-swell indices, and soil properties. The mineralogy classes, soil series, and (physiographic provinces, parent materials) examined include smectitic -- Jackland and Waxpool (Triassic, diabase), Iredell (Piedmont, hornblende); vermiculitic -- Kelly (Triassic, thermal shale); kaolinitic -- Cecil (Piedmont, granite gneiss), Davidson (Triassic, diabase); and mixed -- Carbo and Frederick (Valley and Ridge, limestone), Craven and Peawick (Coastal Plain, fluvial and marine sediments), and Mayodan and Creedmoor (Triassic, sandstones). Three sites in each of the 12 map units were described and major horizons sampled for physical, chemical, and mineralogical laboratory analysis. An expansive soil rating system, termed the Expansive Soil Index (ESI), was developed using the soil properties best correlated with shrink-swell potential. The sum of swelling 2:1 minerals, swell index, liquid limit, and CEC gave expansive soil potential ratings (ESI) for each soil series. The higher the ESI, the greater the shrink-swell potential. Smectite distributions within the soil profiles were investigated. Smectite concentration in the clay fraction increases with depth in soils formed from diabase and thermally altered shale. Smectite weathers to kaolinite and hydroxy-interlayered vermiculite with increasing proximity to the soil surface thus accounting for the observed decrease in smectite toward the soil surface. The highest amount of smectite from the granite gneiss, limestone, sandstones and shales, and Coastal Plain sediments were in the Bt2 horizon where maximum expression of the argillic horizon occurs. Smectite contents decrease away (upwards and downwards) from the maximum in the Bt2 horizon. A satellite study focused on locating and quantifying the variability within five map units in the Culpeper (Triassic) Basin in northern Virginia. Variability of the shrink-swell indices and related properties are high in all map units. Dissimilar inclusions could adversely affect foundations if a home is sited on both moderate and high shrink-swell soils. Although there is extreme variability in the map units, the variability occurs within the delineations of each map unit. Each delineation within an individual map unit contains similar levels of variability. / Ph. D.

Because of the dynamic nature of shrinking and swelling of soils that are classified as Vertisols, partitioning of rainfall into infiltration and runoff in a Vertic watershed is more temporally and spatially unique than in most other watersheds. Hydrology models that account for realistic representation of crack dynamics are rarely used because the spatial and temporal patterns of cracking across a catena and under different land uses are poorly understood. The objectives of the study were to 1) determine if variability in soil cracking on a Vertisol catena, having the same soil and land cover, could be explained by shrink-swell potential of the soil and changes in soil water content; 2) characterize the temporal and spatial variability of the shrinkage of a Vertisol under different land uses; and 3) determine the relationship between specific volume and water content of soils, particularly between saturation and field capacity. The research was conducted in Vertisol catenae of the Houston Black and Heiden soil series. The catenae were located within the USDA-ARS Grassland, Soil and Water Research Laboratory, Riesel Texas. Soil samples were taken to characterize the general properties of the soils. In situ bi-weekly measurements of vertical soil movements and soil water contents were made over a two-year span. Because shrink-swell potential was high at most landscape positions, soil water content was the primary factor driving the spatial and temporal variability of soil shrinking and swelling. The measured relationship between the amount of soil subsidence and water loss generally agreed with what would be theoretically expected. Maximum soil subsidence was 120 mm in the grazed pasture, 75 mm in the native prairie, and 76 mm in the row cropped field. Shrinkage of the whole soil was not equidimensional, and the study generally indicates more horizontal shrinkage than vertical shrinkage. Laboratory analysis showed an appreciable change in volume of soils between saturation and field capacity, suggests a layer of soil layer can subside up to 4% while drying from saturation to field capacity, which indicates the common laboratory measure of shrink swell potential does not capture the complete shrink-swell behavior of soils.

Plugged and abandoned well pads throughout the Uintah Basin face reclamation challenges due to factors including a harsh climate, invasive species, and high salt loads. Finding ways to alleviate soil sodicity could improve soil reclamation success. Gypsum, sulfur, activated carbon, and Biochar are being applied to improve soil parameters negatively impacted by sodicity, but the direct impact of these amendments on Uintah Basin soils is still largely unknown. The aim of this study was two-fold. (1) Evaluate the effectiveness of gypsum, sulfuric acid, Biochar, activated carbon, and combinations of these amendments in reducing the impact of soil sodicity of the Desilt and Conglomerate soils by measuring amendment impact on percent dispersion, saturated hydraulic conductivity, crust bulk density, infiltration, and crust formation. (2) Compare a crust bulk density method using ImageJ to the clod wax density method and a modified linear extensibility percent equation to the linear extensibility percent equation to assess whether the novel methods can be used to accurately measure and calculate soil crust bulk density and shrink swell potential while reducing human error and analysis time.